The WOPR family protein Ryp1 is a key regulator of gene expression, development, and virulence in the thermally dimorphic fungal pathogen Coccidioides posadasii

M. Alejandra Mandel; Sinem Beyhan; Mark Voorhies; Lisa F. Shubitz; John N. Galgiani; Marc J. Orbach; Anita Sil

doi:10.1371/journal.ppat.1009832

Abstract

Coccidioides spp. are mammalian fungal pathogens endemic to the Southwestern US and other desert regions of Mexico, Central and South America, with the bulk of US infections occurring in California and Arizona. In the soil, Coccidioides grows in a hyphal form that differentiates into 3–5 micron asexual spores (arthroconidia). When arthroconidia are inhaled by mammals they undergo a unique developmental transition from polar hyphal growth to isotropic expansion with multiple rounds of nuclear division, prior to segmentation, forming large spherules filled with endospores. Very little is understood about the molecular basis of spherule formation. Here we characterize the role of the conserved transcription factor Ryp1 in Coccidioides development. We show that Coccidioides Δryp1 mutants have altered colony morphology under hypha-promoting conditions and are unable to form mature spherules under spherule-promoting conditions. We analyze the transcriptional profile of wild-type and Δryp1 mutant cells under hypha- and spherule-promoting conditions, thereby defining a set of hypha- or spherule-enriched transcripts (“morphology-regulated” genes) that are dependent on Ryp1 for their expression. Forty percent of morphology-regulated expression is Ryp1-dependent, indicating that Ryp1 plays a dual role in both hyphal and spherule development. Ryp1-dependent transcripts include key virulence factors such as SOWgp, which encodes the spherule outer wall glycoprotein. Concordant with its role in spherule development, we find that the Δryp1 mutant is completely avirulent in the mouse model of coccidioidomycosis, indicating that Ryp1-dependent pathways are essential for the ability of Coccidioides to cause disease. Vaccination of C57BL/6 mice with live Δryp1 spores does not provide any protection from lethal C. posadasii intranasal infection, consistent with our findings that the Δryp1 mutant fails to make mature spherules and likely does not express key antigens required for effective vaccination. Taken together, this work identifies the first transcription factor that drives mature spherulation and virulence in Coccidioides.

Author summary

Coccidioides species, C. immitis and C. posadasii, are dimorphic fungal pathogens that commonly infect humans in North, Central, and South America, causing the respiratory fungal disease known as Valley Fever. Coccidioides grows as hyphae in the soil and differentiates into unique structures called spherules in the mammalian host. Spherules expand and internally divide to form endospores, which are released to facilitate dissemination of the pathogen within the host. The mechanisms underlying spherule differentiation remain largely unknown. In this study, we generated knockout mutants (Δryp1) of the conserved transcription factor Ryp1 in C. posadasii and characterized its role in spherule formation and virulence. We found that Ryp1 is required for the formation of mature spherules and colonization of mouse lungs. Transcriptional profiling of the Δryp1 mutant and the wild-type strain shows that Ryp1 regulates the expression of a subset of the transcripts that are either upregulated in wild-type spherules or in wild-type hyphae. These findings suggest that Ryp1 has a dual role in regulating morphology and virulence under host conditions as well as regulating genes involved in hyphal growth in the environment.

Citation: Mandel MA, Beyhan S, Voorhies M, Shubitz LF, Galgiani JN, Orbach MJ, et al. (2022) The WOPR family protein Ryp1 is a key regulator of gene expression, development, and virulence in the thermally dimorphic fungal pathogen Coccidioides posadasii. PLoS Pathog 18(4): e1009832. https://doi.org/10.1371/journal.ppat.1009832

Editor: Alex Andrianopoulos, University of Melbourne, AUSTRALIA

Received: July 22, 2021; Accepted: March 1, 2022; Published: April 6, 2022

Copyright: © 2022 Mandel et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are contained within the paper and/or Supporting Information files. For high-throughput sequencing data, the raw data are available at the NCBI Gene Expression Omnibus (GEO) databases under GEO accession GSE178277.

Funding: This work was supported by NIH R01AI137418 and 1U19AI166059 to S.B.; 2R37AI066224, 5R01AI146584, VFR-19-633952, R21AI144855, and MRP-17-454959 and 1U19AI166798 to AS; U54AI65359, and 1R01AI132140 to MJO, LSF and JNG; 1U01AI122275, to LSF and JNG. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Coccidioides spp. are fungal pathogens endemic to California, Arizona, and other desert regions in the Americas [1]. Coccidioides infects otherwise healthy individuals when they inhale spores from the soil. The prevalence of Coccidioides infection rose 9-fold between 1998 and 2011 and has continued to spike, with incidence now at an all-time high [2]. Elegant observational studies of Coccidioides have established its complex life cycle, its disease etiology, and its interaction with the mammalian immune system. The molecular basis of these attributes remains poorly understood, although candidate-based studies have implicated a handful of Coccidioides genes in growth and virulence.

Coccidioides is one of a group of thermally dimorphic fungal pathogens that grow in a sporulating hyphal form in the soil [3]. Upon introduction into the mammalian host, spores undergo differentiation into a parasitic form. Indeed, the defining characteristic of Coccidioides pathogenesis is development of the spore (arthroconidium) into a multicellular structure called the spherule [1,4]. In spherule development, instead of germinating and growing as polar hyphae, the arthroconidium undergoes isotropic enlargement in the mammalian lung with multiple rounds of mitosis prior to formation of internal spores (endospores) resulting in a 60 to 100 μm micron multi-cellular spherule, encased by a spherule outer wall [5]. If a spherule ruptures, endospores are released and can re-initiate the spherule cycle, either locally or following dissemination to other sites within the host. Known virulence genes in Coccidioides are expressed concomitant with spherule development [6–9]. One critical unanswered question in Coccidioides biology is the nature of the regulatory network that drives spherulation and virulence.

To identify putative regulators of spherulation, we assessed the role of the conserved fungal transcription factor, Ryp1, which is involved in development of the parasitic form in response to temperature for the other thermally dimorphic fungi. We had previously shown that the WOPR transcription factor Ryp1 is a master regulator that is required for the formation of the host form in the thermally dimorphic fungal pathogen Histoplasma [10,11]. Additionally, Ryp1 is a master regulator of gene expression in response to temperature in Histoplasma, as it is required for the vast majority of the gene expression program at 37°C [10,11]. Investigating the role of Ryp1 in Coccidioides was particularly compelling since Ryp1 orthologs are found throughout the fungal kingdom and are required for key developmental transitions in numerous fungi [12–16]. For example, the Candida albicans ortholog Wor1 regulates cell-type specification by driving the switch from “white” cells to “opaque” cells [17,18].

Here we deleted the RYP1 gene in Coccidioides posadasii and showed that the resultant mutant is unable to undergo mature spherulation. We used RNA-seq to show that Ryp1 has a role in gene expression in both spherules and hyphae. The Δryp1 mutant failed to express the normal complement of hypha-enriched transcripts, and most notably, the immature Δryp1 spherules expressed only a subset of normal spherule-enriched genes, consistent with their inability to fully differentiate. The Δryp1 mutant was completely avirulent in the mouse model of Coccidioides infection, demonstrating that the Ryp1 transcription factor is a key regulator that links the ability to form mature spherules and the expression of critical virulence traits.

Results

Deletion of Coccidioides posadasii RYP1

To define potential regulators of Coccidioides parasitic phase development, we identified the C. posadasii ortholog of the conserved regulator of fungal development, RYP1, based on similarity to the H. capsulatum Ryp1 protein [10]. Ryp1 is a member of the WOPR family of transcription factors. The C. posadasii strain Silveira RYP1 gene, CPSG_00528, encodes a 404 amino acid protein that is 68% similar to the 487 amino acid H. capsulatum RYP1 protein [10], with 84% identity and 91% similarity over the N-terminal half of the protein including the WOPRa and WOPRb regions and the putative nuclear localization signal. To determine the role of RYP1 in the C. posadasii life cycle, gene replacement strains were made via Agrobacterium transformation [19,20]. Transformed lines were screened for homologous recombination of a deletion construct where the E. coli hphB gene was cloned between 1.2 kb 5’ and 3’ flanking regions of the RYP1 gene, using an approach similar to that used for creation of the Δcps1 strain [19]. A total of 40 transformed lines were generated between two transformations with 42% being RYP1 gene replacements as analyzed by DNA hybridization for the first set of transformants and PCR for the second set. The RYP1 gene replacement strains (Δryp1) have a single insertion of the 1.4 kb E. coli hphB expression cassette replacing the full 1.3 kb RYP1 coding region. Phenotypic analyses of Δryp1 mutants were performed using several independent transformed strains and compared to the parental Silveira strain (WT) and one strain (1563.1) where the RYP1 deletion construct had integrated randomly in the genome leaving the RYP1 gene intact.

Ryp1 is required for Coccidioides development

