Biosynthetic gene clusters and the evolution of fungal chemodiversity

Covering: up to 2019Fungi produce a remarkable diversity of secondary metabolites: small, bioactive molecules not required for growth but which are essential to their ecological interactions with other organisms. Genes that participate in the same secondary metabolic pathway typically reside next to each other in fungal genomes and form biosynthetic gene clusters (BGCs). By synthesizing state-of-the-art knowledge on the evolution of BGCs in fungi, we propose that fungal chemodiversity stems from three molecular evolutionary processes involving BGCs: functional divergence, horizontal transfer, and de novo assembly. We provide examples of how these processes have contributed to the generation of fungal chemodiversity, discuss their relative importance, and outline major, outstanding questions in the field.

The Aspergillus fumigatus transcription factor RglT is important for gliotoxin biosynthesis and self-protection, and virulence

Aspergillus fumigatus is an opportunistic fungal pathogen that secretes an array of immune-modulatory molecules, including secondary metabolites (SMs), which contribute to enhancing fungal fitness and growth within the mammalian host. Gliotoxin (GT) is a SM that interferes with the function and recruitment of innate immune cells, which are essential for eliminating A. fumigatus during invasive infections. We identified a C6 Zn cluster-type transcription factor (TF), subsequently named RglT, important for A. fumigatus oxidative stress resistance, GT biosynthesis and self-protection. RglT regulates the expression of several gli genes of the GT biosynthetic gene cluster, including the oxidoreductase-encoding gene gliT, by directly binding to their respective promoter regions. Subsequently, RglT was shown to be important for virulence in a chemotherapeutic murine model of invasive pulmonary aspergillosis (IPA). Homologues of RglT and GliT are present in eurotiomycete and sordariomycete fungi, including the non-GT-producing fungus A. nidulans, where a conservation of function was described. Phylogenetically informed model testing led to an evolutionary scenario in which the GliT-based resistance mechanism is ancestral and RglT-mediated regulation of GliT occurred subsequently. In conclusion, this work describes the function of a previously uncharacterised TF in oxidative stress resistance, GT biosynthesis and self-protection in both GT-producing and non-producing Aspergillus species.

Genomic and Phenotypic Heterogeneity of Clinical Isolates of the Human Pathogens Aspergillus fumigatus, Aspergillus lentulus, and Aspergillus fumigatiaffinis

Fungal pathogens are a global threat to human health. For example, fungi from the genus Aspergillus cause a spectrum of diseases collectively known as aspergillosis. Most of the >200,000 life-threatening aspergillosis infections per year worldwide are caused by Aspergillus fumigatus. Recently, molecular typing techniques have revealed that aspergillosis can also be caused by organisms that are phenotypically similar to A. fumigatus but genetically distinct, such as Aspergillus lentulus and Aspergillus fumigatiaffinis. Importantly, some of these so-called cryptic species are thought to exhibit different virulence and drug susceptibility profiles than A. fumigatus, however, our understanding of their biology and pathogenic potential has been stymied by the lack of genome sequences and phenotypic profiling of multiple clinical strains. To fill this gap, we phenotypically characterized the virulence and drug susceptibility of 15 clinical strains of A. fumigatus, A. lentulus, and A. fumigatiaffinis from Spain and sequenced their genomes. We found heterogeneity in drug susceptibility across species and strains. We further found heterogeneity in virulence within each species but no significant differences in the virulence profiles between the three species. Genes known to influence drug susceptibility (cyp51A and fks1) vary in paralog number and sequence among these species and strains and correlate with differences in drug susceptibility. Similarly, genes known to be important for virulence in A. fumigatus showed variability in number of paralogs across strains and across species. Characterization of the genomic similarities and differences of clinical strains of A. lentulus, A. fumigatiaffinis, and A. fumigatus that vary in disease-relevant traits will advance our understanding of the variance in pathogenicity between Aspergillus species and strains that are collectively responsible for the vast majority of aspergillosis infections in humans.

Functional Characterization of Clinical Isolates of the Opportunistic Fungal Pathogen Aspergillus nidulans

