CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

The wtf meiotic driver gene family has unexpectedly persisted for over 100 million years

Meiotic drivers are selfish elements that bias their own transmission into more than half of the viable progeny produced by a driver+/driver- heterozygote. Meiotic drivers are thought to exist for relatively short evolutionary timespans because a driver gene or gene family is often found in a single species or in a group of very closely related species. Additionally, drivers are generally considered doomed to extinction when they spread to fixation or when suppressors arise. In this study, we examine the evolutionary history of the wtf meiotic drivers first discovered in the fission yeast Schizosaccharomyces pombe. We identify homologous genes in three other fission yeast species, S. octosporus, S. osmophilus, and S. cryophilus, which are estimated to have diverged over 100 million years ago from the S. pombe lineage. Synteny evidence supports that wtf genes were present in the common ancestor of these four species. Moreover, the ancestral genes were likely drivers as wtf genes in S. octosporus cause meiotic drive. Our findings indicate that meiotic drive systems can be maintained for long evolutionary timespans.

Diverse mating phenotypes impact the spread of wtf meiotic drivers in Schizosaccharomyces pombe

Meiotic drivers are genetic elements that break Mendel's law of segregation to be transmitted into more than half of the offspring produced by a heterozygote. The success of a driver relies on outcrossing (mating between individuals from distinct lineages) because drivers gain their advantage in heterozygotes. It is, therefore, curious that Schizosaccharomyces pombe, a species reported to rarely outcross, harbors many meiotic drivers. To address this paradox, we measured mating phenotypes in S. pombe natural isolates. We found that the propensity for cells from distinct clonal lineages to mate varies between natural isolates and can be affected both by cell density and by the available sexual partners. Additionally, we found that the observed levels of preferential mating between cells from the same clonal lineage can slow, but not prevent, the spread of a wtf meiotic driver in the absence of additional fitness costs linked to the driver. These analyses reveal parameters critical to understanding the evolution of S. pombe and help explain the success of meiotic drivers in this species.

A Novel Resistance Pathway for Calcineurin Inhibitors in the Human-Pathogenic Mucorales Mucor circinelloides

Mucormycosis is an emerging lethal fungal infection in immunocompromised patients. Mucor circinelloides is a causal agent of mucormycosis and serves as a model system to understand genetics in Mucorales. Calcineurin is a conserved virulence factor in many pathogenic fungi, and calcineurin inhibition or deletion of the calcineurin regulatory subunit (CnbR) in Mucor results in a shift from hyphal to yeast growth. We analyzed 36 calcineurin inhibitor-resistant or bypass mutants that exhibited hyphal growth in the presence of calcineurin inhibitors or in the yeast-locked cnbRΔ mutant background without carrying any mutations in known calcineurin components. We found that a majority of the mutants had altered sequence in a gene, named here bycA (bypass of calcineurin). bycA encodes an amino acid permease. We verified that both the bycAΔ single mutant and the bycAΔ cnbRΔ double mutant are resistant to calcineurin inhibitor FK506, thereby demonstrating a novel mechanism of resistance against calcineurin inhibitors. We also found that the level of expression of bycA was significantly higher in the wild-type strain treated with FK506 and in the cnbRΔ mutants but was significantly lower in the wild-type strain without FK506 treatment. These findings suggest that bycA is a negative regulator of hyphal growth and/or a positive regulator of yeast growth in Mucor and that calcineurin suppresses expression of the bycA gene at the mRNA level to promote hyphal growth. BycA is involved in the Mucor hypha-yeast transition as our data demonstrate positive correlations among bycA expression, protein kinase A activity, and Mucor yeast growth. Also, calcineurin, independently of its role in morphogenesis, contributes to virulence traits, including phagosome maturation blockade, host cell damages, and proangiogenic growth factor induction during interactions with hosts.IMPORTANCEMucor is intrinsically resistant to most known antifungals, which makes mucormycosis treatment challenging. Calcineurin is a serine/threonine phosphatase that is widely conserved across eukaryotes. When calcineurin function is inhibited in Mucor, growth shifts to a less virulent yeast growth form, which makes calcineurin an attractive target for development of new antifungal drugs. Previously, we identified two distinct mechanisms through which Mucor can become resistant to calcineurin inhibitors involving Mendelian mutations in the gene for FKBP12, including mechanisms corresponding to calcineurin A or B subunits and epimutations silencing the FKBP12 gene. Here, we identified a third novel mechanism where loss-of-function mutations in the amino acid permease corresponding to the bycA gene contribute to resistance against calcineurin inhibitors. When calcineurin activity is absent, BycA can activate protein kinase A (PKA) to promote yeast growth via a cAMP-independent pathway. Our data also show that calcineurin activity contributes to host-pathogen interactions primarily in the pathogenesis of Mucor.