Deletion of RYP1 resulted in several in vitro phenotypes (Fig 1). The Δryp1 arthroconidia were delayed in germination, with visible colonies appearing after 7 days at 25°C on 2X GYE media while the progenitor WT and the ectopic transformed strains gave visible colonies three days after plating. The radial growth of Δryp1 colonies was also retarded relative to WT or ectopic transformed strains (Fig 1B) and Δryp1 colonies failed to grow to the edge of petri dishes (Fig 1A). The most dramatic phenotype was observed under in vitro spherulation conditions: the Δryp1 mutant failed to differentiate into mature spherules. Instead, the arthroconidia remained as thin-walled structures at all time points and by 72 hours had formed short hyphal-like segmented structures (Fig 2A and 2B). In addition, although Δryp1 mutants had reduced spore viability compared to WT and ectopic transformed strains, their viability was independent of temperature (Fig 2C).

Download:

Fig 1. Growth morphology of Δryp1 mutants.

(A) C. posadasii Δryp1 strains (1563.4 and 1563.7) exhibit a growth-limiting phenotype, in comparison to WT and an ectopic transformant. (B) Plugs of 6 mm were transferred from freshly growing plates to 2X GYE plates and incubated at room temperature, to measure colony radial growth. Colony diameters were measured on days 8, 14 and 20. Four independent Δryp1 strains (1563.4, 1563.6, 1563.7 and 1563.14), as well as the WT parental strain and an ectopic transformed strain (1563.1) were tested.

https://doi.org/10.1371/journal.ppat.1009832.g001

Download:

Fig 2. C. posadasii Δryp1 does not produce normal spherules in vitro.

(A, B) Spherules of WT, a Δryp1 mutant and an ectopic transformed strain, grown in Converse medium at 37°C with 20% CO₂ and aliquots analyzed during in vitro spherulation. (A) Images at 24 h show that while the WT and ectopic strains are rounding up to form early spherules, Δryp1 does not form round structures. At 96 h, the WT and ectopic strains produced mature spherules, while the Δryp1 exhibited limited polar hyphal growth. (B) Images stained with lactophenol blue at 120 h show mature spherules produced by WT and polar hyphal growth by Δryp1. Black bars represent 10 μm. (C) Percentage of spore viabilities for WT and ectopic strains, and Δryp1 mutants (1563.7 and 1563.10) were determined based on the growth on 2X GYE plates at 37°C and 25°C.

https://doi.org/10.1371/journal.ppat.1009832.g002

Ryp1 regulates distinct sets of genes in spherules and hyphae

Our previous studies revealed that H. capsulatum Ryp1 is required for yeast-phase growth and is responsible for the expression of the majority of yeast-phase specific genes [10]. Given that C. posadasii Δryp1 mutants are defective in spherulation, we postulated that C. posadasii Ryp1 is also responsible for transcriptional regulation of morphology-associated genes. To this end, we performed RNA-seq to transcriptionally profile WT and Δryp1 strains of C. posadasii grown under either spherulation or hyphal conditions. Under our experimental conditions, 9711 transcripts were observable (transcripts per million, TPM, ≥ 10) in at least one sample. Of these, 4175 transcripts (about 43% of all detected transcripts) were significantly differentially expressed (Fig 3) with at least a two-fold change in at least one of the following three comparisons: (1) a comparison of WT spherules and WT hyphae (WT^Sph/WT^Hy) identified 2985 differentially expressed transcripts that we refer to as “morphology-regulated” genes; (2) a comparison of WT vs Δryp1 under spherule conditions identified 1742 Ryp1-dependent transcripts; and (3) a comparison of WT vs. Δryp1 under hyphal conditions identified 925 Ryp1-dependent transcripts (Fig 3B). Not all of these Ryp1-dependent transcripts in comparisons 2 and 3 are morphology-regulated genes, and in fact, 1157 morphology-regulated transcripts (39% of the 2985 WT^Sph/WT^Hy comparison) were dependent on Ryp1 either under spherulation or hyphal conditions. It is likely that some of the morphology-regulated genes are differentially expressed due to the differences in spherule and hyphae growth conditions (differences in temperature, CO₂ levels and growth medium). To further explore the 4175 differentially expressed transcripts, we grouped them in 17 classes by expression pattern (Fig 3A and S1 Table). Ryp1-dependent transcription under spherulation conditions is correlated with WT morphology-regulated transcription: spherule-enriched genes are more likely to be Ryp1-induced (class 2 in Fig 3A) while hypha-enriched genes are more likely to be Ryp1-repressed in spherules (Fig 3C and class 4 in Fig 3A). In contrast, Ryp1-dependent transcription under hyphal conditions is not correlated with morphology-regulated transcription, as can be seen from the shape of the scatter plot (Fig 3D).

Download:

Fig 3. Ryp1 regulates distinct sets of genes in spherules and hyphae.

(A) Heatmap showing differential expression of the 4175 genes described in the text. Genes are grouped by expression pattern (significantly up, down, or neutral in each of the three contrasts) with different classes labeled in the color bar to the right of the heatmap. The position of genes of interest in the heatmap are indicated to the right of the color bar. (B) Venn diagram of genes observed to be morphology-regulated (differentially expressed in WT^Sph/WT^Hy), Ryp1-dependent in spherulation conditions (differentially expressed in WT^Sph/Δryp1^Sph), or Ryp1-dependent in hyphal conditions (differentially expressed in WT^Hy/Δryp1^Hy). (C and D) Scatter plots of the differentially expressed genes comparing the WT spherule-enrichment (right) or hyphal-enrichment (left) to ryp1-depletion (top) or ryp1-enrichment (bottom) in (C) spherules or (D) hyphae. Genes in the scatter plots are colored by expression pattern, as in the color bar of (A).

https://doi.org/10.1371/journal.ppat.1009832.g003

Ryp1 is required for the expression of a set of spherule-enriched transcripts

Early subtractive cDNA hybridization experiments identified four spherule-specific genes (ALD1, OPS1, PSP1, and MDR1) [21]. Of these, we found that ALD1 is spherule-enriched and Ryp1-independent (S1 Table, class 1) while the remaining three are both spherule-enriched and Ryp1-dependent at 37°C (S1 Table, class 2). Altogether, 406 (25%) of the 1593 transcripts that are more highly expressed in wild-type spherules versus hyphae are dependent on Ryp1 for their expression. This set includes SOWgp, the best characterized Coccidioides virulence factor [22], and CPSG_05265 (CIMG_00509), a previously noted spherule enriched gene [6] at the boundary of a region of introgression between C. immitis and C. posadasii [23].

Intriguingly, about 5% (84 of 1593) of the spherule-enriched transcripts are more highly expressed in the Δryp1 mutant than WT under spherulation conditions, (S1 Table, class 3), indicating their expression is normally repressed by Ryp1 in wild-type spherules. This set includes the sulfur assimilation pathway (MET3, MET14, MET16, and MET10) required to reduce intracellular sulfate to sulfide. In contrast, we found that the genes acting in opposition to sulfur assimilation (both the sulfite oxidase SOX1 and the sulfite exporter SSU1) are spherule-enriched and Ryp1-dependent. These data are consistent with Ryp1 regulation promoting sulfur efflux versus influx.

Ryp1 regulates the expression of a set of hypha-enriched transcripts

Analysis of the hypha-enriched transcripts showed that about 25% of this set (351 out of 1392) are Ryp1-repressed when cells are grown under spherulation conditions (S1 Table, class 4). Notably, this set includes the APSES transcription factor STU1, which regulates hyphal morphology in Histoplasma [24,25] and has conserved hyphal enrichment in C. immitis [8] and all of the thermally dimorphic pezizomycetes for which RNA-seq data is available [25]. Additional genes in this set consistent with hyphal growth include 5 septins (CDC10, CDC12, CDC3, ASPA, and ASPE), which are cytoskeletal proteins involved in cell cycle regulation [26]. Intriguingly, 33 of the 351 genes in this set are also Ryp1-dependent under hyphal growth conditions. In particular, the Coccidioides ortholog of the OefC transcription factor, which confers a “fluffy” phenotype in Aspergillus nidulans when overexpressed [27], exhibits this Ryp1-dependent regulation in hyphae. Taken together, our results suggest that Ryp1 has a dual role in Coccidioides development, regulating the expression of subsets of both spherule- and hypha-enriched genes.

Ryp1 regulates expression of morphology- or temperature-independent transcripts

