Gene expression in human fungal pathogen Coccidioides immitis changes as arthroconidia differentiate into spherules and mature
- Suganya Viriyakosol†1,
- Akul Singhania†2,
- Joshua Fierer1, 2, 3,
- Jonathan Goldberg4,
- Theo N Kirkland1, 2, 3Email author and
- Christopher H Woelk†1, 2
© Viriyakosol et al.; licensee BioMed Central Ltd. 2013
Received: 14 January 2013
Accepted: 20 May 2013
Published: 28 May 2013
Coccidioides immitis is a dimorphic fungus that causes disease in mammals, including human beings. It grows as a mycelium containing arthroconidia in the soil and in the host arthroconidia differentiates into a unique structure called a spherule. We used a custom open reading frame oligonucleotide microarray to compare the transcriptome of C. immitis mycelia with early (day 2) and late stage (day 8) spherules grown in vitro. All hybridizations were done in quadruplicate and stringent criteria were used to identify significantly differentially expressed genes.
22% of C. immitis genes were differentially expressed in either day 2 or day 8 spherules compared to mycelia, and about 12% of genes were differentially expressed comparing the two spherule time points. Oxireductases, including an extracellular superoxide dismutase, were upregulated in spherules and they may be important for defense against oxidative stress. Many signal transduction molecules, including pleckstrin domain proteins, protein kinases and transcription factors were downregulated in day 2 spherules. Several genes involved in sulfur metabolism were downregulated in day 8 spherules compared to day 2 spherules. Transcription of amylase and α (1,3) glucan synthase was upregulated in spherules; these genes have been found to be important for differentiation to yeast in Histoplasma. There were two homologs of 4-hydroxyphenylpyruvate dioxygenase (4-HPPD); transcription of one was up- and the other downregulated. We tested the effect of a 4-HPPD inhibitor, nitisinone, on mycelial and spherule growth and found that it inhibited mycelial but not spherule growth.
Transcription of many genes was differentially expressed in the process of arthroconidia to spherule conversion and spherule maturation, as would be expected given the magnitude of the morphologic change. The transcription profile of early stage (day 2) spherules was different than late stage (day 8) endosporulating spherules. In addition, very few genes that are important for spore to yeast conversion in other dimorphic fungi are differentially expressed in C. immitis mycelia and spherules suggesting that dimorphic fungi may have evolved different mechanisms to differentiate from mycelia to tissue invasive forms.
KeywordsCoccidioides Fungus Spherules Mycelia Gene expression Transcriptome
There are two closely related species of Coccidioides: immitis and posadasii that are found in different geographic areas . They belong to the order Onygenales and are members of the phylum Ascomycota. Both Coccidioides species are indigenous to the New World where they grow as molds in the alkaline desert soils, primarily in North America, but also in scattered desert areas in South America . These organisms are pathogenic for mammals (including humans), and it is estimated that there are ~150,000 infections annually in the US, primarily in the southwestern region . The soil form of these fungi is a mold that produces infectious spores (arthroconidia) that can become airborne if the soil is disturbed. Arthroconidia are ~ 4 micron in diameter and when inhaled into the lung they can cause pneumonia. Inside the host, under the influence of temperature and partial pressure of CO2, the organism undergoes a remarkable transformation into spherules, the pathognomonic tissue form of Coccidioides. Arthroconidia first round up and then they start to enlarge and transform. As they grow their cytoplasm undergoes internal segmentation to produce hundreds of endospores that are released when a spherule ruptures [4, 5]. These endospores in turn enlarge into spherules and replication continues until the host immune response controls the process or the host dies. The two species of Coccidioides have the same life cycle and there is no known difference in the clinical disease caused by infection with the two species.
The natural history of coccidioidomycosis is very variable. About 60% of infections are mild and go undiagnosed, but in Arizona (a hyper-endemic region) coccidioidomycosis is a leading cause of symptomatic pneumonia . Most of those infections resolve spontaneously but they can leave residual solitary granulomas or occasionally thin-walled cavities. In about 5% of cases infection does not remain confined to the lung and spreads to extra-pulmonary sites. This spread can be an overwhelming, life threatening process, or it can manifest as isolated skin, bone, joint, or meningeal infections. The last is uniformly fatal without treatment. Most people with dissemination suffer from prolonged and debilitating infections that are difficult to treat . People who are immunosuppressed, either by disease or iatrogenically, are at high risk for dissemination but the majority of disseminated cases occur in previously healthy individuals with no known immunological defects .
As with all the dimorphic pathogenic fungi (Blastomyces dermatitidis, Histoplasma capsulatum, Paracoccidioides brasiliensis and Sporothrix schenckii) the pathogenic form in tissue looks completely different form the saprobic mycelial form found in the environment. In coccidioidomycosis spherule formation is required for pathogenicity , as exemplified by two mutant strains [10, 11]. Both of these mutants fail to form spherules and although mice can be infected, the infection resolves spontaneously.
One fundamental but poorly understood issue that is relevant to all the dimorphic pathogenic fungi is how they differentiate from a mold (i.e., arthroconidia in mycelia) to the pathogenic form (i.e., spherules). It is possible to induce spherule formation in vitro by incubating arthroconidia at an elevated temperature (42°C) in a 14% CO2 atmosphere in a medium designed to promote the growth of spherules (Converse media) . We chose to study gene expression in early spherules (day 2 in culture) that have not yet begun to form endospores and late spherules (day 8 in culture) that have formed endospores. The development of early and late spherules has been described [4, 5]. C. immitis spherules do not rupture and release endospores when cultured in Converse media in our hands. We chose to compare gene expression in early and late spherules to mycelial gene expression to see whether gene expression varied as the spherules matured. We analyzed gene expression using a custom C. immitis oligonucleotide array platform constructed to probe the expression of every known and predicted open reading frame (ORF). Our hypothesis was that a large fraction of the genome would be differentially expressed in spherules compared to mycelia. We also hypothesized that many of the genes that are known to be important for mycelial to yeast conversion in other dimorphic pathogenic fungi would also be differentially expressed in the transition to spherules. Microarray gene expression analysis identified a large number of genes differentially expressed between the mycelial and spherule forms of the pathogen. The protein families (PFAM) and gene ontology (GO) terms significantly over-represented in the sets of differentially expressed genes were identified in order to better understand the higher biological processes being affected. Many genes associated with such families or terms were confirmed by real-time quantitative PCR (RT-qPCR). A study of C. immitis gene expression by Whiston et al. using RNA-Seq comparing transcript differences between mycelia and day 4 spherules was recently published and allowed us to compare our results to their results obtained at a time point intermediate in spherule development .
Growth of mycelia and spherules
C. immitis RS strain directly isolated from infected mice was grown on Mycosel agar (3.6% Mycosel agar, 0.5% yeast extract, and 50 μg/ml gentamicin). The animal protocol for infecting mice was approved by the Subcommittee on Animal Studies #07-017. The plates were incubated at 30°C until the mycelia covered the surface of the agar. Arthroconidia were harvested from the plate after 6 weeks of incubation at 25°C by adding 25 ml of saline. The plate was gently scraped using cell scraper and the fluid transferred to a 50 ml tube that was then vigorously mixed for 10 sec and centrifuged at 3000 rpm for 10 min at 4°C. The supernatant containing floating mycelia was discarded. The pellet containing arthroconidia was re-suspended in 10 ml saline and the suspension was passed through a 40 μM nylon cell strainer (BD Falcon). The strained suspension was centrifuged again and the pellet used to produce mycelia and spherules.
To grow mycelia, arthroconidia were washed 2 times with glucose-yeast extract (GYE) media and 2×106 spores/ml were incubated in 250 ml flat-bottom Erlenmeyer flasks (Corning) in 50 ml GYE media. Four flasks were cultured in a 30°C incubator without shaking for 5 days. To grow spherules, arthroconidia were washed 2 times in modified Converse media . The spores were inoculated at 4×106 arthroconidia/ml into a 250 ml baffled Erlenmeyer flask containing 50 ml of modified Converse media. Eight identical flasks were set up and grown on a shaker at 160 rpm, in 14% CO2 at 42°C. Four flasks were harvested 2 days after inoculation and the remaining four flasks after 8 days. The spherules did not rupture and release endospores within that time in this culture system.
Inhibition of growth with nitisinone
Nitisinone, 2-(2-nitro-4-trifluoromethylbenzoyl)-cyclohexane-1, 3 dione, a potent specific inhibitor of 4-HPPD was purchased from Swedish Orphan Biovitrum, Sweden. A stock solution of 30 mg/ml was made in 0.2 M NaOH. Nitisinone was added at several concentrations to glucose yeast extract media (GYE) or modified converse media in the presence of 2×106 spores/ml in a 15 ml round-bottom tissue culture tubes (BD Falcon). The culture was grown as described above for mycelial and spherule growth. The control tubes contained equal amounts of 0.2 M NaOH without Nitisinone. For microscopy, 1% formaldehyde was added to the culture overnight and the tubes were centrifuged 10,000 rpm for 10 min. The pellet was re-suspended in Lactophenol Aniline blue stain (Remel) and examined microscopically.
C. immitis mycelia were harvested by straining the media from four cultures through a 40 μM nylon cell strainer (BD Falcon). The mycelia were picked up from the cell strainer using a sterile disposable loop (BD Falcon) and dropped in a 2 ml ZR BashingBead lysis tube with 0.5 mm beads (Zymoresearch) and 0.5 ml Qiazol reagent (Qiagen). The tubes were arranged in a pre-cooled Tissuelyzer II adapter (Qiagen) and mycelia was disrupted by shaking at 50 Hz for 25 min. Spherules in Converse media were harvested from four 2 day cultures and four 8 day cultures. The cell concentration was determined by counting the spherules in Lactophenol Aniline blue stain. The media was centrifuged at 10,000 rpm for 10 min at 4°C. Qiazol (Qiagen) was added to the cell pellet at 4×106 spherules/ml and 0.5 ml of the mixture added to a 2 ml ZR BashingBead lysis tube with 0.5 mm beads (Zymoresearch). Total RNA was purified from mycelia and spherule samples (4 replicates/condition) using the RNeasy Microarray tissue mini-kit (Qiagen) in a Qiacube machine (Qiagen). If necessary RNA was concentrated or re-purified using RNeasy Minelute Cleanup kit (Qiagen) according to the manufacturer’s protocol. Total RNA concentrations were initially determined using a NanoDrop spectrophotometer (NanoDrop) and RNA quality was assessed using a 2100 Bioanalyzer (Agilent). All samples had a RNA integrity number greater than 7.
Microarray design and hybridization
Known and predicted ORFs from the C. immitis genome (RS strain) were previously identified using sequence data available at the Broad Institute . This information was supplied to Roche Nimblegen in order to manufacture a custom oligonucleotide array consisting of 68,927 probes (Nimblegen custom array OID30589). Probes were 60 nucleotides in length and the expression of the majority of genes was assayed using 7 different probes printed in duplicate. The expression of small genes was assayed with fewer probes. Twelve custom microarrays fit on a single slide such that all the samples in this study (4 × mycelia, 4 × day 2 spherule, and 4 × day 8 spherule) could be assayed for gene expression in a single experiment to eliminate technical batch effects. Ten μg of total RNA at a concentration greater than 1 μg/ml from each sample was used for microarray hybridization. Total RNA was converted to cDNA, labeled with dye, and hybridized to the microarray by the VA San Diego Gene Chip Microarray Core according to the Nimblegen protocol. All C. immitis genes are referred to by their locus tag and further information about these genes can be found at the Coccidioides group database at the Broad Institute http://www.broadinstitute.org/annotation/genome/coccidioides_group/MultiHome.html. FungiDB (http://fungidb.org/fungidb/) was also used for annotation because it has more informative gene names for many genes.
Microarray data analysis
Quality control analysis and normalization of microarray gene expression data were performed as previously described . Briefly, several quality control assessments (e.g., boxplots and volcano plots) were applied to assess microarray data quality. Unsupervised clustering was also performed using the web-based tool ANAIS  to determine if samples clustered as expected based on the expression of genes in each sample. All arrays passed quality control filters and no outliers were found. Differentially expressed probes were identified between mycelium, day 2 spherule and day 8 spherule conditions using a one-way ANOVA and the Tukey post hoc test implemented in GeneSpring GX version 11.5 (Agilent Technologies Inc.). The false discovery rate (FDR) associated with multiple tests was corrected for using the Benjamini-Hochberg method . In a conservative approach, a gene was only identified as differentially expressed if all probes for that gene had a fold change greater than 2 or less than −2 and an ANOVA p-value (Tukey and FDR corrected) less than 0.05. Fold changes were calculated for each gene that passed this filter by averaging across the seven probes. The overlap in up and downregulated genes between day 2 and day 8 spherules compared to mycelia was visualized in Venn diagrams using BxToolBox (http://bioinforx.com/). The ANOVA analysis was also used to identify genes that were differentially expressed between day 2 and day 8 spherules where positive fold changes are indicative of greater expression at day 8 compared to day 2 and negative fold changes suggest decreased expression. Gene expression data are available at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) under accession number GSE44225.
PFAM and GO analysis
PFAM enrichment was determined using a tool at the Broad Institute http://www.broadinstitute.org/annotation/genome/coccidioides_group/BatchSelect.html?target=GeneEnrichment.html. This tool looks for over-representation of PFAMs in up- or downregulated genes using a hypergeometric test, and only PFAMs with an FDR-corrected p-value <0.05 were considered significant.
GO terms were assigned to C. immitis genes by reciprocal homology searches at the protein level against the Saccharomyces cerevisiae proteome using BLAST (Additional file 1: Table S1). UniProt IDs were obtained using the C. posadasii homologs of C. immitis genes because many more C. posadasii genes have UniProt IDs. The Biological Networks Gene Ontology (BiNGO) plugin (version 2.441)  for Cytoscape (version 2.8.3) was used to identify those GO terms related to biological processes that were over-represented for differentially expressed genes identified between each of the three comparison groups (mycelia vs. day 2 spherule, mycelia vs. day 8 spherule, day 8 vs. day 2 spherule). BiNGO preserves the hierarchical relationship between GO terms. Significance was assessed with a hypergeometric test and only GO terms with an FDR-corrected p-value <0.05 were considered significant.
C. immitis protein kinases were identified and classified by orthology with the curated Trichophyton rubrum kinome . Non-orthologous kinases were identified and classified by searching the proteome with a protein kinase HMM built from an alignment of Dictyostelium kinases  followed by a BLAST against the curated kinase database (http://kinase.com/) . Kinase abbreviations are provided in Additional file 2: Table S3. Signal peptides in the proteins coded for in the C. immitis genome were identified using artificial neural networks implemented in SignalP version 4.0 .
RT-qPCR confirmation of gene expression
Microarray gene expression was confirmed by RT-qPCR for 24 genes. Three highly expressed genes with low standard deviation across the 12 samples were selected as normalizers (CIMG_01599, CIMG_10083 and CIMG_12902). SYBR® Green primers were designed using Primer Express version 3.0 (Applied Biosystems Inc.) and obtained from Integrated DNA Technologies, Inc. (Coralville, IA). Reverse primers were designed to span a splice site in the same region of the gene probed by the microarray. Duplicate RT-qPCR reactions were performed for each sample using 50 ng of reverse transcribed RNA per reaction. Fold changes were calculated as described previously using the 2-∆∆CT method  implemented in the DataAssist software version 3.0 (ABI), and significance was determined using one-way ANOVA in the R statistical package (version 2.13.2).
Results and discussion
Genes differentially expressed in mycelia and spherules
A recent paper by Whiston et al. assessed transcription in C. immitis and C. posadasii mycelia and day 4 spherules by RNA-seq . We have compared our results to theirs. The two studies used different methods for assessing changes in gene expression. We used microarray technology to estimate transcript abundance while Whiston et al. used RNA-seq to estimate transcript abundance . The literature suggests that these methods should yield comparable results . Despite this difference in methodology, we confirmed the upregulation of 25% of the genes that Whiston found to be upregulated in spherules. Conversely, 43% of genes that we have found to be upregulated in day 2 and day 8 spherules were also upregulated in day 4 spherules in the Whiston study (Additional file 5: Figure S2). Despite the differences in the two studies many of our conclusions are similar (see below).
We know from previous experiments that some genes are overexpressed in spherules compared to mycelia. Some of these genes, such as the spherule outer wall glycoprotein (CIMG_04613)  and the parasitic-phase specific protein PSP-1 (CIMG_05758)  were up regulated more than four fold in spherules in this experiment (Additional file 4: Table S2). Other genes, such as the metalloproteinase Mep1 (CIMG_06703), which has been found to be expressed at high levels in endosporulating spherules in C. posadasii was not found to be over-expressed in this experiment . We also examined the expression level of the Mep1 gene by RT-qPCR and found that its expression was slightly downregulated in spherules compared to mycelia, rather than upregulated as previously reported (see below). Whiston et al. also examined the expression of this gene and found that it was upregulated in C. posadasii spherules but not C. immitis spherules .
Confirmation of differential expression by RT-qPCR
Twenty-four differentially expressed genes as detected by microarray analysis were selected for confirmation by RT-qPCR (Figure 3). Genes were selected for RT-qPCR confirmation of gene expression based on the magnitude of fold change (up- or downregulation) between mycelia and day 2 spherules, mycelia and day 8 spherules, and day 2 and day 8 spherules, and their identification in the PFAM or GO analysis. The significant differential expression (p < 0.05, t-test) of each of these 24 genes was confirmed for at least one of the three comparison groups. In the majority of cases, RT-qPCR analysis detected larger differences in expression than microarray analysis. This is probably because RT-qPCR has a greater dynamic range than microarray does .
Day 2 and day 8 spherules vs mycelia
PFAM functional enrichment
# Specified Genes
Specified Gene %
Whole Genome %
Corrected P Value
Day 2 Spherules Upregulated
Thioesterase superfamily (4HBT)
Short chain dehydrogenase
Day 8 Spherules Upregulated
Day 8 vs 2 Spherules Upregulated
Day 2 Spherules Downregulated
Day 8 Spherules Downregulated
Day 8 vs 2 Spherules Downregulated
The short chain dehydrogenases family was also upregulated in day 2 spherules (maximum upregulation 10.27 fold, CIMG_09765). This family of enzymes catalyzes oxidation/reduction reactions of alcohols and cyclic compounds. Up regulation of this family of enzymes seems plausible given the shift in growth conditions from air which contains less than 0.05% CO2 (mycelial growth) to 14% CO2 (spherule growth) which mimics the shift in oxidation/reduction potential that occurs when the organism grows in the mammalian host. The aldol-keto reductase family was also significantly enriched in both day 2 and day 8 spherules. These genes were also found to be upregulated in spherules by Whiston et al. . This protein family plays a role in reducing oxidative stress  and may be required for resistance to the oxidative burst in mammals or to mitochondrial generated reactive oxygen species. C. immitis spherules are more resistant to oxidative killing in vitro than Aspergillus fumigatus spores .
The major facilitator super family (MFS-1) that was enriched in day 8 spherules is an important transporter of small molecules and includes a number of transporters for uptake as well as efflux . Two highly upregulated genes are sugar transporters (CIMG_03001 and CIMG_08310). A nicotinic acid transporter (CIMG_06071) and a sialate transporter (CIMG_03956) were also upregulated. The nicotinic acid transporter is presumably involved in NAD metabolism ; we have been unable to find a role for the sialate transporter in fungi in the literature.
The pleckstrin domain occurs in a wide range of proteins involved in intracellular signaling or as constituents of the cytoskeleton. Pleckstrin domain transcripts were downregulated in day 2 spherules; in fact, one pleckstrin domain gene is the most downregulated of all the day 2 genes (CIMG_07982, -53.53 fold). The downregulated pleckstrin domain containing genes may be required for polar mycelial growth but not isotropic spherule growth. One downregulated gene in this family is the anucleate primary sterigmata protein A (CIMG_06141, -4.93), which is critical for movement of nuclei into spores on the sterigmata of A. nidulans. This gene may well be required for arthroconidia formation in C. immitis mycelia but not endospore formation in spherules.
A significant proportion of proteins containing SH3 domains were downregulated in day 2 spherules. SH3 protein families include some protein kinases, phosphoinositol 3 kinases, Ras GTPase activating proteins, and the guanine nucleotide exchange factors cdc24 and cdc25. Two of these genes, CIMG_04361 and CIMG_04531, were downregulated in day 2 spherules. CIMG_04531 is annotated as a polarized growth protein, and is highly homologous to cytoskeleton assembly proteins in many fungi. CIMG_02193 is cytoskeletal protein SLA1 and it is downregulated (−4.61 fold change) in day 2 spherules. Perhaps these proteins predispose to polar mycelial growth rather than isotropic spherule growth.
Modulated protein kinases in day 2 and day 8 spherules
C. immitis annotation
Serine/threonine-protein kinase; meiosis induction protein kinase
Kinase domain containing protein
Protein kinase domain containing protein
Protein kinase domain containing protein; serine/threonine protein kinase 24
serine/threonine protein kinase ssp1
serine/threonine protein kinase
Kinase domain containing protein; checkpoint kinase
Protein kinase, sky 1
serine/threonine protein kinase
Ran1-like protein kinase
cytokenesis protein sepH
serine/threonine protein kinase ATG1
MAP kinase kinase kinase SskB
serine/threonine protein kinase
Kinase domain containing protein
Kinase domain containing protein
Response regulator receiver RIM15p
Serine threonine protein kinase
Kinase domain containing protein
Protein kinase domain containing protein
Protein kinase domain containing protein
RWD domain protein
C2H2 zinc finger domain containing proteins were downregulated in day 2 spherules. Most of the proteins containing this domain are transcription factors and the zinc finger is involved in DNA binding [43, 44]. Some of the genes that were downregulated include the transcription factors CIMG_04642 (−9.24, FlbC), CIMG_03725 (−5.06, zinc finger transcription factor PacC) and CIMG_06050 (−3.06, transcription factor steA). In fact, steA is a negative regulator of transcription in Aspergillus, so downregulation of this gene probably results in upregulated transcription of some genes . Ste12 is the Saccharomyces homolog of this gene . Ste12 is involved in the mating response and is involved in the up- and downregulation of many genes. Eight of these 19 C2H2 zinc finger genes were also found to be downregulated in spherules by Whiston et al. .
Day 8 spherule/day 2 spherule comparison
Several of the gene families that were downregulated in day 2 spherules were upregulated in the day 8 spherules (Table 1). Examples are the Ras GTPase activating proteins, the guanine nucleotide exchange factors cdc24 and cdc25. 13 of 19 of the kinases downregulated in the day 2 spherules had returned to mycelial levels in day 8 spherules (Table 2). Two genes in this family were downregulated in both day 2 and day 8 spherules: CIMG_00940 (−5.28 fold in day 2 spherules and −10.04 in day 8 spherules both compared to mycelia) and CIMG_04103 (−3.97 fold in day 2 spherules and −6.75 in day 8 spherules both compared to mycelia). CIMG_00940 was also found to be downregulated in spherules by Whiston et al. . CIMG_00940 is a Swe1 kinase and CIMG_04103 is a STE/STE11/CDC15 kinase. Both of these genes are involved in regulation of mitosis [47–49]. One function of Wee kinases in S. cerevisiae is to prevent small cells from entering mitosis ; endospores are very small so downregulation of this gene may be important for endospore division. A function of the CDC15 kinases is to bind the spindle pole body and facilitate exit from mitosis . There is no obvious reason why this kinase should be downregulated in the internally dividing spherule.
Gene ontology analysis
Mycelia/day 8 spherules comparison
The gene ontology (GO) terms significantly over-represented in lists of differentially expressed genes identified between mycelia and either day 2 or day 8 spherules, as well as between day 2 and day 8 spherules, were identified using BiNGO . There were no GO terms that survived FDR correction between mycelia and day 2 spherules but a large number of significant terms were identified between mycelia and day 8 spherules (Additional file 6: Figure S3). The most significant enriched GO term was “small molecule metabolic process” (corrected p = 0.004). Thirty-one members of this heterogeneous set of genes were upregulated and 75 were downregulated. Twelve of the downregulated genes coded for nucleotide synthesis or DNA replication. For example, a homeobox domain-containing protein was downregulated −8.68 fold (CIMG_09071); thymidylate synthase was down −3.57 fold (CIMG_08646); cell division control protein Cdc6 was down −3.05 fold (CIMG_07523) and DNA topoisomerase 2 was down −3.09 fold (CIMG_02836). This suggests that the rate of DNA synthesis is slower in the day 8 spherules than in mycelia. 10 genes coding for amino acid synthesis were downregulated as well. This suggests that not only is DNA synthesis relatively slow compared to mycelia but protein synthesis is too. Other genes involved in vitamin synthesis and energy generation were also downregulated. This is consistent with the notion that day 8 spherules have produced their endospores. Rupturing and releasing endospores should not be a metabolically expensive process. The observation that MFS-1 sugar transporters are upregulated suggests that that the low metabolic needs may not be universal.
The most strongly upregulated genes in day 8 spherules with the GO term “small metabolic process” included glutamate decarboxylase (21.47), three ABC transporters and parasitic phase specific protein-1 (6.66) previously described by Delgado . The PSP-1 gene is also upregulated in day 2 spherules and in day 4 spherules as reported by Whiston et al. . PSP-1 contains a RTA-1 domain, which is involved in resistance to 7-aminocholesterol . This family of proteins has multiple membrane spanning domains and is thought to be involved in binding 7-aminocholesterol and related substances and preventing toxicity. They are not thought to be efflux pumps .
A group of genes assigned the GO term “carbohydrate metabolic processes” was also enriched in the day 8 spherules dataset. 15 genes were upregulated and 17 genes were downregulated. The upregulated genes included polysaccharide deacetylase (CIMG_02628, 34.82) and 1,4 (α)-amylase (CIMG_03529, 2.70). The most striking downregulated gene in this group is calmodulin (CIMG_04786, -10.38). Two other genes coding for calmodulin (CIMG_02413 and CIMG_08162) are not differentially expressed in day 8 spherules. We looked for differential expression of six calmodulin-dependent kinases and found that they were not up- or downregulated. Calmodulin-dependent signaling enzymes are required for growth of Cryptococcus neoformans at 37°C  and for conversion of Sporothrix shenkii from mold to yeast .
Another enriched GO term in day 8 spherules was “oxidation-reduction processes” (corrected p-value = 0.012). In this group of genes, twenty-seven genes were upregulated in spherules and 27 were downregulated in spherules compared to mycelia. The five most upregulated genes were aldol-keto reductases (maximum 30.2 fold), alcohol dehydrogenases (maximum 11.65 fold) and nitrate reductase (8.06). The aldol-keto reductases were also identified in the PFAM functional enrichment in day 2 and day 8 spherules (Table 1). All these responses make sense in terms of the difference in growth conditions of mycelia (grown in air containing less than 0.05% CO2) and spherules (grown in 14% CO2 in air). These responses also seem reasonable in terms of the spherule living in the high CO2 environment of the mammalian host compared to the mycelia living in the soil. Other upregulated genes include Y20 (CIMG_04756), which was upregulated 11.18 fold in day 2 spherules and 6.81 fold in day 8 spherules. This gene was also upregulated in spherules in the Whiston study . This gene codes for a flavodoxin that plays a role in the cellular responses to oxidative stress . A homolog of this gene is highly upregulated in the yeast phase of P. brasiliensis. Presumably this protein is protecting the fungus against oxidative attack by the mammalian host. Additional genes coding for response to oxidative stress that were upregulated in day 8 spherules include a Cu/Zn super-oxide dismutase (CIMG_06994, 5.88) and a catalase. CIMG_06994 was up regulated 6.51 fold in day 2 spherules too. This gene is highly homologous to the extracellular SOD3 (e = 4 × 10-50) recently identified as a secreted protein in the Histoplasma yeast phase . Gene deletion experiments have shown that this gene is important for defense against oxidative stress . SOD3 is a secreted protein in H. capsulatum; CIMG_06944 has a predicted signal sequence suggesting that it is a secreted protein too. Extracellular SOD may be more protective against mammalian oxidative stress as suggested by Youseff . Presumably the C. immitis homolog of SOD3 is up regulated to protect the spherule against oxidative stress in the host. C. immitis also contains genes highly homologous to A. fumigatus SOD2 and SOD4 but neither of those is up- or downregulated.
Several NAD or NADPH dependent oxireductases were downregulated 3–8 fold in day 8 spherules. An NADPH oxidase was downregulated 3.48 fold. This enzyme makes reactive oxygen intermediates in fungi just as it does in mammalian phagocytes . Reactive oxygen intermediates are important for apical growth and formation of spores in filamentous fungi [57, 58]. It is possible that the NADP oxidase may interfere with isotropic spherule growth and differentiation.
Day 8/day 2 spherule comparison
Several GO terms were significantly over-represented (FDR corrected p-value <0.05) in the list of genes differentially expressed between day 8 and day 2 spherules. The most significant enriched GO term was “sulfur compound metabolic process” (corrected p-value = 0.046). Sixteen genes were downregulated in this data set and one was upregulated. We see downregulation of 5′-methylthioadenosine phosphorylase (CIMG_01361, -7.45 fold), phosphoadenosine phosphosulfate reductase (CIMG_00456, -4.65 fold), two genes coding for adenylyl-sulfate kinase (CIMG_00454, -4.22 fold and CIMG_06918, -2.65 fold) and sulfite reductase (CIMG_00269, -2.94 fold) in day 8 spherules. Two of these genes were upregulated in day 2 spherules compared to mycelia. All these genes are involved in accumulating sulfide. This suggests that C. immitis spherules have no difficulty accumulating enough sulfur for their needs as they mature.
Upregulated or downregulated genes in day 2, day 4 and day 8 spherules
We identified 153 genes that were upregulated more than two fold in day 2 spherules, day 8 spherules and day 4 spherules in the Whiston study . 140 genes were downregulated more than two fold in all three spherule samples. 57% of the upregulated genes and 42% of the downregulated genes had no function annotation (Additional file 7: Table S4). Many of these unannotated genes were highly differentially expressed suggesting that they may be important for spherule development. One upregulated gene that has not been discussed is opsin-1 (CIMG_08103, maximal upregulation 25.72). This gene is closely related to the bacterial rhodopsin gene coding for G protein coupled receptors; its function in fungi has not been determined . Another gene that was upregulated at all three spherule time points was the sulfite transporter Ssu1 (CIMG_05899, maximum upregulation 6.37). Downregulated genes of interest include several glucosidases, glucanases and a chitosanase. Two septin genes are downregulated in spherules. Septin genes are important regulators of the cytoskeleton and play a role in determining cell shape . Why these genes are downregulated is unclear since the spherule is undergoing extensive cellular remodeling. Perhaps septins are required for polar growth and other regulators are needed for isotropic spherule growth. Further analysis of the relatively small group of genes that are consistently up- or downregulated throughout spherule development may be useful for understanding the pathogenic phase of this organism.
Comparison of results to those obtained in other pathogenic fungi
The dimorphic pathogenic fungi are phylogenetically closely related  so it is reasonable to ask if genes important for conidium to yeast transformation in those pathogens are also important for arthroconidia to spherule transformation in Coccidioides. One H. capsulatum gene that is required for mycelium to yeast transformation is the alpha amylase (AMY1) gene. This gene affects the amount of α (1,3) glucan in the yeast cell wall, and that carbohydrate appears to be important for the ability of the yeast to survive in macrophages and in mice . This gene has a nearly identical homolog in C. immitis, CIMG_03142, that was upregulated 3.6 fold in day 2 spherules and 3.39 fold in day 8 spherules. Whiston et al. also found it to be upregulated in spherules . Another H. capsulatum gene that is required for yeast formation is α glucan synthase (AGS1) gene . This enzyme catalyzes the production of α (1,3) glucan in the cell wall that obscures the β (1,3) glucan and prevents activation of innate immunity via the dectin-1 receptor . C. immitis has an AGS1 gene (CIMG_13256) that was upregulated in the day 8 spherule (2.48 fold) but not day 2 spherules. Whiston et al. found this gene to be upregulated 1.93 fold in spherules compared to mycelia . There is no literature describing the relative amounts of α (1,3) glucan and β (1,3) glucan in C. immitis mycelia or spherules. We know, however, that there is enough exposed β (1,3) glucan in Coccidioides spherules to stimulate macrophages to produce cytokines via dectin-1 .
Two genes coding for transcription factors, Ryp2 and Ryp3, have been found to be essential for conversion from filaments to yeast in H. capsulatum. These genes are overexpressed in the yeast phase of H. capsulatum. C. immitis has nearly identical homologs of these genes but they were not overexpressed in either day 2 or day 8 spherules, suggesting that they may not be required for the transformation from mycelium to spherule.
Gene disruption experiments in B. dermatitidis have shown that a histidine kinase, DRK1, is required for the transformation from filaments to yeast . It is not clear from the literature whether or not this gene is overexpressed in the B. dermatitidis yeast phase. C. immitis has a very closely related homolog of this gene (CIMG_04512) but it was not up or down regulated in day 2 or day 8 spherules. In another dimorphic pathogenic fungus, S. schenckii, the calcium/calmodulin kinase I gene (SSMK1) was found to be required for formation of yeast . There are two genes in C. immitis that are highly homologous to the S. schenckii SSMK1 gene; neither one of these was up- or downregulated in day 2 or day 8 spherules.
A number of studies have been done studying the transcriptome of P. brasiliensis[66, 67]. One study identified the 4-HPPD gene to be required for P. brasiliensis conidia to convert to yeast . They found that the 4-HPPD gene expression was upregulated in the yeast form and that a biochemical inhibitor of this enzyme, nitisinone, inhibited mycelium conversion to yeast. 4-HPPD (E.C. 1.13.1127) is an enzyme that converts 4-hydroxyphenylpyruvate to homogentisate that is involved in the synthesis of tyrosine, phenylalanine, and ubiqinone (KEGG, whttp://www.genome.jp/keg).
Conversion from the arthroconidia phase to the parasitic spherule phase in C. immitis requires major transcriptional reprogramming with 22% of the entire genome being differentially expressed between the two conditions. Further, gene expression within spherules is dynamic with 12% of the entire genome being differentially expressed as they mature from day 2 to day 8. It is evident from the transcriptional profile at day 2 compared to mycelia that differentiation of C. immitis is associated with the regulation of specific genes. For example, a number of genes were downregulated during mycelia to spherule conversion including transcriptional repressors (genes encoding zinc finger proteins), pleckstrin domain containing genes, and genes coding for proteins with SH3 signaling domains. Additionally, twenty-four protein kinase genes homologous to S. cerevisiae genes coding for sexual or meiotic function or mitosis or filamentous growth are downregulated and may play a role in arthroconidia differentiation to spherules. About 75% of the protein kinase genes return to mycelial levels of expression in 8 day spherules, suggesting they may be important in arthroconidia to spherule differentiation but not in spherule maturation. Some genes are persistently upregulated or downregulated in spherules at both time points. These include some genes that have previously been shown to be important for yeast development in H. capsulatum such as amylase gene AMY-1. An enzyme such as extra-cellular superoxide dismutase (SOD3),that provides protection against oxidative stress was also upregulated in both day 2 and day 8 spherules . In contrast to Paracoccidioides brasiliensis, where inhibition of the enzyme 4-HPPD inhibits conversion of the mold to the yeast form , inhibition of the enzyme 4-HPPD inhibits mycelial growth in C. immitis but has no effect in arthroconidia to spherule conversion in C. immitis. Arthroconidia to spherule conversion in C. immitis is a complex process requiring modulation of a large number of genes.
False discovery rate
Glucose yeast extract
Open reading frame
We thank Daniel Neafsey and Diego Martinez (Broad Institute, Cambridge MA, 02141) for help retrieving data, reviewing this manuscript, and for providing us with the Trichophyton rubrum kinome ahead of publication, respectively. We also thank Gerald Manning (Salk Institute, San Diego, CA) for providing us with updated protein kinase classifications ahead of publication at http://kinase.com/. This project was supported by funds received from the State of California, Department of Public Health, Award No. 09–11657 (T.N.K.), a grant from the Academic Senate of the University of California (T.N.K.) and The Research Service of the Veterans Administration (J.F.). This work was also performed with the support of the Genomics Core at the UCSD Center for AIDS Research funded through the National Institutes of Health (AI36214)(C.W.). This material is based upon work supported in part by the Department of Veterans Affairs, Veterans Health Administration, Office of Research and Development.
- Neafsey DE, Barker BM, Sharpton TJ, Stajich JE, Park DJ, Whiston E, Hung CY, McMahan C, White J, Sykes S, et al: Population genomic sequencing of Coccidioides fungi reveals recent hybridization and transposon control. Genome Res. 2010, 20 (7): 938-946. 10.1101/gr.103911.109.PubMedPubMed CentralView Article
- Fisher MC, Koenig GL, White TJ, San-Blas G, Negroni R, Alvarez IG, Wanke B, Taylor JW: Biogeographic range expansion into South America by Coccidioides immitis mirrors New world patterns of human migration. Proc Natl Acad Sci USA. 2001, 98 (8): 4558-4562. 10.1073/pnas.071406098.PubMedPubMed CentralView Article
- CDC: Increase in coccidioidomycosis - California, 2000–2007. MMWR Morb Mortal Wkly Rep. 2009, 58 (5): 105-109.
- Cole GT, Hung CY: The parasitic cell wall of Coccidioides immitis. Med Mycoll. 2001, 39 (Suppl 1): 31-40.View Article
- Cole GT, Sun SH: Arthroconidium-spherule-endospore transformation in Coccidioides immitis. Dimorphism. Edited by: Szaniszlo P. 1985, New York: Plenum, 281-333.View Article
- Chang DC, Anderson S, Wannemuehler K, Engelthaler DM, Erhart L, Sunenshine RH, Burwell LA, Park BJ: Testing for coccidioidomycosis among patients with community-acquired pneumonia. Emerg Infect Dis. 2008, 14 (7): 1053-1059. 10.3201/eid1407.070832.PubMedPubMed CentralView Article
- Chiller TM, Galgiani JN, Stevens DA: Coccidioidomycosis. Infect Dis Clin North Am. 2003, 17 (1): 41-57. 10.1016/S0891-5520(02)00040-5. viiiPubMedView Article
- Kirkland TN, Fierer J: Coccidioidomycosis: a reemerging infectious disease. Emerg Infect Dis. 1996, 2 (3): 192-199. 10.3201/eid0203.960305.PubMedPubMed CentralView Article
- Smith CE, Pappagianis D, Levine HB, Saito M: Human coccidioidomycosis. Bacteriol Rev. 1961, 25: 310-320.PubMedPubMed Central
- Kirkland TN, Fierer J: Inbred mouse strains differ in resistance to lethal Coccidioides immitis infection. Infect Immun. 1983, 40: 912-917.PubMedPubMed Central
- Xue J, Chen X, Selby D, Hung CY, Yu JJ, Cole GT: A genetically engineered live attenuated vaccine of Coccidioides posadasii protects BALB/c mice against coccidioidomycosis. Infect Immun. 2009, 77 (8): 3196-3208. 10.1128/IAI.00459-09.PubMedPubMed CentralView Article
- Converse JL: Effect of surface active agents on endosporulation of Coccidioides immitis in a chemically defined medium. J Bacteriology. 1957, 74: 106-107.
- Whiston E, Zhang Wise H, Sharpton TJ, Jui G, Cole GT, Taylor JW: Comparative transcriptomics of the saprobic and parasitic growth phases in Coccidioides spp. PLoS One. 2012, 7 (7): e41034-10.1371/journal.pone.0041034.PubMedPubMed CentralView Article
- Sharpton TJ, Stajich JE, Rounsley SD, Gardner MJ, Wortman JR, Jordar VS, Maiti R, Kodira CD, Neafsey DE, Zeng Q, et al: Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives. Genome Res. 2009, 19 (10): 1722-1731. 10.1101/gr.087551.108.PubMedPubMed CentralView Article
- Woelk CH, Beliakova-Bethell N, Goicoechea M, Zhao Y, Du P, Rought SE, Lozach J, Perez-Santiago J, Richman DD, Smith DM, et al: Gene expression before HAART initiation predicts HIV-infected individuals at risk of poor CD4+ T-cell recovery. AIDS. 2010, 24 (2): 217-222. 10.1097/QAD.0b013e328334f1f0.PubMedPubMed CentralView Article
- Simon A, Biot E: ANAIS: analysis of NimbleGen arrays interface. Bioinformatics. 2010, 26 (19): 2468-2469. 10.1093/bioinformatics/btq410.PubMedView Article
- Klipper-Aurbach Y, Wasserman M, Braunspiegel-Weintrob N, Borstein D, Peleg S, Assa S, Karp M, Benjamini Y, Hochberg Y, Laron Z: Mathematical formulae for the prediction of the residual beta cell function during the first two years of disease in children and adolescents with insulin-dependent diabetes mellitus. Med Hypotheses. 1995, 45 (5): 486-490. 10.1016/0306-9877(95)90228-7.PubMedView Article
- Maere S, Heymans K, Kuiper M: BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005, 21 (16): 3448-3449. 10.1093/bioinformatics/bti551.PubMedView Article
- Martinez DA, Oliver BG, Gräser Y, Goldberg JM, Li W, Martinez-Rossi NM, Monod M, Shelest E, Barton RC, Birch E, et al: Comparative genome analysis of Trichophyton rubrum and related dermatophytes reveals candidate genes involved in infection. MBio. 2012, 3 (5): e00259-12. 10.1128/mBio.00259-12. Print 2012.PubMedPubMed CentralView Article
- Goldberg JM, Manning G, Liu A, Fey P, Pilcher KE, Xu Y, Smith JL: The dictyostelium kinome–analysis of the protein kinases from a simple model organism. PLoS Genet. 2006, 2 (3): e38-10.1371/journal.pgen.0020038.PubMedPubMed CentralView Article
- Manning G, Whyte DB, Martinez R, Hunter T, Sudarsanam S: The protein kinase complement of the human genome. Science. 2002, 298 (5600): 1912-1934. 10.1126/science.1075762.PubMedView Article
- Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8 (10): 785-786. 10.1038/nmeth.1701.PubMedView Article
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView Article
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.PubMedPubMed CentralView Article
- Hung C, Ampel NM, Christian L, Seshan KR, Cole GT: A major cell surface antigen of Coccidioides immitis which elicits both humoral and cellular immune responses. Infect Immun. 2000, 68 (2): 584-593. 10.1128/IAI.68.2.584-593.2000.PubMedPubMed CentralView Article
- Delgado N, Hung CY, Tarcha E, Gardner MJ, Cole GT: Profiling gene expression in Coccidioides posadasii. Med Mycol. 2004, 42 (1): 59-71. 10.1080/1369378031000156890.PubMedView Article
- Hung CY, Seshan KR, Yu JJ, Schaller R, Xue J, Basrur V, Gardner MJ, Cole GT: A metalloproteinase of Coccidioides posadasii contributes to evasion of host detection. Infect Immun. 2005, 73 (10): 6689-6703. 10.1128/IAI.73.10.6689-6703.2005.PubMedPubMed CentralView Article
- Rajeevan MS, Vernon SD, Taysavang N, Unger ER: Validation of array-based gene expression profiles by real-time (kinetic) RT-PCR. JMD. 2001, 3 (1): 26-31.PubMedPubMed CentralView Article
- Smotrys JE, Linder ME: Palmitoylation of intracellular signaling proteins: regulation and function. Annu Rev Biochem. 2004, 73: 559-587. 10.1146/annurev.biochem.73.011303.073954.PubMedView Article
- Wheat RW, Tritschler C, Conant NF, Lowe EP: Comparison of Coccidioides immitis arthrospore, mycelium, and spherule cell walls, and influence of growth medium on mycelial cell wall composition. Infect Immun. 1977, 17 (1): 91-97.PubMedPubMed Central
- Chang Q, Petrash JM: Disruption of aldo-keto reductase genes leads to elevated markers of oxidative stress and inositol auxotrophy in Saccharomyces cerevisiae. Biochim Biophys Acta. 2008, 1783 (2): 237-245. 10.1016/j.bbamcr.2007.08.008.PubMedPubMed CentralView Article
- Margolis DA, Viriyakosol S, Fierer J, Kirkland TN: The role of reactive oxygen intermediates in experimental coccidioidomycois in mice. BMC Microbiol. 2011, 11: 71-10.1186/1471-2180-11-71.PubMedPubMed CentralView Article
- Gaur M, Puri N, Manoharlal R, Rai V, Mukhopadhayay G, Choudhury D, Prasad R: MFS transportome of the human pathogenic yeast Candida albicans. BMC Genomics. 2008, 9: 579-10.1186/1471-2164-9-579.PubMedPubMed CentralView Article
- Bogan KL, Evans C, Belenky P, Song P, Burant CF, Kennedy R, Brenner C: Identification of Isn1 and Sdt1 as glucose- and vitamin-regulated nicotinamide mononucleotide and nicotinic acid mononucleotide [corrected] 5′-nucleotidases responsible for production of nicotinamide riboside and nicotinic acid riboside. J Biol Chem. 2009, 284 (50): 34861-34869. 10.1074/jbc.M109.056689.PubMedPubMed CentralView Article
- Fischer R, Timberlake WE: Aspergillus nidulans apsA (anucleate primary sterigmata) encodes a coiled-coil protein required for nuclear positioning and completion of asexual development. J Cell Biol. 1995, 128 (4): 485-498. 10.1083/jcb.128.4.485.PubMedView Article
- Mayer BJ, Baltimore D: Signalling through SH2 and SH3 domains. Trends Cell Biol. 1993, 3 (1): 8-13. 10.1016/0962-8924(93)90194-6.PubMedView Article
- Fraser JA, Stajich JE, Tarcha EJ, Cole GT, Inglis DO, Sil A, Heitman J: Evolution of the mating type locus: insights gained from the dimorphic primary fungal pathogens Histoplasma capsulatum, Coccidioides immitis, and Coccidioides posadasii. Eukaryot Cell. 2007, 6 (4): 622-629. 10.1128/EC.00018-07.PubMedPubMed CentralView Article
- Bowman BH, Taylor JW, White TJ: Molecular evolution of the fungi: human pathogens. Mol Biol Evol. 1992, 9 (5): 893-904.PubMed
- Menacho-Marquez M, Perez-Valle J, Arino J, Gadea J, Murguia JR: Gcn2p regulates a G1/S cell cycle checkpoint in response to DNA damage. Cell Cycle. 2007, 6 (18): 2302-2305. 10.4161/cc.6.18.4668.PubMedView Article
- Siebel CW, Feng L, Guthrie C, Fu XD: Conservation in budding yeast of a kinase specific for SR splicing factors. Proc Natl Acad Sci USA. 1999, 96 (10): 5440-5445. 10.1073/pnas.96.10.5440.PubMedPubMed CentralView Article
- Stajich JE, Wilke SK, Ahren D, Au CH, Birren BW, Borodovsky M, Burns C, Canback B, Casselton LA, Cheng CK, et al: Insights into evolution of multicellular fungi from the assembled chromosomes of the mushroom Coprinopsis cinerea (Coprinus cinereus). Proc Natl Acad Sci USA. 2010, 107 (26): 11889-11894. 10.1073/pnas.1003391107.PubMedPubMed CentralView Article
- Desjardins CA, Champion MD, Holder JW, Muszewska A, Goldberg J, Bailao AM, Brigido MM, Ferreira ME, Garcia AM, Grynberg M, et al: Comparative genomic analysis of human fungal pathogens causing paracoccidioidomycosis. PLoS Genet. 2011, 7 (10): e1002345-10.1371/journal.pgen.1002345.PubMedPubMed CentralView Article
- Krishna SS, Majumdar I, Grishin NV: Structural classification of zinc fingers: survey and summary. Nucleic Acids Res. 2003, 31 (2): 532-550. 10.1093/nar/gkg161.PubMedPubMed CentralView Article
- Christy B, Nathans D: DNA binding site of the growth factor-inducible protein Zif268. Proc Natl Acad Sci USA. 1989, 86 (22): 8737-8741. 10.1073/pnas.86.22.8737.PubMedPubMed CentralView Article
- Vallim MA, Miller KY, Miller BL: Aspergillus SteA (sterile12-like) is a homeodomain-C2/H2-Zn + 2 finger transcription factor required for sexual reproduction. Mol Microbiol. 2000, 36 (2): 290-301. 10.1046/j.1365-2958.2000.01874.x.PubMedView Article
- Roberts RL, Fink GR: Elements of a single MAP kinase cascade in Saccharomyces cerevisiae mediate two developmental programs in the same cell type: mating and invasive growth. Genes Dev. 1994, 8 (24): 2974-2985. 10.1101/gad.8.24.2974.PubMedView Article
- Wittenberg C, La Valle R: Cell-cycle-regulatory elements and the control of cell differentiation in the budding yeast. Bioessays. 2003, 25 (9): 856-867. 10.1002/bies.10327.PubMedView Article
- La Valle R, Wittenberg C: A role for the Swe1 checkpoint kinase during filamentous growth of Saccharomyces cerevisiae. Genetics. 2001, 158 (2): 549-562.PubMedPubMed Central
- Konig C, Maekawa H, Schiebel E: Mutual regulation of cyclin-dependent kinase and the mitotic exit network. J Cell Biol. 2010, 188 (3): 351-368. 10.1083/jcb.200911128.PubMedPubMed CentralView Article
- Nurse P: Genetic control of cell size at cell division in yeast. Nature. 1975, 256 (5518): 547-551. 10.1038/256547a0.PubMedView Article
- Soustre I, Letourneux Y, Karst F: Characterization of the Saccharomyces cerevisiae RTA1 gene involved in 7-aminocholesterol resistance. Curr Genet. 1996, 30 (2): 121-125. 10.1007/s002940050110.PubMedView Article
- Brown SM, Campbell LT, Lodge JK: Cryptococcus neoformans, a fungus under stress. Curr Opin Microbiol. 2007, 10 (4): 320-325. 10.1016/j.mib.2007.05.014.PubMedPubMed CentralView Article
- Rodriguez-Caban J, Gonzalez-Velazquez W, Perez-Sanchez L, Gonzalez-Mendez R, Valle NR: Calcium/calmodulin kinase1 and its relation to thermotolerance and HSP90 in Sporothrix schenckii: an RNAi and yeast two-hybrid study. BMC Microbiol. 2011, 11: 162-10.1186/1471-2180-11-162.PubMedPubMed CentralView Article
- Daher BS, Venancio EJ, de Freitas SM, Bao SN, Vianney PV, Andrade RV, Dantas AS, Soares CM, Silva-Pereira I, Felipe MS: The highly expressed yeast gene pby20 from Paracoccidioides brasiliensis encodes a flavodoxin-like protein. Fungal Genet Biol. 2005, 42 (5): 434-443. 10.1016/j.fgb.2005.01.007.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedPubMed CentralView Article
- Semighini CP, Harris SD: Regulation of apical dominance in Aspergillus nidulans hyphae by reactive oxygen species. Genetics. 2008, 179 (4): 1919-1932. 10.1534/genetics.108.089318.PubMedPubMed CentralView Article
- Scott B, Eaton CJ: Role of reactive oxygen species in fungal cellular differentiations. Curr Opin Microbiol. 2008, 11 (6): 488-493. 10.1016/j.mib.2008.10.008.PubMedView Article
- Krishnan A, Almen MS, Fredriksson R, Schioth HB: The origin of GPCRs: identification of mammalian like Rhodopsin, adhesion, glutamate and frizzled GPCRs in fungi. PLoS One. 2012, 7 (1): e29817-10.1371/journal.pone.0029817.PubMedPubMed CentralView Article
- Gladfelter AS: Guides to the final frontier of the cytoskeleton: septins in the filamentous fungi. Curr Opin Microbiol. 2010, 13 (6): 720-726. 10.1016/j.mib.2010.09.012.PubMedView Article
- Bowman BH, White TJ, Taylor JW: Human pathogeneic fungi and their close nonpathogenic relatives. Mol Phylogenet Evol. 1996, 6 (1): 89-96. 10.1006/mpev.1996.0061.PubMedView Article
- 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. Epub 2006 Oct 13PubMedView Article
- Viriyakosol S, Fierer J, Brown GD, Kirkland TN: Innate immunity to the pathogenic fungus Coccidioides posadasii is dependent on TLR2 and dectin-1. Infect Immun. 2005, 73 (3):
- 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 USA. 2008, 105 (38): 14573-14578. 10.1073/pnas.0806221105.PubMedPubMed CentralView Article
- Nemecek JC, Wuthrich M, Klein BS: Global control of dimorphism and virulence in fungi. Science. 2006, 312 (5773): 583-588. 10.1126/science.1124105.PubMedView Article
- Nunes LR, Costa De Oliveira R, Leite DB, Da Silva VS, Dos Reis Marques E, Da Silva Ferreira ME, Ribeiro DC, De Souza Bernardes LA, Goldman MH, Puccia R, et al: Transcriptome analysis of Paracoccidioides brasiliensis cells undergoing mycelium-to-yeast transition. Eukaryot Cell. 2005, 4 (12): 2115-2128. 10.1128/EC.4.12.2115-2128.2005.PubMedPubMed CentralView Article
- Monteiro JP, Clemons KV, Mirels LF, Coller JA, Wu TD, Shankar J, Lopes CR, Stevens DA: Genomic DNA microarray comparison of gene expression patterns in Paracoccidioides brasiliensis mycelia and yeasts in vitro. Microbiology. 2009, 155 (Pt 8): 2795-2808.PubMedPubMed CentralView Article
- Moran GR: 4-Hydroxyphenylpyruvate dioxygenase. Arch Biochem Biophys. 2005, 433 (1): 117-128. 10.1016/j.abb.2004.08.015.PubMedView Article
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.