Global deceleration of gene evolution following recent genome hybridizations in fungi

Extensive Copy Number Variation Explains Genome Size Variation in the Unicellular Zygnematophycean Alga, Closterium peracerosum-strigosum-littorale Complex

Genome sizes are known to vary within and among closely related species, but the knowledge about genomic factors contributing to the variation and their impacts on gene functions is limited to only a small number of species. This study identified a more than 2-fold heritable genome size variation among the unicellular Zygnematophycean alga, Closterium peracerosum-strigosum-littorale (C. psl.) complex, based on short-read sequencing analysis of 22 natural strains and F1 segregation analysis. Six de novo assembled genomes revealed that genome size variation is largely attributable to genome-wide copy number variation (CNV) among strains rather than mating type-linked genomic regions or specific repeat sequences such as rDNA. Notably, about 30% of genes showed CNV even between strains that can mate with each other. Transcriptome and gene ontology analysis demonstrated that CNV is distributed nonrandomly in terms of gene functions, such that CNV was more often observed in the gene set with stage-specific expression. Furthermore, in about 30% of these genes with CNV, the expression level does not increase proportionally with the gene copy number, suggesting presence of dosage compensation, which was overrepresented in genes involved in basic biological functions, such as translation. Nonrandom patterns in gene duplications and corresponding expression changes in terms of gene functions may contribute to maintaining the high level of CNV associated with extensive genome size variation in the C. psl. complex, despite its possible detrimental effects.

Draft genome sequences of strains CBS6241 and CBS6242 of the basidiomycetous yeast Filobasidium floriforme

The Tremellomycetes are a species-rich group within the basidiomycete fungi; however, most analyses of this group to date have focused on pathogenic Cryptococcus species within the order Tremellales. Recent genome-assisted studies of other Tremellomycetes have identified interesting features with respect to biotechnological applications as well as the evolution of genes involved in mating and sexual development. Here, we report genome sequences of two strains of Filobasidium floriforme, a species from the order Filobasidiales, which branches basally to the Tremellales, Trichosporonales, and Holtermanniales. The assembled genomes of strains CBS6241 and CBS6242 are 27.4 Mb and 26.4 Mb in size, respectively, with 8314 and 7695 predicted protein-coding genes. Overall sequence identity at nucleic acid level between the strains is 97%. Among the predicted genes are pheromone precursor and pheromone receptor genes as well as two genes encoding homedomain (HD) transcription factors, which are predicted to be part of the mating type (MAT) locus. Sequence analysis indicates that CBS6241 and CBS6242 carry different alleles for both the pheromone/receptor genes as well as the HD transcription factors. Orthology inference identified 1482 orthogroups exclusively found in F. floriforme, some of which were involved in carbohydrate transport and metabolism. Subsequent CAZyme repertoire characterization identified 267 and 247 enzymes for CBS6241 and CBS6242, respectively, the second highest number of CAZymes among the analyzed Tremellomycete species. In addition, F. floriforme contains five CAZymes absent in other species and several plant-cell-wall degrading CAZymes with the highest copy number in Tremellomycota, indicating the biotechnological potential of this species.

Transcriptional profile of the industrial hybrid Saccharomyces pastorianus reveals temperature-dependent allele expression bias and preferential orthologous protein assemblies

Saccharomyces pastorianus is a natural yeast evolved from different hybridization events between the mesophilic S. cerevisiae and the cold-tolerant S. eubayanus. This complex aneuploid hybrid carries multiple copies of the parental alleles alongside specific hybrid genes and encodes for multiple protein isoforms which impart novel phenotypes, such as the strong ability to ferment at low temperature. These characteristics lead to agonistic competition for substrates and a plethora of biochemical activities, resulting in a unique cellular metabolism. Here, we investigated the transcriptional signature of the different orthologous alleles in S. pastorianus during temperature shifts. We identified temperature-dependent media-independent genes and showed that 35% has their regulation dependent on extracellular leucine uptake, suggesting an interplay between leucine metabolism and temperature response. The analysis of the expression of ortholog parental alleles unveiled that the majority of the genes expresses preferentially one parental allele over the other and that S. eubayanus-like alleles are significantly over-represented among the genes involved in the cold acclimatization. The presence of functionally redundant parental alleles may impact on the nature of protein complexes established in the hybrid, where both parental alleles are competing. Our expression data indicate that the majority of the protein complexes investigated in the hybrid are likely to be either exclusively chimeric or unispecific and that the redundancy is discouraged, a scenario that fits well with the gene balance hypothesis. This study offers the first overview of the transcriptional pattern of S. pastorianus and provides a rationalization for its unique industrial traits at the expression level.

