Evolutionary Persistence of DNA Methylation for Millions of Years after Ancient Loss of a De Novo Methyltransferase

Cytosine methylation of DNA is a widespread modification of DNA that plays numerous critical roles. In the yeast Cryptococcus neoformans, CG methylation occurs in transposon-rich repeats and requires the DNA methyltransferase Dnmt5. We show that Dnmt5 displays exquisite maintenance-type specificity in vitro and in vivo and utilizes similar in vivo cofactors as the metazoan maintenance methylase Dnmt1. Remarkably, phylogenetic and functional analysis revealed that the ancestral species lost the gene for a de novo methylase, DnmtX, between 50-150 mya. We examined how methylation has persisted since the ancient loss of DnmtX. Experimental and comparative studies reveal efficient replication of methylation patterns in C. neoformans, rare stochastic methylation loss and gain events, and the action of natural selection. We propose that an epigenome has been propagated for >50 million years through a process analogous to Darwinian evolution of the genome.

Complete Genome Sequences for Two Talaromyces marneffei Clinical Isolates from Northern and Southern Vietnam

Talaromyces marneffei is a thermally dimorphic fungus endemic in China and Southeast Asia that causes fatal infections in immunocompromised individuals, particularly in patients with advanced HIV disease. Here, we report the complete genome sequences of two clinical isolates from northern and southern Vietnam.

Cryptococcus neoformans Recovered From Olive Trees (Olea europaea) in Turkey Reveal Allopatry With African and South American Lineages

Cryptococcus species are life-threatening human fungal pathogens that cause cryptococcal meningoencephalitis in both immunocompromised and healthy hosts. The natural environmental niches of Cryptococcus include pigeon (Columba livia) guano, soil, and a variety of tree species such as Eucalyptus camaldulensis, Ceratonia siliqua, Platanus orientalis, and Pinus spp. Genetic and genomic studies of extensive sample collections have provided insights into the population distribution and composition of different Cryptococcus species in geographic regions around the world. However, few such studies examined Cryptococcus in Turkey. We sampled 388 Olea europaea (olive) and 132 E. camaldulensis trees from seven locations in coastal and inland areas of the Aegean region of Anatolian Turkey in September 2016 to investigate the distribution and genetic diversity present in the natural Cryptococcus population. We isolated 84 Cryptococcus neoformans strains (83 MATα and 1 MAT a) and 3 Cryptococcus deneoformans strains (all MATα) from 87 (22.4% of surveyed) O. europaea trees; a total of 32 C. neoformans strains were isolated from 32 (24.2%) of the E. camaldulensis trees, all of which were MATα. A statistically significant difference was observed in the frequency of C. neoformans isolation between coastal and inland areas (P < 0.05). Interestingly, the MAT a C. neoformans isolate was fertile in laboratory crosses with VNI and VNB MATα tester strains and produced robust hyphae, basidia, and basidiospores, thus suggesting potential sexual reproduction in the natural population. Sequencing analyses of the URA5 gene identified at least five different genotypes among the isolates. Population genetics and genomic analyses revealed that most of the isolates in Turkey belong to the VNBII lineage of C. neoformans, which is predominantly found in southern Africa; these isolates are part of a distinct minor clade within VNBII that includes several isolates from Zambia and Brazil. Our study provides insights into the geographic distribution of different C. neoformans lineages in the Mediterranean region and highlights the need for wider geographic sampling to gain a better understanding of the natural habitats, migration, epidemiology, and evolution of this important human fungal pathogen.

Genetic and Genomic Analyses Reveal Boundaries between Species Closely Related to Cryptococcus Pathogens

Speciation is a central mechanism of biological diversification. While speciation is well studied in plants and animals, in comparison, relatively little is known about speciation in fungi. One fungal model is the Cryptococcus genus, which is best known for the pathogenic Cryptococcus neoformans/Cryptococcus gattii species complex that causes >200,000 new human infections annually. Elucidation of how these species evolved into important human-pathogenic species remains challenging and can be advanced by studying the most closely related nonpathogenic species, Cryptococcus amylolentus and Tsuchiyaea wingfieldii However, these species have only four known isolates, and available data were insufficient to determine species boundaries within this group. By analyzing full-length chromosome assemblies, we reappraised the phylogenetic relationships of the four available strains, confirmed the genetic separation of C. amylolentus and T. wingfieldii (now Cryptococcus wingfieldii), and revealed an additional cryptic species, for which the name Cryptococcus floricola is proposed. The genomes of the three species are ∼6% divergent and exhibit significant chromosomal rearrangements, including inversions and a reciprocal translocation that involved intercentromeric ectopic recombination, which together likely impose significant barriers to genetic exchange. Using genetic crosses, we show that while C. wingfieldii cannot interbreed with any of the other strains, C. floricola can still undergo sexual reproduction with C. amylolentus However, most of the resulting spores were inviable or sterile or showed reduced recombination during meiosis, indicating that intrinsic postzygotic barriers had been established. Our study and genomic data will foster additional studies addressing fungal speciation and transitions between nonpathogenic and pathogenic Cryptococcus lineages.IMPORTANCE The evolutionary drivers of speciation are critical to our understanding of how new pathogens arise from nonpathogenic lineages and adapt to new environments. Here we focus on the Cryptococcus amylolentus species complex, a nonpathogenic fungal lineage closely related to the human-pathogenic Cryptococcus neoformans/Cryptococcus gattii complex. Using genetic and genomic analyses, we reexamined the species boundaries of four available isolates within the C. amylolentus complex and revealed three genetically isolated species. Their genomes are ∼6% divergent and exhibit chromosome rearrangements, including translocations and small-scale inversions. Although two of the species (C. amylolentus and newly described C. floricola) were still able to interbreed, the resulting hybrid progeny were usually inviable or sterile, indicating that barriers to reproduction had already been established. These results advance our understanding of speciation in fungi and highlight the power of genomics in assisting our ability to correctly identify and discriminate fungal species.

Coordinated host-pathogen transcriptional dynamics revealed using sorted subpopulations and single macrophages infected with Candida albicans

The outcome of fungal infections depends on interactions with innate immune cells. Within a population of macrophages encountering Candida albicans, there are distinct host-pathogen trajectories; however, little is known about the molecular heterogeneity that governs these fates. Here we developed an experimental system to separate interaction stages and single macrophage cells infected with C. albicans from uninfected cells and assessed transcriptional variability in the host and fungus. Macrophages displayed an initial up-regulation of pathways involved in phagocytosis and proinflammatory response after C. albicans exposure that declined during later time points. Phagocytosed C. albicans shifted expression programs to survive the nutrient poor phagosome and remodeled the cell wall. The transcriptomes of single infected macrophages and phagocytosed C. albicans displayed a tightly coordinated shift in gene expression co-stages and revealed expression bimodality and differential splicing that may drive infection outcome. This work establishes an approach for studying host-pathogen trajectories to resolve heterogeneity in dynamic populations.

The outcomes of the interactions between individual host cells and pathogens are heterogeneous. Here, the authors assess transcriptional variability in both host and pathogen during infection of macrophages with the fungus Candida albicans, using sorted subpopulations and single macrophages.

Pathways of Pathogenicity: Transcriptional Stages of Germination in the Fatal Fungal Pathogen Rhizopus delemar

