Abstract
Fungal infections are an increasing cause of morbidity and mortality. Current antifungal drugs are limited in spectrum, few new drugs are in development, and resistance is an increasing issue. Drug synergy can enhance available drugs and extend their lifetime, however, few synergistic combinations are in clinical use and mechanistic data on how combinations work is lacking. The multifunctional glycoprotein lactoferrin (LF) acts synergistically with amphotericin B (AMB) in a range of fungal species. Whole LF binds and sequesters iron, and LF can also be digested enzymatically to produce cationic peptides with distinct antimicrobial functions. To understand how LF synergizes AMB, we previously undertook a transcriptomic analysis in Saccharomyces and found a paradoxical down-regulation of iron and stress response, suggesting stress pathway interference was dysregulating an appropriate response, resulting in cell death. To extend this to a fungal pathogen, we here perform the same analysis in Cryptococcus neoformans. While both fungi responded to AMB in a similar way, the addition of LF produced remarkably contrasting results, with the Cryptococcus transcriptome enriched for processes relating to cellular stress, up-regulation of endoplasmic-reticulum-associated protein degradation (ERAD), stress granule disassembly and protein folding, endoplasmic reticulum-Golgi-vacuole trafficking and autophagy, suggesting an overall disruption of protein and lipid biosynthesis. These studies demonstrate that the mechanism of LF-mediated synergy is species-specific, possibly due to differences in the way LF peptides are generated, bind to and enter cells and act on intracellular targets, illustrating how very different cellular processes can underlie what appears to be a similar phenotypic response.
Introduction
Pathogenic Cryptococcus species belonging to the Cryptococcus neoformans and Cryptococcus gattii species complexes can cause invasive disease in both healthy and immunocompromised people (). Anti-cryptococcal therapy involves induction with amphotericin B (AMB) and 5-flucytosine (5-FC), followed by maintenance on fluconazole (FLC) (). However, mortality remains high, AMB and 5-FC are toxic, and 5-FC is not registered in many countries where the burden of cryptococcosis is highest (; ; ).
Although it is widely recognized that there is an urgent need for new antifungal therapies, development is hindered by high risk and cost combined with relatively limited revenue (). An alternative approach that aims to get more life out of existing drugs is drug synergy, where an antifungal is combined with a second drug or agent. The AMB-5FC combination used for cryptococcosis is a model example of improved antifungal efficacy, however, synergistic combinations are not widely employed in medical mycology ().
Previously, we tested a range of commonly used antifungal drugs with various iron chelating agents to determine whether iron limitation might enhance antifungal potency. Synergistic combinations were rare, however, the combination of AMB with lactoferrin (LF), a mammalian glycoprotein, produced a synergistic response in Saccharomyces cerevisiae and various Cryptococcus species, which had also been reported in Candida and Aspergillus (; ). We further demonstrated that AMB-LF synergy was not due to iron starvation, as the addition of iron rescued the yeast cells from inhibition by LF but not from AMB-LF (). To understand the mechanism of synergy, we performed comparative analyses of the transcriptomic responses to AMB and AMB-LF in S. cerevisiae. Paradoxically, we found that while AMB alone increased the expression of stress- and metal-related genes, the AMB-LF combination caused metal and stress-response transcripts to decrease. We hypothesized that AMB-LF synergy was mediated by dysregulation of metal homeostasis and disrupted stress responses, possibly via stress pathway interference ().
To follow up this analysis in the context of a pathogenic yeast, we here report a companion transcriptomic analysis of AMB-LF synergy in C. neoformans strain H99. We show that, in contrast to S. cerevisiae, synergy is mediated by an enhanced stress response in Cryptococcus that includes disruption of endoplasmic reticulum-related processes. This study highlights the species-specific nature of responses to toxic agents by different fungal species, even when these appear to generate similar outcomes.
Materials and Methods
Strains, Culture and Agents
Cryptococcus neoformans strain H99 was cultured in yeast nitrogen broth (YNB) at 37°C with shaking at 180 rpm. Stock solutions of 1,600 μg/mL AMB (Sigma-Aldrich, United States) and 5,120 μg/mL LF (MP Biomedical, United States) were made in MilliQ water.
Time Kill-Curves and RNA Extraction
Time-kill curves for AMB only, LF only and AMB-LF treatment were performed to determine ID20 (inhibition of cell growth by 20%) for RNA extraction and comparative transcriptomic analyses [as performed in ]. Synergistic fractional inhibitory concentrations (FIC) of AMB-LF for C. neoformans [determined in ] were used. ID20 was determined to be 1 h for AMB and 50 min for AMB-LF treatment. At these timepoints LF treatment did not affect growth (Supplementary Figure S1). RNA was extracted from the following cultures: (i) AMB treatment; (ii) growth control for (i); (iii) AMB-LF treatment; and (iv) growth control for (iii). Three independent biological replicates were performed for each treatment.
RNA was extracted from freeze dried cells using bead beating and the Qiagen RNeasy mini kit as detailed in . RNA was sequenced at the Ramaciotti Centre for Genomics using an Illumina HiSeq 2000 to generate 100 bp paired-end reads.
Processing of RNA-Seq Data
Data analysis was performed as in our previous study () but with updated software. Briefly, reads were mapped using the HISAT tool (; version 2.1.0) to the C. neoformans var. grubii H99 reference genome (RefSeq accession: GCF_000149245.1). Reads counting was performed with the featureCount function of Subreads (; ; version 3.5.3) with the edgeR library (; version 3.24.3) and significant differentially expressed genes were defined as those with adjusted p-value < 0.05. RUVSeq (; version 1.16.1) was used to remove unwanted variation. After testing with zero to six unwanted factors with RUVSeq, the number of unwanted factors that maximized the number of significantly differentially expressed genes in total for each fungal species was used.
Protein Orthologs and Gene Ontology Annotation
One-to-one protein orthologs from S. cerevisiae and C. neoformans were defined as reciprocal best hits from BLASTp searches (E-value < 1 × 10–6). The log fold-change of the unique ortholog pairs was graphed, a linear trend line was fitted, and Pearson’s correlation scores were calculated. Gene Ontology (GO) annotations were obtained from Uniprot (; release 2019_02), FungiDB () and QuickGO (; downloaded 10th April 2019). Additional GO annotations were obtained using Blast2GO (; version 5.25) and the Diamond sequence similar search tool (; version 0.9.24) was used to compare the C. neoformans proteome against the Uniprot/Swiss-Prot database (; release 2019_02). Further annotations were obtained by transferring GO terms from orthologs of other fungal species using OrthoMCL ().
Identifying Significantly Differentially Expressed Genes Unique to AMB or AMB-LF Treatments
Sets analysis was used to identify two mutually exclusive subsets of genes with significant differential expression unique to either AMB treatment or AMB-LF treatment. Each of these was then sub-divided to include genes with positive or negative log fold-change compared to the untreated control. The number of genes in each subset was represented as an UpSet plot (). Each of the four subsets of genes was independently analyzed for enriched GO biological process terms using the methods described below.
Gene Ontology Enrichment Analysis and Self-Organizing Maps
The “GOstats” R package and Fisher’s exact test (; version 2.48.0) were used for all GO enrichment analyses. Enriched GO terms were defined as those with Benjamini–Hochberg adjusted p-values < 0.05. For the global GO enrichment analysis, genes with increased or decreased expression were independently analyzed for enriched GO terms. Self-organizing map (SOM) analysis was used to identify clusters of significantly differentially expressed genes (adjusted p-values < 0.05) with similar expression profiles across the different treatments. The read counts in log counts per million were scaled to a mean of 0 and standard deviation of 1 across all samples per gene. Five-by-five SOM clusters were calculated using the R “kohonen” library (; version 3.0.8). Co-expressed genes from each SOM cluster were analyzed for enriched GO terms. To perform exploratory data analyses on the co-expression relationships and annotations of significantly differentially expressed genes (adjusted p-values < 0.05) in each SOM cluster, GO terms with moderate confidence of unadjusted p-values < 0.01 in each SOM cluster were analyzed in addition to the those that passed the stringent threshold of adjusted p-values < 0.05.
Codes and Data Availability
Cryptococcus neoformans H99 RNA-Seq data generated from this study have been deposited in the NCBI Gene Expression Omnibus (GEO) () and are accessible through GEO Series accession number GSE130375. S. cerevisiae S288C RNA-Seq data were previously described in and are accessible through GEO Series accession number GSE80357. The scripts for all the above analyses are available from Github: https://github.com/IgnatiusPang/Fungal_Drug_Synergy under the GNU Lesser General Public License v3.0.
Results
Differential Expression of C. neoformans Genes in Response to AMB and AMB-LF Treatment
We chose ID20 as the point to undertake transcriptome analysis, where cells should be responding to the primary effects of drug-related stress (). C. neoformans cultures were treated with AMB and AMB-LF and harvested at their ID20, (60 min for AMB and 50 min for AMB-LF) along with their matched controls (Supplementary Figure S1). Initial QC analysis found one biological replicate of AMB-LF treatment was an outlier and it was excluded from subsequent analyses. In final analysis, the transcriptome of Cryptococcus following AMB treatment had 1625 up-regulated and 1957 down-regulated genes, relative to its matched control. Following AMB-LF treatment, 1599 genes were up-regulated while 1860 genes were down-regulated, compared to the matched control (Supplementary Table S1).
AMB-LF Synergy Induces a Different Response in C. neoformans H99 and S. cerevisiae S288C
We previously performed a transcriptomic study of the response of S. cerevisiae to drug synergy induced by the addition of LF to AMB () (Supplementary Table S2). To compare this with the response by C. neoformans, genes that were differentially expressed following AMB and AMB-LF treatment (Supplementary Table S3) were analyzed for functional enrichment by GO () (Supplementary Tables S4, S5). In C. neoformans, the responses were similar across both treatments, with GO terms related to energy production, response to oxidative stress, carbohydrate metabolism, nucleic acid metabolism, protein processing and metal ion transport induced, and terms related to microtubule organization and lipid biosynthesis repressed. This contrasted with the situation in S. cerevisiae, where synergy between AMB and LF resulted in stress and metal homeostasis responses that were quite different to the effect of AMB treatment alone ().
Comparison of the response to AMB treatment found the majority of enriched biological processes that were common to both species had the same direction of regulation (Figure 1A). Enrichments related to energy generation, metal regulation and stress responses, amino acid biosynthesis and protein catabolism were all induced, while lipid metabolism was repressed. Functional enrichments that were absent from C. neoformans but were present in S. cerevisiae in response to AMB included the induction of autophagy, actin polymerization, iron homeostasis and copper transport, cell wall integrity MAPK cascade and trehalose metabolism in response to stress. Down-regulated enrichments that were observed only in S. cerevisiae included ergosterol biosynthesis, nucleobase metabolism, protein translation and folding-related processes.
FIGURE 1
In contrast to AMB, the transcriptional response to AMB-LF was drastically different in C. neoformans and S. cerevisiae. Functional enrichments related to energy generation processes, protein processing, response to stress and pyridine nucleotide metabolism were up-regulated in C. neoformans but down-regulated in S. cerevisiae, while amino acid metabolism was down-regulated in C. neoformans but up-regulated in S. cerevisiae (Figure 1A). Biological processes that were enriched in C. neoformans only included the up-regulation of apoptosis, transition metal transport, protein catabolism in the vacuole, stress granule assembly and hydrogen peroxide metabolism related to stress responses, and the down-regulation of microtubule organization. Enrichments that were absent in C. neoformans and specific to S. cerevisiae in response to AMB-LF were the induction of iron-sulfur cluster assembly, nucleobase metabolism and translation-related processes, and the repression of ER-associated processes including protein folding and ubiquitination, protein membrane localization, response to misfolded proteins and ER stress.
To verify that responses to AMB treatment were similar while responses to AMB-LF treatment differed between C. neoformans and S. cerevisiae, the log fold-change of genes expressed in C. neoformans and S. cerevisiae were plotted against one another and the correlation between genes that were orthologous between the two species were analyzed. The log fold-changes of orthologous genes following AMB treatment plotted in a positively correlated manner (ρ = 0.38, p < 0.05) suggesting similar responses (Figure 1B, left panel),while those in response to AMB-LF treatment were negatively correlated (ρ = −0.13, p-value < 0.05), demonstrating an overall opposing response by the two species to synergistic treatment (Figure 1B, right panel).
Comparison of Gene Ontology Enrichments in C. neoformans Transcriptome Data Produced Following AMB and AMB-LF Treatment
While there were many similarities in the regulation of functional enrichments observed between AMB and AMB-LF treatments in C. neoformans, there were some GO enrichments that were unique to AMB-LF synergy (Figure 1A). These included down-regulation of ergosterol biosynthesis and aromatic amino acid family biosynthetic processes, and up-regulation of protein folding, stress granule disassembly and apoptotic processes. Cellular amino acid biosynthetic processes were up-regulated in response to AMB but down-regulated by AMB-LF. The functional enrichment that was seen only in AMB treatment was stress response to metal ion, and this was induced and involved genes encoding ABC transporters, a bile acid transporter, a zinc transporter and sulfide reductase (Supplementary Table S4).
To further tease out differences between AMB and AMB-LF, genes that were differentially expressed and present only in the response to AMB or to AMB-LF were separately analyzed for GO enrichments and are presented visually in an UpSets plot (Figure 2 and Supplementary Tables S3, S4). The sets analysis found a majority of differentially expressed genes to be shared between the AMB and AMB-LF treatment and identified transcripts that were unique to each (Figure 2A). Each group containing unique transcripts was analyzed to identify enriched biological process GO terms and their level of enrichment (Figure 2B). Functionally enriched GO terms that encompassed AMB-LF-specific genes mainly involved down-regulated processes related to protein synthesis (ribosome biogenesis and translation), amino acid biosynthesis (nitrogen compound and S-adenosyl methionine metabolism), and nucleic acid metabolism (RNA transport and response to DNA damage and repair). Up-regulated enrichments were mainly involved in actin polymerization and organization, membrane fusion (endocytosis, SNARE complex assembly, vesicle, and vacuolar fusion) and vesicle mediated transport, as well as ER and protein related processes (protein folding and proteolysis) and cell wall biosynthesis (glucan and melanin synthesis). Genes that were differentially regulated only in response to AMB had enrichments related to the down-regulation of RNA processing, protein folding and ER-to-Golgi transport (Figure 2). There were no biological enrichments identified among the AMB treatment up-regulated genes (Supplementary Table S4).
FIGURE 2
SOMs Analysis Indicates AMB-LF Disrupts Iron Transmembrane Transport in C. neoformans
To investigate whether differentially expressed genes were co-expressed, we grouped transcripts with similar expression patterns using SOMs and analyzed these for the enrichment of biological process GO terms (Figure 3 and Supplementary Tables S6, S7). From this, we focused on clusters that included responses related to metal homeostasis, stress and ER-associated functions (Figure 3B), which we had found to differ between AMB and AMB-LF treatment (Figures 1, 2).
FIGURE 3
Self-organizing maps analysis indicated a general up-regulation of genes in enrichments related to metal homeostasis for both AMB and AMB-LF treatments. Enrichments related to copper and metal ion transport were seen for both treatments, with the genes making up the enrichments encoding iron, zinc and copper transporters (Figure 3B). There was down-regulation of genes related to iron ion transmembrane transport in AMB-LF that was not seen with AMB alone. These genes were co-expressed with protein and ergosterol synthesis in cluster 10 of the AMB-LF SOM (Figure 3B), suggesting a disruption in the synthesis of iron transport proteins.
AMB-LF Synergy Increases Global Cellular Stress in C. neoformans and This Is Associated With Golgi and ER-Related Processes
For both AMB and AMB-LF treatment, analysis of co-expressed enrichments found stress-related responses to be induced and enriched across multiple different SOMs clusters. For example, in the AMB response terms related to free radical detoxification, oxidative and nitrosative stress and glutathione metabolism were enriched in clusters 10, 14, 15, 19, and 24, and similar terms were enriched in clusters 1, 7, 11, 16, and 17 in the AMB-LF response (Figure 3B). This suggests that stress is generalized across the cell and is not associated with any particular biological process.
Similar to previous studies, our analysis indicated that AMB induces responses linked to oxidative and nitrosative stress in C. neoformans (). ER-dependent peroxisome organization, where peroxisomes form in response to hydrogen peroxide production from fatty oxidation, was also induced in response to AMB (cluster 19) along with autophagy (cluster 20) and apoptotic processes (cluster 15).
Oxidative and nitrosative stress, glutathione metabolism and autophagy were likewise induced in response to AMB-LF (clusters 1, 7, 11, 16, 17, and 22), however, these and other stress responses were co-enriched with processes related to protein synthesis and transit through the ER (Figure 3B). Induction of the enrichment related to stress granules (cluster 3), which are assemblies of untranslated messenger ribonucleoproteins (), suggests stress may be resulting from disrupted protein translation. ER stress was further suggested in cluster 3 by induced enrichments related to topologically incorrect proteins and the ERAD (ER-associated protein degradation) response, which is activated by the accumulation of misfolded proteins (). Across other SOMs clusters, enrichments related to Golgi-related protein processing and trafficking to the vacuole, such as Golgi vesicle docking (cluster 1), glycoprotein biosynthesis (cluster 7), endosome vacuole transport (cluster 16), and protein catabolism in the vacuole (cluster 11) were induced and associated with stress responses. The overall relation of stress responses across different clusters and biological processes indicates that AMB-LF treatment induces substantially more stress than AMB alone. Although enrichments related to the induction of apoptosis were not evident in the AMB-LF SOMs data, apoptosis was up-regulated in the functional enrichment profile for this treatment overall (Figure 1), further supporting the notion that the addition of LF to AMB substantially enhances cellular stress.
Addition of LF to AMB Disrupts Biosynthetic, Folding, Sorting and Transport Processes Undertaken by the ER and Golgi and Induces Protein Trafficking
ER- and Golgi-related processes were enriched across different SOMs clusters following both AMB and AMB-LF treatment (Figure 3). For AMB treatment, these enrichments likely reflect a down-regulation in protein sorting and trafficking through the secretory pathway, and include localization of protein to the ER (cluster 11), ER processes such as glycerolipid biosynthesis, protein folding and signal peptide processing for translocation (cluster 11, 12, and 22), and vesicle mediated transport to the Golgi (cluster 16). This AMB-mediated repression of protein biosynthesis and sorting processes in C. neoformans is consistent with the transcriptome responses to AMB seen for S. cerevisiae and C. albicans (; ; ). Other enrichments associated with proteasomal protein degradation and localization to the vacuole membrane were also down-regulated in clusters 11 and 12.
Up- and down-regulation of ER- and Golgi-related functional enrichments were observed in response to AMB-LF treatment (Figure 3B). Similar to AMB treatment, ER- and Golgi-related processes were repressed, including synthesis of lipids, peptides, glycosylphosphatidylinositol (GPI) anchors and proteins involved in iron transport that require GPI modifications, which were all co-expressed in cluster 10. The exception was glycoprotein synthesis, which was induced in cluster 7. Unlike the response to AMB treatment, however, protein trafficking was induced and SOMs enrichments reflected three routes of transport from the ER and Golgi that were co-enriched with stress responses. ER-to-Golgi transport through the secretory pathway was suggested by enrichments involving localization of proteins to the ER and Golgi vesicle docking in clusters 1 and 3. The enrichment relating to endosome-to-vacuole transport in cluster 16 suggests Golgi-to-vacuole trafficking of proteins for storage or degradation. ER-to-cytosol transport-related enrichments were co-enriched in cluster 3 and reflect ERAD activation, which catalyzes the transport of misfolded proteins between the ER and the cytosol for proteolysis (). The difference in regulation of ER and Golgi functions seen between AMB and AMB-LF treatments suggests that the addition of LF to AMB disrupts cellular protein and lipid synthesis, with misfolded proteins inducing sorting and trafficking.
Discussion
The current study follows our analysis of AMB-LF synergy in S. cerevisiae, which was undertaken to allow a detailed analysis of synergy using the extensive genetic and bioinformatic resources that are available for yeast. We chose ID20 for both studies, where growth had been inhibited by 20% compared to the untreated control, and hypothesized that synergy would be mediated by similar mechanisms that could be used to find novel drug targets that are conserved across fungal species (). However, to our surprise we have found synergy to occur via two different and opposing mechanisms, with dysregulation of the metal ion response and down-regulation of stress in Saccharomyces, and up-regulation of cellular stress in Cryptococcus. For the former, we suggested that the independent assaults by AMB and LF resulted in stress pathway interference, preventing the S. cerevisiae cells from sensing and responding to toxic stress in an appropriate way (). For Cryptococcus, we propose that the combined agents result in excessive oxidative and nitrosative damage and an up-regulation of autophagy that results from increased ERAD and ER-Golgi-vacuole trafficking and disrupted protein and lipid biosynthesis. A model of the responses of both species to AMB and AMB-LF that summarizes the data overall and the results shown in Figures 1–3 is presented in Figure 4.
FIGURE 4
Amphotericin B has a broad spectrum of activity against yeasts and molds. AMB binds membrane ergosterol, and it is widely accepted that antifungal activity is mediated by the disruption of cellular integrity (). However, recent studies in yeast indicate that AMB can also enter cells and uses autophagy-dependent transport to the vacuole, which it then damages and disrupts (). AMB also induces oxidative and nitrosative stress (; ), and killing in C. neoformans involves induction of an oxidative burst (). Our transcriptomic data for AMB alone are in good agreement with these studies, with induction of ROS and the overall stress response and an up-regulation of autophagy and vacuole-mediated protein catabolism seen in both Cryptococcus and Saccharomyces. In addition, AMB caused an up-regulation of genes related to metal ion transport and energy production and down-regulation of cell growth, protein and ribosome production and ergosterol biosynthesis in both species (Figures 4A,C), which are strong indicators of cellular inhibition and response to membrane stress ().
Like AMB, LF appears to bind to and disrupt cell membranes, which together with the ability to sequester iron is important in its antimicrobial activity. Unbiased ‘omics approaches have not been used to study the fungal response to LF, and at the FIC levels used here LF does not affect the transcriptome (). However, analysis of the apoptotic response of Saccharomyces to LF found inhibition requires de novo protein synthesis and energy, with LF causing mitochondrial dysfunction resulting in ROS accumulation (). LF targeting of ER and Golgi processes has not been assessed in current literature. However, our finding of up-regulated terms related to ERAD and vesicle trafficking with autophagy, which were co-expressed and linked to the up-regulation of misfolded proteins and proteolysis and the accumulation of stress granules, together with the down-regulation of ergosterol, GPI anchor and iron transporter synthesis (Figure 4B), suggest disrupted protein and lipid biosynthesis are likely to be an important component of LF-mediated synergy in Cryptococcus.
Also in common with AMB, LF is a natural product with multiple biological functions (). In addition to being antimicrobial, LF has anti-tumor and anti-oxidant properties and it can both stimulate and repress the inflammatory response to infection (), presumably by differential binding to receptors on target cell surfaces and subsequent uptake (). There is emerging evidence that following binding, LF can enter microbial cells and interact directly with intracellular targets (). Furthermore, proteolysis results in the production of cationic peptides, in particular lactoferricin and lactoferrampin, that have even greater antimicrobial activity than intact LF (). Thus, LF activity toward fungi is likely to be complex and multifactorial: secreted fungal enzymes may digest LF to a greater or lesser extent, and whole LF along with lactoferricin, lactoferrampin and potentially other peptides may interact with variable receptors on the fungal cell membrane and disrupt a range of intracellular processes. Any of these processes may differ between Saccharomyces and Cryptococcus and account for their opposing responses to AMB-LF synergy.
The ‘omics approach used here provides a whole-of-organism view of the response to the drug combination (), and clustering transcripts by SOMs identified co-expressed genes that enable identification of stress-response pathways (). However, the transcriptome is an indirect assessment of cellular changes, and validation studies are required to determine if the assignment of biological terms based on differential gene expression translates to actual cellular processes. With the pressing need to find new ways to treat fungal infections, understanding how synergy works will be valuable for both augmenting existing therapies and for developing novel, effective antifungal strategies.
Conclusion
We conclude that a phenotypically similar outcome to drug synergy is mediated by distinctly different pathways in Saccharomyces and Cryptococcus. The current study indicates that the combined assaults of AMB and LF produce overwhelming stress to the C. neoformans cell due to a disruption of protein and lipid biosynthesis that results in the up-regulation of ERAD, trafficking from the ER to the Golgi to the vacuole, with autophagy and apoptosis. The current study suggests LF or LF-derived peptides have potential to augment AMB and lower the required dose, thereby both improving efficacy and reducing toxicity.
Statements
Data availability statement
The datasets generated for this study can be found in the NCBI Gene Expression Omnibus (GEO) accession number GSE130375.
Author contributions
DC, MW, and SC conceived the study, obtained grant funding, oversaw the work, and assisted with the data analysis. Y-WL undertook all laboratory work and analysis with supervision from DC, LC, and SC. CP undertook all bioinformatic analyses with supervision from MW. Y-WL, CP, and DC drafted the manuscript. MW, SC, and LC commented and edited the drafts. All authors read and approved the submitted version of the manuscript.
Funding
This work was financially supported by the National Health and Medical Research Council of Australia, grant # APP1021267.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02195/full#supplementary-material
References
1
AbassiM.BoulwareD. R.RheinJ. (2015). Cryptococcal meningitis: diagnosis and management update.Curr. Trop. Med. Rep.290–99. 10.1007/s40475-015-0046-y
2
Acosta-ZaldívarM.AndrésM. T.RegoA.PereiraC. S.FierroJ. F.Côrte-RealM. (2016). Human lactoferrin triggers a mitochondrial- and caspase-dependent regulated cell death in Saccharomyces cerevisiae.Apoptosis21163–173. 10.1007/s10495-015-1199-9
3
AgarwalA. K.RogersP. D.BaersonS. R.JacobM. R.BarkerK. S.ClearyJ. D.et al (2003). Genome-wide expression profiling of the response to polyene, pyrimidine, azole, and echinocandin antifungal agents in Saccharomyces cerevisiae.J. Biol. Chem.27834998–35015. 10.1074/jbc.m306291200
4
BasenkoE. Y.PulmanJ. A.ShanmugasundramA.HarbO. S.CrouchK.StarnsD.et al (2018). FungiDB: an integrated bioinformatic resource for fungi and oomycetes.J. Fungi4:E39. 10.3390/jof4010039
5
BenjaminiY.HochbergY. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing.J. R. Stat. Soc. Series B Methodol.57289–300. 10.1111/j.2517-6161.1995.tb02031.x
6
BinnsD.DimmerE.HuntleyR.BarrellD.O’DonovanC.ApweilerR. (2009). QuickGO: a web-based tool for gene ontology searching.Bioinformatics253045–3046. 10.1093/bioinformatics/btp536
7
BrownG. D.DenningD. W.GowN. A. R.LevitzS. M.NeteaM. G.WhiteT. C. (2012). Hidden killers: human fungal infections.Sci. Transl. Med.4:165rv13. 10.1126/scitranslmed.3004404
8
BuchfinkB.XieC.HusonD. H. (2015). Fast and sensitive protein alignment using diamond.Nat. Methods1259–60. 10.1038/nmeth.3176
9
DayJ. N.ChauT. T. H.WolbersM.MaiP. P.DungN. T.MaiN. H.et al (2013). Combination antifungal therapy for cryptococcal meningitis.N. Engl. J. Med.3681291–1302.
10
Drago-SerranoM.Campos-RodríguezR.CarreroJ.De La GarzaM. (2017). Lactoferrin: balancing ups and downs of inflammation due to microbial infections.Int. J. Mol. Sci.18:E501. 10.3390/ijms18030501
11
EdgarR.DomrachevM.LashA. E. (2002). Gene expression omnibus: NCBI gene expression and hybridization array data repository.Nucleic Acids Res.30207–210. 10.1093/nar/30.1.207
12
FalconS.GentlemanR. (2007). Using GOstats to test gene lists for GO term association.Bioinformatics23257–258. 10.1093/bioinformatics/btl567
13
FernandesK. E.CarterD. A. (2017). The antifungal activity of lactoferrin and its derived peptides: mechanisms of action and synergy with drugs against fungal pathogens.Front. Microbiol.8:2. 10.3389/fmicb.2017.00002
14
FerreiraG. F.Baltazar LdeM.SantosJ. R.MonteiroA. S.FragaL. A.Resende-StoianoffM. A.et al (2013). The role of oxidative and nitrosative bursts caused by azoles and amphotericin B against the fungal pathogen Cryptococcus gattii.J. Antimicrob. Chemother.681801–1811. 10.1093/jac/dkt114
15
FischerS.BrunkB. P.ChenF.GaoX.HarbO. S.IodiceJ. B.et al (2011). Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups.Curr. Protoc. Bioinformatics121–19.
16
FronteraL. S.MoyanoS.QuassolloG.Lanfredi-RangelA.RópoloA. S.TouzM. C. (2018). Lactoferrin and lactoferricin endocytosis halt Giardia cell growth and prevent infective cyst production.Sci. Rep.8:18020. 10.1038/s41598-018-36563-1
17
GotzS.Garcia-GomezmJ. M.TerolJ.WilliamsT. D.NagarajS. H.NuedaM. J.et al (2008). High-throughput functional annotation and data mining with the Blast2GO suite.Nucleic Acids Res.363420–3435. 10.1093/nar/gkn176
18
HudsonN. J.DalrympleB. P.ReverterA. (2012). Beyond differential expression: the quest for causal mutations and effector molecules.BMC Genomics13:356. 10.1186/1471-2164-13-356
19
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements.Nat. Methods12357–360. 10.1038/nmeth.3317
20
LaiY. W.CampbellL. T.WilkinsM. R.PangC. N.ChenS.CarterD. A. (2016). Synergy and antagonism between iron chelators and antifungal drugs in Cryptococcus.Int. J. Antimicrob. Agents48388–394. 10.1016/j.ijantimicag.2016.06.012
21
LexA.GehlenborgN.StrobeltH.VuillemotR.PfisterH. (2014). UpSet: visualization of intersecting sets.IEEE Trans. Vis. Comput. Graph.201983–1992. 10.1109/TVCG.2014.2346248
22
LiaoY.SmythG. K.ShiW. (2014). Featurecounts: an efficient general purpose program for assigning sequence reads to genomic features.Bioinformatics30923–930. 10.1093/bioinformatics/btt656
23
LiuT. T.LeeR. E.BarkerK. S.LeeR. E.WeiL.HomayouniR.et al (2005). Genome-wide expression profiling of the response to azole, polyene, echinocandin, and pyrimidine antifungal agents in Candida albicans.Antimicrob. Agents Chemother.492226–2236. 10.1128/aac.49.6.2226-2236.2005
24
MackS. G.TurnerR. L.DwyerD. J. (2018). Achieving a predictive understanding of antimicrobial stress physiology through systems biology.Trends Microbiol.26296–312. 10.1016/j.tim.2018.02.004
25
Mesa-ArangoA. C.ScorzoniL.ZaragozaO. (2012). It only takes one to do many jobs: amphotericin B as antifungal and immunomodulatory drug.Front. Microbiol.3:286. 10.3389/fmicb.2012.00286
26
PangC. N.LaiY. W.CampbellL. T.ChenS. C.CarterD. A.WilkinsM. R. (2017). Transcriptome and network analyses in Saccharomyces cerevisiae reveal that amphotericin B and lactoferrin synergy disrupt metal homeostasis and stress response.Sci. Rep.7:40232. 10.1038/srep40232
27
PerfectJ. R.DismukesW. E.DromerF.GoldmanD. L.GraybillJ. R.HamillR. J.et al (2010). Clinical practice guidelines for the management of cryptococcal disease: 2010 update by the Infectious diseases society of America.Clin. Infect. Dis.50291–322. 10.1086/649858
28
ProtterD. S. W.ParkerR. (2016). Principles and properties of stress granules.Trends Cell Biol.26668–679. 10.1016/j.tcb.2016.05.004
29
R Core Team (2019). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
30
RissoD.NgaiJ.SpeedT. P.DudoitS. (2014). Normalization of RNA-seq data using factor analysis of control genes or samples.Nat. Biotechnol.32896–902. 10.1038/nbt.2931
31
Sangalli-LeiteF.ScorzoniL.Mesa-ArangoA. C.CasasC.HerreroE.GianinniM. J.et al (2011). Amphotericin B mediates killing in Cryptococcus neoformans through the induction of a strong oxidative burst.Microbes Infect.13457–467. 10.1016/j.micinf.2011.01.015
32
Santos-GandelmanJ.Machado-SilvaA. (2019). Drug development for cryptococcosis treatment: what can patents tell us?Mem. Inst. Oswaldo Cruz114:e180391. 10.1590/0074-02760180391
33
SinhaM.KaushikS.KaurP.SharmaS.SinghT. P. (2013). Antimicrobial lactoferrin peptides: the hidden players in the protective function of amultifunctional protein.Int. J. Pept.20131–12. 10.1155/2013/390230
34
StolzA.WolfD. H. (2010). Endoplasmic reticulum associated protein degradation: a chaperone assisted journey to hell.Biochim. Biophys. Acta Mol. Cell Res.1803694–705. 10.1016/j.bbamcr.2010.02.005
35
SuzukiY. A.LopezV.LönnerdalB. (2005). Mammalian lactoferrin receptors: structure and function.Cell. Mol. Life Sci.622560–2575. 10.1007/s00018-005-5371-1
36
UniProt Consortium (2019). UniProt: a worldwide hub of protein knowledge.Nucleic Acids Res.47D506–D515. 10.1093/nar/gky1049
37
UpadhyaR.CampbellL. T.DonlinM. J.AuroraR.LodgeJ. K. (2013). Global transcriptome profile of Cryptococcus neoformans during exposure to hydrogen peroxide induced oxidative stress.PLoS One8:e55110. 10.1371/journal.pone.0055110
38
VenkateshM. P.RongL. (2008). Human recombinant lactoferrin acts synergistically with antimicrobials commonly used in neonatal practice against coagulase-negative staphylococci and Candida albicans causing neonatal sepsis.J. Med. Microbiol.571113–1121. 10.1099/jmm.0.2008/001263-0
39
WehrensR.KruisselbrinkJ. (2018). Flexible Self-Organizing Maps in kohonen 3.0.J. Stat. Softw.871–18.
40
YoshiokaM.YamadaK.YamaguchiY.OgitaA.FujitaK.-I.TanakaT. (2016). The fungicidal activity of amphotericin B requires autophagy-dependent targeting to the vacuole under a nutrient-starved condition in Saccharomyces cerevisiae.Microbiology162848–854. 10.1099/mic.0.000269
41
ZaremberK. A.CruzA. R.HuangC.-Y.GallinJ. I. (2009). Antifungal activities of natural and synthetic iron chelators alone and in combination with azole and polyene antibiotics against Aspergillus fumigatus.Antimicrob. Agents Chemother.532654–2656. 10.1128/AAC.01547-08
42
ZhangL.ZhangY.ZhouY.AnS.ZhouY.ChengJ. (2002). Response of gene expression in Saccharomyces cerevisiae to amphotericin B and nystatin measured by microarrays.J. Antimicrob. Chemother.49905–915. 10.1093/jac/dkf001
Summary
Keywords
drug synergy, RNA-Seq, Cryptococcus neoformans, lactoferrin, amphotericin B, cellular stress response, Saccharomyces
Citation
Lai Y-W, Pang CNI, Campbell LT, Chen SCA, Wilkins MR and Carter DA (2019) Different Pathways Mediate Amphotericin-Lactoferrin Drug Synergy in Cryptococcus and Saccharomyces. Front. Microbiol. 10:2195. doi: 10.3389/fmicb.2019.02195
Received
14 June 2019
Accepted
06 September 2019
Published
01 October 2019
Volume
10 - 2019
Edited by
Teun Boekhout, Westerdijk Fungal Biodiversity Institute, Netherlands
Reviewed by
Won Hee Jung, Chung-Ang University, South Korea; Marilene Henning Vainstein, Federal University of Rio Grande do Sul, Brazil
Updates
Copyright
© 2019 Lai, Pang, Campbell, Chen, Wilkins and Carter.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dee A. Carter, dee.carter@sydney.edu.au
†These authors have contributed equally to this work
This article was submitted to Fungi and Their Interactions, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.