Aspergillus nidulans is an opportunistic fungal pathogen in patients with immunodeficiency, and virulence of A. nidulans isolates has mainly been studied in the context of chronic granulomatous disease (CGD), with characterization of clinical isolates obtained from non-CGD patients remaining elusive. This study therefore carried out a detailed biological characterization of two A. nidulans clinical isolates (CIs), obtained from a patient with breast carcinoma and pneumonia and from a patient with cystic fibrosis that underwent lung transplantation, and compared them to the reference, nonclinical FGSC A4 strain. Both CIs presented increased growth in comparison to that of the reference strain in the presence of physiologically relevant carbon sources. Metabolomic analyses showed that the three strains are metabolically very different from each other in these carbon sources. Furthermore, the CIs were highly susceptible to cell wall-perturbing agents but not to other physiologically relevant stresses. Genome analyses identified several frameshift variants in genes encoding cell wall integrity (CWI) signaling components. Significant differences in CWI signaling were confirmed by Western blotting among the three strains. In vivo virulence studies using several different models revealed that strain MO80069 had significantly higher virulence in hosts with impaired neutrophil function than the other strains. In summary, this study presents detailed biological characterization of two A. nidulanssensu stricto clinical isolates. Just as in Aspergillus fumigatus, strain heterogeneity exists in A. nidulans clinical strains that can define virulence traits. Further studies are required to fully characterize A. nidulans strain-specific virulence traits and pathogenicity.IMPORTANCE Immunocompromised patients are susceptible to infections with opportunistic filamentous fungi from the genus Aspergillus Although A. fumigatus is the main etiological agent of Aspergillus species-related infections, other species, such as A. nidulans, are prevalent in a condition-specific manner. A. nidulans is a predominant infective agent in patients suffering from chronic granulomatous disease (CGD). A. nidulans isolates have mainly been studied in the context of CGD although infection with A. nidulans also occurs in non-CGD patients. This study carried out a detailed biological characterization of two non-CGD A. nidulans clinical isolates and compared the results to those with a reference strain. Phenotypic, metabolomic, and genomic analyses highlight fundamental differences in carbon source utilization, stress responses, and maintenance of cell wall integrity among the strains. One clinical strain had increased virulence in models with impaired neutrophil function. Just as in A. fumigatus, strain heterogeneity exists in A. nidulans clinical strains that can define virulence traits.

Into the wild: new yeast genomes from natural environments and new tools for their analysis

Genomic studies of yeasts from the wild have increased considerably in the past few years. This revolution has been fueled by advances in high-throughput sequencing technologies and a better understanding of yeast ecology and phylogeography, especially for biotechnologically important species. The present review aims to first introduce new bioinformatic tools available for the generation and analysis of yeast genomes. We also assess the accumulated genomic data of wild isolates of industrially relevant species, such as Saccharomyces spp., which provide unique opportunities to further investigate the domestication processes associated with the fermentation industry and opportunistic pathogenesis. The availability of genome sequences of other less conventional yeasts obtained from the wild has also increased substantially, including representatives of the phyla Ascomycota (e.g. Hanseniaspora) and Basidiomycota (e.g. Phaffia). Here, we review salient examples of both fundamental and applied research that demonstrate the importance of continuing to sequence and analyze genomes of wild yeasts.

Gliotoxin, a Known Virulence Factor in the Major Human Pathogen Aspergillus fumigatus, Is Also Biosynthesized by Its Nonpathogenic Relative Aspergillus fischeri