Rhizopus delemar is an invasive fungal pathogen responsible for the frequently fatal disease mucormycosis. Germination, a crucial mechanism by which infectious spores of Rhizopus delemar cause disease, is a key developmental process that transforms the dormant spore state into a vegetative one. The molecular mechanisms that underpin this transformation may be key to controlling mucormycosis; however, the regulation of germination remains poorly understood. This study describes the phenotypic and transcriptional changes that take place over the course of germination. This process is characterized by four distinct stages: dormancy, isotropic swelling, germ tube emergence, and hyphal growth. Dormant spores are shown to be transcriptionally unique, expressing a subset of transcripts absent in later developmental stages. A large shift in the expression profile is prompted by the initiation of germination, with genes involved in respiration, chitin, cytoskeleton, and actin regulation appearing to be important for this transition. A period of transcriptional consistency can be seen throughout isotropic swelling, before the transcriptional landscape shifts again at the onset of hyphal growth. This study provides a greater understanding of the regulation of germination and highlights processes involved in transforming Rhizopus delemar from a single-cellular to multicellular organism.IMPORTANCE Germination is key to the growth of many organisms, including fungal spores. Mucormycete spores exist abundantly within the environment and germinate to form hyphae. These spores are capable of infecting immunocompromised individuals, causing the disease mucormycosis. Germination from spore to hyphae within patients leads to angioinvasion, tissue necrosis, and often fatal infections. This study advances our understanding of how spore germination occurs in the mucormycetes, identifying processes we may be able to inhibit to help prevent or treat mucormycosis.

Genomic characterization of recurrent mold infections in thoracic transplant recipients

Invasive mold disease in thoracic organ transplant recipients is a well-recognized complication, but the long-term persistence of molds within the human body and evasion of host defenses has not been well-described. We present 2 cases of invasive mold disease (Verruconis gallopava and Aspergillus fumigatus) in thoracic transplant recipients who had the same mold cultured years prior to the invasive disease presentation. The paired isolates from the index and recurrent infections in both patients were compared using whole-genome sequencing to determine if the same strain of mold caused both the index and recurrent infections. In Case 1, the isolates were found to be of the same strain indicating that the initial colonizing isolate identified pre-transplant eventually caused invasive mold disease post-transplant while in Case 2, the 2 isolates were not of the same strain. These results demonstrate the distinct possibility of molds both persisting within the human body for years prior to invasive mold disease or the long-term risk of recurrent, persistent infection with more than one strain. Further studies of long-term molecular epidemiology of IMD and risk factors for mold persistence in transplant recipients are encouraged.

Titan cells formation in Cryptococcus neoformans is finely tuned by environmental conditions and modulated by positive and negative genetic regulators

The pathogenic fungus Cryptococcus neoformans exhibits morphological changes in cell size during lung infection, producing both typical size 5 to 7 μm cells and large titan cells (> 10 μm and up to 100 μm). We found and optimized in vitro conditions that produce titan cells in order to identify the ancestry of titan cells, the environmental determinants, and the key gene regulators of titan cell formation. Titan cells generated in vitro harbor the main characteristics of titan cells produced in vivo including their large cell size (>10 μm), polyploidy with a single nucleus, large vacuole, dense capsule, and thick cell wall. Here we show titan cells derived from the enlargement of progenitor cells in the population independent of yeast growth rate. Change in the incubation medium, hypoxia, nutrient starvation and low pH were the main factors that trigger titan cell formation, while quorum sensing factors like the initial inoculum concentration, pantothenic acid, and the quorum sensing peptide Qsp1p also impacted titan cell formation. Inhibition of ergosterol, protein and nucleic acid biosynthesis altered titan cell formation, as did serum, phospholipids and anti-capsular antibodies in our settings. We explored genetic factors important for titan cell formation using three approaches. Using H99-derivative strains with natural genetic differences, we showed that titan cell formation was dependent on LMP1 and SGF29 genes. By screening a gene deletion collection, we also confirmed that GPR4/5-RIM101, and CAC1 genes were required to generate titan cells and that the PKR1, TSP2, USV101 genes negatively regulated titan cell formation. Furthermore, analysis of spontaneous Pkr1 loss-of-function clinical isolates confirmed the important role of the Pkr1 protein as a negative regulator of titan cell formation. Through development of a standardized and robust in vitro assay, our results provide new insights into titan cell biogenesis with the identification of multiple important factors/pathways.

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.

Population genomics and the evolution of virulence in the fungal pathogen Cryptococcus neoformans

Cryptococcus neoformans is an opportunistic fungal pathogen that causes approximately 625,000 deaths per year from nervous system infections. Here, we leveraged a unique, genetically diverse population of C. neoformans from sub-Saharan Africa, commonly isolated from mopane trees, to determine how selective pressures in the environment coincidentally adapted C. neoformans for human virulence. Genome sequencing and phylogenetic analysis of 387 isolates, representing the global VNI and African VNB lineages, highlighted a deep, nonrecombining split in VNB (herein, VNBI and VNBII). VNBII was enriched for clinical samples relative to VNBI, while phenotypic profiling of 183 isolates demonstrated that VNBI isolates were significantly more resistant to oxidative stress and more heavily melanized than VNBII isolates. Lack of melanization in both lineages was associated with loss-of-function mutations in the BZP4 transcription factor. A genome-wide association study across all VNB isolates revealed sequence differences between clinical and environmental isolates in virulence factors and stress response genes. Inositol transporters and catabolism genes, which process sugars present in plants and the human nervous system, were identified as targets of selection in all three lineages. Further phylogenetic and population genomic analyses revealed extensive loss of genetic diversity in VNBI, suggestive of a history of population bottlenecks, along with unique evolutionary trajectories for mating type loci. These data highlight the complex evolutionary interplay between adaptation to natural environments and opportunistic infections, and that selection on specific pathways may predispose isolates to human virulence.

Advances in Cryptococcus genomics: insights into the evolution of pathogenesis

Cryptococcus species are the causative agents of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals. Initial work on the molecular epidemiology of this fungal pathogen utilized genotyping approaches to describe the genetic diversity and biogeography of two species, Cryptococcus neoformans and Cryptococcus gattii. Whole genome sequencing of representatives of both species resulted in reference assemblies enabling a wide array of downstream studies and genomic resources. With the increasing availability of whole genome sequencing, both species have now had hundreds of individual isolates sequenced, providing fine-scale insight into the evolution and diversification of Cryptococcus and allowing for the first genome-wide association studies to identify genetic variants associated with human virulence. Sequencing has also begun to examine the microevolution of isolates during prolonged infection and to identify variants specific to outbreak lineages, highlighting the potential role of hyper-mutation in evolving within short time scales. We can anticipate that further advances in sequencing technology and sequencing microbial genomes at scale, including metagenomics approaches, will continue to refine our view of how the evolution of Cryptococcus drives its success as a pathogen.

Novel taxa of thermally dimorphic systemic pathogens in the Ajellomycetaceae (Onygenales)

Recent discoveries of novel systemic fungal pathogens with thermally dimorphic yeast-like phases have challenged the current taxonomy of the Ajellomycetaceae, a family currently comprising the genera Blastomyces, Emmonsia, Emmonsiellopsis, Helicocarpus, Histoplasma, Lacazia and Paracoccidioides. Our morphological, phylogenetic and phylogenomic analyses demonstrated species relationships and their specific phenotypes, clarified generic boundaries and provided the first annotated genome assemblies to support the description of two new species. A new genus, Emergomyces, accommodates Emmonsia pasteuriana as type species, and the new species Emergomyces africanus, the aetiological agent of case series of disseminated infections in South Africa. Both species produce small yeast cells that bud at a narrow base at 37°C and lack adiaspores, classically associated with the genus Emmonsia. Another novel dimorphic pathogen, producing broad-based budding cells at 37°C and occurring outside North America, proved to belong to the genus Blastomyces, and is described as Blastomyces percursus.

Genomic innovations linked to infection strategies across emerging pathogenic chytrid fungi