5-fluorocytosine resistance is associated with hypermutation and alterations in capsule biosynthesis in Cryptococcus

Patients infected with the fungal pathogen Cryptococcus are most effectively treated with a combination of 5-fluorocytosine (5FC) and amphotericin B. 5FC acts as a prodrug, which is converted into toxic 5-fluorouracil (5FU) upon uptake into fungal cells. However, the pathogen frequently develops resistance through unclear mechanisms. Here we show that resistance to 5FC in Cryptococcus deuterogattii is acquired more frequently in isolates with defects in DNA mismatch repair that confer an elevated mutation rate. We use whole genome sequencing of 16 independent isolates to identify mutations associated with 5FC resistance in vitro. We find mutations in known resistance genes (FUR1 and FCY2) and in a gene UXS1, previously shown to encode an enzyme that converts UDP-glucuronic acid to UDP-xylose for capsule biosynthesis, but not known to play a role in 5FC metabolism. Mutations in UXS1 lead to accumulation of UDP-glucuronic acid and alterations in nucleotide metabolism, which appear to suppress toxicity of both 5FC and its toxic derivative 5FU.

The authors show that resistance to the antifungal 5-fluorocytosine in Cryptococcus deuterogattii is acquired more frequently in isolates with elevated mutation rate, and is associated with alterations in capsule biosynthesis and nucleotide metabolism.

Roles for Stress Response and Cell Wall Biosynthesis Pathways in Caspofungin Tolerance in Cryptococcus neoformans

Limited antifungal diversity and availability are growing problems for the treatment of fungal infections in the face of increasing drug resistance. The echinocandins, one of the newest classes of antifungal drugs, inhibit production of a crucial cell wall component. However, these compounds do not effectively inhibit the growth of the opportunistic fungal pathogen Cryptococcus neoformans, despite potent inhibition of the target enzyme in vitro Therefore, we performed a forward genetic screen to identify cellular processes that mediate the relative tolerance of this organism to the echinocandin drug caspofungin. Through these studies, we identified 14 genetic mutants that enhance caspofungin antifungal activity. Rather than directly affecting caspofungin antifungal activity, these mutations seem to prevent the activation of various stress-induced compensatory cellular processes. For example, the pfa4Δ mutant has defects in the palmitoylation and localization of many of its target proteins, including the Ras1 GTPase and the Chs3 chitin synthase, which are both required for caspofungin tolerance. Similarly, we have confirmed the link between caspofungin treatment and calcineurin signaling in this organism, but we suggest a deeper mechanism in which caspofungin tolerance is mediated by multiple pathways downstream of calcineurin function. In summary, we describe here several pathways in C. neoformans that contribute to the complex caspofungin tolerance phenotype in this organism.

Gastrointestinal microbiota alteration induced by Mucor circinelloides in a murine model

Mucor circinelloides is a pathogenic fungus and etiologic agent of mucormycosis. In 2013, cases of gastrointestinal illness after yogurt consumption were reported to the US FDA, and the producer found that its products were contaminated with Mucor. A previous study found that the Mucor strain isolated from an open contaminated yogurt exhibited virulence in a murine systemic infection model and showed that this strain is capable of surviving passage through the gastrointestinal tract of mice. In this study, we isolated another Mucor strain from an unopened yogurt that is closely related but distinct from the first Mucor strain and subsequently examined if Mucor alters the gut microbiota in a murine host model. DNA extracted from a ten-day course of stool samples was used to analyze the microbiota in the gastrointestinal tracts of mice exposed via ingestion of Mucor spores. The bacterial 16S rRNA gene and fungal ITS1 sequences obtained were used to identify taxa of each kingdom. Linear regressions revealed that there are changes in bacterial and fungal abundance in the gastrointestinal tracts of mice which ingested Mucor. Furthermore, we found an increased abundance of the bacterial genus Bacteroides and a decreased abundance of the bacteria Akkermansia muciniphila in the gastrointestinal tracts of exposed mice. Measurements of abundances show shifts in relative levels of multiple bacterial and fungal taxa between mouse groups. These findings suggest that exposure of the gastrointestinal tract to Mucor can alter the microbiota and, more importantly, illustrate an interaction between the intestinal mycobiota and bacteriota. In addition, Mucor was able to induce increased permeability in epithelial cell monolayers in vitro, which might be indicative of unstable intestinal barriers. Understanding how the gut microbiota is shaped is important to understand the basis of potential methods of treatment for gastrointestinal illness. How the gut microbiota changes in response to exposure, even by pathogens not considered to be causative agents of food-borne illness, may be important to how commercial food producers prevent and respond to contamination of products aimed at the public. This study provides evidence that the fungal microbiota, though understudied, may play an important role in diseases of the human gut.