In Histoplasma, Ryp1 directly interacts with two other transcriptional regulators, Ryp2 and Ryp3, to regulate yeast-phase growth [11]. The orthologs of RYP2 and RYP3 in Coccidioides are not differentially regulated in wild-type spherules and hyphae; however, RYP2 displays a Ryp1-dependent expression pattern under spherulation conditions (S1 Table, class 6). In addition, 282 transcripts show a similar morphology-independent but Ryp1-dependent expression pattern under spherulation conditions, including the cytosolic catalase CATP, the tyrosinase TYR2, as well as ~40 genes involved in primary carbohydrate, nucleotide, and amino acid metabolism. An additional 75 transcripts display a morphology-independent but Ryp1-dependent expression pattern under both spherule- and hypha-promoting conditions (S1 Table, class 7). Among these genes are MEP1, a protease that is thought to contribute to immune evasion via degradation of SOWgp [28]; urease (CPSG_08438), a known virulence factor of Coccidioides [29]; UAZ1, a uricase that catalyzes the first step in the breakdown of uric acid; and CTR2, a copper transporter observed to be spherule-enriched in C. immitis.

In addition to the Ryp1-dependent transcripts, there are 432 transcripts that are equivalently expressed in WT spherules and hyphae but are upregulated in Δryp1 mutants compared to wild-type under spherulation conditions (S1 Table, class 8), suggesting that Ryp1 represses their expression. This set includes LAE1, a master regulator of secondary metabolism [30]; orthologs of the Aspergillus developmental regulators NsdD (NSD4) and FlbC (FBC1); the MAPK pathway components STE11, PBS2, and SDP1; and the histidine kinase PhkB (PHK2). In Histoplasma, PHKB shares a Ryp1-associated divergent promoter with the histidine kinase PHKA, and both genes are Ryp1-induced in Histoplasma yeast [11]. In Coccidioides, these genes likewise share a divergent promoter, and have correlated expression profiles, but the differential expression of PHKA is too modest to pass our 2-fold change criterion.

RYP1 deletion mutants are avirulent

The critical role of Ryp1 in spherule formation and gene regulation suggested that it may be either reduced in virulence, or avirulent, in vivo. To test this, we assessed the virulence of the Δryp1 strain using the mouse model of coccidioidomycosis in C57BL/6 mice. Twelve mice were intranasally infected with 50 or 1,000 viable arthroconidia of Δryp1 strain 1563.7 and disease progression was compared to infection with 50 WT viable arthroconidia or 4,400 Δcps1 viable arthroconidia which were previously shown to be avirulent [19]. Fungal burden in the lungs and spleen as well as survival was monitored for 28 days post-infection. Additionally, two mice from each group were sacrificed after 14 days for histopathological studies. The WT-infected mice all became moribund and were euthanized between day 12 and 19, while all Δryp1 and Δcps1-infected mice survived to day 28. Lung and spleen homogenates did not yield Coccidioides colony-forming units (CFUs) from infections with the Δryp1 and Δcps1 mice, consistent with an inability of these mutant strains to survive in vivo. The histopathology did not reveal Coccidioides organisms in any of the mice. Two animals displayed one mild, focal, inflammatory lesion with no observable fungal cells. In a follow-up experiment, four additional Δryp1 strains were screened for virulence along with an ectopic transformed line (Fig 4). As in the first experiment, Δryp1-infected mice, which received between 525 and 1133 viable arthroconidia, survived for the full duration of the 28-day study (Fig 4). In contrast, seven of eight mice infected with 52 viable spores of the ectopically integrated RYP1 deletion construct or infected with 49 WT viable spores died by day 17. One mouse in each of these two groups survived for 28 days but had significant pulmonary disease and dissemination and would have ultimately succumbed to the infection. As we observed previously, there was no growth of Coccidioides from the lungs or spleens of the Δryp1-infected mice. The mean lung CFUs of the mice infected with ectopic strain 1563.1 were 2.87 x 10⁷ (range 1.24 x 10⁵–4.7 x 10⁷), and for those infected with WT were 1.31 x 10⁷ (range 1.06 x 10⁵–4.8 x 10⁷), with dissemination to the spleens in all mice as reflected by spleen CFUs.

Download:

Fig 4. Infection of mice with C. posadasii Δryp1.

Survival results for C57BL/6 mice (N = 8 mice/group) infected intranasally with four independent Δryp1 strains, 1563.7 (circle), 1563.10 (square), 1563.14 (inverted triangle) and 1563.15 (triangle); one strain with the ryp1 deletion cassette integrated ectopically into WT, 1563.1 (half-filled circle); and WT (half-filled square). The Δryp1 strains were infected at between 525 and 1133 viable arthroconidia, while 50 viable arthroconidia of the ectopic stain and WT were used.

https://doi.org/10.1371/journal.ppat.1009832.g004

Δryp1 mutants do not provide protection against WT infection

Due to the failure of Δryp1 to cause disease or persist in the lungs of C57BL/6 mice, the Δryp1 strain was tested to determine whether it provides protective immunity against subsequent WT infection, as is the case for the Δcps1 strain. The Δcps1 vaccine strain and the Δryp1 strain 1563.7 were used to vaccinate C57BL/6 mice either intraperitoneally (IP) or subcutaneously (SC) as described previously [19]. Mice were vaccinated, boosted after two weeks, challenged four weeks later with 90 WT arthroconidia, and sacrificed two weeks later. While Δcps1 provided dramatic protection by either SC or IP vaccination, resulting in a 5-log lower lung fungal burden than mice vaccinated with a control adjuvant, vaccination with Δryp1 resulted in no reduction in lung fungal burden and thus, no protection (Fig 5). Mice in the control group had a mean lung fungal burden of 5.3 x 10⁶ CFU, which was similar to that of the Δryp1 SC-vaccinated mice (P<0.05) which had mean lung fungal burdens of 4.1 x 10⁶ CFU and the Δryp1 IP-vaccinated mice with an average of 6.4 x 10⁶ CFU per lung. In contrast, the Δcps1-vaccinated mice had mean lung fungal burdens of 341 CFU for the IP-vaccinated group and 84 CFU for the SC-vaccinated group. Furthermore, when whole spleens were plated for fungal growth, those from Δcps1-vaccinated mice were almost fully free of fungi, with only one spleen of a single Δcps1 IP-vaccinated mouse having a small area of growth, likely from a single spherule, appearing seven days after plating, while the control mice and Δryp1-vaccinated mice all grew Coccidioides by three days after plating. Thus it is clear that the Δryp1 strain does not provide protective immunity against Coccidioides infection.

Download:

Fig 5. Vaccination of mice with C. posadasii Δryp1 mutant.

Protection results for C57BL/6 mice vaccinated with either Δcps1, Δryp1 or a S. cerevisiae culture supernatant (Control). C57BL/6 mice (N = 8 mice/group) were injected either by intraperitoneal (IP) or subcutaneous (SC) injection twice 2 weeks apart with Δcps1, Δryp1, or Control. Vaccinated mice were challenged with a lethal dose of WT (90 spores) and total lung CFUs were determined 14 days post-infection. Δcps1 significantly protects mice from C. posadasii infection compared to Control (P<0.004 both comparisons) while Δryp1 shows no protection (P = 0.99 vs. control). Bar equals geometric mean.

https://doi.org/10.1371/journal.ppat.1009832.g005

Discussion

Spherulation is a critical step of pathogenesis in Coccidioides infection, allowing the fungus to adapt to growth within the host and producing endospores which allow expansion and dissemination of the fungus. Here we identified the first transcription factor that is essential for the process of Coccidioides spherulation by interrogating the function of Ryp1, a conserved transcriptional regulator that plays a critical role in fungal development and virulence in many fungi, including Histoplasma, a close relative of Coccidioides. Our results showed that as in Histoplasma, Ryp1 is required for development of the host parasitic form; namely spherules in Coccidioides and yeast-phase cells in Histoplasma. Additionally, in Coccidioides, we demonstrated that Ryp1 is absolutely required for virulence in the mouse model of infection, consistent with its critical role in spherulation. The Δryp1 strain failed to differentiate into spherules in vitro and failed to persist in mice, producing no symptoms of disease or viable propagules at 28 days post-infection, even when exposed to 10-20X the level of spores that results in a lethal infection for WT strains. Interestingly, in contrast to Histoplasma, where Ryp1 is required for the differential expression of the vast majority of yeast-enriched transcripts, there were a significant number of spherule-enriched genes whose expression was Ryp1-independent, indicating that other transcription factors likely function in parallel to Ryp1 to regulate spherule development. Additionally, Coccidioides Ryp1 is required for normal rates of radial growth of hyphal colonies and regulates the expression of a subset of hypha-enriched transcripts. Thus Coccidioides Ryp1 has a dual role in both the environmental form and the host form of the fungus where it regulates gene expression and development.