To understand the evolutionary pathways that lead to emerging infections of vertebrates, here we explore the genomic innovations that allow free-living chytrid fungi to adapt to and colonize amphibian hosts. Sequencing and comparing the genomes of two pathogenic species of Batrachochytrium to those of close saprophytic relatives reveals that pathogenicity is associated with remarkable expansions of protease and cell wall gene families, while divergent infection strategies are linked to radiations of lineage-specific gene families. By comparing the host–pathogen response to infection for both pathogens, we illuminate the traits that underpin a strikingly different immune response within a shared host species. Our results show that, despite commonalities that promote infection, specific gene-family radiations contribute to distinct infection strategies. The breadth and evolutionary novelty of candidate virulence factors that we discover underscores the urgent need to halt the advance of pathogenic chytrids and prevent incipient loss of biodiversity.

Batrachochytrium dendrobatidis and B. salamandrivorans are both important pathogens of amphibians, but they differ in their host ranges, infection strategies, and host immune responses. Here, Farrer and colleagues compare their genomes and transcriptomes to identify the genetic basis of these differences.

Comparative Genomics Integrated with Association Analysis Identifies Candidate Effector Genes Corresponding to Lr20 in Phenotype-Paired Puccinia triticina Isolates from Australia

Leaf rust is one of the most common and damaging diseases of wheat, and is caused by an obligate biotrophic basidiomycete, Puccinia triticina (Pt). In the present study, 20 Pt isolates from Australia, comprising 10 phenotype-matched pairs with contrasting pathogenicity for Lr20, were analyzed using whole genome sequencing. Compared to the reference genome of the American Pt isolate 1-1 BBBD Race 1, an average of 404,690 single nucleotide polymorphisms (SNPs) per isolate was found and the proportion of heterozygous SNPs was above 87% in the majority of the isolates, demonstrating a high level of polymorphism and a high rate of heterozygosity. From the genome-wide SNPs, a phylogenetic tree was inferred, which consisted of a large clade of 15 isolates representing diverse presumed clonal lineages including 14 closely related isolates and the more diverged isolate 670028, and a small clade of five isolates characterized by lower heterozygosity level. Principle component analysis detected three distinct clusters, corresponding exactly to the two major subsets of the small clade and the large clade comprising all 15 isolates without further separation of isolate 670028. While genome-wide association analysis identified 302 genes harboring at least one SNP associated with Lr20 virulence (p < 0.05), a Wilcoxon rank sum test revealed that 36 and 68 genes had significant (p < 0.05) and marginally significant (p < 0.1) differences in the counts of non-synonymous mutations between Lr20 avirulent and virulent groups, respectively. Twenty of these genes were predicted to have a signal peptide without a transmembrane segment, and hence identified as candidate effector genes corresponding to Lr20. SNP analysis also implicated the potential involvement of epigenetics and small RNA in Pt pathogenicity. Future studies are thus warranted to investigate the biological functions of the candidate effectors as well as the gene regulation mechanisms at epigenetic and post-transcription levels. Our study is the first to integrate phenotype-genotype association with effector prediction in Pt genomes, an approach that may circumvent some of the technical difficulties in working with obligate rust fungi and accelerate avirulence gene identification.

Transcriptional analysis of mating and pre-infection stages of the anther smut, Microbotryum lychnidis-dioicae

Microbotryum lychnidis-dioicae is an obligate biotrophic parasite of the wildflower species Silene latifolia. This dikaryotic fungus, commonly known as an anther smut, requires that haploid, yeast-like sporidia of opposite mating types fuse and differentiate into dikaryotic hyphae that penetrate host tissue as part of the fungal life cycle. Mating occurs under conditions of cool temperatures and limited nutrients. Further development requires host cues or chemical mimics, including a variety of lipids, e.g. phytols. To identify global changes in transcription associated with developmental shifts, RNA-Seq was conducted at several in vitro stages of fungal propagation, i.e. haploid cells grown independently on rich and nutrient-limited media, mated cells on nutrient-limited media as well as a time course of such mated cells exposed to phytol. Comparison of haploid cells grown under rich and nutrient-limited conditions identified classes of genes probably associated with general nutrient availability, including components of the RNAi machinery. Some gene enrichment patterns comparing the nutrient-limited and mated transcriptomes suggested gene expression changes associated with the mating programme (e.g. homeodomain binding proteins, secreted proteins, proteins unique to M. lychnidis-dioicae¸ multicopper oxidases and RhoGEFs). Analysis for phytol treatment compared with mated cells alone allowed identification of genes likely to be involved in the dikaryotic switch (e.g. oligopeptide transporters). Gene categories of particular note in all three conditions included those in the major facilitator superfamily, proteins containing PFAM domains of the secretory lipase family as well as proteins predicted to be secreted, many of which have the hallmarks of fungal effectors with potential roles in pathogenicity.

Paracoccidioides spp. catalases and their role in antioxidant defense against host defense responses

Dimorphic human pathogenic fungi interact with host effector cells resisting their microbicidal mechanisms. Yeast cells are able of surviving within the tough environment of the phagolysosome by expressing an antioxidant defense system that provides protection against host-derived reactive oxygen species (ROS). This includes the production of catalases (CATs). Here we identified and analyzed the role of CAT isoforms in Paracoccidioides, the etiological agent of paracoccidioidomycosis. Firstly, we found that one of these isoforms was absent in the closely related dimorphic pathogen Coccidioides and dermatophytes, but all of them were conserved in Paracoccidioides, Histoplasma and Blastomyces species. We probed the contribution of CATs in Paracoccidioides by determining the gene expression levels of each isoform through quantitative RT-qPCR, in both the yeast and mycelia phases, and during the morphological switch (transition and germination), as well as in response to oxidative agents and during interaction with neutrophils. PbCATP was preferentially expressed in the pathogenic yeast phase, and was associated to the response against exogenous H₂O₂. Therefore, we created and analyzed the virulence defects of a knockdown strain for this isoform, and found that CATP protects yeast cells from H₂O₂ generated in vitro and is relevant during lung infection. On the other hand, CATA and CATB seem to contribute to ROS homeostasis in Paracoccidioides cells, during endogenous oxidative stress. CAT isoforms in Paracoccidioides might be coordinately regulated during development and dimorphism, and differentially expressed in response to different stresses to control ROS homeostasis during the infectious process, contributing to the virulence of Paracoccidioides.

Microevolutionary traits and comparative population genomics of the emerging pathogenic fungus Cryptococcus gattii

Emerging fungal pathogens cause an expanding burden of disease across the animal kingdom, including a rise in morbidity and mortality in humans. Yet, we currently have only a limited repertoire of available therapeutic interventions. A greater understanding of the mechanisms of fungal virulence and of the emergence of hypervirulence within species is therefore needed for new treatments and mitigation efforts. For example, over the past decade, an unusual lineage of Cryptococcus gattii, which was first detected on Vancouver Island, has spread to the Canadian mainland and the Pacific Northwest infecting otherwise healthy individuals. The molecular changes that led to the development of this hypervirulent cryptococcal lineage remain unclear. To explore this, we traced the history of similar microevolutionary events that can lead to changes in host range and pathogenicity. Here, we detail fine-resolution mapping of genetic differences between two highly related Cryptococcus gattii VGIIc isolates that differ in their virulence traits (phagocytosis, vomocytosis, macrophage death, mitochondrial tubularization and intracellular proliferation). We identified a small number of single site variants within coding regions that potentially contribute to variations in virulence. We then extended our methods across multiple lineages of C. gattii to study how selection is acting on key virulence genes within different lineages.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Simultaneous Emergence of Multidrug-Resistant Candida auris on 3 Continents Confirmed by Whole-Genome Sequencing and Epidemiological Analyses