Broad antifungal resistance mediated by RNAi-dependent epimutation in the basal human fungal pathogen Mucor circinelloides

Mucormycosis—an emergent, deadly fungal infection—is difficult to treat, in part because the causative species demonstrate broad clinical antifungal resistance. However, the mechanisms underlying drug resistance in these infections remain poorly understood. Our previous work demonstrated that one major agent of mucormycosis, Mucor circinelloides, can develop resistance to the antifungal agents FK506 and rapamycin through a novel, transient RNA interference-dependent mechanism known as epimutation. Epimutations silence the drug target gene and are selected by drug exposure; the target gene is re-expressed and sensitivity is restored following passage without drug. This silencing process involves generation of small RNA (sRNA) against the target gene via core RNAi pathway proteins. To further elucidate the role of epimutation in the broad antifungal resistance of Mucor, epimutants were isolated that confer resistance to another antifungal agent, 5-fluoroorotic acid (5-FOA). We identified epimutant strains that exhibit resistance to 5-FOA without mutations in PyrF or PyrG, enzymes which convert 5-FOA into the active toxic form. Using sRNA hybridization as well as sRNA library analysis, we demonstrate that these epimutants harbor sRNA against either pyrF or pyrG, and further show that this sRNA is lost after reversion to drug sensitivity. We conclude that epimutation is a mechanism capable of targeting multiple genes, enabling Mucor to develop resistance to a variety of antifungal agents. Elucidation of the role of RNAi in epimutation affords a fuller understanding of mucormycosis. Furthermore, it improves our understanding of fungal pathogenesis and adaptation to stresses, including the evolution of drug resistance.

The emerging infection mucormycosis causes high mortality in part because the major causative fungi, including Mucor circinelloides, are resistant to most clinically available antifungal drugs. We previously discovered an RNA interference-based resistance mechanism, epimutation, through which M. circinelloides develops transient resistance to the antifungal agent FK506 by altering endogenous RNA expression. We further characterize this novel mechanism by isolating epimutations in two genes that confer resistance to another antifungal agent, 5-fluoroorotic acid. Thus, we demonstrate epimutation can induce resistance to multiple antifungals by targeting a variety of genes. These results reveal epimutation plays a broad role enabling rapid and reversible fungal responses to environmental stresses, including drug exposure, and controlling antifungal drug resistance and RNA expression. As resistance to antifungals emerges, a deeper understanding of the causative mechanisms is crucial for improving treatment.

Whole-Genome Analysis Illustrates Global Clonal Population Structure of the Ubiquitous Dermatophyte Pathogen Trichophyton rubrum

Dermatophytes include fungal species that infect humans, as well as those that also infect other animals or only grow in the environment. The dermatophyte species Trichophyton rubrum is a frequent cause of skin infection in immunocompetent individuals. While members of the T. rubrum species complex have been further categorized based on various morphologies, their population structure and ability to undergo sexual reproduction are not well understood. In this study, we analyze a large set of T. rubrum and T. interdigitale isolates to examine mating types, evidence of mating, and genetic variation. We find that nearly all isolates of T. rubrum are of a single mating type, and that incubation with T. rubrum "morphotype" megninii isolates of the other mating type failed to induce sexual development. While the region around the mating type locus is characterized by a higher frequency of SNPs compared to other genomic regions, we find that the population is remarkably clonal, with highly conserved gene content, low levels of variation, and little evidence of recombination. These results support a model of recent transition to asexual growth when this species specialized to growth on human hosts.

Natural mismatch repair mutations mediate phenotypic diversity and drug resistance in Cryptococcus deuterogattii