Interspecific hybridization as a driver of fungal evolution and adaptation

Cross-species gene transfer is often associated with bacteria, which have evolved several mechanisms that facilitate horizontal DNA exchange. However, the increased availability of whole-genome sequences has revealed that fungal species also exchange DNA, leading to intertwined lineages, blurred species boundaries or even novel species. In contrast to prokaryotes, fungal DNA exchange originates from interspecific hybridization, where two genomes are merged into a single, often highly unstable, polyploid genome that evolves rapidly into stabler derivatives. The resulting hybrids can display novel combinations of genetic and phenotypic variation that enhance fitness and allow colonization of new niches. Interspecific hybridization led to the emergence of important pathogens of humans and plants (for example, various Candida and 'powdery mildew' species, respectively) and industrially important yeasts, such as Saccharomyces hybrids that are important in the production of cold-fermented lagers or cold-cellared Belgian ales. In this Review, we discuss the genetic processes and evolutionary implications of fungal interspecific hybridization and highlight some of the best-studied examples. In addition, we explain how hybrids can be used to study molecular mechanisms underlying evolution, adaptation and speciation, and serve as a route towards development of new variants for industrial applications.

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families

Motivation

Reconstruction of gene copy number evolution is an essential approach for understanding how complex biological systems have been organized. Although various models have been proposed for gene copy number evolution, existing evolutionary models have not appropriately addressed the fact that different gene families can have very different gene gain/loss rates.

Results

In this study, we developed Mirage (MIxtuRe model for Ancestral Genome Estimation), which allows different gene families to have flexible gene gain/loss rates. Mirage can use three models for formulating heterogeneous evolution among gene families: the discretized Γ model, probability distribution-free model and pattern mixture (PM) model. Simulation analysis showed that Mirage can accurately estimate heterogeneous gene gain/loss rates and reconstruct gene-content evolutionary history. Application to empirical datasets demonstrated that the PM model fits genome data from various taxonomic groups better than the other heterogeneous models. Using Mirage, we revealed that metabolic function-related gene families displayed frequent gene gains and losses in all taxa investigated.

Availability and implementation

The source code of Mirage is freely available at https://github.com/fukunagatsu/Mirage.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

The Interspecific Fungal Hybrid Verticillium longisporum Displays Subgenome-Specific Gene Expression

Hybridization is an important evolutionary mechanism that can enable organisms to adapt to environmental challenges. It has previously been shown that the fungal allodiploid species Verticillium longisporum, the causal agent of verticillium stem striping in rapeseed, originated from at least three independent hybridization events between two haploid Verticillium species. To reveal the impact of genome duplication as a consequence of hybridization, we studied the genome and transcriptome dynamics upon two independent V. longisporum hybridization events, represented by the hybrid lineages “A1/D1” and “A1/D3.” We show that V. longisporum genomes are characterized by extensive chromosomal rearrangements, including between parental chromosomal sets. V. longisporum hybrids display signs of evolutionary dynamics that are typically associated with the aftermath of allodiploidization, such as haploidization and more relaxed gene evolution. The expression patterns of the two subgenomes within the two hybrid lineages are more similar than those of the shared A1 parent between the two lineages, showing that the expression patterns of the parental genomes homogenized within a lineage. However, as genes that display differential parental expression in planta do not typically display the same pattern in vitro, we conclude that subgenome-specific responses occur in both lineages. Overall, our study uncovers genomic and transcriptomic plasticity during the evolution of the filamentous fungal hybrid V. longisporum and illustrates its adaptive potential.