BACKGROUND

Candida auris, a multidrug-resistant yeast that causes invasive infections, was first described in 2009 in Japan and has since been reported from several countries.

METHODS

To understand the global emergence and epidemiology of C. auris, we obtained isolates from 54 patients with C. auris infection from Pakistan, India, South Africa, and Venezuela during 2012-2015 and the type specimen from Japan. Patient information was available for 41 of the isolates. We conducted antifungal susceptibility testing and whole-genome sequencing (WGS).

RESULTS

Available clinical information revealed that 41% of patients had diabetes mellitus, 51% had undergone recent surgery, 73% had a central venous catheter, and 41% were receiving systemic antifungal therapy when C. auris was isolated. The median time from admission to infection was 19 days (interquartile range, 9-36 days), 61% of patients had bloodstream infection, and 59% died. Using stringent break points, 93% of isolates were resistant to fluconazole, 35% to amphotericin B, and 7% to echinocandins; 41% were resistant to 2 antifungal classes and 4% were resistant to 3 classes. WGS demonstrated that isolates were grouped into unique clades by geographic region. Clades were separated by thousands of single-nucleotide polymorphisms, but within each clade isolates were clonal. Different mutations in ERG11 were associated with azole resistance in each geographic clade.

CONCLUSIONS

C. auris is an emerging healthcare-associated pathogen associated with high mortality. Treatment options are limited, due to antifungal resistance. WGS analysis suggests nearly simultaneous, and recent, independent emergence of different clonal populations on 3 continents. Risk factors and transmission mechanisms need to be elucidated to guide control measures.

A Multistate Toggle Switch Defines Fungal Cell Fates and Is Regulated by Synergistic Genetic Cues

Heritable epigenetic changes underlie the ability of cells to differentiate into distinct cell types. Here, we demonstrate that the fungal pathogen Candida tropicalis exhibits multipotency, undergoing stochastic and reversible switching between three cellular states. The three cell states exhibit unique cellular morphologies, growth rates, and global gene expression profiles. Genetic analysis identified six transcription factors that play key roles in regulating cell differentiation. In particular, we show that forced expression of Wor1 or Efg1 transcription factors can be used to manipulate transitions between all three cell states. A model for tristability is proposed in which Wor1 and Efg1 are self-activating but mutually antagonistic transcription factors, thereby forming a symmetrical self-activating toggle switch. We explicitly test this model and show that ectopic expression of WOR1 can induce white-to-hybrid-to-opaque switching, whereas ectopic expression of EFG1 drives switching in the opposite direction, from opaque-to-hybrid-to-white cell states. We also address the stability of induced cell states and demonstrate that stable differentiation events require ectopic gene expression in combination with chromatin-based cues. These studies therefore experimentally test a model of multistate stability and demonstrate that transcriptional circuits act synergistically with chromatin-based changes to drive cell state transitions. We also establish close mechanistic parallels between phenotypic switching in unicellular fungi and cell fate decisions during stem cell reprogramming.

Genome Diversity, Recombination, and Virulence across the Major Lineages of Paracoccidioides

The Paracoccidioides genus includes two species of thermally dimorphic fungi that cause paracoccidioidomycosis, a neglected health-threatening human systemic mycosis endemic to Latin America. To examine the genome evolution and the diversity of Paracoccidioides spp., we conducted whole-genome sequencing of 31 isolates representing the phylogenetic, geographic, and ecological breadth of the genus. These samples included clinical, environmental and laboratory reference strains of the S1, PS2, PS3, and PS4 lineages of P. brasiliensis and also isolates of Paracoccidioides lutzii species. We completed the first annotated genome assemblies for the PS3 and PS4 lineages and found that gene order was highly conserved across the major lineages, with only a few chromosomal rearrangements. Comparing whole-genome assemblies of the major lineages with single-nucleotide polymorphisms (SNPs) predicted from the remaining 26 isolates, we identified a deep split of the S1 lineage into two clades we named S1a and S1b. We found evidence for greater genetic exchange between the S1b lineage and all other lineages; this may reflect the broad geographic range of S1b, which is often sympatric with the remaining, largely geographically isolated lineages. In addition, we found evidence of positive selection for the GP43 and PGA1 antigen genes and genes coding for other secreted proteins and proteases and lineage-specific loss-of-function mutations in cell wall and protease genes; these together may contribute to virulence and host immune response variation among natural isolates of Paracoccidioides spp. These insights into the recent evolutionary events highlight important differences between the lineages that could impact the distribution, pathogenicity, and ecology of Paracoccidioides. IMPORTANCE Characterization of genetic differences between lineages of the dimorphic human-pathogenic fungus Paracoccidioides can identify changes linked to important phenotypes and guide the development of new diagnostics and treatments. In this article, we compared genomes of 31 diverse isolates representing the major lineages of Paracoccidioides spp. and completed the first annotated genome sequences for the PS3 and PS4 lineages. We analyzed the population structure and characterized the genetic diversity among the lineages of Paracoccidioides, including a deep split of S1 into two lineages (S1a and S1b), and differentiated S1b, associated with most clinical cases, as the more highly recombining and diverse lineage. In addition, we found patterns of positive selection in surface proteins and secreted enzymes among the lineages, suggesting diversifying mechanisms of pathogenicity and adaptation across this species complex. These genetic differences suggest associations with the geographic range, pathogenicity, and ecological niches of Paracoccidioides lineages.

Lower prevalence but similar fitness in a parasitic fungus at higher radiation levels near Chernobyl

Nuclear disasters at Chernobyl and Fukushima provide examples of effects of acute ionizing radiation on mutations that can affect the fitness and distribution of species. Here, we investigated the prevalence of Microbotryum lychnidis-dioicae, a pollinator-transmitted fungal pathogen of plants causing anther-smut disease in Chernobyl, its viability, fertility and karyotype variation, and the accumulation of nonsynonymous mutations in its genome. We collected diseased flowers of Silene latifolia from locations ranging by more than two orders of magnitude in background radiation, from 0.05 to 21.03 μGy/h. Disease prevalence decreased significantly with increasing radiation level, possibly due to lower pollinator abundance and altered pollinator behaviour. Viability and fertility, measured as the budding rate of haploid sporidia following meiosis from the diploid teliospores, did not vary with increasing radiation levels and neither did karyotype overall structure and level of chromosomal size heterozygosity. We sequenced the genomes of twelve samples from Chernobyl and of four samples collected from uncontaminated areas and analysed alignments of 6068 predicted genes, corresponding to 1.04 × 10(7)  base pairs. We found no dose-dependent differences in substitution rates (neither dN, dS, nor dN/dS). Thus, we found no significant evidence of increased deleterious mutation rates at higher levels of background radiation in this plant pathogen. We even found lower levels of nonsynonymous substitution rates in contaminated areas compared to control regions, suggesting that purifying selection was stronger in contaminated than uncontaminated areas. We briefly discuss the possibilities for a mechanistic basis of radio resistance in this nonmelanized fungus.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Dimorphic fungal pathogens including Blastomyces are the cause of major fungal diseases in North and South America. The genus Emmonsia includes species infecting small mammals as well as a newly emerging pathogenic species recently reported in HIV-positive patients in South Africa. Here, we synthesize both genome sequencing of four isolates of Blastomyces and two species of Emmonsia as well as deep sequencing of Blastomyces RNA to draw major new insights into the evolution of this group and the pathogen response to infection. We investigate the trajectory of genome evolution of this group, characterizing the phylogenetic relationships of these species, a remarkable genome expansion that formed large isochore-like regions of low GC content in Blastomyces, and variation of gene content, related to host interaction, among the dimorphic fungal pathogens. Using RNA-Seq, we profile the response of Blastomyces to macrophage and mouse pulmonary infection, identifying key pathways and novel virulence factors. The identification of key fungal genes involved in adaptation to the host suggests targets for further study and therapeutic intervention in Blastomyces and related dimorphic fungal pathogens.