Pathogenic microbes confront an evolutionary conflict between the pressure to maintain genome stability and the need to adapt to mounting external stresses. Bacteria often respond with elevated mutation rates, but little evidence exists of stable eukaryotic hypermutators in nature. Whole genome resequencing of the human fungal pathogen Cryptococcus deuterogattii identified an outbreak lineage characterized by a nonsense mutation in the mismatch repair component MSH2. This defect results in a moderate mutation rate increase in typical genes, and a larger increase in genes containing homopolymer runs. This allows facile inactivation of genes with coding homopolymer runs including FRR1, which encodes the target of the immunosuppresive antifungal drugs FK506 and rapamycin. Our study identifies a eukaryotic hypermutator lineage spread over two continents and suggests that pathogenic eukaryotic microbes may experience similar selection pressures on mutation rate as bacterial pathogens, particularly during long periods of clonal growth or while expanding into new environments.

As humans, we often think of genetic mutations as being bad. Over the past several decades we have seen health warnings issued on a variety of environmental exposures, from cigarettes to tanning beds, and with good reason because they cause mutations. For multicellular organisms like humans, these mutations are strongly associated with cancer. But in bacteria, this is not true. In fact, the rate at which mutations occur sometimes increases to help bacteria cope with stressful environments.

Unlike bacteria, humans are eukaryotes – the name given to organisms whose cells contain different compartments separated by membranes, such as the nucleus of the cell. For years, we have assumed that eukaryotic microbes, like fungi and parasites, act more like humans than like bacteria because work in budding yeast (another eukaryote) has suggested this to be the case. However, recent work in disease-causing fungi has shown that, much like bacteria, elevated mutation rates may help them to respond to stress. This could also enable fungi to become resistant to drugs used to treat fungal infections.

Cryptococcus deuterogattii is a fungus that causes human diseases including meningoencephalitis and a lung infection called pulmonary cryptococcosis. An ongoing outbreak of the fungus began in the Pacific Northwest of Canada in the late 1990s and emerged in the United States in 2006/2007. Among isolates closely related to those fungi causing the outbreak, three were found that appear to have a specific mutation in their DNA mismatch repair pathway, meaning that they may also experience a higher mutation rate. These strains are also less able to cause disease than others.

Billmyre et al. now demonstrate experimentally that all three isolates have a specific DNA mismatch repair defect, and show that these fungi experience elevated mutation rates, resulting in what is known as a hypermutator state. Furthermore, whole genome sequencing and phylogenetic analysis showed that these hypermutator strains are derived from the outbreak-causing fungi, and that their reduced ability to cause disease is likely a result of accumulating mutations and the loss of the ability to grow at the higher temperatures found in the human body.

Fungal infections are difficult to treat, in part because there are a limited number of available drugs. Elevated mutation rates will likely increase how often and how rapidly fungi develop resistance to these drugs. Understanding how commonly fungi exhibit a hypermutator state that could impact the development of drug resistance will therefore be important for treating patients with fungal infections, which account for millions of infections and hundreds of thousands of deaths annually worldwide.

Proteogenomics produces comprehensive and highly accurate protein-coding gene annotation in a complete genome assembly of Malassezia sympodialis

Complete and accurate genome assembly and annotation is a crucial foundation for comparative and functional genomics. Despite this, few complete eukaryotic genomes are available, and genome annotation remains a major challenge. Here, we present a complete genome assembly of the skin commensal yeast Malassezia sympodialis and demonstrate how proteogenomics can substantially improve gene annotation. Through long-read DNA sequencing, we obtained a gap-free genome assembly for M. sympodialis (ATCC 42132), comprising eight nuclear and one mitochondrial chromosome. We also sequenced and assembled four M. sympodialis clinical isolates, and showed their value for understanding Malassezia reproduction by confirming four alternative allele combinations at the two mating-type loci. Importantly, we demonstrated how proteomics data could be readily integrated with transcriptomics data in standard annotation tools. This increased the number of annotated protein-coding genes by 14% (from 3612 to 4113), compared to using transcriptomics evidence alone. Manual curation further increased the number of protein-coding genes by 9% (to 4493). All of these genes have RNA-seq evidence and 87% were confirmed by proteomics. The M. sympodialis genome assembly and annotation presented here is at a quality yet achieved only for a few eukaryotic organisms, and constitutes an important reference for future host-microbe interaction studies.

Importance of Resolving Fungal Nomenclature: the Case of Multiple Pathogenic Species in the Cryptococcus Genus

