Histoplasma yeast and mycelial transcriptomes reveal pathogenic-phase and lineage-specific gene expression profiles
- Jessica A Edwards†1, 2,
- Chenxi Chen†3,
- Megan M Kemski1, 2,
- Jinnan Hu3,
- Thomas K Mitchell3 and
- Chad A Rappleye1, 2Email author
© Edwards et al.; licensee BioMed Central Ltd. 2013
Received: 23 August 2013
Accepted: 7 October 2013
Published: 10 October 2013
The dimorphic fungus Histoplasma capsulatum causes respiratory and systemic disease in mammalian hosts by expression of factors that enable survival within phagocytic cells of the immune system. Histoplasma’s dimorphism is distinguished by growth either as avirulent mycelia or as pathogenic yeast. Geographically distinct strains of Histoplasma differ in their relative virulence in mammalian hosts and in production of and requirement for specific virulence factors. The close similarity in the genome sequences of these diverse strains suggests that phenotypic variations result from differences in gene expression rather than gene content. To provide insight into how the transcriptional program translates into morphological variation and the pathogenic lifestyle, we compared the transcriptional profile of the pathogenic yeast phase and the non-pathogenic mycelial phase of two clinical isolates of Histoplasma.
To overcome inaccuracies in ab initio genome annotation of the Histoplasma genome, we used RNA-seq methodology to generate gene structure models based on experimental evidence. Quantitative analyses of the sequencing reads revealed 6% to 9% of genes are differentially regulated between the two phases. RNA-seq-based mRNA quantitation was strongly correlated with gene expression levels determined by quantitative RT-PCR. Comparison of the yeast-phase transcriptomes between strains showed 7.6% of all genes have lineage-specific expression differences including genes contributing, or potentially related, to pathogenesis. GFP-transcriptional fusions and their introduction into both strain backgrounds revealed that the difference in transcriptional activity of individual genes reflects both variations in the cis- and trans-acting factors between Histoplasma strains.
Comparison of the yeast and mycelial transcriptomes highlights genes encoding virulence factors as well as those involved in protein glycosylation, alternative metabolism, lipid remodeling, and cell wall glycanases that may contribute to Histoplasma pathogenesis. These studies lay an essential foundation for understanding how gene expression variations contribute to the strain- and phase-specific virulence differences of Histoplasma.
Pulmonary infections with the dimorphic fungal pathogen Histoplasma capsulatum constitute one of the most common respiratory mycoses, affecting both immunocompromised as well as immunocompetent individuals [1, 2]. Histoplasma is acquired by inhalation of infectious conidia, which are produced by the environmental mycelial form of the fungus. The dimorphism of Histoplasma is evident within the mammalian host where the elevated temperature causes differentiation of the conidia into the virulent yeast form. The yeasts infect, survive, and replicate within lung alveolar macrophages. The extent of histoplasmosis disease results from the net contributions of initial inoculum size, the inherent virulence of the strain, and the immune status of the host .
The differentiation of Histoplasma into yeasts and expression of the yeast-phase transcriptional program are necessary for virulence. Histoplasma cells genetically or chemically prevented from transitioning into yeast are avirulent, highlighting the essentiality of differentiation to the pathogenic phase [4–7]. However, it is most likely the expression of yeast-phase-specific genes, rather than the morphology itself, that contributes to Histoplasma virulence. Accordingly, most Histoplasma virulence factors identified to date are restricted to pathogenic-phase yeast cells [8–12].
The H. capsulatum species is constituted of several geographically and phylogenetically distinct groups. Two clades, the North American clade 2 (NAm2; e.g., clinical isolate G217B) and the Panamanian clade (Pan; e.g., clinical isolate G186A) [13, 14], typify the diversity among Histoplasma strains both at the genomic and phenotypic levels. Both strains are virulent, although in murine models of histoplasmosis, G217B infection results in higher organ fungal burdens and increased lethality compared to G186A [15–17]. The G217B genome is roughly 30% larger than the G186A genome (41.0 megabases vs. 30.4 megabases, respectively). Most of the excess DNA in G217B is located in intergenic, repetitive DNA. Both genomes are predicted by in silico analyses of the genome sequence to encode between 9,000 and 10,000 genes (http://www.genome.wustl.edu; http://www.broadinstitute.org). The only differences in gene content determined to date, are the high-affinity iron transport genes, FET3 and FTR1, which are found in the G186A genome but not the G217B genome .
The close similarity in gene content and in coding sequences among Histoplasma strains suggests that phenotypic differences likely result from differences in gene expression rather than variations in gene content. Two examples of known virulence factors clearly illustrate this. First, most Histoplasma lineages, including G186A, have cell walls containing α-glucan and rely on production of this polysaccharide for disease establishment . In the G186A strain, deletion of the α-glucan synthase gene AGS1 causes attenuation due to exposure of immunostimulatory cell wall β-glucans that are normally masked by α-glucan . In contrast to G186A, the cell walls of many North American isolates, represented by G217B, lack α-glucan. Although G186A and G217B have nearly identical α-glucan synthase (AGS1) coding sequences, production α-glucanis absent in G217B during yeast-phase growth. In G217B, the AGS1 promoter is interrupted by a 2.7-kb insertion of repetitive DNA elements that alters expression levels of the synthase. Despite this, G217B remains virulent, suggesting that this strain uses an alternative mechanism to circumvent the need for α-glucan . The second example of expression-based phenotypic differences between Histoplasma strains is the YPS3 gene, which encodes a yeast phase specific factor related to the Blastomyces Bad1 protein [20–22]. The genomes of both G186A and G217B contain the YPS3 gene, but only G217B yeasts produce the Yps3 protein, which contributes to G217B virulence [23, 24]. Similar to α-glucan, differential production of Yps3 appears to result from transcriptional regulation since placing the YPS3 gene under control of an ectopic promoter in the G186A background is sufficient to restore Yps3 protein production .
To better understand how gene expression differences between Histoplasma strains translate into phenotypic differences including yeast virulence, a more complete examination of gene expression profiles is needed. Past technologies for defining the identity and quantity of all transcripts expressed by an organism have included both hybridization-based (e.g., microarrays) and sequence-based (e.g., Sanger sequencing of cDNA or EST libraries) approaches . Hybridization-based approaches have been used to study differences in Histoplasma expression between non-pathogenic (mycelial) phase and pathogenic (yeast) phases [26, 27] as well as during nitrosative stress . Both of these studies were limited to intra-strain expression differences. Recently, a microarray-based analysis of two strains documenting mycelial, yeast, and conidial gene expression was determined . The analytical power of microarrays, however, is limited since they are highly dependent on the accuracy of the predicted gene sets. Inaccurate gene structures derived from ab initio predictions misses or mispredicts genes resulting in errors in the annotations and subsequent microarray data. In contrast, next-generation sequencing-based transcriptome determination defines genes directly from experimentally derived mRNA sequence evidence. Furthermore, there is no upper limit to the expression level with the number of mapped reads being highly correlated to actual gene expression level .
In this study, we use next-generation sequencing (i.e., RNA-seq) of Histoplasma G186A and G217B yeast and mycelial mRNAs to profile the respective pathogenic and non-pathogenic-phase transcriptomes and to identify interstrain pathogenic-phase expression differences that may contribute to variations in virulence. Furthermore, we show for a subset of differentially expressed genes, that differences in expression result from both cis- and trans-acting factors that affect promoter activity in the different genetic backgrounds. These findings will improve our understanding of the mechanisms underpinning morphological, biochemical, and virulence differences among strains of Histoplasma capsulatum.
Determination of the G186A and G217B Histoplasma transcriptomes using RNA-seq
To provide an experimental-evidence-based annotation of the Histoplasma genome, we used RNA-seq methodology to construct gene models with transcriptional support. For a more comprehensive gene definition, we sequenced the mRNAs from pathogenic-phase (yeast) and non-pathogenic-phase (mycelia) Histoplasma cells, the two distinguishing lifestyles of this dimorphic fungal pathogen. In addition, RNAs from two strains of Histoplasma (G186A and G217B), representing the Pan and NAm2 clades, respectively, were analyzed. These two strains are clinical isolates of Histoplasma, have complete genome sequences (http://www.genome.wustl.edu; http://www.broadinstitute.org), and are the two strains for which molecular genetic methodologies have been established [8, 9, 15, 30–34]. For yeast-phase RNA samples, strains were grown to late-exponential phase (approximately 72 hours) in conditions approximating the mammalian host environment (growth at 37°C in 5% CO2 / 95% air). The growth of two cultures for biological replicate RNA samples was nearly identical as measured by yeast culture turbidity (data not shown). For mycelia RNA samples, strains were grown in liquid medium at 25°C in normal air without shaking until sufficient mycelial biomass formed (approximately 3 weeks). The RNA integrity number (RIN)  for all RNA samples used was ≥ 8.5. Next generation sequencing (Illumina RNA-seq) of the eight mRNA libraries (two biological replicates for each condition) yielded approximately 118 million paired-end reads for derivation of the Histoplasma set of expressed genes (38.6 million and 79.5 million reads for yeasts and mycelia, respectively).
Experimental-evidence-based annotation of the Histoplasma genome
Due to the high similarity of the genomes of G186A and G217B and the increased complexity of the G217B genome from the large amount of repetitive DNA, we used G186A gene structures to inform construction of the G217B genes. Short reads in the four G217B mRNA libraries (2 yeast-phase and 2 mycelial-phase) were aligned to the G186A reference genome using Tophat with relaxed parameters to account for the nucleotide variation between strains. We allowed for 6 mismatches in the 75-bp G217B short reads, providing at least 92% sequence identity, which is similar to the identity of known orthologous genes between strains (identity ranges from 93% to 99%). Using these parameters, 72% of the processed reads from each G217B library were matched to the G186A reference genome with high confidence and these were used to derive the base G217B gene set. As some genes unique to the G217B genome or only expressed by G217B cells would be missed, the reference-based alignment and de novo transcript assembly pipeline (Figure 1) was then applied to the remaining high quality G217B reads that were not directly matched to the G186A transcriptome. This identified an additional 62 G217B genes, which were added to the G217B models to derive the final G217B gene set of 9004 genes.
G186A gene statistics
Middle 90% range
Gene length (bp)
333 - 4857
Exons per gene
1 - 7
Exon length (bp)
Introns per gene
1 - 6
Intron length (bp)
54 - 237
Gene density (per kb)
RNA-seq improved gene annotation for G186A
To determine the false positive rate, we also compared the exon structure accuracy by calculating the exon coverage rate by mRNA short reads (Figure 3B). An exon with 99% or more coverage was defined as perfect support, 70-99% coverage as fair support, and less than 70% coverage as poor support. 95.2% of exons defined by our optimized gene model pipeline show perfect support by mRNA short reads, compared to 89.8% of exons defined by in the ab initio predictions. Only 0.8% of exons in our gene models had poor support compared to 2.8% of exons in the in ab initio gene prediction set. Thus, the gene models derived from our optimized pipeline are more supported by experimental evidence and thus have more accurate structures.
As further validation of the gene structures defined by RNA-seq, we sequenced the mRNAs for some genes with discrepant structures between the RNA-seq-derived gene set and the ab initio predictions (Figure 3C-E). For each selected gene, a cDNA spanning all exons was generated by RT-PCR and the amplicon was sequenced to provide nucleotide-level validation of the gene structure. We resolved differences in gene structure predictions for the MFS5 gene, which encodes a predicted major facilitator superfamily membrane transporter and two hypothetical genes (i.e., genes without recognizable functions), HYP12 and HYP13. The sequence of the MFS5 cDNA shows the MFS5 gene is composed of 6 exons, which completely matches the RNA-seq-derived structure (Figure 3C). The ab inito MFS5 gene prediction has 5 exons, inaccurately missing the intron between exons 2 and 3. For the HYP12 gene, 5 exons were correctly defined by RNA-seq, but the ab initio predictions missed exon 2, instead including an abnormally large 451 base-pair intron (Figure 3D). For the HYP13 gene, RNA-seq accurately defined the gene with 2 exons whereas ab initio mispredicted the location of the second exon and added a third exon further downstream (Figure 3E). Together, these data demonstrate the greater experimental support and the improved accuracy in gene structure definition from the RNA-seq based transcriptome compared to the ab initio predicted gene models.
Quantitative gene expression profiling
The transcription profile and relative gene expression levels for genes expressed by G186A and G217B Histoplasma cells were determined by counting the number of matching RNA-seq reads from each strain. Relative expression levels were calculated using the Cufflinks algorithm . To enable cross-species and cross-phase comparisons, gene expression levels were normalized and calculated as Fragments Per Kilobase of exon per Million fragments mapped (FPKM; ). Analysis of the FPKM values for the biological replicates of each phase for each Histoplasma strain shows the replicates are highly similar (Additional file 1: Figure S1). Thus, the mean FPKM between replicate libraries was used for calculation of the fold change in gene expression levels.
G186A and G217B gene expression statistics
0 – 13514
0 – 15819
0 – 6727
0 – 17942
Genes comprising 50% of total transcripts
Examination of genes whose expression is upregulated in the yeast phase highlights features potentially linked to the virulence that characterizes this phase. The set of upregulated yeast-phase genes includes CBP1, SOD3, CATB, CFP4, CFP8, ENG1, TRL1, all well established yeast-phase expressed genes in both strains [31, 45], many of which have been linked to virulence [9, 10, 12]. Consistently, CBP1 is one of the most highly expressed genes (FPKM values at least 8000) and one of the most differentially expressed (at least 1000-fold increased expression in yeast compared to mycelia). The virulence-promoting AGS1 and YPS3 genes are also enriched in yeast, but only in G186A and G217B, respectively.
The list of genes with significant yeast-phase enriched expression (at least 5-fold compared to mycelia) is presented in Additional file 4: Table S3. In total 275 genes were upregulated in yeasts with 43 genes upregulated in both strains (100 were upregulated only in G186A yeasts compared to G186A mycelia and 132 were upregulated only in G217B yeasts compared to G217B mycelia). In G186A, the most differentially expressed gene is HC186_02213 (1735-fold induced in yeast), a gene of unknown function. In G217B, the gene most upregulated in yeast is CBP1. In both strains, the GNT1 gene, which encodes an N-acetylglucosaminyltransferase involved in N-linked protein glycosylation, is consistently highly upregulated (53-fold and 78-fold in G186A and G217B, respectively). G217B yeasts also upregulate an α-mannosyltransferase (MNN2; 104-fold).
Genes upregulated in yeasts compared to mycelia suggests different metabolism between the phases. Genes encoding enzymes involved in coenzyme A (CoA) synthesis (2-dehydropantoate 2-reductase (PAN5) and pantetheine-phosphate adenylyltransferase (CAB4)) are also significantly upregulated in yeasts despite pantothenate in the growth medium. Various membrane transporters are also more highly expressed in yeasts compared to mycelia including MFS- and ABC-type transporters and proteins putatively transporting metabolites such as amino acids (3 transporters including the general amino acid permease GAP1), zinc (a zinc transporter in G186A yeasts and ZRT2 in G217B yeasts), sulfate, and phosphate. Yeasts of both strains also upregulate the ATG1 kinase.
Yeasts are also characterized by increased expression of various transcriptional control proteins, which may underlie the differing expression profiles of yeasts and mycelia. In G186A and G217B, a subunit of TFIIE (encoded by an ortholog of TFA2) is upregulated 32-fold and 73-fold, respectively. G186A yeasts and G217B yeasts compared to their respective mycelia also have higher expression of Zn-finger transcription factors (two in G186A and one in G217B). Both strains upregulate expression of the gene encoding a subunit of DNA Polymerase II, which may indicate a greater DNA synthesis capacity, is required for faster cell cycles in yeast compared to mycelia.
Although a few inferences can be made from the yeast-phase regulon, the vast majority of yeast-phase regulated genes (80% of G186A and 73% of G217B) encode hypothetical proteins with no ortholog with known biochemical or molecular function. This highlights how little is currently understood about the gene expression profiles that provide for a pathogenic lifestyle compared to non-pathogenic growth.
Identification of strain-specific, pathogenic-phase gene expression profiles
The complete list of genes unique to G186A or with enhanced expression in G186A is detailed in Additional file 5: Table S4. Genes structurally unique to G186A or only expressed by G186A yeasts comprise most of the enriched kinase/phosphatase functions but also include FET3, FTR1, an ARN-family siderophore transporter (ARN2), and a predicted quinone oxidoreductase (ZTA1). Among the genes with enriched expression in G186A yeasts compared to G217B are those encoding a glycosyl transferase (ATG26) and a chitinase (CHI1), two phospholipases (PLB1 and PLD1), a glyoxylate reductase (GOR1), and an N-acetyltransferase (ARD1).
In contrast, G217B yeast has unique or enriched expression of a different set of genes that are detailed in Additional file 6: Table S5. Expressed genes unique to the G217B genome, and for which clear orthologous proteins could be identified, encode a putative beta-glucanase (TOS1) and two kinases potentially involved in cellular signaling (PKH2 and GSK3). A v-SNARE protein is also uniquely encoded in the G217B genome by the YKT6 gene, which may indicate differential trafficking of vesicles in this strain. Genes of the siderophore biosynthesis gene cluster (i.e., SID3, SID4, OXR1, MFS1, and NPS1);  are preferentially expressed by G217B yeasts compared to G186A yeasts even though both genomes have this cluster. Interestingly, genes encoding histones (H1, H2A, H2B, H3, and H4) all show increased expression in G217B yeasts. Genes encoding a copper transporter (CTR3) and one of three carbonic anhydrases (CAH1) are 68-fold and 30-fold more expressed, respectively, in G217B than G186A. The genes for a ras-GTPase activating protein (IQG1), a MAP kinase kinase kinase (SSK2), a Ca2+/calcineurin-dependent transcription factor (CRZ1), and a histone deactylase (HOS1) show increased expression in G217B. G217B-enriched genes also include factors linked to Histoplasma virulence, namely Yps3 (nearly 300-fold), but also the secreted superoxide dismutase (Sod3; ) and the secreted catalase (CatB; ).
Quantitative RT-PCR was used to confirm the expression differences between strains for a subset of 41 differentially expressed and similarly expressed genes. Both genes encoding proteins with recognizable functions as well as genes without known functions (HYP genes) were included. In addition the CBP1, SOD3, CATB, and YPS3 genes were included to determine if expression differences in these known virulence factors [9, 10, 12, 23] correlated with differences in strain virulence. The relative fold-change in expression between strains determined by FPKM analysis was highly correlated with differences determined by qRT-PCR (R2 = 0.88; Figure 7B) providing validation of the FPKM analysis. Genes significantly more expressed by G186A included HYP18 (234.5-fold), ENV9 (68.0-fold), a putative alkaline phosphatase (PHO1; 43.9-fold), MFS5 (39.2-fold), HYP8 (11.0-fold), TRL1 (10.6-fold), HYP7 (9.8-fold), SNF3 (9.6-fold), HYP17 (9.4-fold), ACS1 (8.2-fold), HYP3 (6.9-fold), and KIN2 (4.9-fold). Genes significantly more expressed by G217B included YPS3 (1019.3-fold), HYP16 (58.4-fold), HYP15 (34.4-fold), HYP1 (23.4-fold), SID4 (17.4-fold), H2B (16.0-fold), CAH1 (15.7-fold), SOD3 (15.7-fold), CTR3 (13.9-fold), OXR1 (12.3-fold), HYP10 (7.8-fold), and CATB (6-fold) (Figure 7C). Regarding known virulence factors, YPS3, SOD3, and CATB were more highly expressed by G217B while CBP1 was slightly more expressed in G186A.
Factors determining inter-strain variation in expression of genes
For genes more highly expressed in G217B, we tested whether the promoter sequence was responsible for the variation in gene expression. Three genes more highly expressed in G217B were selected for investigation: CTR3, SOD3, and YPS3. For CTR3, 1.45 kb of upstream sequence was sufficient to cause expression of the gfp reporter (Figure 8B). The CTR3 promoter nucleotide sequences are 92% identical between G217B and G186A. GFP reporter fluorescence was significantly higher in the G217B background (2.1 to 2.9-fold) regardless of the strain from which the promoter sequence originated. This indicates that the enhanced expression of CTR3 in G217B is independent of polymorphisms in the promoter but highly dependent on the genetic background (i.e., trans-acting factors that differ between the strains). In contrast to CTR3, cis-acting factors (i.e., promoter sequences) are the major determinants of interstrain expression differences for SOD3 and YPS3. The G217B sequence of the SOD3 promoter (1.9 kb) provided higher transcription of the GFP reporter than the G186A SOD3 promoter in both strain backgrounds (1.8-fold and 1.9-fold; Figure 8C). For regulation of YPS3 expression, the G186A sequence of the YPS3 promoter (1.9 kb of upstream sequence) could not drive transcription of the GFP reporter in either strain background whereas the G217B sequence of the promoter enabled transcription of the reporter gene (Figure 8D). Thus, transcription of YPS3 is controlled primarily by the promoter sequence, although transcription of the reporter was 2.4-fold higher in the G217B background than the G186A background suggesting the contribution of some additional trans-acting factors operating in G217B.
To determine whether cis- or trans-acting factors determined expression differences of genes more highly expressed in G186A, promoter fusions were created using sequences upstream of the AGS1 (1.9 kb), MFS5 (1.9 kb), and ENV9 (1.9 kb) genes. Previously we showed that AGS1 transcription was significantly attenuated in G217B due to insertional disruption of the promoter  indicating cis-control of AGS1. Removal of the disrupting sequences from the G217B AGS1 promoter restored transcriptional activity in both G186A and G217B backgrounds indicating deficient AGS1 expression in G217B is primarily due to the disruption in the promoter (Figure 8E). For MFS5 promoter-gfp fusions, transcription of the reporter gene was strongly influenced by the strain background with both G186A and G217B MFS5 promoter sequences driving 2.0-fold to 3.0-fold more transcription in the G186A background (Figure 8F). Thus, MFS5 transcriptional differences arise predominantly from differences in the strain background rather than promoter polymorphisms. Transcriptional gfp fusions to the ENV9 promoter (up to 1.9 kb) failed to produce any GFP fluorescence in transformants irrespective of the genetic background or which promoter sequence was used. We suspect that the lack of GFP protein production results from unknown alternative translational start sites, which could shift gfp translation out of frame. As an alternative means of testing for transcriptional activity of the different ENV9 promoters without requiring reporter protein synthesis, we used quantitative RT-PCR of the gfp transgene. For ENV9, both the promoter and background influenced the expression of the transgene with increased expression driven by the G186A promoter. This indicates that ENV9 expression is strongly influenced by promoter sequence, but the strongest expression was in the G186A background indicating trans-acting factors also contribute to ENV9 expression differences (Figure 8G). Together, these data provide examples that show gene expression differences among Histoplasma strains that are determined by cis acting factors (i.e., promoter sequence) and/or variations in trans acting factors (e.g., translation factor abundance or activity) that exist between strain backgrounds.
While genome sequencing provides an excellent starting point for characterization and analysis of Histoplasma capsulatum’s nearly 10,000 genes, accurate gene structure determination requires experimental evidence. Use of tiling microarrays with isolated RNAs improved gene structure definition in Histoplasma. In this study, we used RNA-seq to inform gene models for two Histoplasma clinical isolates, G186A and G217B, representing divergent phylogenetic clades. In addition to spliced-alignment of the mRNAs to a references genome, we used de novo assembly of the mRNAs at multiple points of the annotation pipeline to further refine annotations with experimental evidence. The assembled mRNA sequences generated more precise gene models than those derived from ab initio predictions, and improved the evidence-based resolution of the gene structures to the nucleotide level.
In addition, RNA-seq enabled a comparative gene expression approach to identify phase-specific as well as strain-specific gene profiles. As Histoplasma cells exist as yeasts during mammalian infection, we directed our efforts at identification of genes preferentially expressed by these virulent cells compared to mycelia that are unable to establish disease. Overall, we found that 6% to 9% of the genome is differentially expressed between these two phases. Microarray-based studies found 5%  to 19%  of genes had phase-dependent regulation. The dramatically increased percentage of regulated genes in Inglis et al., likely results from the 3-fold differential expression criteria they used, whereas we required a more stringent threshold of 5-fold differential expression. By RNA-seq, we found that growth as yeast cells results in induced expression of 3% to 5% of the genes. Assuming that increased expression confers the characteristics required for pathogenic-phase growth, 300 to 400 genes could potentially contribute to Histoplasma virulence.
Genes preferentially expressed by yeast-phase cells of both strains suggest some functions that characterize pathogenic-phase growth. Yeast cells upregulate gene expression of enzymes for glycosylation of proteins (Gnt1 and an α-mannosyltransferase). These yeast-phase-expressed glycoysltransferases suggest that yeasts, but not mycelia, rely upon glycosylation of extracellular proteins as they transit through the ER and Golgi. In support of this, Holbrook et al. found that proteins from yeast culture filtrates, but not mycelia culture filtrates, are heavily glycosylated . Yeast cells also upregulate enzymes involved in CoA synthesis (e.g., Pan5 and Cab4) that may indicate yeast have additional CoA/acetyl-CoA need for carbon and lipid metabolism pathways. Yeast cells have increased expression of ATG1, a kinase that is involved in regulating autophagy in response to environmental signals [48–50]. As the cDNA libraries were prepared from yeasts and mycelia grown in identical growth medium, the yeast-phase upregulated genes presumably represent gene expressions regulated by lifestyle rather than nutrient availability.
The phylogenetic groups represented by G186A and G217B strains are notably dissimilar in many of the yeast-phase-enriched genes. Only 43 of the 275 genes differentially expressed by yeasts compared to mycelia are common between the two strains. This may not be surprising as G186A and G217B have been shown to have different relative virulence in murine models [15–17] and utilize different virulence factors to mediate their pathogenesis [15, 19, 23]. Gene expression differences between strains indicate that the different Histoplasma lineages are more dissimilar than previously assumed. For example, expression of a high affinity copper transporter (CTR3) and a carbonic anhydrase (CAH1) are substantially increased in G217B while an MFS-family transporter (MFS5) shows 100-fold more expression in G186A. Our data indicate that the differential regulation of genes between strains results from differences in cis-acting sequence polymorphisms in the promoter of genes as well as trans-acting factors that reflect the different background of the strains.
Our expression profiling highlights differences in iron acquisition between G186A and G217B. We found that the SID1 gene, which is involved in siderophore synthesis, is only yeast-phase enriched in G217B contrary to the microarray analysis . Both G186A and G217B express ARN-family siderophore transporters. In G217B, one family member has enriched expression in yeasts compared to mycelia. On the other hand, G186A yeasts express an ARN-family transporter unique to the G186A genome and the G186A-specific iron oxidase and permease encoding genes FET3 and FTR1. This data indicates G217B and G186A differ in the mechanisms for iron acquisition.
The G186A-enriched gene set suggests G186A yeasts differ from G217B yeasts in lipid metabolism. G186A has increased expression of two phospholipases, Plb1 and Pld1, as well as the acyl chain desaturase Ole1, which has been linked to Histoplasma thermotolerance [51, 52]. Together, these factors may suggest that G186A is better able to remodel lipids to adapt to thermal stresses encountered during infection. G186A also preferentially expresses a number of kinases (a SKY1-related kinase and a HOG1-related kinase) that may help it sense and respond to conditions related to infection of phagocytes.
The cell walls of G186A and G217B yeasts are biochemically different from each other. Most notable is the α-glucan difference between strains. Our expression analysis highlights additional glycan-related enzymes that have differential expression between the strains. G186A yeasts express higher levels of the CHI1 gene encoding a chitinase and a glucosyl transferase (ATG26). On the other hand, G217B yeasts have increased expression of the TOS1 gene and a putative exo-glucanase gene (EXG2). The Tos1 factor is a putative β-1,3-glucanase [53, 54] that responds to cell stress . Cell wall modification by Tos1 and/or Exg2 may be essential for modifying the yeast cell wall of G217B, which lacks the virulence-promoting α-glucan of G186A yeasts. Thus, the surface composition differences between these two lineages may extend beyond the α-glucan content and this may directly affect how G186A and G217B yeasts interface with host cells.
Compared to G186A, G217B cells express higher levels of oxidative stress defense genes. G217B yeasts show higher transcription of the extracellular catalase (CATB) and the extracellular superoxide dismutase (SOD3) genes, the products of which are required for Histoplasma virulence [10, 12]. This may indicate G217B yeasts may survive the phagocyte oxidative burst better than G186A yeasts. Together these data suggest G217B may rely on enhanced defense mechanisms rather than evasion of phagocyte detection through modification of the glycan composition of the cell wall.
Although inferences and hypotheses based on differentially expressed genes can be made, the majority of differentially expressed genes encode proteins with no known function, whether comparing pathogenic and non-pathogenic phases or comparing between strains. This underscores how little we understand the biology of Histoplasma and the need for functional studies. The experimental evidence-based transcriptomes established in this study will provide an important framework for identifying genes that underlie the pathogenesis differences between phases of this dimorphic fungal pathogen.
In this study, we used mRNA sequencing (i.e., RNA-seq) to refine gene models for two Histoplasma clinical isolates, G186A and G217B, representing divergent phylogenetic clades with different virulence. Quantitation of the transcriptional profiles identified phase-specific and strain-specific expression differences that correlate with differences in fungal virulence. Depending on the strain, between 6% and 9% of genes are differentially regulated between the virulent yeast and avirulent mycelial phases. As the gene content is equivalent between yeasts and mycelia of a given strain, these findings underscore the fact that Histoplasma pathogenesis is primarily a function of gene expression differences between the dimorphic phases.
Histoplasma strains and growth conditions
The wild-type Histoplasma strains used were the clinical isolates G186A (ATCC 26029) and G217B (ATCC 26032). Uracil auxotroph strains for transformation with URA5-based plasmids were the ura5-deletion strains OSU1 and WU15, which were derived from the G186A and G217B wild-type isolates, respectively [8, 56]. Yeast and mycelial-phase fungal cells were cultured in Histoplasma-macrophage media (HMM) . For growth of uracil auxotrophs, HMM was supplemented with uracil (100 μg/ml). For growth on plates, HMM was solidified with agarose and supplemented with 25 μM FeSO4. To maintain yeast-phase morphology, yeasts were cultured at 37°C with agitation. Mycelia cultures were grown at 25°C.
Yeasts grown to late exponential phase (approximately 72 hours) were collected by centrifugation, resuspended in RNAlater (Ambion), and frozen at -80°C. To collect mycelia, hyphae were separated from the culture media by filtration through Whatman #5 filter paper and the mycelial cells placed in RNAlater and frozen at -80°C. Total RNA was isolated using the RiboPure-Yeast Kit (Ambion) using mechanical disruption of fungal cells and purification of RNA from the lysate on an RNA-binding column. RNA quality was assessed with the Bioanalyzer platform (Agilent). Two biological replicate cultures of yeast and mycelia were prepared for transcriptome libraries.
Library preparation and transcriptome sequencing
Library preparation and sequencing of mRNAs by RNA-seq were performed at the Molecular and Cellular Imaging Center at the Ohio Agriculture Research and Development Center at Ohio State University. cDNA libraries were prepared using the TruSeq RNA sample preparation kit (Illumina). Briefly, mRNA was isolated from total RNA using poly-A capture. The mRNAs were fragmented by cation treatment with heat and then converted to cDNA by reverse transcription and second strand synthesis. cDNA ends were repaired and adenylated to facilitate ligation of indexed adaptors. Following ligation of adaptors, 15 cycles of PCR were performed to enrich for cDNA fragments with adaptors on both ends of the molecule. Indexed libraries were pooled and sequenced using the Illumina GAII platform to generate paired-end reads. Library reads were deconvolved post-sequencing by virtue of the adaptor index on each cDNA molecule.
Gene modeling and annotation
The spliced alignment tool Tophat  was used to align short reads from G186A and G217B mRNAs to the to the G186A reference genome (Broad Institute). Strict parameters were used (0 mismatches and no gaps) for mapping G186A reads to generate a highly accurate alignment result. Alignments were used to indicate the transcript regions and identify splice junctions, which were subsequently used as hints to derive gene structure models with Augustus [37, 38]. Augustus was also used to produce a gene prediction set for additional sensitivity. In parallel, Inchworm  was used for de novo transcript assembly of RNA-seq short reads and previously sequenced EST reads. BestORF (Molquest package, Softberry) was used to identify open reading frames in the de novo transcript assembly. The Tophat/Augustus gene models were refined by the de novo transcript-based gene structure evidence using the PASA algorithm . The three data sets were integrated by using the spliced-alignment model as a base and adding in unique genes (based on genomic locations) from the other two data sets. Manual inspection and refinement of the gene structures was included to divide likely gene fusions (genes with unusually large introns, i.e., introns > 350 base pairs) or to combine potentially split genes (genes separated by unusually small intergenic distances, i.e., < 500 base pairs). Repetitive genes were identified in the final gene set by BLAST search against genome sequences and identification of those genes with two or more matches to the genome with at least 50% coverage or an E-value less than 10-30. To identify splicing signals, 15 bp of both ends of all introns in the G186A gene models were extracted and fed into the MEME motif finder . Functional annotations of genes were assigned to the G186A and G217 gene sets using Blast2GO , KAAS , and BLAST searches of NCBI protein databases to identify homologous genes and/or protein functions. Reciprocal top-hit BLAST was used to assign orthologous identities between strains.
Gene expression analysis of RNA-seq data
For RNA-seq-based expression analysis, the normalized gene expression (FPKM) for each defined gene was calculated by Cufflinks and Cuffdiff . For these analyses, repetitive genes or genes with extremely low expression in all libraries (FPKM less than 0.1) were excluded. For cross-species gene expression comparisons, the G186A gene set was used as the reference gene models for counting reads. G217B mRNA reads were matched to the G186A gene set by allowing for 6 mismatches (which maintains 92% nucleotide sequence identity between strains). G217B mRNA reads that did not match any G186A gene model were extracted and assembled de novo into transcripts using Inchworm  and the genes designated as G217B unique genes. The common genes between strains were then compared using Cuffdiff to identify those genes with significantly different expression (at least five fold, q value < 0.01) between G186A and G217B. To avoid artificially high ratios of expression due to very low expression in one strain, all FPKM values less than 0.5 were set to 0.5 before ratios of expression were calculated. The final set of genes with differential expression was subsequently filtered to remove those genes with low expression (FPKM less than 0.5 in both strains). The differentially expressed gene sets were then combined with the structurally unique genes of each strain to generate the final set of genes with strain-specific expression.
Endpoint and quantitative RT-PCR
Three micrograms of total RNA was reverse transcribed using SuperScript III reverse transcriptase (Invitrogen) and 15-mer Oligo (dT) primers and genomic DNA removed by DNAse treatment. For endpoint RT-PCR, reaction mixes included 0.5 μM gene-specific primers, 0.2 mM dNTPs, and 1:10 dilution of the reverse-transcribed RNA. Reactions using RNA in the absence of reverse transcription were performed to verify the lack of genomic DNA in RNA isolations. For quantitative PCR, reverse-transcribed RNA templates were used at a 1:10 concentrations in a PCR master mix with SYBR green (Bio-Rad) and 0.5 μM each gene-specific primer (Additional file 7: Table S6). PCR products were amplified for 30 cycles at 94°C for 10 seconds, 52°C -55°C for 15 seconds, and 68°C for 30 seconds using a realplex2 thermal cycler (Eppendorf). Cycle threshold (CT values) were calculated with the realplex software (v2.2) using the CalQplex algorithm (Eppendorf). Transcript levels were normalized to the TEF3 gene. Relative fold changes in gene expression between strains were calculated using the ΔΔC T method . For determination of ENV9 promoter activity, RNA was isolated from nine independent transformants containing the P ENV9 -gfp transcriptional reporter fusion and the RNA was reverse transcribed as above. Transcriptional activity of the reporter was determined by qPCR of the gfp gene using gfp-specific primers and the relative gfp mRNA abundance compared to TEF1 after normalization of all C T values to the TEF3 gene .
Analysis of promoter activity
Promoter activity was determined by creation of transcriptional fusions to a gfp reporter gene. Promoters were amplified with Phusion High-Fidelity Polymerase (NEB) and cloned into plasmids pCR623 or pCR628 that contain the gfp reporter. Putative promoter regions encompassed sequences upstream of the start codon from G186A and G217B for the TEF1 (661 bp; pCR640 and pCR639), YPS3 (1873 bp; pMK43 and pMK42), SOD3 (1951 bp; pMK49 and pMK48), CTR3 (1454 bp; pMK33 and pMK32), AGS1 (1916 bp; pCR637 and pCR635), MFS5 (1943 bp; pMK51 and pMK52), and ENV9 (1961 bp; pMK47 and pMK46) genes. Promoter constructs were sequenced and then transformed into Histoplasma OSU1 or WU15 using Agrobacterium tumefaciens-mediated transformation . Ura+ transformants were selected by plating on HMM media with 10 μg/mL tetracycline at 37°C. Individual transformant colonies were picked and spotted onto HMM medium and the GFP fluorescence of individual spots quantified using an AlphaImager UV transillumination system (CellBiosciences: . GFP fluorescence measurements were all normalized to GFP expression levels of transformants with GFP under control of the TEF1 promoter.
Short reads of Histoplasma mRNAs generated by Illumina sequencing have been submitted to the NCBI Sequence Read Archive as accession numbers SRX332607 (G186A yeast), SRX332751 (G186A mycelia), SRX332749 (G217B yeast), and SRX332752 (G217B mycelia). Assembled transcripts and FPKM expression values are available at http://microbiology.osu.edu/RappleyeHistoplasma and the Additional files.
We thank Dr. Tea Meulia and Dr. Asela Wijeratne at the Molecular and Cellular Imaging core at Ohio State University-OARDC for library preparation and sequencing. The Broad Institute (MIT) and the Genome Institute (Washington University) provided publically available G186A and G217B reference genome sequences. This work was supported by research grant AI083335 from the National Institutes of Health and by grant 0865450D from the American Heart Association to Chad A. Rappleye.
- Kauffman CA: Histoplasmosis. Clin Chest Med. 2009, 30 (2): 217-225. 10.1016/j.ccm.2009.02.002. vView ArticlePubMedGoogle Scholar
- Pfaller MA, Diekema DJ: Epidemiology of invasive mycoses in North America. Crit Rev Microbiol. 2010, 36 (1): 1-53. 10.3109/10408410903241444.View ArticlePubMedGoogle Scholar
- Rippon JW: Histoplasmosis (Histoplasmosis casulati). Medical Mycology: the Pathogenic Fungi and the Pathogenic Actinomycetes. 1988, Philadelphia: W. B. Saunders Co, 381-423. 3Google Scholar
- Medoff G, Sacco M, Maresca B, Schlessinger D, Painter A, Kobayashi GS, Carratu L: Irreversible block of the mycelial-to-yeast phase transition of Histoplasma capsulatum. Science. 1986, 231 (4737): 476-479. 10.1126/science.3001938.View ArticlePubMedGoogle Scholar
- Nemecek JC, Wuthrich M, Klein BS: Global control of dimorphism and virulence in fungi. Science. 2006, 312 (5773): 583-588. 10.1126/science.1124105.View ArticlePubMedGoogle Scholar
- Webster RH, Sil A: Conserved factors Ryp2 and Ryp3 control cell morphology and infectious spore formation in the fungal pathogen Histoplasma capsulatum. Proc Natl Acad Sci U S A. 2008, 105 (38): 14573-14578. 10.1073/pnas.0806221105.PubMed CentralView ArticlePubMedGoogle Scholar
- Nguyen VQ, Sil A: Temperature-induced switch to the pathogenic yeast form of Histoplasma capsulatum requires Ryp1, a conserved transcriptional regulator. Proc Natl Acad Sci U S A. 2008, 105 (12): 4880-4885. 10.1073/pnas.0710448105.PubMed CentralView ArticlePubMedGoogle Scholar
- Rappleye CA, Engle JT, Goldman WE: RNA interference in Histoplasma capsulatum demonstrates a role for alpha-(1,3)-glucan in virulence. Mol Microbiol. 2004, 53 (1): 153-165. 10.1111/j.1365-2958.2004.04131.x.View ArticlePubMedGoogle Scholar
- Sebghati TS, Engle JT, Goldman WE: Intracellular parasitism by Histoplasma capsulatum: fungal virulence and calcium dependence. Science. 2000, 290 (5495): 1368-1372. 10.1126/science.290.5495.1368.View ArticlePubMedGoogle Scholar
- Youseff BH, Holbrook ED, Smolnycki KA, Rappleye CA: Extracellular superoxide dismutase protects Histoplasma yeast cells from host-derived oxidative stress. PLoS Pathog. 2012, 8 (5): e1002713-10.1371/journal.ppat.1002713.PubMed CentralView ArticlePubMedGoogle Scholar
- Keath EJ, Painter AA, Kobayashi GS, Medoff G: Variable expression of a yeast-phase-specific gene in Histoplasma capsulatum strains differing in thermotolerance and virulence. Infect Immun. 1989, 57 (5): 1384-1390.PubMed CentralPubMedGoogle Scholar
- Holbrook ED, Smolnycki KA, Youseff BH, Rappleye CA: Redundant Catalases Detoxify Phagocyte Reactive Oxygen and Facilitate Histoplasma capsulatum Pathogenesis. Infect Immun. 2013, 81 (7): 2334-2346. 10.1128/IAI.00173-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Kasuga T, Taylor JW, White TJ: Phylogenetic relationships of varieties and geographical groups of the human pathogenic fungus Histoplasma capsulatum Darling. J Clin Microbiol. 1999, 37 (3): 653-663.PubMed CentralPubMedGoogle Scholar
- Kasuga T, White TJ, Koenig G, McEwen J, Restrepo A, Castaneda E, Da Silva Lacaz C, Heins-Vaccari EM, De Freitas RS, Zancope-Oliveira RM: Phylogeography of the fungal pathogen Histoplasma capsulatum. Mol Ecol. 2003, 12 (12): 3383-3401. 10.1046/j.1365-294X.2003.01995.x.View ArticlePubMedGoogle Scholar
- Edwards JA, Alore EA, Rappleye CA: The Yeast-Phase Virulence Requirement for alpha-Glucan Synthase Differs among Histoplasma capsulatum Chemotypes. Eukaryot Cell. 2011, 10 (1): 87-97. 10.1128/EC.00214-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Mayfield JA, Rine J: The genetic basis of variation in susceptibility to infection with Histoplasma capsulatum in the mouse. Genes Immun. 2007, 8 (6): 468-474. 10.1038/sj.gene.6364411.View ArticlePubMedGoogle Scholar
- Tewari RP, Berkhout FJ: Comparative pathogenicity of albino and brown types of Histoplasma capsulatum for mice. J Infect Dis. 1972, 125 (5): 504-508. 10.1093/infdis/125.5.504.View ArticlePubMedGoogle Scholar
- Hilty J: George Smulian A. 2011, Histoplasma capsulatum utilizes siderophores for intracellular iron acquisition in macrophages. Med Mycol: Newman SLGoogle Scholar
- Rappleye CA, Eissenberg LG, Goldman WE: Histoplasma capsulatum alpha-(1,3)-glucan blocks innate immune recognition by the beta-glucan receptor. Proc Natl Acad Sci USA. 2007, 104 (4): 1366-1370. 10.1073/pnas.0609848104.PubMed CentralView ArticlePubMedGoogle Scholar
- Rooney PJ, Klein BS: Sequence elements necessary for transcriptional activation of BAD1 in the yeast phase of Blastomyces dermatitidis. Eukaryot Cell. 2004, 3 (3): 785-794. 10.1128/EC.3.3.785-794.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Bohse ML, Woods JP: Expression and interstrain variability of the YPS3 gene of Histoplasma capsulatum. Eukaryot Cell. 2007, 6 (4): 609-615. 10.1128/EC.00010-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Bohse ML, Woods JP: Surface localization of the Yps3p protein of Histoplasma capsulatum. Eukaryot Cell. 2005, 4 (4): 685-693. 10.1128/EC.4.4.685-693.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Bohse ML, Woods JP: RNA interference-mediated silencing of the YPS3 gene of Histoplasma capsulatum reveals virulence defects. Infect Immun. 2007, 75 (6): 2811-2817. 10.1128/IAI.00304-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Brandhorst T, Wuthrich M, Finkel-Jimenez B, Klein B: A C-terminal EGF-like domain governs BAD1 localization to the yeast surface and fungal adherence to phagocytes, but is dispensable in immune modulation and pathogenicity of Blastomyces dermatitidis. Mol Microbiol. 2003, 48 (1): 53-65. 10.1046/j.1365-2958.2003.03415.x.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.PubMed CentralView ArticlePubMedGoogle Scholar
- Hwang L, Hocking-Murray D, Bahrami AK, Andersson M, Rine J, Sil A: Identifying phase-specific genes in the fungal pathogen Histoplasma capsulatum using a genomic shotgun microarray. Mol Biol Cell. 2003, 14 (6): 2314-2326. 10.1091/mbc.E03-01-0027.PubMed CentralView ArticlePubMedGoogle Scholar
- Voorhies M, Foo CK, Sil A: Experimental annotation of the human pathogen Histoplasma capsulatum transcribed regions using high-resolution tiling arrays. BMC Microbiol. 2011, 11: 216-10.1186/1471-2180-11-216.PubMed CentralView ArticlePubMedGoogle Scholar
- Nittler MP, Hocking-Murray D, Foo CK, Sil A: Identification of Histoplasma capsulatum transcripts induced in response to reactive nitrogen species. Mol Biol Cell. 2005, 16 (10): 4792-4813. 10.1091/mbc.E05-05-0434.PubMed CentralView ArticlePubMedGoogle Scholar
- Inglis DO, Voorhies M, Hocking Murray DR, Sil A: Comparative transcriptomics of infectious spores from the fungal pathogen Histoplasma capsulatum reveals a core set of transcripts that specify infectious and pathogenic states. Eukaryot Cell. 2013, 12 (6): 828-852. 10.1128/EC.00069-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Magrini V, Goldman WE: Molecular mycology: a genetic toolbox for Histoplasma capsulatum. Trends Microbiol. 2001, 9 (11): 541-546. 10.1016/S0966-842X(01)02204-1.View ArticlePubMedGoogle Scholar
- Kugler S, Young B, Miller VL, Goldman WE: Monitoring phase-specific gene expression in Histoplasma capsulatum with telomeric GFP fusion plasmids. Cell Microbiol. 2000, 2 (6): 537-547. 10.1046/j.1462-5822.2000.00078.x.View ArticlePubMedGoogle Scholar
- Retallack DM, Deepe GS, Woods JP: Applying in vivo expression technology (IVET) to the fungal pathogen Histoplasma capsulatum. Microb Pathog. 2000, 28 (3): 169-182. 10.1006/mpat.1999.0337.View ArticlePubMedGoogle Scholar
- Woods JP, Heinecke EL, Goldman WE: Electrotransformation and expression of bacterial genes encoding hygromycin phosphotransferase and beta-galactosidase in the pathogenic fungus Histoplasma capsulatum. Infect Immun. 1998, 66 (4): 1697-1707.PubMed CentralPubMedGoogle Scholar
- Woods JP, Retallack DM, Heinecke EL, Goldman WE: Rare homologous gene targeting in Histoplasma capsulatum: disruption of the URA5 Hc gene by allelic replacement. J Bacteriol. 1998, 180 (19): 5135-5143.PubMed CentralPubMedGoogle Scholar
- Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, Lightfoot S, Menzel W, Granzow M, Ragg T: The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006, 7: 3-10.1186/1471-2199-7-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- Stanke M, Diekhans M, Baertsch R, Haussler D: Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008, 24 (5): 637-644. 10.1093/bioinformatics/btn013.View ArticlePubMedGoogle Scholar
- Stanke M, Waack S: Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003, 19 (Suppl 2): ii215-ii225. 10.1093/bioinformatics/btg1080.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652. 10.1038/nbt.1883.PubMed CentralView ArticlePubMedGoogle Scholar
- Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK, Hannick LI, Maiti R, Ronning CM, Rusch DB, Town CD: Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003, 31 (19): 5654-5666. 10.1093/nar/gkg770.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34 (Web Server issue): W369-W373.PubMed CentralView ArticlePubMedGoogle Scholar
- Shapiro MB, Senapathy P: RNA splice junctions of different classes of eukaryotes: sequence statistics and functional implications in gene expression. Nucleic Acids Res. 1987, 15 (17): 7155-7174. 10.1093/nar/15.17.7155.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-515. 10.1038/nbt.1621.PubMed CentralView ArticlePubMedGoogle Scholar
- Garber M, Grabherr MG, Guttman M, Trapnell C: Computational methods for transcriptome annotation and quantification using RNA-seq. Nat Methods. 2011, 8 (6): 469-477. 10.1038/nmeth.1613.View ArticlePubMedGoogle Scholar
- Holbrook ED, Edwards JA, Youseff BH, Rappleye CA: Definition of the extracellular proteome of pathogenic-phase Histoplasma capsulatum. J Proteome Res. 2011, 10 (4): 1929-1943. 10.1021/pr1011697.PubMed CentralView ArticlePubMedGoogle Scholar
- Hwang LH, Mayfield JA, Rine J, Sil A: Histoplasma requires SID1, a member of an iron-regulated siderophore gene cluster, for host colonization. PLoS Pathog. 2008, 4 (4): e1000044-10.1371/journal.ppat.1000044.PubMed CentralView ArticlePubMedGoogle Scholar
- Bookout AL, Mangelsdorf DJ: Quantitative real-time PCR protocol for analysis of nuclear receptor signaling pathways. Nucl Recept Signal. 2003, 1: e012-PubMed CentralView ArticlePubMedGoogle Scholar
- Chang YY, Neufeld TP: An Atg1/Atg13 complex with multiple roles in TOR-mediated autophagy regulation. Mol Biol Cell. 2009, 20 (7): 2004-2014. 10.1091/mbc.E08-12-1250.PubMed CentralView ArticlePubMedGoogle Scholar
- Kamada Y, Sekito T, Ohsumi Y: Autophagy in yeast: a TOR-mediated response to nutrient starvation. Curr Top Microbiol Immunol. 2004, 279: 73-84. 10.1007/978-3-642-18930-2_5.PubMedGoogle Scholar
- Stipanuk MH: Macroautophagy and its role in nutrient homeostasis. Nutr Rev. 2009, 67 (12): 677-689. 10.1111/j.1753-4887.2009.00252.x.View ArticlePubMedGoogle Scholar
- Gargano S, Di Lallo G, Kobayashi GS, Maresca B: A temperature-sensitive strain of Histoplasma capsulatum has an altered delta 9-fatty acid desaturase gene. Lipids. 1995, 30 (10): 899-906. 10.1007/BF02537480.View ArticlePubMedGoogle Scholar
- Carratu L, Franceschelli S, Pardini CL, Kobayashi GS, Horvath I, Vigh L, Maresca B: Membrane lipid perturbation modifies the set point of the temperature of heat shock response in yeast. Proc Natl Acad Sci USA. 1996, 93 (9): 3870-3875. 10.1073/pnas.93.9.3870.PubMed CentralView ArticlePubMedGoogle Scholar
- Steczkiewicz K, Knizewski L, Rychlewski L, Ginalski K: TOS1 is circularly permuted 1,3-beta-glucanase. Cell Cycle. 2010, 9 (1): 201-204. 10.4161/cc.9.1.10510.View ArticlePubMedGoogle Scholar
- Yin QY, de Groot PW, Dekker HL, de Jong L, Klis FM, de Koster CG: Comprehensive proteomic analysis of Saccharomyces cerevisiae cell walls: identification of proteins covalently attached via glycosylphosphatidylinositol remnants or mild alkali-sensitive linkages. J Biol Chem. 2005, 280 (21): 20894-20901. 10.1074/jbc.M500334200.View ArticlePubMedGoogle Scholar
- Heilmann CJ, Sorgo AG, Klis FM: News from the fungal front: wall proteome dynamics and host-pathogen interplay. PLoS Pathog. 2012, 8 (12): e1003050-10.1371/journal.ppat.1003050.PubMed CentralView ArticlePubMedGoogle Scholar
- Marion CL, Rappleye CA, Engle JT, Goldman WE: An alpha-(1,4)-amylase is essential for alpha-(1,3)-glucan production and virulence in Histoplasma capsulatum. Mol Microbiol. 2006, 62 (4): 970-983. 10.1111/j.1365-2958.2006.05436.x.View ArticlePubMedGoogle Scholar
- Worsham PL, Goldman WE: Quantitative plating of Histoplasma capsulatum without addition of conditioned medium or siderophores. J Med Vet Mycol. 1988, 26 (3): 137-143. 10.1080/02681218880000211.View ArticlePubMedGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35 (Web Server issue): W182-W185.PubMed CentralView ArticlePubMedGoogle Scholar
- Gavrilin MA, Bouakl IJ, Knatz NL, Duncan MD, Hall MW, Gunn JS, Wewers MD: Internalization and phagosome escape required for Francisella to induce human monocyte IL-1beta processing and release. Proc Natl Acad Sci USA. 2006, 103 (1): 141-146. 10.1073/pnas.0504271103.PubMed CentralView ArticlePubMedGoogle Scholar
- Zemska O, Rappleye CA: Agrobacterium-mediated insertional mutagenesis in Histoplasma capsulatum. Methods Mol Biol. 2012, 845: 51-66. 10.1007/978-1-61779-539-8_4.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.