Sex and parasites: genomic and transcriptomic analysis of Microbotryum lychnidis-dioicae, the biotrophic and plant-castrating anther smut fungus

Background

The genus Microbotryum includes plant pathogenic fungi afflicting a wide variety of hosts with anther smut disease. Microbotryum lychnidis-dioicae infects Silene latifolia and replaces host pollen with fungal spores, exhibiting biotrophy and necrosis associated with altering plant development.

Results

We determined the haploid genome sequence for M. lychnidis-dioicae and analyzed whole transcriptome data from plant infections and other stages of the fungal lifecycle, revealing the inventory and expression level of genes that facilitate pathogenic growth. Compared to related fungi, an expanded number of major facilitator superfamily transporters and secretory lipases were detected; lipase gene expression was found to be altered by exposure to lipid compounds, which signaled a switch to dikaryotic, pathogenic growth. In addition, while enzymes to digest cellulose, xylan, xyloglucan, and highly substituted forms of pectin were absent, along with depletion of peroxidases and superoxide dismutases that protect the fungus from oxidative stress, the repertoire of glycosyltransferases and of enzymes that could manipulate host development has expanded. A total of 14 % of the genome was categorized as repetitive sequences. Transposable elements have accumulated in mating-type chromosomal regions and were also associated across the genome with gene clusters of small secreted proteins, which may mediate host interactions.

Conclusions

The unique absence of enzyme classes for plant cell wall degradation and maintenance of enzymes that break down components of pollen tubes and flowers provides a striking example of biotrophic host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1660-8) contains supplementary material, which is available to authorized users.

Contrasting host-pathogen interactions and genome evolution in two generalist and specialist microsporidian pathogens of mosquitoes

Obligate intracellular pathogens depend on their host for growth yet must also evade detection by host defenses. Here we investigate host adaptation in two Microsporidia, the specialist Edhazardia aedis and the generalist Vavraia culicis, pathogens of disease vector mosquitoes. Genomic analysis and deep RNA-Seq across infection time courses reveal fundamental differences between these pathogens. E. aedis retains enhanced cell surface modification and signalling capacity, upregulating protein trafficking and secretion dynamically during infection. V. culicis is less dependent on its host for basic metabolites and retains a subset of spliceosomal components, with a transcriptome broadly focused on growth and replication. Transcriptional profiling of mosquito immune responses reveals that response to infection by E. aedis differs dramatically depending on the mode of infection, and that antimicrobial defensins may play a general role in mosquito defense against Microsporidia. This analysis illuminates fundamentally different evolutionary paths and host interplay of specialist and generalist pathogens.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.

Comparative genomics of Australian isolates of the wheat stem rust pathogen Puccinia graminis f. sp. tritici reveals extensive polymorphism in candidate effector genes

The wheat stem rust fungus Puccinia graminis f. sp. tritici (Pgt) is one of the most destructive pathogens of wheat. In this study, a draft genome was built for a founder Australian Pgt isolate of pathotype (pt.) 21-0 (collected in 1954) by next generation DNA sequencing. A combination of reference-based assembly using the genome of the previously sequenced American Pgt isolate CDL 75-36-700-3 (p7a) and de novo assembly were performed resulting in a 92 Mbp reference genome for Pgt isolate 21-0. Approximately 13 Mbp of de novo assembled sequence in this genome is not present in the p7a reference assembly. This novel sequence is not specific to 21-0 as it is also present in three other Pgt rust isolates of independent origin. The new reference genome was subsequently used to build a pan-genome based on five Australian Pgt isolates. Transcriptomes from germinated urediniospores and haustoria were separately assembled for pt. 21-0 and comparison of gene expression profiles showed differential expression in ∼10% of the genes each in germinated spores and haustoria. A total of 1,924 secreted proteins were predicted from the 21-0 transcriptome, of which 520 were classified as haustorial secreted proteins (HSPs). Comparison of 21-0 with two presumed clonal field derivatives of this lineage (collected in 1982 and 1984) that had evolved virulence on four additional resistance genes (Sr5, Sr11, Sr27, SrSatu) identified mutations in 25 HSP effector candidates. Some of these mutations could explain their novel virulence phenotypes.

Genetic and phenotypic intra-species variation in Candida albicans

Candida albicans is a commensal fungus of the human gastrointestinal tract and a prevalent opportunistic pathogen. To examine diversity within this species, extensive genomic and phenotypic analyses were performed on 21 clinical C. albicans isolates. Genomic variation was evident in the form of polymorphisms, copy number variations, chromosomal inversions, subtelomeric hypervariation, loss of heterozygosity (LOH), and whole or partial chromosome aneuploidies. All 21 strains were diploid, although karyotypic changes were present in eight of the 21 isolates, with multiple strains being trisomic for Chromosome 4 or Chromosome 7. Aneuploid strains exhibited a general fitness defect relative to euploid strains when grown under replete conditions. All strains were also heterozygous, yet multiple, distinct LOH tracts were present in each isolate. Higher overall levels of genome heterozygosity correlated with faster growth rates, consistent with increased overall fitness. Genes with the highest rates of amino acid substitutions included many cell wall proteins, implicating fast evolving changes in cell adhesion and host interactions. One clinical isolate, P94015, presented several striking properties including a novel cellular phenotype, an inability to filament, drug resistance, and decreased virulence. Several of these properties were shown to be due to a homozygous nonsense mutation in the EFG1 gene. Furthermore, loss of EFG1 function resulted in increased fitness of P94015 in a commensal model of infection. Our analysis therefore reveals intra-species genetic and phenotypic differences in C. albicans and delineates a natural mutation that alters the balance between commensalism and pathogenicity.

Genome update of the dimorphic human pathogenic fungi causing paracoccidioidomycosis