The evolving species concepts used for yeasts: from phenotypes and genomes to speciation networks

Here we review how evolving species concepts have been applied to understand yeast diversity. Initially, a phenotypic species concept was utilized taking into consideration morphological aspects of colonies and cells, and growth profiles. Later the biological species concept was added, which applied data from mating experiments. Biophysical measurements of DNA similarity between isolates were an early measure that became more broadly applied with the advent of sequencing technology, leading to a sequence-based species concept using comparisons of parts of the ribosomal DNA. At present phylogenetic species concepts that employ sequence data of rDNA and other genes are universally applied in fungal taxonomy, including yeasts, because various studies revealed a relatively good correlation between the biological species concept and sequence divergence. The application of genome information is becoming increasingly common, and we strongly recommend the use of complete, rather than draft genomes to improve our understanding of species and their genome and genetic dynamics. Complete genomes allow in-depth comparisons on the evolvability of genomes and, consequently, of the species to which they belong. Hybridization seems a relatively common phenomenon and has been observed in all major fungal lineages that contain yeasts. Note that hybrids may greatly differ in their post-hybridization development. Future in-depth studies, initially using some model species or complexes may shift the traditional species concept as isolated clusters of genetically compatible isolates to a cohesive speciation network in which such clusters are interconnected by genetic processes, such as hybridization.

Multidrug-resistant Trichosporon species: underestimated fungal pathogens posing imminent threats in clinical settings

Species of Trichosporon and related genera are widely used in biotechnology and, hence, many species have their genome sequenced. Importantly, yeasts of the genus Trichosporon have been increasingly identified as a cause of life-threatening invasive trichosporonosis (IT) in humans and are associated with an exceptionally high mortality rate. Trichosporon spp. are intrinsically resistant to frontline antifungal agents, which accounts for numerous reports of therapeutic failure when echinocandins are used to treat IT. Moreover, these fungi have low sensitivity to polyenes and azoles and, therefore, are potentially regarded as multidrug-resistant pathogens. However, despite the clinical importance of Trichosporon spp., our understanding of their antifungal resistance mechanisms is quite limited. Furthermore, antifungal susceptibility testing is not standardized, and there is a lack of interpretive epidemiological cut-off values for minimal inhibitory concentrations to distinguish non-wild type Trichosporon isolates. The route of infection remains obscure and detailed clinical and environmental studies are required to determine whether the Trichosporon infections are endogenous or exogenous in nature. Although our knowledge on effective IT treatments is rather limited and future randomized clinical trials are required to identify the best antifungal agent, the current paradigm advocates the use of voriconazole, removal of central venous catheters and recovery from neutropenia.

Evolutionary and genomic comparisons of hybrid uninucleate and nonhybrid Rhizoctonia fungi

The basidiomycetous fungal genus, Rhizoctonia, can cause severe damage to many plants and is composed of multinucleate, binucleate, and uninucleate species differing in pathogenicity. Here we generated chromosome-scale genome assemblies of the three nuclear types of Rhizoctonia isolates. The genomic comparisons revealed that the uninucleate JN strain likely arose by somatic hybridization of two binucleate isolates, and maintained a diploid nucleus. Homeolog gene pairs in the JN genome have experienced both decelerated or accelerated evolution. Homeolog expression dominance occurred between JN subgenomes, in which differentially expressed genes show potentially less evolutionary constraint than the genes without. Analysis of mating-type genes suggested that Rhizoctonia maintains the ancestral tetrapolarity of the Basidiomycota. Long terminal repeat-retrotransposons displayed a reciprocal correlation with the chromosomal GC content in the three chromosome-scale genomes. The more aggressive multinucleate XN strain had more genes encoding enzymes for host cell wall decomposition. These findings demonstrate some evolutionary changes of a recently derived hybrid and in multiple nuclear types of Rhizoctonia.

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Description of four Apiotrichum and two Cutaneotrichosporon species isolated from guano samples from bat-inhabited caves in Japan