Cryptococcosis is a major fungal disease caused by members of the Cryptococcus gattii and Cryptococcus neoformans species complexes. After more than 15 years of molecular genetic and phenotypic studies and much debate, a proposal for a taxonomic revision was made. The two varieties within C. neoformans were raised to species level, and the same was done for five genotypes within C. gattii. In a recent perspective (K. J. Kwon-Chung et al., mSphere 2:e00357-16, 2017, https://doi.org/10.1128/mSphere.00357-16), it was argued that this taxonomic proposal was premature and without consensus in the community. Although the authors of the perspective recognized the existence of genetic diversity, they preferred the use of the informal nomenclature "C. neoformans species complex" and "C. gattii species complex." Here we highlight the advantage of recognizing these seven species, as ignoring these species will impede deciphering further biologically and clinically relevant differences between them, which may in turn delay future clinical advances.

Elucidation of the calcineurin-Crz1 stress response transcriptional network in the human fungal pathogen Cryptococcus neoformans

Calcineurin is a highly conserved Ca²⁺/calmodulin-dependent serine/threonine-specific protein phosphatase that orchestrates cellular Ca²⁺ signaling responses. In Cryptococcus neoformans, calcineurin is activated by multiple stresses including high temperature, and is essential for stress adaptation and virulence. The transcription factor Crz1 is a major calcineurin effector in Saccharomyces cerevisiae and other fungi. Calcineurin dephosphorylates Crz1, thereby enabling Crz1 nuclear translocation and transcription of target genes. Here we show that loss of Crz1 confers phenotypes intermediate between wild-type and calcineurin mutants, and demonstrate that deletion of the calcineurin docking domain results in the inability of Crz1 to translocate into the nucleus under thermal stress. RNA-sequencing revealed 102 genes that are regulated in a calcineurin-Crz1-dependent manner at 37°C. The majority of genes were down-regulated in cna1Δ and crz1Δ mutants, indicating these genes are normally activated by the calcineurin-Crz1 pathway at high temperature. About 58% of calcineurin-Crz1 target genes have unknown functions, while genes with known or predicted functions are involved in cell wall remodeling, calcium transport, and pheromone production. We identified three calcineurin-dependent response element motifs within the promoter regions of calcineurin-Crz1 target genes, and show that Crz1 binding to target gene promoters is increased upon thermal stress in a calcineurin-dependent fashion. Additionally, we found a large set of genes independently regulated by calcineurin, and Crz1 regulates 59 genes independently of calcineurin. Given the intermediate crz1Δ mutant phenotype, and our recent evidence for a calcineurin regulatory network impacting mRNA in P-bodies and stress granules independently of Crz1, calcineurin likely acts on factors beyond Crz1 that govern mRNA expression/stability to operate a branched transcriptional/post-transcriptional stress response network necessary for fungal virulence. Taken together, our findings reveal the core calcineurin-Crz1 stress response cascade is maintained from ascomycetes to a pathogenic basidiomycete fungus, but its output in C. neoformans appears to be adapted to promote fungal virulence.

The ubquitiously conserved serine/threonine-specific protein phosphatase calcineurin is crucial for virulence of several opportunistic human fungal pathogens including Candida albicans, Aspergillus fumigatus, and Cryptococcus neoformans. We show that Crz1 acts downstream of calcineurin, to 1) govern expression of genes involved in cell wall integrity, and calcium and small molecule transport, and 2) contribute to stress survival and virulence of C. neoformans. Our studies reveal that calcineurin also controls mRNA expression levels of other genes independently of Crz1. We propose that calcineurin operates in a branched signal transduction cascade controlling targets at both the transcriptional and post-transcriptional levels.

Identification of small RNAs in extracellular vesicles from the commensal yeast Malassezia sympodialis

Malassezia is the dominant fungus in the human skin mycobiome and is associated with common skin disorders including atopic eczema (AE)/dermatitis. Recently, it was found that Malassezia sympodialis secretes nanosized exosome-like vesicles, designated MalaEx, that carry allergens and can induce inflammatory cytokine responses. Extracellular vesicles from different cell-types including fungi have been found to deliver functional RNAs to recipient cells. In this study we assessed the presence of small RNAs in MalaEx and addressed if the levels of these RNAs differ when M. sympodialis is cultured at normal human skin pH versus the elevated pH present on the skin of patients with AE. The total number and the protein concentration of the released MalaEx harvested after 48 h culture did not differ significantly between the two pH conditions nor did the size of the vesicles. From small RNA sequence data, we identified a set of reads with well-defined start and stop positions, in a length range of 16 to 22 nucleotides consistently present in the MalaEx. The levels of small RNAs were not significantly differentially expressed between the two different pH conditions indicating that they are not influenced by the elevated pH level observed on the AE skin.