Paracoccidiodomycosis (PCM) is a clinically important fungal disease that can acquire serious systemic forms and is caused by the thermodimorphic fungal Paracoccidioides spp. PCM is a tropical disease that is endemic in Latin America, where up to ten million people are infected; 80% of reported cases occur in Brazil, followed by Colombia and Venezuela. To enable genomic studies and to better characterize the pathogenesis of this dimorphic fungus, two reference strains of P. brasiliensis (Pb03, Pb18) and one strain of P. lutzii (Pb01) were sequenced [1]. While the initial draft assemblies were accurate in large scale structure and had high overall base quality, the sequences had frequent small scale defects such as poor quality stretches, unknown bases (N's), and artifactual deletions or nucleotide duplications, all of which caused larger scale errors in predicted gene structures. Since assembly consensus errors can now be addressed using next generation sequencing (NGS) in combination with recent methods allowing systematic assembly improvement, we re-sequenced the three reference strains of Paracoccidioides spp. using Illumina technology. We utilized the high sequencing depth to re-evaluate and improve the original assemblies generated from Sanger sequence reads, and obtained more complete and accurate reference assemblies. The new assemblies led to improved transcript predictions for the vast majority of genes of these reference strains, and often substantially corrected gene structures. These include several genes that are central to virulence or expressed during the pathogenic yeast stage in Paracoccidioides and other fungi, such as HSP90, RYP1-3, BAD1, catalase B, alpha-1,3-glucan synthase and the beta glucan synthase target gene FKS1. The improvement and validation of these reference sequences will now allow more accurate genome-based analyses. To our knowledge, this is one of the first reports of a fully automated and quality-assessed upgrade of a genome assembly and annotation for a non-model fungus.

The fungal genus Paracoccidioides is the causal agent of paracoccidioidomycosis (PCM), a neglected tropical disease that is endemic in several countries of South America. Paracoccidioides is a pathogenic dimorphic fungus that is capable of converting to a virulent yeast form after inhalation by the host. Therefore the molecular biology of the switch to the yeast phase is of particular interest for understanding the virulence of this and other human pathogenic fungi, and ultimately for reducing the morbidity and mortality caused by such fungal infections. We here present the strategy and methods we used to update and improve accuracy of three reference genome sequences of Paracoccidioides spp. utilizing state-of-the-art Illumina re-sequencing, assembly improvement, re-annotation, and quality assessment. The resulting improved genome resource should be of wide use not solely for advancing research on the genetics and molecular biology of Paracoccidioides and the closely related pathogenic species Histoplasma and Blastomyces, but also for fungal diagnostics based on sequencing or molecular assays, characterizing rapidly changing proteins that may be involved in virulence, SNP-based population analyses and other tasks that require high sequence accuracy. The genome update and underlying strategy and methods also serve as a proof of principle that could encourage similar improvements of other draft genomes.

Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3–5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.

Fungi on the skin: dermatophytes and Malassezia

Several human skin diseases and disorders are associated with two groups of fungi, the dermatophytes and Malassezia. Although these skin-related problems are not generally life threatening, they are among the most common diseases and disorders of mankind. These fungi are phylogenetically divergent, with the dermatophytes within the Ascomycota and Malassezia within Basidiomycota. Genome analysis indicates that the adaptations to the skin environment are different in these two groups of fungi. Malassezia are dependent on host lipids and secrete lipases and phospholipases that likely release host fatty acids. The dermatophytes encode multiple enzymes with potential roles in modulating host interactions: polyketide synthases, nonribosomal peptide synthetases, LysM, proteases, kinases, and pseudokinases. These two fungal groups have maximized their interactions with the host using two very different mechanisms.

Ubiquitin-mediated response to microsporidia and virus infection in C. elegans

Microsporidia comprise a phylum of over 1400 species of obligate intracellular pathogens that can infect almost all animals, but little is known about the host response to these parasites. Here we use the whole-animal host C. elegans to show an in vivo role for ubiquitin-mediated response to the microsporidian species Nematocida parisii, as well to the Orsay virus, another natural intracellular pathogen of C. elegans. We analyze gene expression of C. elegans in response to N. parisii, and find that it is similar to response to viral infection. Notably, we find an upregulation of SCF ubiquitin ligase components, such as the cullin ortholog cul-6, which we show is important for ubiquitin targeting of N. parisii cells in the intestine. We show that ubiquitylation components, the proteasome, and the autophagy pathway are all important for defense against N. parisii infection. We also find that SCF ligase components like cul-6 promote defense against viral infection, where they have a more robust role than against N. parisii infection. This difference may be due to suppression of the host ubiquitylation system by N. parisii: when N. parisii is crippled by anti-microsporidia drugs, the host can more effectively target pathogen cells for ubiquitylation. Intriguingly, inhibition of the ubiquitin-proteasome system (UPS) increases expression of infection-upregulated SCF ligase components, indicating that a trigger for transcriptional response to intracellular infection by N. parisii and virus may be perturbation of the UPS. Altogether, our results demonstrate an in vivo role for ubiquitin-mediated defense against microsporidian and viral infections in C. elegans.

Standardized metadata for human pathogen/vector genomic sequences

High throughput sequencing has accelerated the determination of genome sequences for thousands of human infectious disease pathogens and dozens of their vectors. The scale and scope of these data are enabling genotype-phenotype association studies to identify genetic determinants of pathogen virulence and drug/insecticide resistance, and phylogenetic studies to track the origin and spread of disease outbreaks. To maximize the utility of genomic sequences for these purposes, it is essential that metadata about the pathogen/vector isolate characteristics be collected and made available in organized, clear, and consistent formats. Here we report the development of the GSCID/BRC Project and Sample Application Standard, developed by representatives of the Genome Sequencing Centers for Infectious Diseases (GSCIDs), the Bioinformatics Resource Centers (BRCs) for Infectious Diseases, and the U.S. National Institute of Allergy and Infectious Diseases (NIAID), part of the National Institutes of Health (NIH), informed by interactions with numerous collaborating scientists. It includes mapping to terms from other data standards initiatives, including the Genomic Standards Consortium's minimal information (MIxS) and NCBI's BioSample/BioProjects checklists and the Ontology for Biomedical Investigations (OBI). The standard includes data fields about characteristics of the organism or environmental source of the specimen, spatial-temporal information about the specimen isolation event, phenotypic characteristics of the pathogen/vector isolated, and project leadership and support. By modeling metadata fields into an ontology-based semantic framework and reusing existing ontologies and minimum information checklists, the application standard can be extended to support additional project-specific data fields and integrated with other data represented with comparable standards. The use of this metadata standard by all ongoing and future GSCID sequencing projects will provide a consistent representation of these data in the BRC resources and other repositories that leverage these data, allowing investigators to identify relevant genomic sequences and perform comparative genomics analyses that are both statistically meaningful and biologically relevant.

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.

Using transcription of six Puccinia triticina races to identify the effective secretome during infection of wheat

Wheat leaf rust, caused by the basidiomycete Puccinia triticina, can cause yield losses of up to 20% in wheat producing regions. During infection, the fungus forms haustoria that secrete proteins into the plant cell and effect changes in plant transcription, metabolism, and defense. It is hypothesized that new races emerge as a result of overcoming plant resistance via changes in the secreted effector proteins. To understand gene expression during infection and find genetic differences associated with races, RNA from wheat leaves infected with six different rust races, at 6 days post inoculation, was sequenced using Illumina. As P. triticina is an obligate biotroph, RNA from both the host and fungi were present and separated by alignment to the P. triticina genome and a wheat EST reference. A total of 222,571 rust contigs were assembled from 165 million reads. An examination of the resulting contigs revealed 532 predicted secreted proteins among the transcripts. Of these, 456 were found in all races. Fifteen genes were found with amino acid changes, corresponding to putative avirulence effectors potentially recognized by 11 different leaf rust resistance (Lr) genes. Twelve of the potential avirulence effectors have no homology to known genes. One gene had significant similarity to cerato-platanin, a known fungal elicitor, and another showed similarity to fungal tyrosinase, an enzyme involved in melanin synthesis. Temporal expression profiles were developed for these genes by qRT-PCR and show that the genes expression patterns were consistent between races from infection initiation to just prior to spore eruption.

One fungus, one name: defining the genus Fusarium in a scientifically robust way that preserves longstanding use

In this letter, we advocate recognizing the genus Fusarium as the sole name for a group that includes virtually all Fusarium species of importance in plant pathology, mycotoxicology, medicine, and basic research. This phylogenetically guided circumscription will free scientists from any obligation to use other genus names, including teleomorphs, for species nested within this clade, and preserve the application of the name Fusarium in the way it has been used for almost a century. Due to recent changes in the International Code of Nomenclature for algae, fungi, and plants, this is an urgent matter that requires community attention. The alternative is to break the longstanding concept of Fusarium into nine or more genera, and remove important taxa such as those in the F. solani species complex from the genus, a move we believe is unnecessary. Here we present taxonomic and nomenclatural proposals that will preserve established research connections and facilitate communication within and between research communities, and at the same time support strong scientific principles and good taxonomic practice.

Isolation of Blastomyces dermatitidis yeast from lung tissue during murine infection for in vivo transcriptional profiling

Blastomyces dermatitidis belongs to a group of thermally dimorphic fungi that grow as sporulating mold in the soil and convert to pathogenic yeast in the lung following inhalation of spores. Knowledge about the molecular events important for fungal adaptation and survival in the host remains limited. The development of high-throughput analytic tools such as RNA sequencing (RNA-Seq) has potential to provide novel insight on fungal pathogenesis especially if applied in vivo during infection. However, in vivo transcriptional profiling is hindered by the low abundance of fungal cells relative to mammalian tissue and difficulty in isolating fungal cells from the tissues they infect. For the purpose of obtaining B. dermatitidis RNA for in vivo transcriptional analysis by RNA-Seq, we developed a simple technique for isolating yeast from murine lung tissue. Using a two-step approach of filtration and centrifugation following lysis of murine lung cells, 91% of yeast cells causing infection were isolated from lung tissue. B. dermatitidis recovered from the lung yielded high-quality RNA with minimal murine contamination and was suitable for RNA-Seq. Approximately 87% of the sequencing reads obtained from the recovered yeast aligned with the B. dermatitidis genome. This was similar to 93% alignment for yeast grown in vitro. The use of near-freezing temperature along with short ex vivo time minimized transcriptional changes that would have otherwise occurred with higher temperature or longer processing time. In conclusion, we have developed a technique that recovers the majority of yeast causing pulmonary infection and yields high-quality fungal RNA with minimal contamination by mammalian RNA.

Sequencing and characterization of the complete mitochondrial genomes of three Pneumocystis species provide new insights into divergence between human and rodent Pneumocystis

Pneumocystis jirovecii is an important opportunistic pathogen associated with AIDS and other immunodeficient conditions. Currently, very little is known about its nuclear and mitochondrial genomes. In this study, we sequenced the complete mitochondrial genome (mtDNA) of this organism and its closely related species Pneumocystis carinii and Pneumocystis murina by a combination of sequencing technologies. Our study shows that P. carinii and P. murina mtDNA share a nearly identical number and order of genes in a linear configuration, whereas P. jirovecii has a circular mtDNA containing nearly the same set of genes but in a different order. Detailed studies of the mtDNA terminal structures of P. murina and P. carinii suggest a unique replication mechanism for linear mtDNA. Phylogenetic analysis supports a close association of Pneumocystis species with Taphrina, Saitoella, and Schizosaccharomyces, and divergence within Pneumocystis species, with P. murina and P. carinii being more closely related to each other than either is to P. jirovecii. Comparative analysis of four complete P. jirovecii mtDNA sequences in this study and previously reported mtDNA sequences for diagnosing and genotyping suggests that the current diagnostic and typing methods can be improved using the complete mtDNA data. The availability of the complete P. jirovecii mtDNA also opens the possibility of identifying new therapeutic targets.

Conserved loci of leaf and stem rust fungi of wheat share synteny interrupted by lineage-specific influx of repeat elements

Background

Wheat leaf rust (Puccinia triticina Eriks; Pt) and stem rust fungi (P. graminis f.sp. tritici; Pgt) are significant economic pathogens having similar host ranges and life cycles, but different alternate hosts. The Pt genome, currently estimated at 135 Mb, is significantly larger than Pgt, at 88 Mb, but the reason for the expansion is unknown. Three genomic loci of Pt conserved proteins were characterized to gain insight into gene content, genome complexity and expansion.

Results

A bacterial artificial chromosome (BAC) library was made from P. triticina race 1, BBBD and probed with Pt homologs of genes encoding two predicted Pgt secreted effectors and a DNA marker mapping to a region of avirulence. Three BACs, 103 Kb, 112 Kb, and 166 Kb, were sequenced, assembled, and open reading frames were identified. Orthologous genes were identified in Pgt and local conservation of gene order (microsynteny) was observed. Pairwise protein identities ranged from 26 to 99%. One Pt BAC, containing a RAD18 ortholog, shares syntenic regions with two Pgt scaffolds, which could represent both haplotypes of Pgt. Gene sequence is diverged between the species as well as within the two haplotypes. In all three BAC clones, gene order is locally conserved, however, gene shuffling has occurred relative to Pgt. These regions are further diverged by differing insertion loci of LTR-retrotransposon, Gypsy, Copia, Mutator, and Harbinger mobile elements. Uncharacterized Pt open reading frames were also found; these proteins are high in lysine and similar to multiple proteins in Pgt.

Conclusions

The three Pt loci are conserved in gene order, with a range of gene sequence divergence. Conservation of predicted haustoria expressed secreted protein genes between Pt and Pgt is extended to the more distant poplar rust, Melampsora larici-populina. The loci also reveal that genome expansion in Pt is in part due to higher occurrence of repeat-elements in this species.

Microsporidian genome analysis reveals evolutionary strategies for obligate intracellular growth

Microsporidia comprise a large phylum of obligate intracellular eukaryotes that are fungal-related parasites responsible for widespread disease, and here we address questions about microsporidia biology and evolution. We sequenced three microsporidian genomes from two species, Nematocida parisii and Nematocida sp1, which are natural pathogens of Caenorhabditis nematodes and provide model systems for studying microsporidian pathogenesis. We performed deep sequencing of transcripts from a time course of N. parisii infection. Examination of pathogen gene expression revealed compact transcripts and a dramatic takeover of host cells by Nematocida. We also performed phylogenomic analyses of Nematocida and other microsporidian genomes to refine microsporidian phylogeny and identify evolutionary events of gene loss, acquisition, and modification. In particular, we found that all microsporidia lost the tumor-suppressor gene retinoblastoma, which we speculate could accelerate the parasite cell cycle and increase the mutation rate. We also found that microsporidia acquired transporters that could import nucleosides to fuel rapid growth. In addition, microsporidian hexokinases gained secretion signal sequences, and in a functional assay these were sufficient to export proteins out of the cell; thus hexokinase may be targeted into the host cell to reprogram it toward biosynthesis. Similar molecular changes appear during formation of cancer cells and may be evolutionary strategies adopted independently by microsporidia to proliferate rapidly within host cells. Finally, analysis of genome polymorphisms revealed evidence for a sexual cycle that may provide genetic diversity to alleviate problems caused by clonal growth. Together these events may explain the emergence and success of these diverse intracellular parasites.

Approaches to Fungal Genome Annotation

Fungal genome annotation is the starting point for analysis of genome content. This generally involves the application of diverse methods to identify features on a genome assembly such as protein-coding and non-coding genes, repeats and transposable elements, and pseudogenes. Here we describe tools and methods leveraged for eukaryotic genome annotation with a focus on the annotation of fungal nuclear and mitochondrial genomes. We highlight the application of the latest technologies and tools to improve the quality of predicted gene sets. The Broad Institute eukaryotic genome annotation pipeline is described as one example of how such methods and tools are integrated into a sequencing center's production genome annotation environment.

Comparative genomic analysis of human fungal pathogens causing paracoccidioidomycosis

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

Comparative genomics yields insights into niche adaptation of plant vascular wilt pathogens

The vascular wilt fungi Verticillium dahliae and V. albo-atrum infect over 200 plant species, causing billions of dollars in annual crop losses. The characteristic wilt symptoms are a result of colonization and proliferation of the pathogens in the xylem vessels, which undergo fluctuations in osmolarity. To gain insights into the mechanisms that confer the organisms' pathogenicity and enable them to proliferate in the unique ecological niche of the plant vascular system, we sequenced the genomes of V. dahliae and V. albo-atrum and compared them to each other, and to the genome of Fusarium oxysporum, another fungal wilt pathogen. Our analyses identified a set of proteins that are shared among all three wilt pathogens, and present in few other fungal species. One of these is a homolog of a bacterial glucosyltransferase that synthesizes virulence-related osmoregulated periplasmic glucans in bacteria. Pathogenicity tests of the corresponding V. dahliae glucosyltransferase gene deletion mutants indicate that the gene is required for full virulence in the Australian tobacco species Nicotiana benthamiana. Compared to other fungi, the two sequenced Verticillium genomes encode more pectin-degrading enzymes and other carbohydrate-active enzymes, suggesting an extraordinary capacity to degrade plant pectin barricades. The high level of synteny between the two Verticillium assemblies highlighted four flexible genomic islands in V. dahliae that are enriched for transposable elements, and contain duplicated genes and genes that are important in signaling/transcriptional regulation and iron/lipid metabolism. Coupled with an enhanced capacity to degrade plant materials, these genomic islands may contribute to the expanded genetic diversity and virulence of V. dahliae, the primary causal agent of Verticillium wilts. Significantly, our study reveals insights into the genetic mechanisms of niche adaptation of fungal wilt pathogens, advances our understanding of the evolution and development of their pathogenesis, and sheds light on potential avenues for the development of novel disease management strategies to combat destructive wilt diseases.

SREB, a GATA transcription factor that directs disparate fates in Blastomyces dermatitidis including morphogenesis and siderophore biosynthesis

Blastomyces dermatitidis belongs to a group of human pathogenic fungi that exhibit thermal dimorphism. At 22 degrees C, these fungi grow as mold that produce conidia or infectious particles, whereas at 37 degrees C they convert to budding yeast. The ability to switch between these forms is essential for virulence in mammals and may enable these organisms to survive in the soil. To identify genes that regulate this phase transition, we used Agrobacterium tumefaciens to mutagenize B. dermatitidis conidia and screened transformants for defects in morphogenesis. We found that the GATA transcription factor SREB governs multiple fates in B. dermatitidis: phase transition from yeast to mold, cell growth at 22 degrees C, and biosynthesis of siderophores under iron-replete conditions. Insertional and null mutants fail to convert to mold, do not accumulate significant biomass at 22 degrees C, and are unable to suppress siderophore biosynthesis under iron-replete conditions. The defect in morphogenesis in the SREB mutant was independent of exogenous iron concentration, suggesting that SREB promotes the phase transition by altering the expression of genes that are unrelated to siderophore biosynthesis. Using bioinformatic and gene expression analyses, we identified candidate genes with upstream GATA sites whose expression is altered in the null mutant that may be direct or indirect targets of SREB and promote the phase transition. We conclude that SREB functions as a transcription factor that promotes morphogenesis and regulates siderophore biosynthesis. To our knowledge, this is the first gene identified that promotes the conversion from yeast to mold in the dimorphic fungi, and may shed light on environmental persistence of these pathogens.

Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium

Fusarium species are among the most important phytopathogenic and toxigenic fungi. To understand the molecular underpinnings of pathogenicity in the genus Fusarium, we compared the genomes of three phenotypically diverse species: Fusarium graminearum, Fusarium verticillioides and Fusarium oxysporum f. sp. lycopersici. Our analysis revealed lineage-specific (LS) genomic regions in F. oxysporum that include four entire chromosomes and account for more than one-quarter of the genome. LS regions are rich in transposons and genes with distinct evolutionary profiles but related to pathogenicity, indicative of horizontal acquisition. Experimentally, we demonstrate the transfer of two LS chromosomes between strains of F. oxysporum, converting a non-pathogenic strain into a pathogen. Transfer of LS chromosomes between otherwise genetically isolated strains explains the polyphyletic origin of host specificity and the emergence of new pathogenic lineages in F. oxysporum. These findings put the evolution of fungal pathogenicity into a new perspective.

The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization

We sequenced and annotated the genome of the filamentous fungus Fusarium graminearum, a major pathogen of cultivated cereals. Very few repetitive sequences were detected, and the process of repeat-induced point mutation, in which duplicated sequences are subject to extensive mutation, may partially account for the reduced repeat content and apparent low number of paralogous (ancestrally duplicated) genes. A second strain of F. graminearum contained more than 10,000 single-nucleotide polymorphisms, which were frequently located near telomeres and within other discrete chromosomal segments. Many highly polymorphic regions contained sets of genes implicated in plant-fungus interactions and were unusually divergent, with higher rates of recombination. These regions of genome innovation may result from selection due to interactions of F. graminearum with its plant hosts.

Comparative proteomics of extracellular proteinsin vitro andin planta from the pathogenic fungusFusarium graminearum

High-throughput MS/MS was used to identify proteins secreted by Fusarium graminearum (Gibberella zeae) during growth on 13 media in vitro and in planta during infection of wheat heads. In vitro secreted proteins were collected from the culture filtrates, and in planta proteins were collected by vacuum infiltration. A total of 289 proteins (229 in vitro and 120 in planta) were identified with high statistical confidence. Forty-nine of the in planta proteins were not found in any of the in vitro conditions. The majority (91-100%) of the in vitro proteins had predicted signal peptides, but only 56% of the in planta proteins. At least 13 of the nonsecreted proteins found only in planta were single-copy housekeeping enzymes, including enolase, triose phosphate isomerase, phosphoglucomutase, calmodulin, aconitase, and malate dehydrogenase. The presence of these proteins in the in planta but not in vitro secretome might indicate that significant fungal lysis occurs during pathogenesis. On the other hand, several of the proteins lacking signal peptides that were found in planta have been reported to be potent immunogens secreted by animal pathogenic fungi, and therefore could be important in the interaction between F. graminearum and its host plants.

DNA sequence of human chromosome 17 and analysis of rearrangement in the human lineage

Chromosome 17 is unusual among the human chromosomes in many respects. It is the largest human autosome with orthology to only a single mouse chromosome, mapping entirely to the distal half of mouse chromosome 11. Chromosome 17 is rich in protein-coding genes, having the second highest gene density in the genome. It is also enriched in segmental duplications, ranking third in density among the autosomes. Here we report a finished sequence for human chromosome 17, as well as a structural comparison with the finished sequence for mouse chromosome 11, the first finished mouse chromosome. Comparison of the orthologous regions reveals striking differences. In contrast to the typical pattern seen in mammalian evolution, the human sequence has undergone extensive intrachromosomal rearrangement, whereas the mouse sequence has been remarkably stable. Moreover, although the human sequence has a high density of segmental duplication, the mouse sequence has a very low density. Notably, these segmental duplications correspond closely to the sites of structural rearrangement, demonstrating a link between duplication and rearrangement. Examination of the main classes of duplicated segments provides insight into the dynamics underlying expansion of chromosome-specific, low-copy repeats in the human genome.

Analysis of the DNA sequence and duplication history of human chromosome 15

Here we present a finished sequence of human chromosome 15, together with a high-quality gene catalogue. As chromosome 15 is one of seven human chromosomes with a high rate of segmental duplication, we have carried out a detailed analysis of the duplication structure of the chromosome. Segmental duplications in chromosome 15 are largely clustered in two regions, on proximal and distal 15q; the proximal region is notable because recombination among the segmental duplications can result in deletions causing Prader-Willi and Angelman syndromes. Sequence analysis shows that the proximal and distal regions of 15q share extensive ancient similarity. Using a simple approach, we have been able to reconstruct many of the events by which the current duplication structure arose. We find that most of the intrachromosomal duplications seem to share a common ancestry. Finally, we demonstrate that some remaining gaps in the genome sequence are probably due to structural polymorphisms between haplotypes; this may explain a significant fraction of the gaps remaining in the human genome.