Aspergillus fumigatus is a major opportunistic human pathogen. Multiple traits contribute to A. fumigatus pathogenicity, including its ability to produce specific secondary metabolites, such as gliotoxin. Gliotoxin is known to inhibit the host immune response, and genetic mutants that inactivate gliotoxin biosynthesis (or secondary metabolism in general) attenuate A. fumigatus virulence. The genome of Aspergillus fischeri, a very close nonpathogenic relative of A. fumigatus, contains a biosynthetic gene cluster that is homologous to the A. fumigatus gliotoxin cluster. However, A. fischeri is not known to produce gliotoxin. To gain further insight into the similarities and differences between the major pathogen A. fumigatus and the nonpathogen A. fischeri, we examined whether A. fischeri strain NRRL 181 biosynthesizes gliotoxin and whether the production of secondary metabolites influences the virulence profile of A. fischeri We found that A. fischeri biosynthesizes gliotoxin under the same conditions as A. fumigatus However, whereas loss of laeA, a master regulator of secondary metabolite production (including gliotoxin biosynthesis), has previously been shown to reduce A. fumigatus virulence, we found that laeA loss (and loss of secondary metabolite production) in A. fischeri does not influence its virulence. These results suggest that LaeA-regulated secondary metabolites are virulence factors in the genomic and phenotypic background of the major pathogen A. fumigatus but are much less important in the background of the nonpathogen A. fischeri Understanding the observed spectrum of pathogenicity across closely related pathogenic and nonpathogenic Aspergillus species will require detailed characterization of their biological, chemical, and genomic similarities and differences.IMPORTANCEAspergillus fumigatus is a major opportunistic fungal pathogen of humans, but most of its close relatives are nonpathogenic. Why is that so? This important, yet largely unanswered, question can be addressed by examining how A. fumigatus and its close nonpathogenic relatives are similar or different with respect to virulence-associated traits. We investigated whether Aspergillus fischeri, a nonpathogenic close relative of A. fumigatus, can produce gliotoxin, a mycotoxin known to contribute to A. fumigatus virulence. We discovered that the nonpathogenic A. fischeri produces gliotoxin under the same conditions as those of the major pathogen A. fumigatus However, we also discovered that, in contrast to what has previously been observed in A. fumigatus, the loss of secondary metabolite production in A. fischeri does not alter its virulence. Our results are consistent with the "cards of virulence" model of opportunistic fungal disease, in which the ability to cause disease stems from the combination ("hand") of virulence factors ("cards") but not from individual factors per se.

Rapid Phenotypic and Metabolomic Domestication of Wild Penicillium Molds on Cheese

Fermented foods provide novel ecological opportunities for natural populations of microbes to evolve through successive recolonization of resource-rich substrates. Comparative genomic data have reconstructed the evolutionary histories of microbes adapted to food environments, but experimental studies directly demonstrating the process of domestication are lacking for most fermented food microbes. Here, we show that during adaptation to cheese, phenotypic and metabolomic traits of wild Penicillium molds rapidly change to produce domesticated phenotypes with properties similar to those of the industrial cultures used to make Camembert and other bloomy rind cheeses. Over a period of just a few weeks, populations of wild Penicillium strains serially passaged on cheese had reduced pigment, spore, and mycotoxin production. Domesticated strains also had a striking change in volatile metabolite production, shifting from production of earthy or musty volatile compounds (e.g., geosmin) to fatty and cheesy volatiles (e.g., 2-nonanone, 2-undecanone). RNA sequencing demonstrated a significant decrease in expression of 356 genes in domesticated strains, with an enrichment of many secondary metabolite production pathways in these downregulated genes. By manipulating the presence of neighboring microbial species and overall resource availability, we demonstrate that the limited competition and high nutrient availability of the cheese environment promote rapid trait evolution of Penicillium molds.IMPORTANCE Industrial cultures of filamentous fungi are used to add unique aesthetics and flavors to cheeses and other microbial foods. How these microbes adapted to live in food environments is generally unknown as most microbial domestication is unintentional. Our work demonstrates that wild molds closely related to the starter culture Penicillium camemberti can readily lose traits and quickly shift toward producing desirable aroma compounds. In addition to experimentally demonstrating a putative domestication pathway for P. camemberti, our work suggests that wild Penicillium isolates could be rapidly domesticated to produce new flavors and aesthetics in fermented foods.

Treehouse: a user-friendly application to obtain subtrees from large phylogenies

Objective

Phylogenetic trees that contain hundreds to thousands of taxa are now routinely generated. Retrieving the relationships among a subset of taxa in these large phylogenies can be a challenging or time-consuming task. Addressing this challenge requires the development of tools that facilitate the easy retrieval of subtrees from any user-specified set of taxa in a given phylogeny.

Results

We developed treehouse, an open source tool that enables the retrieval of any subtree from a given large phylogeny. With a three-step workflow, treehouse successfully allows a user to obtain a subtree from any phylogeny. Treehouse can help researchers to explore the relationships among any set of taxa from across the tree of life. Treehouse is implemented as a shiny application in the R programming language. Treehouse software and usage instructions are publicly available at https://github.com/JLSteenwyk/treehouse.