The contrast between Histoplasma and Coccidioides Ryp1 should be interpreted in the light of the evolutionary relationship between Histoplasma and Coccidioides—although both are thermally dimorphic fungal pathogens in the order Onygenales, Histoplasma is a member of the family Ajellomycetaceae, which also includes the thermally dimorphic pathogens Paracoccidioides, Blastomyces, and Emergomyces, consistent with a single origin of dimorphism in this family. In contrast, Coccidioides is the only known dimorphic member of the family Onygenaceae and has a unique dimorphic pattern, forming endosporulating spherules during host invasion vs. the more typical hyphal to yeast phase transition. This evolutionary relationship suggests that dimorphism of Histoplasma and Coccidioides evolved separately, and that some regulatory mechanisms may be shared whereas others are distinct. In Histoplasma, since Ryp1 is required for the vast majority of yeast-enriched expression, identifying Ryp-dependent genes by expression analysis of WT vs ryp1 mutants did not further refine the set of yeast-enriched transcripts defined by comparing the gene expression program of wild-type Histoplasma yeast and hyphae [10]. In contrast, only 25% of spherule-enriched transcripts were dependent on Coccidioides Ryp1, thereby suggesting that these genes may be key to spherule development and virulence, both of which are defective in the absence of Ryp1. Further identification of key effectors of Coccidioides spherulation and virulence could come from future studies identifying direct targets of Ryp1, which was useful in identifying virulence factors in Histoplasma [11].

Transcriptome differences between Coccidioides spherules and hyphae have previously been profiled using RNA-seq in both C. posadasii and C. immitis [6,8,31]. In our dataset, about 30% of the detectable transcripts were differentially regulated between C. posadasii spherules and hyphae. When analyzed by our pipeline and significance criteria, 31% of our spherule-enriched transcripts were spherule enriched in the young spherule data of Carlin et al [8] and 48% in the C. posadasii data of Whiston et al [6]. Likewise, 59% of our hypha-enriched transcripts were hypha enriched in Carlin et al and 54% in Whiston et al. While there is not a complete overlap with the other studies, likely due to differences in growth conditions and stage of spherule development, most of the previously observed highly regulated genes are also differentially regulated in our dataset. For example, SOWgp, the best characterized virulence factor of Coccidioides [22], is spherule-enriched in all datasets, and we observe its expression to be Ryp1-dependent. Among the shared hypha-enriched genes in all datasets, STU1, an APSES transcription factor that is a regulator of hyphal growth in Histoplasma [24], is Ryp1-repressed in spherules, underscoring the importance of Ryp1 for the expression of morphology-specific genes.

In this study, we also explored the role of Ryp1 in pathogenesis using the murine model of coccidioidomycosis. Our results show that Ryp1 is essential for in vivo growth and colonization; however, it does not confer protection to the host, unlike the previously characterized Δcps1 mutant. Intriguingly, the Δcps1 mutant is able to convert to spherules but is impaired in development of mature spherules and endospores [19], whereas the morphology of the Δryp1 mutant in vivo is not yet known. Given its morphology defect in vitro, we predict that the Δryp1 mutant may not be able to convert to spherules in vivo, and thus may not display spherule-specific antigens that are required to elicit a protective host response. Consequently, comparison of the Δryp1 and Δcps1 mutants may narrow the search for the molecules responsible for the protective effect of the Δcps1 mutant. Furthermore, while the Δryp1 mutant may not be a good vaccine candidate, Ryp1 and other proteins that are required for spherule growth may be effective targets for the development of therapeutics that inhibit the transition from arthroconidia to spherules in vivo. Additionally, elucidation of the Ryp1 regulon as described here identifies key downstream effectors that are essential for the development and virulence of the host form of Coccidioides, and thus also serve as valuable diagnostic and therapeutic targets.

Materials and methods

Ethics statement

All murine infections and handling were carried out at animal BSL3 facility and all procedures with mice were approved by the University of Arizona Institutional Animal Care and Use Committee (Protocol number 14–562). Animals were housed and cared for according to PHS standards. Mice were euthanized by overdose of inhalant anesthetics (ketamine (80 mg/kg)/xylazine (8 mg/kg)) in a rodent anesthesia chamber.

Strains and growth conditions

Wild-type Coccidioides posadasii (strain Silveira, ATCC 28868) was cultured on 2X GYE medium (2% glucose, 1% yeast extract, and 1.5% agar) at room temperature (approximately 24°C). Mutant strains were selected and maintained on 2X GYE media supplemented with 50 mg/ml of hygromycin, also at room temperature as described previously [19]. Arthroconidia of WT and mutant strains were harvested from 4 week-old cultures, using sterile water by the mini-stir bar method described previously [32] and stored in sterile water at 4°C. Spore numbers were determined with a hemocytometer and viable counts by serial dilution and plating. All manipulations of viable cultures were performed at biosafety level 3 (BSL3). Spherules were prepared by growth in a modified Converse medium at 38°C with 20% CO₂ [33,34], with shaking at 180 rpm for 48 hours. To measure colonial radial growth and assess colony morphology, 6 mm plugs of actively growing cultures were transferred to 2X GYE agar plates using transfertubes (Transfertube Disposable Harvesters, Spectrum Laboratories,) and grown at room temperature.

Creation of RYP1 gene deletion mutants

A gene replacement cassette for the C. posadasii strain Silveira RYP1 gene (CPSG_00528) was constructed using double-joint PCR [35]. DNA fragments of 1.2kb (5’ 1192 and 3’ 1233) representing the 5’ and 3’ flanking sequences of CPSG_00528 were amplified by PCR with oligonucleotide primer combinations OAM1153/OAM1154, and OAM1155/OAM1156, respectively (S3 Table). The hygromycin resistance gene (E. coli hphB) was amplified with primers OAM597/OAM598. These three amplicons were then recombined by double-joint PCR as described [35]. The extension product was amplified by PCR with the nested primers OAM1159/OAM1160 which contain BamHI restriction sites added at the ends. This product was cloned into pGEM-T Easy (Promega) and the insert in plasmid AM1538 was sequenced to verify the construct. The gene replacement cassette was cloned as a BamHI fragment into AM1145, the binary vector for Agrobacterium tumefaciens transformation [36], creating AM1563. AM1563 was transformed into A. tumefaciens strain AD965. C. posadasii strain Silveira was transformed with AD965/AM1563 as previously described [36] and hygromycin resistant strains selected and passaged three times to get homokaryons, prior to molecular analysis. Gene replacement strains and ectopic insertion transformants were identified by Southern hybridization analysis as described previously [36], or by PCR. For PCR analysis, primers OAM597 and OAM598 were used to confirm the presence of the hphB gene in all transformants. Primers OAM1611 located upstream of the RYP1 locus and OAM431 of the hphB gene were used to detect a homologous integration event of the transforming construct at the 5’ end of RYP1, and primers OAM1612 located downstream of the RYP1 locus along with OAM432 were used to detect homologous integration of the transforming construct at the 3’ end of RYP1. For initial experiments including in vitro growth, Δryp1 strains 1563.1, 1563.4, 1563.7 and 1563.14 were used along with ectopic transformed line 1563.19, where the RYP1 deletion construct was integrated elsewhere in the genome. For mouse virulence and protection studies, Δryp1 strain 1563.7 was used. For further experiments to validate mouse virulence results, an additional set of strains were used; Δryp1 strains 1563.7, 1563.10, 1563.14 and 1563.15 and ectopic transformed strain 1563.1.

RNA-seq library preparation and sequencing

For expression studies, total RNAs were isolated from WT and Δryp1 hyphal and spherule cultures. Duplicate hyphal cultures of each strain were grown by inoculation of 5 x 10⁷ arthroconidia into 100 ml of 2X GYE and then incubated with shaking at 180 rpm for 48 h at 28°C. Spherule cultures of WT and Δryp1 were also grown in duplicate by inoculating 100 ml of modified Converse liquid medium with 3 x 10⁸ arthroconidia and growing them at 38°C with 20% CO₂, with shaking at 180 rpm for 48 h [37]. RNAs were isolated using a modified acid-phenol preparation as described previously [37] with the modification that liquid N₂ grinding was used to initially disrupt the cells. RNAs were resuspended in diethyl pyrocarbonate-treated dH₂O prior to assessment of their quality and concentration with an Agilent Bioanalyzer (Agilent Technologies, Palo Alto, CA).

RNA-seq libraries were prepared as previously described [25]. Briefly, for the isolation of mRNAs, five μg total RNA from each sample was subjected to poly-A purification using Dynabeads Oligo (dT)₂₅(Invitrogen-ThermoFisher). RNA libraries were prepared using NEBNext Ultra Directional RNA Library Prep Kit (NEB). Quality of the total RNA, mRNA and library was confirmed using Bioanalyzer (Agilent). Sequencing was done in-house at UCSF Center for Advanced Technology using Illumina HiSeq2500 instruments. 11 to 19 million single-end 50 bp reads were obtained for each sample.

RNA-seq data analysis