Four new yeast species belonging to the genus Apiotrichum and two new yeast species belonging to Cutaneotrichosporon are described for strains isolated from guano samples from bat-inhabited caves in Japan. In 2005, we reported these isolates as Trichosporon species based on sequence analyses of the D1/D2 domain of large subunit (LSU) rRNA genes according to available basidiomycetous yeast classification criteria; however, to date, they have not been officially published as new species with descriptions. Their phylogenetic positions have been reanalysed based on comparison of internal transcribed spacer (ITS) region sequences (including the 5.8S rRNA gene) and the D1/D2 domain of the LSU rRNA gene with those of known species; we confirmed clear separation from previously described species. Physiological and biochemical properties of the isolates also suggest their distinctiveness. Therefore, we describe Apiotrichum akiyoshidainum (holotype JCM 12595^T), Apiotrichum chiropterorum (JCM 12594^T), Apiotrichum coprophilum (JCM 12596^T), Apiotrichum otae (JCM 12593^T), Cutaneotrichosporon cavernicola (JCM 12590^T) and Cutaneotrichosporon middelhovenii (JCM 12592^T) as new species. C. cavernicola showed particularly distinctive morphology including large inflated anomalous cells on the hyphae and germination from the cells, although clear clamp connections on the hyphae were not confirmed. Further study is needed to elucidate the morph of this species.

Genomic insights into the lifestyles, functional capacities and oleagenicity of members of the fungal family Trichosporonaceae

Trichosporonaceae incorporates six genera of physiologically and ecologically diverse fungi including both human pathogenic taxa as well as yeasts of biotechnological interest, especially those oleagenic taxa that accumulate large amounts of single cell oils (SCOs). Here, we have undertaken comparative genomic analysis of thirty-three members of the family with a view to gain insight into the molecular determinants underlying their lifestyles and niche specializations. Phylogenomic analysis revealed potential misidentification of three strains which could impact subsequent analyses. Evaluation of the predicted proteins coding sequences showed that the free-living members of the family harbour greater numbers of carbohydrate active enzymes (CAZYmes), metallo- and serine peptidases compared to their host-associated counterparts. Phylogenies of selected lipid biosynthetic enzymes encoded in the genomes of the studied strains revealed disparate evolutionary histories for some proteins inconsistent with the core genome phylogeny. However, the documented oleagenic members distinctly cluster based on the constitution of the upstream regulatory regions of genes encoding acetyl-CoA carboxylase (ACC), ATP-citrate synthase (ACS) and isocitrate dehydrogenase [NADP] (ICDH), which are among the major proteins in the lipid biosynthetic pathway of these yeasts, suggesting a possible pattern in the regulation of these genes.

Draft Genome Analysis of Trichosporonales Species That Contribute to the Taxonomy of the Genus Trichosporon and Related Taxa

Many nomenclatural changes, including proposals of new taxa, have been carried out in fungi to adapt to the "One fungus = One name" (1F=1N) principle. In yeasts, while some changes have been made in response to 1F=1N, most have resulted from two other factors: i) an improved understanding of biological diversity due to an increase in number of known species, and ii) progress in the methods for analyzing and evaluating biological diversity. The method for constructing a backbone tree, which is a basal tree used to infer phylogeny, has also progressed from single-gene trees to multi-locus trees and further, to genome trees. This paper describes recent advances related to the contribution of genomic data to taxonomy, using the order Trichosporonales as an example.

Convergent evolution of linked mating-type loci in basidiomycete fungi

