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.

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.

Tracing Genetic Exchange and Biogeography of Cryptococcus neoformans var. grubii at the Global Population Level.

Cryptococcus neoformans var. grubii is the causative agent of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals, typically human immunodeficiency virus/AIDS patients from developing countries. Despite the worldwide emergence of this ubiquitous infection, little is known about the global molecular epidemiology of this fungal pathogen. Here we sequence the genomes of 188 diverse isolates and characterize the major subdivisions, their relative diversity, and the level of genetic exchange between them. While most isolates of C. neoformans var. grubii belong to one of three major lineages (VNI, VNII, and VNB), some haploid isolates show hybrid ancestry including some that appear to have recently interbred, based on the detection of large blocks of each ancestry across each chromosome. Many isolates display evidence of aneuploidy, which was detected for all chromosomes. In diploid isolates of C. neoformans var. grubii (serotype AA) and of hybrids with C. neoformans var. neoformans (serotype AD) such aneuploidies have resulted in loss of heterozygosity, where a chromosomal region is represented by the genotype of only one parental isolate. Phylogenetic and population genomic analyses of isolates from Brazil reveal that the previously "African" VNB lineage occurs naturally in the South American environment. This suggests migration of the VNB lineage between Africa and South America prior to its diversification, supported by finding ancestral recombination events between isolates from different lineages and regions. The results provide evidence of substantial population structure, with all lineages showing multi-continental distributions; demonstrating the highly dispersive nature of this pathogen.

Comparative genomic and transcriptomic analysis of wangiella dermatitidis, a major cause of phaeohyphomycosis and a model black yeast human pathogen.

Black or dark brown (phaeoid) fungi cause cutaneous, subcutaneous, and systemic infections in humans. Black fungi thrive in stressful conditions such as intense light, high radiation, and very low pH. Wangiella (Exophiala) dermatitidis is arguably the most studied phaeoid fungal pathogen of humans. Here, we report our comparative analysis of the genome of W. dermatitidis and the transcriptional response to low pH stress. This revealed that W. dermatitidis has lost the ability to synthesize alpha-glucan, a cell wall compound many pathogenic fungi use to evade the host immune system. In contrast, W. dermatitidis contains a similar profile of chitin synthase genes as related fungi and strongly induces genes involved in cell wall synthesis in response to pH stress. The large portfolio of transporters may provide W. dermatitidis with an enhanced ability to remove harmful products as well as to survive on diverse nutrient sources. The genome encodes three independent pathways for producing melanin, an ability linked to pathogenesis; these are active during pH stress, potentially to produce a barrier to accumulated oxidative damage that might occur under stress conditions. In addition, a full set of fungal light-sensing genes is present, including as part of a carotenoid biosynthesis gene cluster. Finally, we identify a two-gene cluster involved in nucleotide sugar metabolism conserved with a subset of fungi and characterize a horizontal transfer event of this cluster between fungi and algal viruses. This work reveals how W. dermatitidis has adapted to stress and survives in diverse environments, including during human infections.

Analysis of a food-borne fungal pathogen outbreak: virulence and genome of a Mucor circinelloides isolate from yogurt.