Predicted mRNA sequences for Coccidioides posadasii Silveira were downloaded from the Broad Institute on 3/18/2015 (http://www.broadinstitute.org/annotation/genome/coccidioides_group/MultiDownloads.html) [23,38].

Relative abundances (reported as TPM values [39]) and estimated counts (est_counts) of each transcript in each sample were estimated by alignment free comparison of k-mers between the preprocessed reads and mRNA sequences using KALLISTO version 0.46.0 [40]. Further analysis was restricted to transcripts with TPM ≥ 10 in at least one sample.

Differentially expressed genes were identified by comparing replicate means for contrasts of interest using LIMMA version 3.30.8 [41,42]. Genes were considered significantly differentially expressed if they were statistically significant (at 5% FDR) with an effect size of at least 2x (absolute log2 fold change ≥ 1) for a given contrast.

For comparison with previous expression profiling, reads were downloaded for SRA accessions SRP313022 [8] and SRP013936 [6] and analyzed with the above KALLISTO/LIMMA pipeline.

Murine virulence and protection studies

To assess the pathogenicity of Δryp1 mutants, 8-week old female C57BL/6 mice (B6) were anesthetized with ketamine (80 mg/kg)/xylazine (8 mg/kg) and infected intranasally (IN) as previously described [43]. Briefly, 8–10 mice per group were challenged with 50 or 1000 viable spores of the Δryp1 mutant strains. As a control, B6 mice were given 50 WT spores, and in some studies, B6 mice were infected with 50 viable spores of a transformed line with an ectopic insertion of the RYP1 gene deletion construct.

The potential for Δryp1 strains to protect mice from subsequent infection by WT was assessed by performing vaccination experiments as described for Δcps1 [19]. Six-week-old B6 mice were primed with 2.5 x 10⁴ Δryp1 (strain 1563.7) arthroconidia either by intraperitoneal (IP) or subcutaneous (SC) injection of spores in groups of eight mice and two weeks later boosted with the same number of spores. Two other groups of mice were used as a positive control for vaccination by being vaccinated in a similar manner with Δcps1 arthroconidia, at a dose of 5 x 10⁴ spores. As a control, mice were vaccinated with a S. cerevisiae culture supernatant, which was previously used as an adjuvant for a chimeric protein antigen that shows some vaccine protection for Coccidioides [44]. Mice were challenged four weeks later with 90 WT arthroconidia and sacrificed two weeks later for lung fungal burden determination. Limit of detection was < 10 CFU. Spleens were cultured whole to determine dissemination [44].

Other software and libraries

We wrote custom scripts and generated plots in Prism (GraphPad Software, San Diego, CA), R 3.3.3 [45] and Python 2.7.13, using Numpy 1.12.1 [46] and Matplotlib 2.0.0 [47]. Jupyter 4.2.3 [48] notebooks and JavaTreeView 1.1.6r4 [49] were used for interactive data exploration.

Supporting information

S1 Table. Table of Kallisto quantification, limma statistics, and annotations for differentially expressed genes.

Excel-compatible tab-delimited text conforming to JavaTreeView extended CDT format. Each row is a transcript with the UNIQID column giving the Broad Cp Silveira systematic gene name. The NAME column gives manually curated short gene names, transferred from Histoplasma capsulatum G217B (HcG217B, GenBank GCA_017607445.1) augmented with Coccidioides-specific gene names. The next three columns give limma BH-adjusted p-values for differential expression in each of the three contrasts. The next 8 columns give kallisto estimated counts for each sample. The next 8 columns give kallisto normalized abundances (TPMs) for each sample. Cp_GenBank and Cp_anno give the GenBank Cp Silveira accession and annotation respectively. CiRS, CiRS_GenBank, and Ci_anno give the Broad systematic gene name, GenBank accession, and annotation respectively for the InParanoid-mapped C. immitis RS ortholog. HcG217B, HcG217B_anno, and HcG217B_GenBank give the GSC systematic gene name and annotation (from [50]) and GenBank accession for the InParanoid-mapped HcG217B ortholog. Cp_new_GenBank gives the GenBank accession for the corresponding gene in the newly available Cp Silveira assembly and annotation (BioProject PRJNA664774) [51]. Class indicates the differential expression pattern, as referenced in the results and discussion sections. BGCOLOR gives the hex code for the class coloring in Fig 3. GWEIGHT is a place-holder column for JavaTreeView compatibility. The final three columns give the limma fit values for the three contrasts plotted in Fig 3.

https://doi.org/10.1371/journal.ppat.1009832.s001

(TXT)

S2 Table. Table of Kallisto quantification, limma statistics, and annotations for all expressed genes.

Excel-compatible tab-delimited text conforming to JavaTreeView extended CDT format. Columns are exactly as for S1 Table, but rows include all 9711 analyzed genes. The estimated counts in this file are sufficient to recapitulate the limma analysis.

https://doi.org/10.1371/journal.ppat.1009832.s002

(TXT)

S3 Table. Table of oligonucleotides used for Δryp1 mutant creation and analysis.

https://doi.org/10.1371/journal.ppat.1009832.s003

(DOCX)

Acknowledgments

We thank members of the Beyhan, Galgiani, Orbach, and Sil labs for helpful discussions and Cadence Flohrschutz for technical assistance with Coccidioides growth and sporulation studies.

References

1. Johnson L, Gaab EM, Sanchez J, Bui PQ, Nobile CJ, Hoyer KK, et al. Valley fever: danger lurking in a dust cloud. Microbes Infect. 2014;16: 591–600. pmid:25038397
- View Article
- PubMed/NCBI
- Google Scholar
2. Centers for Disease Control and Prevention (CDC). Increase in reported coccidioidomycosis—United States, 1998–2011. MMWR Morb Mortal Wkly Rep. 2013;62: 217–221. pmid:23535687
- View Article
- PubMed/NCBI
- Google Scholar
3. Sil A, Andrianopoulos A. Thermally Dimorphic Human Fungal Pathogens—Polyphyletic Pathogens with a Convergent Pathogenicity Trait. Cold Spring Harb Perspect Med. 2014;5: a019794. pmid:25384771
- View Article
- PubMed/NCBI
- Google Scholar
4. Cole GT, Sun SH. Arthroconidium-Spherule-Endospore Transformation in Coccidioides immitis. Fungal Dimorphism. 1985. pp. 281–333.
- View Article
- Google Scholar
5. Cole GT, Hung CY. The parasitic cell wall of Coccidioides immitis. Med Mycol. 2001;39 Suppl 1: 31–40.
- View Article
- Google Scholar
6. 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: e41034. pmid:22911737
- View Article
- PubMed/NCBI
- Google Scholar
7. Johannesson H, Kasuga T, Schaller RA, Good B, Gardner MJ, Townsend JP, et al. Phase-specific gene expression underlying morphological adaptations of the dimorphic human pathogenic fungus, Coccidioides posadasii. Fungal Genet Biol. 2006;43: 545–559. pmid:16697669
- View Article
- PubMed/NCBI
- Google Scholar
8. Carlin AF, Beyhan S, Peña JF, Stajich JE, Viriyakosol S, Fierer J, et al. Transcriptional Analysis of Coccidioides immitis Mycelia and Spherules by RNA Sequencing. Journal of Fungi. 2021. pmid:34067070
- View Article
- PubMed/NCBI
- Google Scholar
9. Viriyakosol S, Singhania A, Fierer J, Goldberg J, Kirkland TN, Woelk CH. Gene expression in human fungal pathogen Coccidioides immitis changes as arthroconidia differentiate into spherules and mature. BMC Microbiol. 2013;13: 1–1. pmid:23286760
- View Article
- PubMed/NCBI
- Google Scholar
10. Nguyen VQ, Sil A. Temperature-induced switch to the pathogenic yeast form of Histoplasma capsulatum requires Ryp1, a conserved transcriptional regulator. Proceedings of the National Academy of Sciences. 2008;105: 4880.
- View Article
- Google Scholar
11. Beyhan S, Gutierrez M, Voorhies M, Sil A. A Temperature-Responsive Network Links Cell Shape and Virulence Traits in a Primary Fungal Pathogen. PLoS Biol. 2013;11: e1001614. pmid:23935449
- View Article
- PubMed/NCBI
- Google Scholar
12. Lohse MB, Zordan RE, Cain CW, Johnson AD. Distinct class of DNA-binding domains is exemplified by a master regulator of phenotypic switching in Candida albicans. Proc Natl Acad Sci U S A. 2010;107: 14105–14110. pmid:20660774
- View Article
- PubMed/NCBI
- Google Scholar
13. Michielse CB, Becker M, Heller J, Moraga J, Collado IG, Tudzynski P. The Botrytis cinerea Reg1 protein, a putative transcriptional regulator, is required for pathogenicity, conidiogenesis, and the production of secondary metabolites. Mol Plant Microbe Interact. 2011;24: 1074–1085. pmid:21635139
- View Article
- PubMed/NCBI
- Google Scholar
14. Jonkers W, Dong Y, Broz K, Kistler HC. The Wor1-like protein Fgp1 regulates pathogenicity, toxin synthesis and reproduction in the phytopathogenic fungus Fusarium graminearum. PLoS Pathog. 2012;8: e1002724. pmid:22693448
- View Article
- PubMed/NCBI
- Google Scholar
15. Brown DW, Busman M, Proctor RH. Fusarium verticillioides SGE1 is required for full virulence and regulates expression of protein effector and secondary metabolite biosynthetic genes. Mol Plant Microbe Interact. 2014;27: 809–823. pmid:24742071
- View Article
- PubMed/NCBI
- Google Scholar
16. Michielse CB, van Wijk R, Reijnen L, Manders EMM, Boas S, Olivain C, et al. The nuclear protein Sge1 of Fusarium oxysporum is required for parasitic growth. PLoS Pathog. 2009;5: e1000637. pmid:19851506
- View Article
- PubMed/NCBI
- Google Scholar
17. Srikantha T, Borneman AR, Daniels KJ, Pujol C, Wu W, Seringhaus MR, et al. TOS9 regulates white-opaque switching in Candida albicans. Eukaryot Cell. 2006;5: 1674–1687. pmid:16950924
- View Article
- PubMed/NCBI
- Google Scholar
18. Zordan RE, Galgoczy DJ, Johnson AD. Epigenetic properties of white–opaque switching in Candida albicans are based on a self-sustaining transcriptional feedback loop. Proceedings of the National Academy of Sciences. 2006;103: 12807. pmid:16899543
- View Article
- PubMed/NCBI
- Google Scholar
19. Narra HP, Shubitz LF, Mandel MA, Trinh HT, Griffin K, Buntzman AS, et al. A Coccidioides posadasii CPS1 Deletion Mutant Is Avirulent and Protects Mice from Lethal Infection. Infect Immun. 2016;84: 3007–3016. pmid:27481239
- View Article
- PubMed/NCBI
- Google Scholar
20. Abuodeh RO, Orbach MJ, Mandel MA, Das A, Galgiani JN. Genetic transformation of Coccidioides immitis facilitated by Agrobacterium tumefaciens. J Infect Dis. 2000;181: 2106–2110. pmid:10837205
- View Article
- PubMed/NCBI
- Google Scholar
21. Delgado N, Hung CY, Tarcha E, Gardner MJ, Cole GT. Profiling gene expression in Coccidioides posadasii. Med Mycol. 2004;42: 59–71. pmid:14982115
- View Article
- PubMed/NCBI
- Google Scholar
22. Hung C-Y, Yu J-J, Seshan KR, Reichard U, Cole GT. A parasitic phase-specific adhesin of Coccidioides immitis contributes to the virulence of this respiratory fungal pathogen. Infect Immun. 2002;70: 3443–3456. pmid:12065484
- View Article
- PubMed/NCBI
- Google Scholar
23. Neafsey DE, Barker BM, Sharpton TJ, Stajich JE, Park DJ, Whiston E, et al. Population genomic sequencing of Coccidioides fungi reveals recent hybridization and transposon control. Genome Res. 2010;20: 938–946. pmid:20516208
- View Article
- PubMed/NCBI
- Google Scholar
24. Longo LVG, Ray SC, Puccia R, Rappleye CA. Characterization of the APSES-family transcriptional regulators of Histoplasma capsulatum. FEMS Yeast Res. 2018;18: 991. pmid:30101348
- View Article
- PubMed/NCBI
- Google Scholar
25. Rodriguez L, Voorhies M, Gilmore S, Beyhan S, Myint A, Sil A. Opposing signaling pathways regulate morphology in response to temperature in the fungal pathogen Histoplasma capsulatum. PLoS Biol. 2019;17: e3000168. pmid:31568523
- View Article
- PubMed/NCBI
- Google Scholar
26. Gladfelter AS. Control of filamentous fungal cell shape by septins and formins. Nat Rev Microbiol. 2006;4: 223–229. pmid:16429163
- View Article
- PubMed/NCBI
- Google Scholar
27. Lee B-Y, Han S-Y, Choi HG, Kim JH, Han K-H, Han D-M. Screening of growth- or development-related genes by using genomic library with inducible promoter in Aspergillus nidulans. J Microbiol. 2005;43: 523–528. pmid:16410769
- View Article
- PubMed/NCBI
- Google Scholar
28. Hung C-Y, Seshan KR, Yu J-J, Schaller R, Xue J, Basrur V, et al. A metalloproteinase of Coccidioides posadasii contributes to evasion of host detection. Infect Immun. 2005;73: 6689–6703. pmid:16177346
- View Article
- PubMed/NCBI
- Google Scholar
29. Mirbod-Donovan F, Schaller R, Hung C-Y, Xue J, Reichard U, Cole GT. Urease produced by Coccidioides posadasii contributes to the virulence of this respiratory pathogen. Infect Immun. 2006;74: 504–515. pmid:16369007
- View Article
- PubMed/NCBI
- Google Scholar
30. Bayram O, Braus GH. Coordination of secondary metabolism and development in fungi: the velvet family of regulatory proteins. FEMS Microbiol Rev. 2012;36: 1–24. pmid:21658084
- View Article
- PubMed/NCBI
- Google Scholar
31. Mead HL, Roe CC, Higgins Keppler EA, Van Dyke MCC, Laux KL, Funke AL, et al. Defining Critical Genes During Spherule Remodeling and Endospore Development in the Fungal Pathogen, Coccidioides posadasii. Front Genet. 2020;11: 483. pmid:32499817
- View Article
- PubMed/NCBI
- Google Scholar
32. Huppert M, Sun SH, Gross AJ. Evaluation of an experimental animal model for testing antifungal substances. Antimicrob Agents Chemother. 1972;1: 367–372. pmid:4670475
- View Article
- PubMed/NCBI
- Google Scholar
33. Converse JL, Besemer AR. Nutrition of the parasitic phase of Coccidioides immitis in a chemically defined liquid medium. J Bacteriol. 1959;78: 231–239. pmid:16561837
- View Article
- PubMed/NCBI
- Google Scholar
34. Orsborn KI, Shubitz LF, Peng T, Kellner EM, Orbach MJ, Haynes PA, et al. Protein expression profiling of Coccidioides posadasii by two-dimensional differential in-gel electrophoresis and evaluation of a newly recognized peroxisomal matrix protein as a recombinant vaccine candidate. Infect Immun. 2006;74: 1865–1872. pmid:16495561
- View Article
- PubMed/NCBI
- Google Scholar
35. Yu J-H, Hamari Z, Han K-H, Seo J-A, Reyes-Domínguez Y, Scazzocchio C. Double-joint PCR: a PCR-based molecular tool for gene manipulations in filamentous fungi. Fungal Genet Biol. 2004;41: 973–981. pmid:15465386
- View Article
- PubMed/NCBI
- Google Scholar
36. Kellner EM, Orsborn KI, Siegel EM, Mandel MA, Orbach MJ, Galgiani JN. Coccidioides posadasii contains a single 1,3-beta-glucan synthase gene that appears to be essential for growth. Eukaryot Cell. 2005;4: 111–120. pmid:15643067
- View Article
- PubMed/NCBI
- Google Scholar
37. Mandel MA, Galgiani JN, Kroken S, Orbach MJ. Coccidioides posadasii contains single chitin synthase genes corresponding to classes I to VII. Fungal Genet Biol. 2006;43: 775–788. pmid:16857399
- View Article
- PubMed/NCBI
- Google Scholar
38. Sharpton TJ, Stajich JE, Rounsley SD, Gardner MJ, Wortman JR, Jordar VS, et al. Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives. Genome Res. 2009;19: 1722–1731. pmid:19717792
- View Article
- PubMed/NCBI
- Google Scholar
39. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12: 323. pmid:21816040
- View Article
- PubMed/NCBI
- Google Scholar
40. Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34: 525–527. pmid:27043002
- View Article
- PubMed/NCBI
- Google Scholar
41. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43: e47. pmid:25605792
- View Article
- PubMed/NCBI
- Google Scholar
42. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3: Article3. pmid:16646809
- View Article
- PubMed/NCBI
- Google Scholar
43. Shubitz L, Peng T, Perrill R, Simons J, Orsborn K, Galgiani JN. Protection of mice against Coccidioides immitis intranasal infection by vaccination with recombinant antigen 2/PRA. Infect Immun. 2002;70: 3287–3289. pmid:12011027
- View Article
- PubMed/NCBI
- Google Scholar
44. Shubitz LF, Yu J-J, Hung C-Y, Kirkland TN, Peng T, Perrill R, et al. Improved protection of mice against lethal respiratory infection with Coccidioides posadasii using two recombinant antigens expressed as a single protein. Vaccine. 2006;24: 5904–5911. pmid:16759762
- View Article
- PubMed/NCBI
- Google Scholar
45. R Core Team. R: A language and environment for statistical computing, v3.3.3. R Foundation for Statistical Computing, Vienna, Austria. 2021 [cited 2021]. Available: https://www.r-project.org/
46. van der Walt S, Colbert SC, Varoquaux G. The NumPy Array: A Structure for Efficient Numerical Computation. Computing in Science Engineering. 2011;13: 22–30.
- View Article
- Google Scholar
47. Hunter. Matplotlib: A 2D Graphics Environment. Computing in Science & Engineering. 2007;9: 90–95.
- View Article
- Google Scholar
48. Perez F, Granger BE. IPython: A System for Interactive Scientific Computing. Computing in Science Engineering. 2007;9: 21–29.
- View Article
- Google Scholar
49. Saldanha AJ. Java Treeview—extensible visualization of microarray data. Bioinformatics. 2004;20: 3246–3248. pmid:15180930
- View Article
- PubMed/NCBI
- Google Scholar
50. Voorhies M, Cohen S, Shea TP, Petrus S, Muñoz JF, Poplawski S, et al. Chromosome-level genome assembly of a human fungal pathogen reveals synteny among geographically distinct dpecies. mBio. 2022 Jan 4;13(1):e0257421. Epub ahead of print. pmid:35089059
- View Article
- PubMed/NCBI
- Google Scholar
51. de Melo Teixeira M, Stajich JE, Sahl JW, Thompson GR, Brem RB, Dubin CA, et al. A chromosomal-level reference genome of the widely utilized laboratory strain “Silveira.” G3 (Bethesda). 2022. Feb 7:jkac031. Epub ahead of print. pmid:35137016
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Johnson L, Gaab EM, Sanchez J, Bui PQ, Nobile CJ, Hoyer KK, et al. Valley fever: danger lurking in a dust cloud. Microbes Infect. 2014;16: 591–600. pmid:25038397
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Centers for Disease Control and Prevention (CDC). Increase in reported coccidioidomycosis—United States, 1998–2011. MMWR Morb Mortal Wkly Rep. 2013;62: 217–221. pmid:23535687
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Sil A, Andrianopoulos A. Thermally Dimorphic Human Fungal Pathogens—Polyphyletic Pathogens with a Convergent Pathogenicity Trait. Cold Spring Harb Perspect Med. 2014;5: a019794. pmid:25384771
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Cole GT, Sun SH. Arthroconidium-Spherule-Endospore Transformation in Coccidioides immitis. Fungal Dimorphism. 1985. pp. 281–333.
View Article
Google Scholar

[14] View Article

[15] Google Scholar

[ref5] 5. Cole GT, Hung CY. The parasitic cell wall of Coccidioides immitis. Med Mycol. 2001;39 Suppl 1: 31–40.
View Article
Google Scholar

[17] View Article

[18] Google Scholar

[ref6] 6. 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: e41034. pmid:22911737
View Article
PubMed/NCBI
Google Scholar

[20] View Article

[21] PubMed/NCBI

[22] Google Scholar

[ref7] 7. Johannesson H, Kasuga T, Schaller RA, Good B, Gardner MJ, Townsend JP, et al. Phase-specific gene expression underlying morphological adaptations of the dimorphic human pathogenic fungus, Coccidioides posadasii. Fungal Genet Biol. 2006;43: 545–559. pmid:16697669
View Article
PubMed/NCBI
Google Scholar

[24] View Article

[25] PubMed/NCBI

[26] Google Scholar

[ref8] 8. Carlin AF, Beyhan S, Peña JF, Stajich JE, Viriyakosol S, Fierer J, et al. Transcriptional Analysis of Coccidioides immitis Mycelia and Spherules by RNA Sequencing. Journal of Fungi. 2021. pmid:34067070
View Article
PubMed/NCBI
Google Scholar

[28] View Article

[29] PubMed/NCBI

[30] Google Scholar

[ref9] 9. Viriyakosol S, Singhania A, Fierer J, Goldberg J, Kirkland TN, Woelk CH. Gene expression in human fungal pathogen Coccidioides immitis changes as arthroconidia differentiate into spherules and mature. BMC Microbiol. 2013;13: 1–1. pmid:23286760
View Article
PubMed/NCBI
Google Scholar

[32] View Article

[33] PubMed/NCBI

[34] Google Scholar

[ref10] 10. Nguyen VQ, Sil A. Temperature-induced switch to the pathogenic yeast form of Histoplasma capsulatum requires Ryp1, a conserved transcriptional regulator. Proceedings of the National Academy of Sciences. 2008;105: 4880.
View Article
Google Scholar

[36] View Article

[37] Google Scholar

[ref11] 11. Beyhan S, Gutierrez M, Voorhies M, Sil A. A Temperature-Responsive Network Links Cell Shape and Virulence Traits in a Primary Fungal Pathogen. PLoS Biol. 2013;11: e1001614. pmid:23935449
View Article
PubMed/NCBI
Google Scholar

[39] View Article

[40] PubMed/NCBI

[41] Google Scholar

[ref12] 12. Lohse MB, Zordan RE, Cain CW, Johnson AD. Distinct class of DNA-binding domains is exemplified by a master regulator of phenotypic switching in Candida albicans. Proc Natl Acad Sci U S A. 2010;107: 14105–14110. pmid:20660774
View Article
PubMed/NCBI
Google Scholar

[43] View Article

[44] PubMed/NCBI

[45] Google Scholar

[ref13] 13. Michielse CB, Becker M, Heller J, Moraga J, Collado IG, Tudzynski P. The Botrytis cinerea Reg1 protein, a putative transcriptional regulator, is required for pathogenicity, conidiogenesis, and the production of secondary metabolites. Mol Plant Microbe Interact. 2011;24: 1074–1085. pmid:21635139
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref14] 14. Jonkers W, Dong Y, Broz K, Kistler HC. The Wor1-like protein Fgp1 regulates pathogenicity, toxin synthesis and reproduction in the phytopathogenic fungus Fusarium graminearum. PLoS Pathog. 2012;8: e1002724. pmid:22693448
View Article
PubMed/NCBI
Google Scholar