Cancer-associated isocitrate dehydrogenase mutations induce mitochondrial DNA instability

A major advance in understanding the progression and prognostic outcome of certain cancers, such as low-grade gliomas, acute myeloid leukaemia, and chondrosarcomas, has been the identification of early-occurring mutations in the NADP⁺-dependent isocitrate dehydrogenase genes IDH1 and IDH2 These mutations result in the production of the onco-metabolite D-2-hydroxyglutarate (2HG), thought to contribute to disease progression. To better understand the mechanisms of 2HG pathophysiology, we introduced the analogous glioma-associated mutations into the NADP⁺isocitrate dehydrogenase genes (IDP1, IDP2, IDP3) in Saccharomyces cerevisiae Intriguingly, expression of the mitochondrial IDP1^R148H mutant allele results in high levels of 2HG production as well as extensive mtDNA loss and respiration defects. We find no evidence for a reactive oxygen-mediated mechanism mediating this mtDNA loss. Instead, we show that 2HG production perturbs the iron sensing mechanisms as indicated by upregulation of the Aft1-controlled iron regulon and a concomitant increase in iron levels. Accordingly, iron chelation, or overexpression of a truncated AFT1 allele that dampens transcription of the iron regulon, suppresses the loss of respirative capacity. Additional suppressing factors include overexpression of the mitochondrial aldehyde dehydrogenase gene ALD5 or disruption of the retrograde response transcription factor RTG1 Furthermore, elevated α-ketoglutarate levels also suppress 2HG-mediated respiration loss; consistent with a mechanism by which 2HG contributes to mtDNA loss by acting as a toxic α-ketoglutarate analog. Our findings provide insight into the mechanisms that may contribute to 2HG oncogenicity in glioma and acute myeloid leukaemia progression, with the promise for innovative diagnostic and prognostic strategies and novel therapeutic modalities.

Gene Network Polymorphism Illuminates Loss and Retention of Novel RNAi Silencing Components in the Cryptococcus Pathogenic Species Complex

RNAi is a ubiquitous pathway that serves central functions throughout eukaryotes, including maintenance of genome stability and repression of transposon expression and movement. However, a number of organisms have lost their RNAi pathways, including the model yeast Saccharomyces cerevisiae, the maize pathogen Ustilago maydis, the human pathogen Cryptococcus deuterogattii, and some human parasite pathogens, suggesting there may be adaptive benefits associated with both retention and loss of RNAi. By comparing the RNAi-deficient genome of the Pacific Northwest Outbreak C. deuterogattii strain R265 with the RNAi-proficient genomes of the Cryptococcus pathogenic species complex, we identified a set of conserved genes that were lost in R265 and all other C. deuterogattii isolates examined. Genetic and molecular analyses reveal several of these lost genes play roles in RNAi pathways. Four novel components were examined further. Znf3 (a zinc finger protein) and Qip1 (a homolog of N. crassa Qip) were found to be essential for RNAi, while Cpr2 (a constitutive pheromone receptor) and Fzc28 (a transcription factor) are involved in sex-induced but not mitosis-induced silencing. Our results demonstrate that the mitotic and sex-induced RNAi pathways rely on the same core components, but sex-induced silencing may be a more specific, highly induced variant that involves additional specialized or regulatory components. Our studies further illustrate how gene network polymorphisms involving known components of key cellular pathways can inform identification of novel elements and suggest that RNAi loss may have been a core event in the speciation of C. deuterogattii and possibly contributed to its pathogenic trajectory.

Genome instability and mutations provoked by transposon movement are counteracted by novel defense mechanisms in organisms as diverse as fungi, plants, and mammals. In the human fungal pathogen Cryptococcus neoformans, an RNAi silencing pathway operates to defend the genome against mobile elements and transgene repeats. RNAi silencing pathways are conserved in the Cryptococcus pathogenic species complex and are mediated by canonical RNAi components. Surprisingly, several of these components are missing from all analyzed C. deuterogattii VGII strains, the molecular type responsible for the North American Pacific Northwest outbreak. To identify novel components of the RNAi pathways, we surveyed the reference genomes of C. deuterogattii, C. gattii, C. neoformans, and C. deneoformans. We identified 14 otherwise conserved genes missing in R265, including the RDP1, AGO1, and DCR1 canonical RNAi components, and focused on four potentially novel RNAi components: ZNF3, QIP1, CPR2, and FZC28. We found that Znf3 and Qip1 are both required for mitotic- and sex-induced silencing, while Cpr2 and Fzc28 contribute to sex-induced but not mitosis-induced silencing. Our studies reveal elements of RNAi pathways that operate to defend the genome during sexual development and vegetative growth and illustrate the power of network polymorphisms to illuminate novel components of biological pathways.