Variation and selection on codon usage bias across an entire subphylum

Variation in synonymous codon usage is abundant across multiple levels of organization: between codons of an amino acid, between genes in a genome, and between genomes of different species. It is now well understood that variation in synonymous codon usage is influenced by mutational bias coupled with both natural selection for translational efficiency and genetic drift, but how these processes shape patterns of codon usage bias across entire lineages remains unexplored. To address this question, we used a rich genomic data set of 327 species that covers nearly one third of the known biodiversity of the budding yeast subphylum Saccharomycotina. We found that, while genome-wide relative synonymous codon usage (RSCU) for all codons was highly correlated with the GC content of the third codon position (GC3), the usage of codons for the amino acids proline, arginine, and glycine was inconsistent with the neutral expectation where mutational bias coupled with genetic drift drive codon usage. Examination between genes' effective numbers of codons and their GC3 contents in individual genomes revealed that nearly a quarter of genes (381,174/1,683,203; 23%), as well as most genomes (308/327; 94%), significantly deviate from the neutral expectation. Finally, by evaluating the imprint of translational selection on codon usage, measured as the degree to which genes' adaptiveness to the tRNA pool were correlated with selective pressure, we show that translational selection is widespread in budding yeast genomes (264/327; 81%). These results suggest that the contribution of translational selection and drift to patterns of synonymous codon usage across budding yeasts varies across codons, genes, and genomes; whereas drift is the primary driver of global codon usage across the subphylum, the codon bias of large numbers of genes in the majority of genomes is influenced by translational selection.

Extensive loss of cell-cycle and DNA repair genes in an ancient lineage of bipolar budding yeasts

Cell-cycle checkpoints and DNA repair processes protect organisms from potentially lethal mutational damage. Compared to other budding yeasts in the subphylum Saccharomycotina, we noticed that a lineage in the genus Hanseniaspora exhibited very high evolutionary rates, low Guanine-Cytosine (GC) content, small genome sizes, and lower gene numbers. To better understand Hanseniaspora evolution, we analyzed 25 genomes, including 11 newly sequenced, representing 18/21 known species in the genus. Our phylogenomic analyses identify two Hanseniaspora lineages, a faster-evolving lineage (FEL), which began diversifying approximately 87 million years ago (mya), and a slower-evolving lineage (SEL), which began diversifying approximately 54 mya. Remarkably, both lineages lost genes associated with the cell cycle and genome integrity, but these losses were greater in the FEL. E.g., all species lost the cell-cycle regulator WHIskey 5 (WHI5), and the FEL lost components of the spindle checkpoint pathway (e.g., Mitotic Arrest-Deficient 1 [MAD1], Mitotic Arrest-Deficient 2 [MAD2]) and DNA-damage-checkpoint pathway (e.g., Mitosis Entry Checkpoint 3 [MEC3], RADiation sensitive 9 [RAD9]). Similarly, both lineages lost genes involved in DNA repair pathways, including the DNA glycosylase gene 3-MethylAdenine DNA Glycosylase 1 (MAG1), which is part of the base-excision repair pathway, and the DNA photolyase gene PHotoreactivation Repair deficient 1 (PHR1), which is involved in pyrimidine dimer repair. Strikingly, the FEL lost 33 additional genes, including polymerases (i.e., POLymerase 4 [POL4] and POL32) and telomere-associated genes (e.g., Repressor/activator site binding protein-Interacting Factor 1 [RIF1], Replication Factor A 3 [RFA3], Cell Division Cycle 13 [CDC13], Pbp1p Binding Protein [PBP2]). Echoing these losses, molecular evolutionary analyses reveal that, compared to the SEL, the FEL stem lineage underwent a burst of accelerated evolution, which resulted in greater mutational loads, homopolymer instabilities, and higher fractions of mutations associated with the common endogenously damaged base, 8-oxoguanine. We conclude that Hanseniaspora is an ancient lineage that has diversified and thrived, despite lacking many otherwise highly conserved cell-cycle and genome integrity genes and pathways, and may represent a novel, to our knowledge, system for studying cellular life without them.