[51] View Article

[52] PubMed/NCBI

[53] Google Scholar

[ref15] 15. Brown DW, Busman M, Proctor RH. Fusarium verticillioides SGE1 is required for full virulence and regulates expression of protein effector and secondary metabolite biosynthetic genes. Mol Plant Microbe Interact. 2014;27: 809–823. pmid:24742071
View Article
PubMed/NCBI
Google Scholar

[55] View Article

[56] PubMed/NCBI

[57] Google Scholar

[ref16] 16. Michielse CB, van Wijk R, Reijnen L, Manders EMM, Boas S, Olivain C, et al. The nuclear protein Sge1 of Fusarium oxysporum is required for parasitic growth. PLoS Pathog. 2009;5: e1000637. pmid:19851506
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. Srikantha T, Borneman AR, Daniels KJ, Pujol C, Wu W, Seringhaus MR, et al. TOS9 regulates white-opaque switching in Candida albicans. Eukaryot Cell. 2006;5: 1674–1687. pmid:16950924
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref18] 18. Zordan RE, Galgoczy DJ, Johnson AD. Epigenetic properties of white–opaque switching in Candida albicans are based on a self-sustaining transcriptional feedback loop. Proceedings of the National Academy of Sciences. 2006;103: 12807. pmid:16899543
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref19] 19. Narra HP, Shubitz LF, Mandel MA, Trinh HT, Griffin K, Buntzman AS, et al. A Coccidioides posadasii CPS1 Deletion Mutant Is Avirulent and Protects Mice from Lethal Infection. Infect Immun. 2016;84: 3007–3016. pmid:27481239
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref20] 20. Abuodeh RO, Orbach MJ, Mandel MA, Das A, Galgiani JN. Genetic transformation of Coccidioides immitis facilitated by Agrobacterium tumefaciens. J Infect Dis. 2000;181: 2106–2110. pmid:10837205
View Article
PubMed/NCBI
Google Scholar