Food-borne pathogens are ongoing problems, and new pathogens are emerging. The impact of fungi, however, is largely underestimated. Recently, commercial yogurts contaminated with Mucor circinelloides were sold, and >200 consumers became ill with nausea, vomiting, and diarrhea. Mucoralean fungi cause the fatal fungal infection mucormycosis, whose incidence has been continuously increasing. In this study, we isolated an M. circinelloides strain from a yogurt container, and multilocus sequence typing identified the strain as Mucor circinelloides f. circinelloides. M. circinelloides f. circinelloides is the most virulent M. circinelloides subspecies and is commonly associated with human infections, whereas M. circinelloides f. lusitanicus and M. circinelloides f. griseocyanus are less common causes of infection. Whole-genome analysis of the yogurt isolate confirmed it as being close to the M. circinelloides f. circinelloides subgroup, with a higher percentage of divergence with the M. circinelloides f. lusitanicus subgroup. In mating assays, the yogurt isolate formed sexual zygospores with the (-) M. circinelloides f. circinelloides tester strain, which is congruent with its sex locus encoding SexP, the (+) mating type sex determinant. The yogurt isolate was virulent in murine and wax moth larva host systems. In a murine gastromucormycosis model, Mucor was recovered from fecal samples of infected mice for up to 10 days, indicating that Mucor can survive transit through the GI tract. In interactions with human immune cells, M. circinelloides f. lusitanicus induced proinflammatory cytokines but M. circinelloides f. circinelloides did not, which may explain the different levels of virulence in mammalian hosts. This study demonstrates that M. circinelloides can spoil food products and cause gastrointestinal illness in consumers and may pose a particular risk to immunocompromised patients. Importance: The U.S. FDA reported that yogurt products were contaminated with M. circinelloides, a mucoralean fungal pathogen, and >200 consumers complained of symptoms, including vomiting, nausea, and diarrhea. The manufacturer voluntarily withdrew the affected yogurt products from the market. Compared to other food-borne pathogens, including bacteria, viruses, and parasites, less focus has been placed on the risk of fungal pathogens. This study evaluates the potential risk from the food-borne fungal pathogen M. circinelloides that was isolated from the contaminated commercial yogurt. We successfully cultured an M. circinelloides isolate and found that the isolate belongs to the species M. circinelloides f. circinelloides, which is often associated with human infections. In murine and insect host models, the isolate was virulent. While information disseminated in the popular press would suggest this fungal contaminant poses little or no risk to consumers, our results show instead that it is capable of causing significant infections in animals.

The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum.

We sequenced and annotated the genomes of four P. vivax strains collected from disparate geographic locations, tripling the number of genome sequences available for this understudied parasite and providing the first genome-wide perspective of global variability in this species. We observe approximately twice as much SNP diversity among these isolates as we do among a comparable collection of isolates of P. falciparum, a malaria-causing parasite that results in higher mortality. This indicates a distinct history of global colonization and/or a more stable demographic history for P. vivax relative to P. falciparum, which is thought to have undergone a recent population bottleneck. The SNP diversity, as well as additional microsatellite and gene family variability, suggests a capacity for greater functional variation in the global population of P. vivax. These findings warrant a deeper survey of variation in P. vivax to equip disease interventions targeting the distinctive biology of this neglected but major pathogen.

A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries.

Genome targeting methods enable cost-effective capture of specific subsets of the genome for sequencing. We present here an automated, highly scalable method for carrying out the Solution Hybrid Selection capture approach that provides a dramatic increase in scale and throughput of sequence-ready libraries produced. Significant process improvements and a series of in-process quality control checkpoints are also added. These process improvements can also be used in a manual version of the protocol.

Comparative functional genomics of the fission yeasts.

The fission yeast clade--comprising Schizosaccharomyces pombe, S. octosporus, S. cryophilus, and S. japonicus--occupies the basal branch of Ascomycete fungi and is an important model of eukaryote biology. A comparative annotation of these genomes identified a near extinction of transposons and the associated innovation of transposon-free centromeres. Expression analysis established that meiotic genes are subject to antisense transcription during vegetative growth, which suggests a mechanism for their tight regulation. In addition, trans-acting regulators control new genes within the context of expanded functional modules for meiosis and stress response. Differences in gene content and regulation also explain why, unlike the budding yeast of Saccharomycotina, fission yeasts cannot use ethanol as a primary carbon source. These analyses elucidate the genome structure and gene regulation of fission yeast and provide tools for investigation across the Schizosaccharomyces clade.

A catalog of reference genomes from the human microbiome.

The human microbiome refers to the community of microorganisms, including prokaryotes, viruses, and microbial eukaryotes, that populate the human body. The National Institutes of Health launched an initiative that focuses on describing the diversity of microbial species that are associated with health and disease. The first phase of this initiative includes the sequencing of hundreds of microbial reference genomes, coupled to metagenomic sequencing from multiple body sites. Here we present results from an initial reference genome sequencing of 178 microbial genomes. From 547,968 predicted polypeptides that correspond to the gene complement of these strains, previously unidentified ("novel") polypeptides that had both unmasked sequence length greater than 100 amino acids and no BLASTP match to any nonreference entry in the nonredundant subset were defined. This analysis resulted in a set of 30,867 polypeptides, of which 29,987 (approximately 97%) were unique. In addition, this set of microbial genomes allows for approximately 40% of random sequences from the microbiome of the gastrointestinal tract to be associated with organisms based on the match criteria used. Insights into pan-genome analysis suggest that we are still far from saturating microbial species genetic data sets. In addition, the associated metrics and standards used by our group for quality assurance are presented.