Nutritional Heterogeneity Among Aspergillus fumigatus Strains Has Consequences for Virulence in a Strain- and Host-Dependent Manner

Acquisition and subsequent metabolism of different carbon and nitrogen sources have been shown to play an important role in virulence attributes of the fungal pathogen Aspergillus fumigatus, such as the secretion of host tissue-damaging proteases and fungal cell wall integrity. We examined the relationship between the metabolic processes of carbon catabolite repression (CCR), nitrogen catabolite repression (NCR) and virulence in a variety of A. fumigatus clinical isolates. A considerable amount of heterogeneity with respect to the degree of CCR and NCR was observed and a positive correlation between NCR and virulence in a neutropenic mouse model of pulmonary aspergillosis (PA) was found. Isolate Afs35 was selected for further analysis and compared to the reference strain A1163, with both strains presenting the same degree of virulence in a neutropenic mouse model of PA. Afs35 metabolome analysis in physiological-relevant carbon sources indicated an accumulation of intracellular sugars that also serve as cell wall polysaccharide precursors. Genome analysis showed an accumulation of missense substitutions in the regulator of protease secretion and in genes encoding enzymes required for cell wall sugar metabolism. Based on these results, the virulence of strains Afs35 and A1163 was assessed in a triamcinolone murine model of PA and found to be significantly different, confirming the known importance of using different mouse models to assess strain-specific pathogenicity. These results highlight the importance of nitrogen metabolism for virulence and provide a detailed example of the heterogeneity that exists between A. fumigatus isolates with consequences for virulence in a strain-specific and host-dependent manner.

Characterizing the Pathogenic, Genomic, and Chemical Traits of Aspergillus fischeri, a Close Relative of the Major Human Fungal Pathogen Aspergillus fumigatus

Aspergillus fischeri is closely related to Aspergillus fumigatus, the major cause of invasive mold infections. Even though A. fischeri is commonly found in diverse environments, including hospitals, it rarely causes invasive disease. Why A. fischeri causes less human disease than A. fumigatus is unclear. A comparison of A. fischeri and A. fumigatus for pathogenic, genomic, and secondary metabolic traits revealed multiple differences in pathogenesis-related phenotypes. We observed that A. fischeri NRRL 181 is less virulent than A. fumigatus strain CEA10 in multiple animal models of disease, grows slower in low-oxygen environments, and is more sensitive to oxidative stress. Strikingly, the observed differences for some traits are of the same order of magnitude as those previously reported between A. fumigatus strains. In contrast, similar to what has previously been reported, the two species exhibit high genomic similarity; ∼90% of the A. fumigatus proteome is conserved in A. fischeri, including 48/49 genes known to be involved in A. fumigatus virulence. However, only 10/33 A. fumigatus biosynthetic gene clusters (BGCs) likely involved in secondary metabolite production are conserved in A. fischeri and only 13/48 A. fischeri BGCs are conserved in A. fumigatus Detailed chemical characterization of A. fischeri cultures grown on multiple substrates identified multiple secondary metabolites, including two new compounds and one never before isolated as a natural product. Additionally, an A. fischeri deletion mutant of laeA, a master regulator of secondary metabolism, produced fewer secondary metabolites and in lower quantities, suggesting that regulation of secondary metabolism is at least partially conserved. These results suggest that the nonpathogenic A. fischeri possesses many of the genes important for A. fumigatus pathogenicity but is divergent with respect to its ability to thrive under host-relevant conditions and its secondary metabolism.IMPORTANCEAspergillus fumigatus is the primary cause of aspergillosis, a devastating ensemble of diseases associated with severe morbidity and mortality worldwide. A. fischeri is a close relative of A. fumigatus but is not generally observed to cause human disease. To gain insights into the underlying causes of this remarkable difference in pathogenicity, we compared two representative strains (one from each species) for a range of pathogenesis-relevant biological and chemical characteristics. We found that disease progression in multiple A. fischeri mouse models was slower and caused less mortality than A. fumigatus Remarkably, the observed differences between A. fischeri and A. fumigatus strains examined here closely resembled those previously described for two commonly studied A. fumigatus strains, AF293 and CEA10. A. fischeri and A. fumigatus exhibited different growth profiles when placed in a range of stress-inducing conditions encountered during infection, such as low levels of oxygen and the presence of chemicals that induce the production of reactive oxygen species. We also found that the vast majority of A. fumigatus genes known to be involved in virulence are conserved in A. fischeri, whereas the two species differ significantly in their secondary metabolic pathways. These similarities and differences that we report here are the first step toward understanding the evolutionary origin of a major fungal pathogen.