[75] View Article

[76] PubMed/NCBI

[77] Google Scholar

[ref21] 21. Delgado N, Hung CY, Tarcha E, Gardner MJ, Cole GT. Profiling gene expression in Coccidioides posadasii. Med Mycol. 2004;42: 59–71. pmid:14982115
View Article
PubMed/NCBI
Google Scholar

[79] View Article

[80] PubMed/NCBI

[81] Google Scholar

[ref22] 22. Hung C-Y, Yu J-J, Seshan KR, Reichard U, Cole GT. A parasitic phase-specific adhesin of Coccidioides immitis contributes to the virulence of this respiratory fungal pathogen. Infect Immun. 2002;70: 3443–3456. pmid:12065484
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref23] 23. Neafsey DE, Barker BM, Sharpton TJ, Stajich JE, Park DJ, Whiston E, et al. Population genomic sequencing of Coccidioides fungi reveals recent hybridization and transposon control. Genome Res. 2010;20: 938–946. pmid:20516208
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref24] 24. Longo LVG, Ray SC, Puccia R, Rappleye CA. Characterization of the APSES-family transcriptional regulators of Histoplasma capsulatum. FEMS Yeast Res. 2018;18: 991. pmid:30101348
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref25] 25. Rodriguez L, Voorhies M, Gilmore S, Beyhan S, Myint A, Sil A. Opposing signaling pathways regulate morphology in response to temperature in the fungal pathogen Histoplasma capsulatum. PLoS Biol. 2019;17: e3000168. pmid:31568523
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref26] 26. Gladfelter AS. Control of filamentous fungal cell shape by septins and formins. Nat Rev Microbiol. 2006;4: 223–229. pmid:16429163
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref27] 27. Lee B-Y, Han S-Y, Choi HG, Kim JH, Han K-H, Han D-M. Screening of growth- or development-related genes by using genomic library with inducible promoter in Aspergillus nidulans. J Microbiol. 2005;43: 523–528. pmid:16410769
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref28] 28. Hung C-Y, Seshan KR, Yu J-J, Schaller R, Xue J, Basrur V, et al. A metalloproteinase of Coccidioides posadasii contributes to evasion of host detection. Infect Immun. 2005;73: 6689–6703. pmid:16177346
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref29] 29. Mirbod-Donovan F, Schaller R, Hung C-Y, Xue J, Reichard U, Cole GT. Urease produced by Coccidioides posadasii contributes to the virulence of this respiratory pathogen. Infect Immun. 2006;74: 504–515. pmid:16369007
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref30] 30. Bayram O, Braus GH. Coordination of secondary metabolism and development in fungi: the velvet family of regulatory proteins. FEMS Microbiol Rev. 2012;36: 1–24. pmid:21658084
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref31] 31. Mead HL, Roe CC, Higgins Keppler EA, Van Dyke MCC, Laux KL, Funke AL, et al. Defining Critical Genes During Spherule Remodeling and Endospore Development in the Fungal Pathogen, Coccidioides posadasii. Front Genet. 2020;11: 483. pmid:32499817
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref32] 32. Huppert M, Sun SH, Gross AJ. Evaluation of an experimental animal model for testing antifungal substances. Antimicrob Agents Chemother. 1972;1: 367–372. pmid:4670475
View Article
PubMed/NCBI
Google Scholar