Unisexual versus bisexual mating in Cryptococcus neoformans: Consequences and biological impacts

Cryptococcus neoformans is an opportunistic human fungal pathogen and can undergo both bisexual and unisexual mating. Despite the fact that one mating type is dispensable for unisexual mating, the two sexual cycles share surprisingly similar features. Both mating cycles are affected by similar environmental factors and regulated by the same pheromone response pathway. Recombination takes place during unisexual reproduction in a fashion similar to bisexual reproduction and can both admix pre-existing genetic diversity and also generate diversity de novo just like bisexual reproduction. These common features may allow the unisexual life cycle to provide phenotypic and genotypic plasticity for the natural Cryptococcus population, which is predominantly α mating type, and to avoid Muller's ratchet. The morphological transition from yeast to hyphal growth during both bisexual and unisexual mating may provide increased opportunities for outcrossing and the ability to forage for nutrients at a distance. The unisexual life cycle is a key evolutionary factor for Cryptococcus as a highly successful global fungal pathogen.

Unisexual reproduction drives meiotic recombination and phenotypic and karyotypic plasticity in Cryptococcus neoformans

In fungi, unisexual reproduction, where sexual development is initiated without the presence of two compatible mating type alleles, has been observed in several species that can also undergo traditional bisexual reproduction, including the important human fungal pathogens Cryptococcus neoformans and Candida albicans. While unisexual reproduction has been well characterized qualitatively, detailed quantifications are still lacking for aspects of this process, such as the frequency of recombination during unisexual reproduction, and how this compares with bisexual reproduction. Here, we analyzed meiotic recombination during α-α unisexual and a-α bisexual reproduction of C. neoformans. We found that meiotic recombination operates in a similar fashion during both modes of sexual reproduction. Specifically, we observed that in α-α unisexual reproduction, the numbers of crossovers along the chromosomes during meiosis, recombination frequencies at specific chromosomal regions, as well as meiotic recombination hot and cold spots, are all similar to those observed during a-α bisexual reproduction. The similarity in meiosis is also reflected by the fact that phenotypic segregation among progeny collected from the two modes of sexual reproduction is also similar, with transgressive segregation being observed in both. Additionally, we found diploid meiotic progeny were also produced at similar frequencies in the two modes of sexual reproduction, and transient chromosomal loss and duplication likely occurs frequently and results in aneuploidy and loss of heterozygosity that can span entire chromosomes. Furthermore, in both α-α unisexual and a-α bisexual reproduction, we observed biased allele inheritance in regions on chromosome 4, suggesting the presence of fragile chromosomal regions that might be vulnerable to mitotic recombination. Interestingly, we also observed a crossover event that occurred within the MAT locus during α-α unisexual reproduction. Our results provide definitive evidence that α-α unisexual reproduction is a meiotic process similar to a-α bisexual reproduction.

Unisexual reproduction has been reported in several fungal species that have been traditionally thought to undergo bisexual reproduction, including major human fungal pathogens such as Cryptococcus neoformans and Candida albicans. While it has been well characterized qualitatively, quantitative description of unisexual reproduction, and detailed comparisons between unisexual and bisexual reproduction, are lacking. Here, by analyzing meiotic progeny generated from both α-α unisexual and a-α bisexual reproduction in C. neoformans, we find that the progeny collected from the two modes of sexual reproduction show similar phenotypic segregation, with transgressive segregation of several phenotypes being observed in both. Additionally, the two modes of sexual reproduction are similar in all the aspects of meiotic recombination that we have examined, providing definitive evidence that α-α unisexual reproduction is a meiotic process that operates similarly as in a-α bisexual reproduction. The ability to undergo both unisexual and bisexual reproduction may provide evolutionary advantages in environments where suitable mating partners are scarce, or where sexual reproduction is favored over asexual reproduction by mixing genetic materials and producing spores that are more tolerant of harsh environments. We discuss the implications of these findings in the context of the evolution of pathogenesis, mating types, and mating systems.

Cryptococcus gattii VGIII isolates causing infections in HIV/AIDS patients in Southern California: identification of the local environmental source as arboreal