Mapping the Fungal Battlefield: Using in situ Chemistry and Deletion Mutants to Monitor Interspecific Chemical Interactions Between Fungi

Fungi grow in competitive environments, and to cope, they have evolved strategies, such as the ability to produce a wide range of secondary metabolites. This begs two related questions. First, how do secondary metabolites influence fungal ecology and interspecific interactions? Second, can these interspecific interactions provide a way to “see” how fungi respond, chemically, within a competitive environment? To evaluate these, and to gain insight into the secondary metabolic arsenal fungi possess, we co-cultured Aspergillus fischeri, a genetically tractable fungus that produces a suite of mycotoxins, with Xylaria cubensis, a fungus that produces the fungistatic compound and FDA-approved drug, griseofulvin. To monitor and characterize fungal chemistry in situ, we used the droplet-liquid microjunction-surface sampling probe (droplet probe). The droplet probe makes a microextraction at defined locations on the surface of the co-culture, followed by analysis of the secondary metabolite profile via liquid chromatography-mass spectrometry. Using this, we mapped and compared the spatial profiles of secondary metabolites from both fungi in monoculture versus co-culture. X. cubensis predominantly biosynthesized griseofulvin and dechlorogriseofulvin in monoculture. In contrast, under co-culture conditions a deadlock was formed between the two fungi, and X. cubensis biosynthesized the same two secondary metabolites, along with dechloro-5′-hydroxygriseofulvin and 5′-hydroxygriseofulvin, all of which have fungistatic properties, as well as mycotoxins like cytochalasin D and cytochalasin C. In contrast, in co-culture, A. fischeri increased the production of the mycotoxins fumitremorgin B and verruculogen, but otherwise remained unchanged relative to its monoculture. To evaluate that secondary metabolites play an important role in defense and territory establishment, we co-cultured A. fischeri lacking the master regulator of secondary metabolism laeA with X. cubensis. We found that the reduced secondary metabolite biosynthesis of the ΔlaeA strain of A. fischeri eliminated the organism’s ability to compete in co-culture and led to its displacement by X. cubensis. These results demonstrate the potential of in situ chemical analysis and deletion mutant approaches for shedding light on the ecological roles of secondary metabolites and how they influence fungal ecological strategies; co-culturing may also stimulate the biosynthesis of secondary metabolites that are not produced in monoculture in the laboratory.

integRATE: a desirability-based data integration framework for the prioritization of candidate genes across heterogeneous omics and its application to preterm birth

BACKGROUND

The integration of high-quality, genome-wide analyses offers a robust approach to elucidating genetic factors involved in complex human diseases. Even though several methods exist to integrate heterogeneous omics data, most biologists still manually select candidate genes by examining the intersection of lists of candidates stemming from analyses of different types of omics data that have been generated by imposing hard (strict) thresholds on quantitative variables, such as P-values and fold changes, increasing the chance of missing potentially important candidates.

METHODS

To better facilitate the unbiased integration of heterogeneous omics data collected from diverse platforms and samples, we propose a desirability function framework for identifying candidate genes with strong evidence across data types as targets for follow-up functional analysis. Our approach is targeted towards disease systems with sparse, heterogeneous omics data, so we tested it on one such pathology: spontaneous preterm birth (sPTB).

RESULTS