Sexual development is a key evolutionary innovation of eukaryotes. In many species, mating involves interaction between compatible mating partners that can undergo cell and nuclear fusion and subsequent steps of development including meiosis. Mating compatibility in fungi is governed by the mating type (MAT) loci. In basidiomycetes, the ancestral state is hypothesized to be tetrapolar, with two genetically unlinked MAT loci containing homeodomain transcription factor genes (HD locus) and pheromone and pheromone receptor genes (P/R locus), respectively. Alleles at both loci must differ between mating partners for completion of sexual development. However, there are also basidiomycetes with bipolar mating systems, which can arise through genomic linkage of the HD and P/R loci. In the order Tremellales, bipolarity is found only in the pathogenic Cryptococcus species. Here, we describe the analysis of MAT loci from 24 species of the Trichosporonales, a sister order to the Tremellales. In all of the species analyzed, the MAT loci are fused and a single HD gene is present in each mating type, similar to the organization in the pathogenic Cryptococci. However, the HD and P/R allele combinations in the Trichosporonales are different from those in the pathogenic Cryptococci. This and the existence of tetrapolar species in the Tremellales suggest that fusion of the HD and P/R loci occurred independently in the Trichosporonales and pathogenic Cryptococci, supporting the hypothesis of convergent evolution towards fused MAT regions, similar to previous findings in other fungal groups. Unlike the fused MAT loci in several other basidiomycete lineages though, the gene content and gene order within the fused MAT loci are highly conserved in the Trichosporonales, and there is no apparent suppression of recombination extending from the MAT loci to adjacent chromosomal regions, suggesting different mechanisms for the evolution of physically linked MAT loci in these groups.

Genome and transcriptome evolve separately in recently hybridized Trichosporon fungi

Genome hybridization is an important evolutionary event that gives rise to species with novel capabilities. However, the merging of distinct genomes also brings together incompatible regulatory networks that must be resolved during the course of evolution. Understanding of the early stages of post-hybridization evolution is particularly important because changes in these stages have long-term evolutionary consequences. Here, via comparative transcriptomic analyses of two closely related, recently hybridized Trichosporon fungi, T. coremiiforme and T. ovoides, and three extant relatives, we show that early post-hybridization evolutionary processes occur separately at the gene sequence and gene expression levels but together contribute to the stabilization of hybrid genome and transcriptome. Our findings also highlight lineage-specific consequences of genome hybridization, revealing that the transcriptional regulatory dynamics in these hybrids responded completely differently to gene loss events: one involving both subgenomes and another that is strictly subgenome-specific.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

A Trichosporonales genome tree based on 27 haploid and three evolutionarily conserved 'natural' hybrid genomes

To construct a backbone tree consisting of basidiomycetous yeasts, draft genome sequences from 25 species of Trichosporonales (Tremellomycetes, Basidiomycota) were generated. In addition to the hybrid genomes of Trichosporon coremiiforme and Trichosporon ovoides that we described previously, we identified an interspecies hybrid genome in Cutaneotrichosporon mucoides (formerly Trichosporon mucoides). This hybrid genome had a gene retention rate of ~55%, and its closest haploid relative was Cutaneotrichosporon dermatis. After constructing the C. mucoides subgenomes, we generated a phylogenetic tree using genome data from the 27 haploid species and the subgenome data from the three hybrid genome species. It was a high-quality tree with 100% bootstrap support for all of the branches. The genome-based tree provided superior resolution compared with previous multi-gene analyses. Although our backbone tree does not include all Trichosporonales genera (e.g. Cryptotrichosporon), it will be valuable for future analyses of genome data. Interest in interspecies hybrid fungal genomes has recently increased because they may provide a basis for new technologies. The three Trichosporonales hybrid genomes described in this study are different from well-characterized hybrid genomes (e.g. those of Saccharomyces pastorianus and Saccharomyces bayanus) because these hybridization events probably occurred in the distant evolutionary past. Hence, they will be useful for studying genome stability following hybridization and speciation events. Copyright © 2017 John Wiley & Sons, Ltd.

The Cell Biology of the Trichosporon-Host Interaction

Fungi of the genus Trichosporon are increasingly recognized as causative agents of superficial and invasive fungal disease in humans. Although most species are considered commensals of the human skin and gastrointestinal tract, these basidiomycetes are an increasing cause of fungal disease among immunocompromised hosts, such as hematological patients and solid organ transplant recipients. The initiation of commensal or pathogenic programs by Trichosporon spp. involves the adaptation to the host microenvironment and its immune system. However, the exact virulence factors activated upon the transition to a pathogenic lifestyle, including the intricate biology of the cell wall, and how these interact with and subvert the host immune responses remain largely unknown. Here, we revisit our current understanding of the virulence attributes of Trichosporon spp., particularly T. asahii, and their interaction with the host immune system, and accommodate this knowledge within novel perspectives on fungal diagnostics and therapeutics.