Ongoing Cryptococcus gattii outbreaks in the Western United States and Canada illustrate the impact of environmental reservoirs and both clonal and recombining propagation in driving emergence and expansion of microbial pathogens. C. gattii comprises four distinct molecular types: VGI, VGII, VGIII, and VGIV, with no evidence of nuclear genetic exchange, indicating these represent distinct species. C. gattii VGII isolates are causing the Pacific Northwest outbreak, whereas VGIII isolates frequently infect HIV/AIDS patients in Southern California. VGI, VGII, and VGIII have been isolated from patients and animals in the Western US, suggesting these molecular types occur in the environment. However, only two environmental isolates of C. gattii have ever been reported from California: CBS7750 (VGII) and WM161 (VGIII). The incongruence of frequent clinical presence and uncommon environmental isolation suggests an unknown C. gattii reservoir in California. Here we report frequent isolation of C. gattii VGIII MATα and MATa isolates and infrequent isolation of VGI MATα from environmental sources in Southern California. VGIII isolates were obtained from soil debris associated with tree species not previously reported as hosts from sites near residences of infected patients. These isolates are fertile under laboratory conditions, produce abundant spores, and are part of both locally and more distantly recombining populations. MLST and whole genome sequence analysis provide compelling evidence that these environmental isolates are the source of human infections. Isolates displayed wide-ranging virulence in macrophage and animal models. When clinical and environmental isolates with indistinguishable MLST profiles were compared, environmental isolates were less virulent. Taken together, our studies reveal an environmental source and risk of C. gattii to HIV/AIDS patients with implications for the >1,000,000 cryptococcal infections occurring annually for which the causative isolate is rarely assigned species status. Thus, the C. gattii global health burden could be more substantial than currently appreciated.

The environmentally-acquired human pathogen C. gattii is responsible for ongoing and expanding outbreaks in the Western United States and Canada. C. gattii comprises four distinct molecular types: VGI, VGII, VGIII, and VGIV. Molecular types VGI, VGII, and VGIII have been isolated from patients and animals throughout the Western US. The Pacific Northwest and Canadian outbreak is primarily caused by C. gattii VGII. VGIII is responsible for ongoing infections in HIV/AIDS patients in Southern California. However, only two environmental C. gattii isolates have ever been identified from the Californian environment: CBS7750 (VGII) and WM161 (VGIII). We sought to collect environmental samples from areas that had confirmed reports of clinical or veterinary infections. Here we report the isolation of C. gattii VGI and VGIII from environmental soil and tree samples. C. gattii isolates were obtained from three novel tree species: Canary Island pine, American sweetgum, and a Pohutukawa tree. Genetic analysis provides robust evidence that these environmental isolates are the source of human infections.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Cryptococcus neoformans var. grubii is a major human pathogen responsible for deadly meningoencephalitis in immunocompromised patients. Here, we report the sequencing and annotation of its genome. Evidence for extensive intron splicing, antisense transcription, non-coding RNAs, and alternative polyadenylation indicates the potential for highly intricate regulation of gene expression in this opportunistic pathogen. In addition, detailed molecular, genetic, and genomic studies were performed to characterize structural features of the genome, including centromeres and origins of replication. Finally, the phenotypic and genome re-sequencing analysis of a collection of isolates of the reference H99 strain resulting from laboratory passage revealed that microevolutionary processes during in vitro culturing of pathogenic fungi can impact virulence.

Photosynthetic microbial fuel cells with positive light response

The current study introduces an aerobic single-chamber photosynthetic microbial fuel cell (PMFC). Evaluation of PMFC performance using naturally growing fresh-water photosynthetic biofilm revealed a weak positive light response, that is, an increase in cell voltage upon illumination. When the PMFC anodes were coated with electrically conductive polymers, the rate of voltage increased and the amplitude of the light response improved significantly. The rapid immediate positive response to light was consistent with a mechanism postulating that the photosynthetic electron-transfer chain is the source of the electrons harvested on the anode surface. This mechanism is fundamentally different from the one exploited in previously designed anaerobic microbial fuel cells (MFCs), sediment MFCs, or anaerobic PMFCs, where the electrons are derived from the respiratory electron-transfer chain. The power densities produced in PMFCs were substantially lower than those that are currently reported for conventional MFC (0.95 mW/m(2) for polyaniline-coated and 1.3 mW/m(2) for polypyrrole-coated anodes). However, the PMFC did not depend on an organic substrate as an energy source and was powered only by light energy. Its operation was CO(2)-neutral and did not require buffers or exogenous electron transfer shuttles.