We developed the software integRATE, which uses desirability functions to rank genes both within and across studies, identifying well-supported candidate genes according to the cumulative weight of biological evidence rather than based on imposition of hard thresholds of key variables. Integrating 10 sPTB omics studies identified both genes in pathways previously suspected to be involved in sPTB as well as novel genes never before linked to this syndrome. integRATE is available as an R package on GitHub ( https://github.com/haleyeidem/integRATE ).

CONCLUSIONS

Desirability-based data integration is a solution most applicable in biological research areas where omics data is especially heterogeneous and sparse, allowing for the prioritization of candidate genes that can be used to inform more targeted downstream functional analyses.

Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum

Budding yeasts (subphylum Saccharomycotina) are found in every biome and are as genetically diverse as plants or animals. To understand budding yeast evolution, we analyzed the genomes of 332 yeast species, including 220 newly sequenced ones, which represent nearly one-third of all known budding yeast diversity. Here, we establish a robust genus-level phylogeny comprising 12 major clades, infer the timescale of diversification from the Devonian period to the present, quantify horizontal gene transfer (HGT), and reconstruct the evolution of 45 metabolic traits and the metabolic toolkit of the budding yeast common ancestor (BYCA). We infer that BYCA was metabolically complex and chronicle the tempo and mode of genomic and phenotypic evolution across the subphylum, which is characterized by very low HGT levels and widespread losses of traits and the genes that control them. More generally, our results argue that reductive evolution is a major mode of evolutionary diversification.

Copy Number Variation in Fungi and Its Implications for Wine Yeast Genetic Diversity and Adaptation

In recent years, copy number (CN) variation has emerged as a new and significant source of genetic polymorphisms contributing to the phenotypic diversity of populations. CN variants are defined as genetic loci that, due to duplication and deletion, vary in their number of copies across individuals in a population. CN variants range in size from 50 base pairs to whole chromosomes, can influence gene activity, and are associated with a wide range of phenotypes in diverse organisms, including the budding yeast Saccharomyces cerevisiae. In this review, we introduce CN variation, discuss the genetic and molecular mechanisms implicated in its generation, how they can contribute to genetic and phenotypic diversity in fungal populations, and consider how CN variants may influence wine yeast adaptation in fermentation-related processes. In particular, we focus on reviewing recent work investigating the contribution of changes in CN of fermentation-related genes in yeast wine strains and offer notable illustrations of such changes, including the high levels of CN variation among the CUP genes, which confer resistance to copper, a metal with fungicidal properties, and the preferential deletion and duplication of the MAL1 and MAL3 loci, respectively, which are responsible for metabolizing maltose and sucrose. Based on the available data, we propose that CN variation is a substantial dimension of yeast genetic diversity that occurs largely independent of single nucleotide polymorphisms. As such, CN variation harbors considerable potential for understanding and manipulating yeast strains in the wine fermentation environment and beyond.

Copy number variation contributes to cryptic genetic variation in outbreak lineages of Cryptococcus gattii from the North American Pacific Northwest

Background

Copy number variants (CNVs) are a class of structural variants (SVs) and are defined as fragments of DNA that are present at variable copy number in comparison with a reference genome. Recent advances in bioinformatics methodologies and sequencing technologies have enabled the high-resolution quantification of genome-wide CNVs. In pathogenic fungi SVs have been shown to alter gene expression, influence host specificity, and drive fungicide resistance, but little attention has focused specifically on CNVs. Using publicly available sequencing data, we identified 90 isolates across 212 Cryptococcus gattii genomes that belong to the VGII subgroups responsible for the recent deadly outbreaks in the North American Pacific Northwest. We generated CNV profiles for each sample to investigate the prevalence and function of CNV in C. gattii.

Results

We identified eight genetic clusters among publicly available Illumina whole genome sequence data from 212 C. gattii isolates through population structure analysis. Three clusters represent the VGIIa, VGIIb, and VGIIc subgroups from the North American Pacific Northwest. CNV was bioinformatically predicted and affected ~300–400 Kilobases (Kb) of the C. gattii VGII subgroup genomes. Sixty-seven loci, encompassing 58 genes, showed highly divergent patterns of copy number variation between VGII subgroups. Analysis of PFam domains within divergent CN variable genes revealed enrichment of protein domains associated with transport, cell wall organization and external encapsulating structure.

Conclusions

CNVs may contribute to pathological and phenotypic differences observed between the C. gattii VGIIa, VGIIb, and VGIIc subpopulations. Genes overlapping with population differentiated CNVs were enriched for several virulence related functional terms. These results uncover novel candidate genes to examine the genetic and functional underpinnings of C. gattii pathogenicity.

Electronic supplementary material

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