[123] View Article

[124] PubMed/NCBI

[125] Google Scholar

[ref33] 33. Converse JL, Besemer AR. Nutrition of the parasitic phase of Coccidioides immitis in a chemically defined liquid medium. J Bacteriol. 1959;78: 231–239. pmid:16561837
View Article
PubMed/NCBI
Google Scholar

[127] View Article

[128] PubMed/NCBI

[129] Google Scholar

[ref34] 34. Orsborn KI, Shubitz LF, Peng T, Kellner EM, Orbach MJ, Haynes PA, et al. Protein expression profiling of Coccidioides posadasii by two-dimensional differential in-gel electrophoresis and evaluation of a newly recognized peroxisomal matrix protein as a recombinant vaccine candidate. Infect Immun. 2006;74: 1865–1872. pmid:16495561
View Article
PubMed/NCBI
Google Scholar

[131] View Article

[132] PubMed/NCBI

[133] Google Scholar

[ref35] 35. Yu J-H, Hamari Z, Han K-H, Seo J-A, Reyes-Domínguez Y, Scazzocchio C. Double-joint PCR: a PCR-based molecular tool for gene manipulations in filamentous fungi. Fungal Genet Biol. 2004;41: 973–981. pmid:15465386
View Article
PubMed/NCBI
Google Scholar

[135] View Article

[136] PubMed/NCBI

[137] Google Scholar

[ref36] 36. Kellner EM, Orsborn KI, Siegel EM, Mandel MA, Orbach MJ, Galgiani JN. Coccidioides posadasii contains a single 1,3-beta-glucan synthase gene that appears to be essential for growth. Eukaryot Cell. 2005;4: 111–120. pmid:15643067
View Article
PubMed/NCBI
Google Scholar

[139] View Article

[140] PubMed/NCBI

[141] Google Scholar

[ref37] 37. Mandel MA, Galgiani JN, Kroken S, Orbach MJ. Coccidioides posadasii contains single chitin synthase genes corresponding to classes I to VII. Fungal Genet Biol. 2006;43: 775–788. pmid:16857399
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref38] 38. Sharpton TJ, Stajich JE, Rounsley SD, Gardner MJ, Wortman JR, Jordar VS, et al. Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives. Genome Res. 2009;19: 1722–1731. pmid:19717792
View Article
PubMed/NCBI
Google Scholar

[147] View Article

[148] PubMed/NCBI

[149] Google Scholar

[ref39] 39. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12: 323. pmid:21816040
View Article
PubMed/NCBI
Google Scholar

[151] View Article

[152] PubMed/NCBI

[153] Google Scholar

[ref40] 40. Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34: 525–527. pmid:27043002
View Article
PubMed/NCBI
Google Scholar

[155] View Article

[156] PubMed/NCBI

[157] Google Scholar

[ref41] 41. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43: e47. pmid:25605792
View Article
PubMed/NCBI
Google Scholar

[159] View Article

[160] PubMed/NCBI

[161] Google Scholar

[ref42] 42. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3: Article3. pmid:16646809
View Article
PubMed/NCBI
Google Scholar

[163] View Article

[164] PubMed/NCBI

[165] Google Scholar

[ref43] 43. Shubitz L, Peng T, Perrill R, Simons J, Orsborn K, Galgiani JN. Protection of mice against Coccidioides immitis intranasal infection by vaccination with recombinant antigen 2/PRA. Infect Immun. 2002;70: 3287–3289. pmid:12011027
View Article
PubMed/NCBI
Google Scholar

[167] View Article

[168] PubMed/NCBI

[169] Google Scholar

[ref44] 44. Shubitz LF, Yu J-J, Hung C-Y, Kirkland TN, Peng T, Perrill R, et al. Improved protection of mice against lethal respiratory infection with Coccidioides posadasii using two recombinant antigens expressed as a single protein. Vaccine. 2006;24: 5904–5911. pmid:16759762
View Article
PubMed/NCBI
Google Scholar

[171] View Article

[172] PubMed/NCBI

[173] Google Scholar

[ref45] 45. R Core Team. R: A language and environment for statistical computing, v3.3.3. R Foundation for Statistical Computing, Vienna, Austria. 2021 [cited 2021]. Available: https://www.r-project.org/

[ref46] 46. van der Walt S, Colbert SC, Varoquaux G. The NumPy Array: A Structure for Efficient Numerical Computation. Computing in Science Engineering. 2011;13: 22–30.
View Article
Google Scholar

[176] View Article

[177] Google Scholar

[ref47] 47. Hunter. Matplotlib: A 2D Graphics Environment. Computing in Science & Engineering. 2007;9: 90–95.
View Article
Google Scholar

[179] View Article

[180] Google Scholar

[ref48] 48. Perez F, Granger BE. IPython: A System for Interactive Scientific Computing. Computing in Science Engineering. 2007;9: 21–29.
View Article
Google Scholar

[182] View Article

[183] Google Scholar

[ref49] 49. Saldanha AJ. Java Treeview—extensible visualization of microarray data. Bioinformatics. 2004;20: 3246–3248. pmid:15180930
View Article
PubMed/NCBI
Google Scholar

[185] View Article

[186] PubMed/NCBI

[187] Google Scholar

[ref50] 50. Voorhies M, Cohen S, Shea TP, Petrus S, Muñoz JF, Poplawski S, et al. Chromosome-level genome assembly of a human fungal pathogen reveals synteny among geographically distinct dpecies. mBio. 2022 Jan 4;13(1):e0257421. Epub ahead of print. pmid:35089059
View Article
PubMed/NCBI
Google Scholar

[189] View Article

[190] PubMed/NCBI

[191] Google Scholar

[ref51] 51. de Melo Teixeira M, Stajich JE, Sahl JW, Thompson GR, Brem RB, Dubin CA, et al. A chromosomal-level reference genome of the widely utilized laboratory strain “Silveira.” G3 (Bethesda). 2022. Feb 7:jkac031. Epub ahead of print. pmid:35137016
View Article
PubMed/NCBI
Google Scholar

[193] View Article

[194] PubMed/NCBI

[195] Google Scholar

Figures

Abstract

Author summary

Introduction

Results

Deletion of Coccidioides posadasii RYP1

Ryp1 is required for Coccidioides development

Ryp1 regulates distinct sets of genes in spherules and hyphae

Ryp1 is required for the expression of a set of spherule-enriched transcripts

Ryp1 regulates the expression of a set of hypha-enriched transcripts

Ryp1 regulates expression of morphology- or temperature-independent transcripts

RYP1 deletion mutants are avirulent

Δryp1 mutants do not provide protection against WT infection

Discussion

Materials and methods

Ethics statement

Strains and growth conditions

Creation of RYP1 gene deletion mutants

RNA-seq library preparation and sequencing

RNA-seq data analysis

Murine virulence and protection studies

Other software and libraries

Supporting information

S1 Table. Table of Kallisto quantification, limma statistics, and annotations for differentially expressed genes.

S2 Table. Table of Kallisto quantification, limma statistics, and annotations for all expressed genes.

S3 Table. Table of oligonucleotides used for Δryp1 mutant creation and analysis.

Acknowledgments

References