Genome-scale phylogeny and contrasting modes of genome evolution in the fungal phylum Ascomycota

Magic-angle spinning NMR spectral editing of polysaccharides in whole cells using the DREAM scheme

Most bacterial, plant and fungal cells possess at their surface a protective layer called the cell wall, conferring strength, plasticity and rigidity to withstand the osmotic pressure. This molecular barrier is crucial for pathogenic microorganisms, as it protects the cell from the local environment and often constitutes the first structural component encountered in the host-pathogen interaction. In pathogenic molds and yeasts, the cell wall constitutes the main target for the development of clinically-relevant antifungal drugs. In the past decade, solid-state NMR has emerged as a powerful analytical technique to investigate the molecular organization of microbial cell walls in the context of intact cells. 13C NMR chemical shift is an exquisite source of information to identify the polysaccharides present in the cell wall, and two-dimensional 13C-13C correlation experiments provide an efficient tool to rapidly access the polysaccharide composition in whole cells. Here we investigate the use of the adiabatic DREAM (for dipolar recoupling enhancement through amplitude modulation) recoupling scheme to improve solid-state NMR analysis of polysaccharides in intact cells. We demonstrate the advantages of two-dimensional 13C-13C experiments using the DREAM recoupling scheme. We report the spectral editing of polysaccharide signals by varying the radio-frequency carrier position. We provide practical considerations for the implementation of DREAM experiments to characterize polysaccharides in whole cells. We demonstrate the approach on intact fungal cells of Neurospora crassa and Aspergillus fumigatus, a model and a pathogenic filamentous fungus, respectively. The approach could be envisioned to efficiently reduce the spectral crowding of more complex cell surfaces, such as cell wall and peptidoglycan in bacteria.

Exploring Saccharomycotina Yeast Ecology Through an Ecological Ontology Framework

Yeasts in the subphylum Saccharomycotina are found across the globe in disparate ecosystems. A major aim of yeast research is to understand the diversity and evolution of ecological traits, such as carbon metabolic breadth, insect association, and cactophily. This includes studying aspects of ecological traits like genetic architecture or association with other phenotypic traits. Genomic resources in the Saccharomycotina have grown rapidly. Ecological data, however, are still limited for many species, especially those only known from species descriptions where usually only a limited number of strains are studied. Moreover, ecological information is recorded in natural language format limiting high throughput computational analysis. To address these limitations, we developed an ontological framework for the analysis of yeast ecology. A total of 1,088 yeast strains were added to the Ontology of Yeast Environments (OYE) and analyzed in a machine-learning framework to connect genotype to ecology. This framework is flexible and can be extended to additional isolates, species, or environmental sequencing data. Widespread adoption of OYE would greatly aid the study of macroecology in the Saccharomycotina subphylum.

Membrane and organelle rearrangement during ascospore formation in budding yeast

SUMMARYIn ascomycete fungi, sexual spores, termed ascospores, are formed after meiosis. Ascospore formation is an unusual cell division in which daughter cells are created within the cytoplasm of the mother cell by de novo generation of membranes that encapsulate each of the haploid chromosome sets created by meiosis. This review describes the molecular events underlying the creation, expansion, and closure of these membranes in the budding yeast, Saccharomyces cerevisiae. Recent advances in our understanding of the regulation of gene expression and the dynamic behavior of different membrane-bound organelles during this process are detailed. While less is known about ascospore formation in other systems, comparison to the distantly related fission yeast suggests that the molecular events will be broadly similar throughout the ascomycetes.

Starships: a new frontier for fungal biology

Transposable elements (TEs) are semiautonomous genetic entities that proliferate in genomes. We recently discovered the Starships, a previously hidden superfamily of giant TEs found in a diverse subphylum of filamentous fungi, the Pezizomycotina. Starships are unlike other eukaryotic TEs because they have evolved mechanisms for both mobilizing entire genes, including those encoding conditionally beneficial phenotypes, and for horizontally transferring between individuals. We argue that Starships have unrivaled capacity to engage their fungal hosts as genetic parasites and mutualists, revealing unexplored terrain for investigating the ecoevolutionary dynamics of TE-eukaryote interactions. We build on existing models of fungal genome evolution by conceptualizing Starships as a distinct genomic compartment whose dynamics profoundly shape fungal biology.

Reticulate evolution: Detection and utility in the phylogenomics era

Phylogenomics has enriched our understanding that the Tree of Life can have network-like or reticulate structures among some taxa and genes. Two non-vertical modes of evolution - hybridization/introgression and horizontal gene transfer - deviate from a strictly bifurcating tree model, causing non-treelike patterns. However, these reticulate processes can produce similar patterns to incomplete lineage sorting or recombination, potentially leading to ambiguity. Here, we present a brief overview of a phylogenomic workflow for inferring organismal histories and compare methods for distinguishing modes of reticulate evolution. We discuss how the timing of coalescent events can help disentangle introgression from incomplete lineage sorting and how horizontal gene transfer events can help determine the relative timing of speciation events. In doing so, we identify pitfalls of certain methods and discuss how to extend their utility across the Tree of Life. Workflows, methods, and future directions discussed herein underscore the need to embrace reticulate evolutionary patterns for understanding the timing and rates of evolutionary events, providing a clearer view of life's history.

Understanding the molecular mechanisms of human diseases: the benefits of fission yeasts

The role of model organisms such as yeasts in life science research is crucial. Although the baker's yeast (Saccharomyces cerevisiae) is the most popular model among yeasts, the contribution of the fission yeasts (Schizosaccharomyces) to life science is also indisputable. Since both types of yeasts share several thousands of common orthologous genes with humans, they provide a simple research platform to investigate many fundamental molecular mechanisms and functions, thereby contributing to the understanding of the background of human diseases. In this review, we would like to highlight the many advantages of fission yeasts over budding yeasts. The usefulness of fission yeasts in virus research is shown as an example, presenting the most important research results related to the Human Immunodeficiency Virus Type 1 (HIV-1) Vpr protein. Besides, the potential role of fission yeasts in the study of prion biology is also discussed. Furthermore, we are keen to promote the uprising model yeast Schizosaccharomyces japonicus, which is a dimorphic species in the fission yeast genus. We propose the hyphal growth of S. japonicus as an unusual opportunity as a model to study the invadopodia of human cancer cells since the two seemingly different cell types can be compared along fundamental features. Here we also collect the latest laboratory protocols and bioinformatics tools for the fission yeasts to highlight the many possibilities available to the research community. In addition, we present several limiting factors that everyone should be aware of when working with yeast models.

Genomic insights into the evolution of secondary metabolism of Escovopsis and its allies, specialized fungal symbionts of fungus-farming ants

The metabolic intimacy of symbiosis often demands the work of specialists. Natural products and defensive secondary metabolites can drive specificity by ensuring infection and propagation across host generations. But in contrast to bacteria, little is known about the diversity and distribution of natural product biosynthetic pathways among fungi and how they evolve to facilitate symbiosis and adaptation to their host environment. In this study, we define the secondary metabolism of Escovopsis and closely related genera, symbionts in the gardens of fungus-farming ants. We ask how the gain and loss of various biosynthetic pathways correspond to divergent lifestyles. Long-read sequencing allowed us to define the chromosomal features of representative Escovopsis strains, revealing highly reduced genomes composed of seven to eight chromosomes. The genomes are highly syntenic with macrosynteny decreasing with increasing phylogenetic distance, while maintaining a high degree of mesosynteny. An ancestral state reconstruction analysis of biosynthetic pathways revealed that, while many secondary metabolites are shared with non-ant-associated Sordariomycetes, 56 pathways are unique to the symbiotic genera. Reflecting adaptation to diverging ant agricultural systems, we observe that the stepwise acquisition of these pathways mirrors the ecological radiations of attine ants and the dynamic recruitment and replacement of their fungal cultivars. As different clades encode characteristic combinations of biosynthetic gene clusters, these delineating profiles provide important insights into the possible mechanisms underlying specificity between these symbionts and their fungal hosts. Collectively, our findings shed light on the evolutionary dynamic nature of secondary metabolism in Escovopsis and its allies, reflecting adaptation of the symbionts to an ancient agricultural system.IMPORTANCEMicrobial symbionts interact with their hosts and competitors through a remarkable array of secondary metabolites and natural products. Here, we highlight the highly streamlined genomic features of attine-associated fungal symbionts. The genomes of Escovopsis species, as well as species from other symbiont genera, many of which are common with the gardens of fungus-growing ants, are defined by seven chromosomes. Despite a high degree of metabolic conservation, we observe some variation in the symbionts' potential to produce secondary metabolites. As the phylogenetic distribution of the encoding biosynthetic gene clusters coincides with attine transitions in agricultural systems, we highlight the likely role of these metabolites in mediating adaptation by a group of highly specialized symbionts.

Unique trajectory of gene family evolution from genomic analysis of nearly all known species in an ancient yeast lineage

Gene gains and losses are a major driver of genome evolution; their precise characterization can provide insights into the origin and diversification of major lineages. Here, we examined gene family evolution of 1,154 genomes from nearly all known species in the medically and technologically important yeast subphylum Saccharomycotina. We found that yeast gene family and genome evolution are distinct from plants, animals, and filamentous ascomycetes and are characterized by small genome sizes and smaller gene numbers but larger gene family sizes. Faster-evolving lineages (FELs) in yeasts experienced significantly higher rates of gene losses-commensurate with a narrowing of metabolic niche breadth-but higher speciation rates than their slower-evolving sister lineages (SELs). Gene families most often lost are those involved in mRNA splicing, carbohydrate metabolism, and cell division and are likely associated with intron loss, metabolic breadth, and non-canonical cell cycle processes. Our results highlight the significant role of gene family contractions in the evolution of yeast metabolism, genome function, and speciation, and suggest that gene family evolutionary trajectories have differed markedly across major eukaryotic lineages.

Exploiting synthetic biology platforms for enhanced biosynthesis of natural products in Yarrowia lipolytica

With the rapid development of synthetic biology, researchers can design, modify, or even synthesize microorganisms de novo, and microorganisms endowed with unnatural functions can be considered "artificial life" and facilitate the development of functional products. Based on this concept, researchers can solve critical problems related to the insufficient supply of natural products, such as low yields, long production cycles, and cumbersome procedures. Due to its superior performance and unique physiological and biochemical characteristics, Yarrowia lipolytica is a favorable chassis cell used for green biomanufacturing by numerous researchers. This paper mainly reviews the development of synthetic biology techniques for Y. lipolytica and summarizes the recent research progress on the synthesis of natural products in Y. lipolytica. This review will promote the continued innovative development of Y. lipolytica by providing theoretical guidance for research on the biosynthesis of natural products.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Regional centromere configuration in the fungal pathogens of the Pneumocystis genus

Centromeres are constricted chromosomal regions that are essential for cell division. In eukaryotes, centromeres display a remarkable architectural and genetic diversity. The basis of centromere-accelerated evolution remains elusive. Here, we focused on Pneumocystis species, a group of mammalian-specific fungal pathogens that form a sister taxon with that of the Schizosaccharomyces pombe, an important genetic model for centromere biology research. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of S. pombe. Using organisms from a short-term in vitro culture or infected animal models and chromatin immunoprecipitation (ChIP)-Seq, we identified CENP-A bound regions in two Pneumocystis species that diverged ~35 million years ago. Each species has a unique short regional centromere (<10 kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. These features suggest an epigenetic specification of centromere function. Analysis of centromeric DNA across multiple Pneumocystis species suggests a vertical transmission at least 100 million years ago. The common ancestry of Pneumocystis and S. pombe centromeres is untraceable at the DNA level, but the overall architectural similarity could be the result of functional constraint for successful chromosomal segregation.IMPORTANCEPneumocystis species offer a suitable genetic system to study centromere evolution in pathogens because of their phylogenetic proximity with the non-pathogenic yeast S. pombe, a popular model for cell biology. We used this system to explore how centromeres have evolved after the divergence of the two clades ~ 460 million years ago. To address this question, we established a protocol combining short-term culture and ChIP-Seq to characterize centromeres in multiple Pneumocystis species. We show that Pneumocystis have short epigenetic centromeres that function differently from those in S. pombe.

Schizosaccharomyces versatilis represents a distinct evolutionary lineage of fission yeast

The fission yeast species Schizosaccharomyces japonicus is currently divided into two varieties-S. japonicus var. japonicus and S. japonicus var. versatilis. Here we examine the var. versatilis isolate CBS5679. The CBS5679 genome shows 88% identity to the reference genome of S. japonicus var. japonicus at the coding sequence level, with phylogenetic analyses suggesting that it has split from the S. japonicus lineage 25 million years ago. The CBS5679 genome contains a reciprocal translocation between chromosomes 1 and 2, together with several large inversions. The products of genes linked to the major translocation are associated with 'metabolism' and 'cellular assembly' ontology terms. We further show that CBS5679 does not generate viable progeny with the reference strain of S. japonicus. Although CBS5679 shares closer similarity to the 'type' strain of var. versatilis as compared to S. japonicus, it is not identical to the type strain, suggesting population structure within var. versatilis. We recommend that the taxonomic status of S. japonicus var. versatilis is raised, with it being treated as a separate species, Schizosaccharomyces versatilis.

Reinstatement of the fission yeast species Schizosaccharomyces versatilis Wickerham et Duprat, a sibling species of Schizosaccharomyces japonicus

Schizosaccharomyces japonicus Yukawa et Maki (1931) and Schizosaccharomyces versatilis Wickerham et Duprat (1945) have been treated as varieties of S. japonicus or as conspecific, based on various approaches including mating trials and nDNA/nDNA optical reassociation studies. However, the type strains of S. japonicus and S. versatilis differ by five substitutions (99.15% identity) and one 1-bp indel in the sequences of the D1/D2 domain of the 26S rRNA gene, and 23 substitutions (96.3% identity) and 31-bp indels in the sequences of internal transcribed spacer (ITS) of rRNA, suggesting that they may not be conspecific. To reassess their taxonomic status, we conducted mating trials and whole-genome analyses. Mating trials using the type strains showed a strong but incomplete prezygotic sterility barrier, yielding interspecies mating products at two orders of magnitude lower efficiency than intraspecies matings. These mating products, which were exclusively allodiploid hybrids, were unable to undergo the haplontic life cycle of the parents. We generated chromosome-level gap-less genome assemblies for both type strains. Whole genome sequences yielded an average nucleotide identity (ANI) of 86.4%, indicating clear separation of S. japonicus and S. versatilis. Based on these findings, we propose the reinstatement of S. versatilis as a distinct species (holotype strain: CBS 103T and ex-types: NRRL Y-1026, NBRC 1607, ATCC 9987, PYCC 7100; Mycobank no.: 847838).

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Natural trait variation across Saccharomycotina species

Among molecular biologists, the group of fungi called Saccharomycotina is famous for its yeasts. These yeasts in turn are famous for what they have in common-genetic, biochemical, and cell-biological characteristics that serve as models for plants and animals. But behind the apparent homogeneity of Saccharomycotina species lie a wealth of differences. In this review, we discuss traits that vary across the Saccharomycotina subphylum. We describe cases of bright pigmentation; a zoo of cell shapes; metabolic specialties; and species with unique rules of gene regulation. We discuss the genetics of this diversity and why it matters, including insights into basic evolutionary principles with relevance across Eukarya.

Taxogenomic placement of Rasporella oleae and Rasporella dianae gen. and spp. nov., two insect associated yeast species

During the course of independent studies in Europe, North America, and Africa, seven yeast strains were isolated from insect frass, decaying wood, tree flux, and olive oil sediment. Phylogenetic analysis of two barcoding DNA regions (internal transcribed spacer and the D1/D2 domain of the LSU rRNA gene) revealed that they belong to two closely related undescribed species distinct from all genera in the family Debaryomycetaceae. For reliable taxonomic placement the genomes of four strains of the two novel species and six type strains of closely related species were sequenced. Orthologous genes from 54 genomes of representatives of the Pichiomycetes and 23 outgroup taxa were concatenated to construct a fully supported phylogenetic tree. Consistent with the assumptions, we found that the two new species belong to a novel genus. In addition, the delimitation of the novel species was supported by genetic distance calculations from average nucleotide identity (ANI) and digital DNA:DNA hybridization (dDDH) values. The physiological characterization of the novel species was generally consistent with their genomic content. All strains had two alleles encoding secretory lipase in either two or three copies depending on the species. However, lipolytic activity was detected only in strains with three copies of the secretory lipase gene. Nevertheless, lipolytic activity might be related to their association with the insect gut. Based on these results, formal descriptions of the new genus Rasporella gen. nov. and of two new species Rasporella dianae sp. nov. (holotype UCDFST 68-643T , MycoBank no.: 850238) and Rasporella oleae sp. nov. (holotype ZIM 2471T , MycoBank no.: 850126) are provided.

Identification of novel proteins regulating lipid droplet biogenesis in filamentous fungi

Lipid droplets (LDs) are storage organelles for neutral lipids which are critical for lipid homeostasis. Current knowledge of fungal LD biogenesis is largely limited to budding yeast, while LD regulation in multinucleated filamentous fungi which exhibit considerable metabolic activity remains unexplored. In this study, 19 LD-associated proteins were identified in the multinucleated species Aspergillus oryzae using a colocalization screening of a previously established enhanced green fluorescent protein (EGFP) fusion library. Functional screening identified 12 lipid droplet-regulating (LDR) proteins whose loss of function resulted in irregular LD biogenesis, particularly in terms of LD number and size. Bioinformatics analysis, targeted mutagenesis, and microscopy revealed four LDR proteins that localize to LD via the putative amphipathic helices (AHs). Further analysis revealed that LdrA with an Opi1 domain is essential for cytoplasmic and nuclear LD biogenesis involving a novel AH. Phylogenetic analysis demonstrated that the patterns of gene evolution were predominantly based on gene duplication. Our study identified a set of novel proteins involved in the regulation of LD biogenesis, providing unique molecular and evolutionary insights into fungal lipid storage.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Genomic and ecological factors shaping specialism and generalism across an entire subphylum

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Paradigms proposed to explain this variation either invoke trade-offs between performance efficiency and breadth or underlying intrinsic or extrinsic factors. We assembled genomic (1,154 yeast strains from 1,049 species), metabolic (quantitative measures of growth of 843 species in 24 conditions), and ecological (environmental ontology of 1,088 species) data from nearly all known species of the ancient fungal subphylum Saccharomycotina to examine niche breadth evolution. We found large interspecific differences in carbon breadth stem from intrinsic differences in genes encoding specific metabolic pathways but no evidence of trade-offs and a limited role of extrinsic ecological factors. These comprehensive data argue that intrinsic factors driving microbial niche breadth variation.

Streamlining N-terminally anchored yeast surface display via structural insights into S. cerevisiae Pir proteins

Surface display co-opts yeast's innate ability to embellish its cell wall with mannoproteins, thus converting the yeast's outer surface into a growing and self-sustaining catalyst. However, the efficient toolbox for converting the enzyme of interest into its surface-displayed isoform is currently lacking, especially if the isoform needs to be anchored to the cell wall near the isoform's N-terminus, e.g., through a short GPI-independent protein anchor. Aiming to advance such N-terminally anchored surface display, we employed in silico and machine-learning strategies to study the 3D structure, function, genomic organisation, and evolution of the Pir protein family, whose members evolved to covalently attach themselves near their N-terminus to the β-1,3-glucan of the cell wall. Through the newly-gained insights, we rationally engineered 14 S. cerevisiae Hsp150 (Pir2)-based fusion proteins. We quantified their performance, uncovering guidelines for efficient yeast surface display while developing a construct that promoted a 2.5-fold more efficient display of a reporter protein than the full-length Hsp150. Moreover, we developed a Pir-tag, i.e., a peptide spanning only 4.5 kDa but promoting as efficient surface display of a reporter protein as the full-length Hsp150. These constructs fortify the existing surface display toolbox, allowing for a prompt and routine refitting of intracellular proteins into their N-terminally anchored isoforms.

Genomic Characterization and Functional Description of Beauveria bassiana Isolates from Latin America

Beauveria bassiana is an entomopathogenic fungus used in agriculture as a biological controller worldwide. Despite being a well-studied organism, there are no genomic studies of B. bassiana isolates from Central American and Caribbean countries. This work characterized the functional potential of eight Neotropical isolates and provided an overview of their genomic characteristics, targeting genes associated with pathogenicity, the production of secondary metabolites, and the identification of CAZYmes as tools for future biotechnological applications. In addition, a comparison between these isolates and reference genomes was performed. Differences were observed according to geographical location and the lineages of the B. bassiana complex to which each isolate belonged.

Schizosaccharomyces lindneri sp. nov., a fission yeast occurring in honey

Two strains of fission yeast were isolated from honey. They differ from the type strain of Schizosaccharomyces octosporus by three substitutions in the D1/D2 domain of the nuclear 26S large subunit ribosomal RNA (rRNA) gene sequence, resulting in a 99.5% identity. In the internal transcribed spacer (ITS) region (consisting of ITS1, 5.8S rDNA, and ITS2), the strains differ from S. octosporus by 16 gaps and 91 substitutions, which is equivalent to an identity of 88.1%. Genome sequencing on one of the new strains revealed that the average nucleotide identity (ANI) between its genome and the reference genome of S. octosporus is 90.43% and there exist major genome rearrangements between the two genomes. Mating analysis revealed that S. octosporus and one of the new strains are completely reproductively separated. A strong prezygotic barrier exists and the few mating products consist of diploid hybrids that do not form recombinant ascospores. In the new strains, asci are either zygotic, arising from conjugation, or they develop without conjugation from asexual cells (azygotic). Compared to the currently recognized Schizosaccharomyces species, the spectrum of nutrients that are assimilated by the new strains is restricted. Of the 43 carbohydrates that were included in the physiological standard tests, only 7 were assimilated. According to the results of the genome sequence analysis, the mating trials, and the phenotypic characterization, the new species Schizosaccharomyces lindneri is described to accommodate the two strains (holotype: CBS 18203^T  and ex-type: MUCL 58363; MycoBank no.: MB 847838).

The Host Adapted Fungal Pathogens of Pneumocystis Genus Utilize Genic Regional Centromeres

Centromeres are genomic regions that coordinate accurate chromosomal segregation during mitosis and meiosis. Yet, despite their essential function, centromeres evolve rapidly across eukaryotes. Centromeres are often the sites of chromosomal breaks which contribute to genome shuffling and promote speciation by inhibiting gene flow. How centromeres form in strongly host-adapted fungal pathogens has yet to be investigated. Here, we characterized the centromere structures in closely related species of mammalian-specific pathogens of the fungal phylum of Ascomycota. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of Schizosaccharomyces pombe. Using organisms from a short-term in vitro culture or infected animal models and ChIP-seq, we identified centromeres in three Pneumocystis species that diverged ~100 million years ago. Each species has a unique short regional centromere (< 10kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. CENP-C, a scaffold protein that links the inner centromere to the kinetochore appears dispensable in one species, suggesting a kinetochore rewiring. Despite the loss of DNA methyltransferases, 5-methylcytosine DNA methylation occurs in these species, though not related to centromere function. These features suggest an epigenetic specification of centromere function.

Functional Clustering of Metabolically Related Genes Is Conserved across Dikarya

Transcriptional regulation is vital for organismal survival, with many layers and mechanisms collaborating to balance gene expression. One layer of this regulation is genome organization, specifically the clustering of functionally related, co-expressed genes along the chromosomes. Spatial organization allows for position effects to stabilize RNA expression and balance transcription, which can be advantageous for a number of reasons, including reductions in stochastic influences between the gene products. The organization of co-regulated gene families into functional clusters occurs extensively in Ascomycota fungi. However, this is less characterized within the related Basidiomycota fungi despite the many uses and applications for the species within this clade. This review will provide insight into the prevalence, purpose, and significance of the clustering of functionally related genes across Dikarya, including foundational studies from Ascomycetes and the current state of our understanding throughout representative Basidiomycete species.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Coding Sequence Insertions in Fungal Genomes are Intrinsically Disordered and can Impart Functionally-Important Properties on the Host Protein

Insertion and deletion mutations (indels) are important mechanisms of generating protein diversity. Indels in coding sequences are under considerable selective pressure to maintain reading frames and to preserve protein function, but once generated, indels provide raw material for the acquisition of new protein properties and functions. We reported recently that coding sequence insertions in the Candida albicans NDU1 protein, a mitochondrial protein involved in the assembly of the NADH:ubiquinone oxidoreductase are imperative for respiration, biofilm formation and pathogenesis. NDU1 inserts are specific to CTG-clade fungi, absent in human ortholog and successfully harnessed as drug targets. Here, we present the first comprehensive report investigating indels and clade-defining insertions (CDIs) in fungal proteomes. We investigated 80 ascomycete proteomes encompassing CTG clade species, the Saccharomycetaceae family, the Aspergillaceae family and the Herpotrichiellaceae (black yeasts) family. We identified over 30,000 insertions, 4,000 CDIs and 2,500 clade-defining deletions (CDDs). Insert sizes range from 1 to over 1,000 residues in length, while maximum deletion length is 19 residues. Inserts are strikingly over-represented in protein kinases, and excluded from structural domains and transmembrane segments. Inserts are predicted to be highly disordered. The amino acid compositions of the inserts are highly depleted in hydrophobic residues and enriched in polar residues. An indel in the Saccharomyces cerevisiae Sth1 protein, the catalytic subunit of the RSC (Remodel the Structure of Chromatin) complex is predicted to be disordered until it forms a ß-strand upon interaction. This interaction performs a vital role in RSC-mediated transcriptional regulation, thereby expanding protein function.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Fungal taxonomy: current status and research agendas for the interdisciplinary and globalisation era

Fungal taxonomy is a fundamental discipline that aims to recognise all fungi and their kinships. Approximately 5% of a practical estimate of 2.2-3.8 million species globally are currently known, and consequently the Fungal Tree of Life (FTOL) is very incompletely reconstructed. With the advances of new technologies, mycology is marching into the interdisciplinary and globalisation era. To make fungal taxonomy relevant, innovative sampling methods and phylogenomics analyses should be performed to reconstruct a much more comprehensive FTOL. In association with this densely sampled FTOL, multiomics will reveal what drives fungal species diversification and how fungal traits evolve to adapt to various environments, while metagenomics will facilitate the understanding and protection of the ecological functions of fungi. A coordinated approach to pursuing these research agendas that includes conceiving of and costing a mission to describe all the fungi on the planet will unlock potential of fungi to support sustainable development of our society.

Metagenomics Shines Light on the Evolution of "Sunscreen" Pigment Metabolism in the Teloschistales (Lichen-Forming Ascomycota)

Fungi produce a vast number of secondary metabolites that shape their interactions with other organisms and the environment. Characterizing the genes underpinning metabolite synthesis is therefore key to understanding fungal evolution and adaptation. Lichenized fungi represent almost one-third of Ascomycota diversity and boast impressive secondary metabolites repertoires. However, most lichen biosynthetic genes have not been linked to their metabolite products. Here we used metagenomic sequencing to survey gene families associated with production of anthraquinones, UV-protectant secondary metabolites present in various fungi, but especially abundant in a diverse order of lichens, the Teloschistales (class Lecanoromycetes, phylum Ascomycota). We successfully assembled 24 new, high-quality lichenized-fungal genomes de novo and combined them with publicly available Lecanoromycetes genomes from taxa with diverse secondary chemistry to produce a whole-genome tree. Secondary metabolite biosynthetic gene cluster (BGC) analysis showed that whilst lichen BGCs are numerous and highly dissimilar, core enzyme genes are generally conserved across taxa. This suggests metabolite diversification occurs via re-shuffling existing enzyme genes with novel accessory genes rather than BGC gains/losses or de novo gene evolution. We identified putative anthraquinone BGCs in our lichen dataset that appear homologous to anthraquinone clusters from non-lichenized fungi, suggesting these genes were present in the common ancestor of the subphylum Pezizomycotina. Finally, we identified unique transporter genes in Teloschistales anthraquinone BGCs that may explain why these metabolites are so abundant and ubiquitous in these lichens. Our results support the importance of metagenomics for understanding the secondary metabolism of non-model fungi such as lichens.

Genomic diversification of the specialized parasite of the fungus-growing ant symbiosis

Fungi shape the diversity of life. Characterizing the evolution of fungi is critical to understanding symbiotic associations across kingdoms. In this study, we investigate the genomic and metabolomic diversity of the genus Escovopsis, a specialized parasite of fungus-growing ant gardens. Based on 25 high-quality draft genomes, we show that Escovopsis forms a monophyletic group arising from a mycoparasitic fungal ancestor 61.82 million years ago (Mya). Across the evolutionary history of fungus-growing ants, the dates of origin of most clades of Escovopsis correspond to the dates of origin of the fungus-growing ants whose gardens they parasitize. We reveal that genome reduction, determined by both genomic sequencing and flow cytometry, is a consistent feature across the genus Escovopsis, largely occurring in coding regions, specifically in the form of gene loss and reductions in copy numbers of genes. All functional gene categories have reduced copy numbers, but resistance and virulence genes maintain functional diversity. Biosynthetic gene clusters (BGCs) contribute to phylogenetic differences among Escovopsis spp., and sister taxa in the Hypocreaceae. The phylogenetic patterns of co-diversification among BGCs are similarly exhibited across mass spectrometry analyses of the metabolomes of Escovopsis and their sister taxa. Taken together, our results indicate that Escovopsis spp. evolved unique genomic repertoires to specialize on the fungus-growing ant-microbe symbiosis.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

Mycology: Metagenomes illuminate evolutionary relationships and reframe symbiotic interactions

An intriguing new study leverages newly generated metagenomes to remap the evolution of the most species-rich clade of fungi, highlighting how some of the most intriguing and visible manifestations of symbioses - lichens - may arise.

Genome-level analyses resolve an ancient lineage of symbiotic ascomycetes

Ascomycota account for about two-thirds of named fungal species.1 Over 98% of known Ascomycota belong to the Pezizomycotina, including many economically important species as well as diverse pathogens, decomposers, and mutualistic symbionts.2 Our understanding of Pezizomycotina evolution has until now been based on sampling traditionally well-defined taxonomic classes.3,4,5 However, considerable diversity exists in undersampled and uncultured, putatively early-diverging lineages, and the effect of these on evolutionary models has seldom been tested. We obtained genomes from 30 putative early-diverging lineages not included in recent phylogenomic analyses and analyzed these together with 451 genomes covering all available ascomycete genera. We show that 22 of these lineages, collectively representing over 600 species, trace back to a single origin that diverged from the common ancestor of Eurotiomycetes and Lecanoromycetes over 300 million years BP. The new clade, which we recognize as a more broadly defined Lichinomycetes, includes lichen and insect symbionts, endophytes, and putative mycorrhizae and encompasses a range of morphologies so disparate that they have recently been placed in six different taxonomic classes. To test for shared hidden features within this group, we analyzed genome content and compared gene repertoires to related groups in Ascomycota. Regardless of their lifestyle, Lichinomycetes have smaller genomes than most filamentous Ascomycota, with reduced arsenals of carbohydrate-degrading enzymes and secondary metabolite gene clusters. Our expanded genome sample resolves the relationships of numerous "orphan" ascomycetes and establishes the independent evolutionary origins of multiple mutualistic lifestyles within a single, morphologically hyperdiverse clade of fungi.

Long-Read Sequencing Analysis Revealed the Impact of Forest Conversion on Soil Fungal Diversity in Limu Mountain, Hainan

Soil fungi are essential to soil microorganisms that play an important role in the ecosystem's soil carbon cycle and mineral nutrient transformation. Understanding the structural characteristics and diversity of soil fungal communities helps understand the health of forest ecosystems. The transition from tropical rainforest to artificial forest greatly impacts the composition and diversity of fungal communities. Hainan Limushan tropical rainforest National Park has a large area of artificial forests. Ecologists have conducted in-depth studies on the succession of animals and plants to regenerate tropical rainforests. There are few reports on the diversity of soil fungi and its influencing factors in the succession of tropical rainforests in Limu Mountain. In this study, 44 soil samples from five different stands were collected in the tropical rainforest of Limushan, Hainan. High-throughput sequencing of rDNA in its region was used to analyze fungal communities and study their α and β diversity. Analysis of variance and multiple regression models was used to analyze soil variables and fungal functional groups to determine the effects of interaction between fungi and environmental factors. A total of 273,996 reads and 1290 operational taxonomic units (OTUs) were obtained, belonging to 418 species, 325 genera, 159 families, eight phyla, 30 classes, and 73 orders. The results showed that the composition of soil fungal communities in the five stands was similar, with ascomycetes accounting for 70.5% and basidiomycetes accounting for 14.7%. α and β diversity analysis showed that soil fungi in Limushan tropical rainforest had high abundance and diversity. Multiple regression analysis between soil variables and functional groups showed that organic matter, TN, TP, TK, and AK were excellent predictors for soil fungi. TP was the strongest predictor in all functional groups except soil saprotroph. Organic matter and total nitrogen were the strongest predictors of soil rot. The transformation from tropical rainforest to artificial forest in Limushan did not change the soil fungal community structure, but the richness and diversity of soil fungi changed. The forest transformation did not lead to decreased soil fungal abundance and diversity. Different vegetation types and soil properties affect the diversity of soil fungal communities. We found that Caribbean pine plantations can improve soil fungal diversity, while long-term Eucalyptus spp. plantations may reduce soil fungal diversity.

Evolution and codon usage bias of mitochondrial and nuclear genomes in Aspergillus section Flavi

The fungal genus Aspergillus contains a diversity of species divided into taxonomic sections of closely related species. Section Flavi contains 33 species, many of industrial, agricultural, or medical relevance. Here, we analyze the mitochondrial genomes (mitogenomes) of 20 Flavi species-including 18 newly assembled mitogenomes-and compare their evolutionary history and codon usage bias patterns to their nuclear counterparts. Codon usage bias refers to variable frequencies of synonymous codons in coding DNA and is shaped by a balance of neutral processes and natural selection. All mitogenomes were circular DNA molecules with highly conserved gene content and order. As expected, genomic content, including GC content, and genome size differed greatly between mitochondrial and nuclear genomes. Phylogenetic analysis based on 14 concatenated mitochondrial genes predicted evolutionary relationships largely consistent with those predicted by a phylogeny constructed from 2,422 nuclear genes. Comparing similarities in interspecies patterns of codon usage bias between mitochondrial and nuclear genomes showed that species grouped differently by patterns of codon usage bias depending on whether analyses were performed using mitochondrial or nuclear relative synonymous usage values. We found that patterns of codon usage bias at gene level are more similar between mitogenomes of different species than the mitogenome and nuclear genome of the same species. Finally, we inferred that, although most genes-both nuclear and mitochondrial-deviated from the neutral expectation for codon usage, mitogenomes were not under translational selection while nuclear genomes were under moderate translational selection. These results contribute to the study of mitochondrial genome evolution in filamentous fungi.

OrthoSNAP: A tree splitting and pruning algorithm for retrieving single-copy orthologs from gene family trees

Molecular evolution studies, such as phylogenomic studies and genome-wide surveys of selection, often rely on gene families of single-copy orthologs (SC-OGs). Large gene families with multiple homologs in 1 or more species-a phenomenon observed among several important families of genes such as transporters and transcription factors-are often ignored because identifying and retrieving SC-OGs nested within them is challenging. To address this issue and increase the number of markers used in molecular evolution studies, we developed OrthoSNAP, a software that uses a phylogenetic framework to simultaneously split gene families into SC-OGs and prune species-specific inparalogs. We term SC-OGs identified by OrthoSNAP as SNAP-OGs because they are identified using a splitting and pruning procedure analogous to snapping branches on a tree. From 415,129 orthologous groups of genes inferred across 7 eukaryotic phylogenomic datasets, we identified 9,821 SC-OGs; using OrthoSNAP on the remaining 405,308 orthologous groups of genes, we identified an additional 10,704 SNAP-OGs. Comparison of SNAP-OGs and SC-OGs revealed that their phylogenetic information content was similar, even in complex datasets that contain a whole-genome duplication, complex patterns of duplication and loss, transcriptome data where each gene typically has multiple transcripts, and contentious branches in the tree of life. OrthoSNAP is useful for increasing the number of markers used in molecular evolution data matrices, a critical step for robustly inferring and exploring the tree of life.

Analysis of plant cell death-inducing proteins of the necrotrophic fungal pathogens Botrytis squamosa and Botrytis elliptica

Fungal plant pathogens secrete proteins that manipulate the host in order to facilitate colonization. Necrotrophs have evolved specialized proteins that actively induce plant cell death by co-opting the programmed cell death machinery of the host. Besides the broad host range pathogen Botrytis cinerea, most other species within the genus Botrytis are restricted to a single host species or a group of closely related hosts. Here, we focused on Botrytis squamosa and B. elliptica, host specific pathogens of onion (Allium cepa) and lily (Lilium spp.), respectively. Despite their occurrence on different hosts, the two fungal species are each other's closest relatives. Therefore, we hypothesize that they share a considerable number of proteins to induce cell death on their respective hosts. In this study, we first confirmed the host-specificity of B. squamosa and B. elliptica. Then we sequenced and assembled high quality genomes. The alignment of these two genomes revealed a high level of synteny with few balanced structural chromosomal arrangements. To assess the cell death-inducing capacity of their secreted proteins, we produced culture filtrates of B. squamosa and B. elliptica that induced cell death responses upon infiltration in host leaves. Protein composition of the culture filtrate was analysed by mass spectrometry, and we identified orthologous proteins that were present in both samples. Subsequently, the expression of the corresponding genes during host infection was compared. RNAseq analysis showed that the majority of the orthogroups of the two sister species display similar expression patterns during infection of their respective host. The analysis of cell death-inducing proteins of B. squamosa and B. elliptica provides insights in the mechanisms used by these two Botrytis species to infect their respective hosts.

Single-chromosome fission yeast models reveal the configuration robustness of a functional genome

In eukaryotic organisms, genetic information is usually carried on multiple chromosomes. Whether and how the number and configuration of chromosomes affect organismal fitness and speciation remain unclear. Here, we have successfully established several single-chromosome fission yeast Schizosaccharomyces pombe strains, in which the three natural chromosomes have been fused into one giant chromosome in different orders. Chromosome fusions accompanied by the deletions of telomeres and centromeres result in the loss of chromosomal interactions and a drastic change of global chromosome organization, but alter gene expression marginally. The single-chromosome strains display little defects in cell morphology, mitosis, genotoxin sensitivity, and meiosis. Crosses between a wild-type strain and a single-chromosome strain or between two single-chromosome strains with different fusion orders suffer defective meiosis and poor spore viability. We conclude that eukaryotic genomes are robust against dramatic chromosomal reconfiguration, and stochastic changes in chromosome number and genome organization during evolution underlie reproductive isolation and speciation.

Beyond Nuclear Ribosomal DNA Sequences: Evolution, Taxonomy, and Closest Known Saprobic Relatives of Powdery Mildew Fungi (Erysiphaceae) Inferred From Their First Comprehensive Genome-Scale Phylogenetic Analyses

Powdery mildew fungi (Erysiphaceae), common obligate biotrophic pathogens of many plants, including important agricultural and horticultural crops, represent a monophyletic lineage within the Ascomycota. Within the Erysiphaceae, molecular phylogenetic relationships and DNA-based species and genera delimitations were up to now mostly based on nuclear ribosomal DNA (nrDNA) phylogenies. This is the first comprehensive genome-scale phylogenetic analysis of this group using 751 single-copy orthologous sequences extracted from 24 selected powdery mildew genomes and 14 additional genomes from Helotiales, the fungal order that includes the Erysiphaceae. Representative genomes of all powdery mildew species with publicly available whole-genome sequencing (WGS) data that were of sufficient quality were included in the analyses. The 24 powdery mildew genomes included in the analysis represented 17 species belonging to eight out of 19 genera recognized within the Erysiphaceae. The epiphytic genera, all but one represented by multiple genomes, belonged each to distinct, well-supported lineages. Three hemiendophytic genera, each represented by a single genome, together formed the hemiendophytic lineage. Out of the 14 other taxa from the Helotiales, Arachnopeziza araneosa, a saprobic species, was the only taxon that grouped together with the 24 genome-sequenced powdery mildew fungi in a monophyletic clade. The close phylogenetic relationship between the Erysiphaceae and Arachnopeziza was revealed earlier by a phylogenomic study of the Leotiomycetes. Further analyses of powdery mildew and Arachnopeziza genomes may discover signatures of the evolutionary processes that have led to obligate biotrophy from a saprobic way of life. A separate phylogeny was produced using the 18S, 5.8S, and 28S nrDNA sequences of the same set of powdery mildew specimens and compared to the genome-scale phylogeny. The nrDNA phylogeny was largely congruent to the phylogeny produced using 751 orthologs. This part of the study has revealed multiple contamination and other quality issues in some powdery mildew genomes. We recommend that the presence of 28S, internal transcribed spacer (ITS), and 18S nrDNA sequences in powdery mildew WGS datasets that are identical to those determined by Sanger sequencing should be used to assess the quality of assemblies, in addition to the commonly used Benchmarking Universal Single-Copy Orthologs (BUSCO) values.

BioKIT: a versatile toolkit for processing and analyzing diverse types of sequence data

Bioinformatic analysis-such as genome assembly quality assessment, alignment summary statistics, relative synonymous codon usage, file format conversion, and processing and analysis-is integrated into diverse disciplines in the biological sciences. Several command-line pieces of software have been developed to conduct some of these individual analyses, but unified toolkits that conduct all these analyses are lacking. To address this gap, we introduce BioKIT, a versatile command line toolkit that has, upon publication, 42 functions, several of which were community-sourced, that conduct routine and novel processing and analysis of genome assemblies, multiple sequence alignments, coding sequences, sequencing data, and more. To demonstrate the utility of BioKIT, we conducted a comprehensive examination of relative synonymous codon usage across 171 fungal genomes that use alternative genetic codes, showed that the novel metric of gene-wise relative synonymous codon usage can accurately estimate gene-wise codon optimization, evaluated the quality and characteristics of 901 eukaryotic genome assemblies, and calculated alignment summary statistics for 10 phylogenomic data matrices. BioKIT will be helpful in facilitating and streamlining sequence analysis workflows. BioKIT is freely available under the MIT license from GitHub (https://github.com/JLSteenwyk/BioKIT), PyPi (https://pypi.org/project/jlsteenwyk-biokit/), and the Anaconda Cloud (https://anaconda.org/jlsteenwyk/jlsteenwyk-biokit). Documentation, user tutorials, and instructions for requesting new features are available online (https://jlsteenwyk.com/BioKIT).

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

Evolution of the human pathogenic lifestyle in fungi

Fungal pathogens cause more than a billion human infections every year, resulting in more than 1.6 million deaths annually. Understanding the natural history and evolutionary ecology of fungi is helping us understand how disease-relevant traits have repeatedly evolved. Different types and mechanisms of genetic variation have contributed to the evolution of fungal pathogenicity and specific genetic differences distinguish pathogens from non-pathogens. Insights into the traits, genetic elements, and genetic and ecological mechanisms that contribute to the evolution of fungal pathogenicity are crucial for developing strategies to both predict emergence of fungal pathogens and develop drugs to combat them.

Insights into the ecology of Schizosaccharomyces species in natural and artificial habitats

The fission yeast genus Schizosaccharomyces contains important model organisms for biological research. In particular, S. pombe is a widely used model eukaryote. So far little is known about the natural and artificial habitats of species in this genus. Finding out where S. pombe and other fission yeast species occur and how they live in their habitats can promote better understanding of their biology. Here we investigate in which substrates S. pombe, S. octosporus, S. osmophilus and S. japonicus are present. To this end about 2100 samples consisting of soil, tree sap fluxes, fresh fruit, dried fruit, honey, cacao beans, molasses and other substrates were analyzed. Effective isolation methods that allow efficient isolation of the above mentioned species were developed. Based on the frequency of isolating different fission yeast species in various substrates and on extensive literature survey, conclusions are drawn on their ecology. The results suggest that the primary habitat of S. pombe and S. octosporus is honeybee honey. Both species were also frequently detected on certain dried fruit like raisins, mango or pineapple to which they could be brought by the honey bees during ripening or during drying. While S. pombe was regularly isolated from grape mash and from fermented raw cacao beans S. octosporus was never isolated from fresh fruit. The main habitat of S. osmophilus seems to be solitary bee beebread. It was rarely isolated from raisins. S. japonicus was mainly found in forest substrates although it occurs on fruit and in fruit fermentations, too.

Fission Yeast Autophagy Machinery

Autophagy is a conserved process that delivers cytoplasmic components to the vacuole/lysosome. It plays important roles in maintaining cellular homeostasis and conferring stress resistance. In the fission yeast Schizosaccharomyces pombe, autophagy is important for cell survival under nutrient depletion and ER stress conditions. Experimental analyses of fission yeast autophagy machinery in the last 10 years have unveiled both similarities and differences in autophagosome biogenesis mechanisms between fission yeast and other model eukaryotes for autophagy research, in particular, the budding yeast Saccharomyces cerevisiae. More recently, selective autophagy pathways that deliver hydrolytic enzymes, the ER, and mitochondria to the vacuole have been discovered in fission yeast, yielding novel insights into how cargo selectivity can be achieved in autophagy. Here, we review the progress made in understanding the autophagy machinery in fission yeast.

Ecological generalism drives hyperdiversity of secondary metabolite gene clusters in xylarialean endophytes

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

OUP accepted manuscript

While thermotolerance is an attractive trait for yeasts used in industrial ethanol production, oxygen requirements of known thermotolerant species are incompatible with process requirements. Analysis of oxygen-sufficient and oxygen-limited chemostat cultures of the facultatively fermentative, thermotolerant species Ogataea parapolymorpha showed its minimum oxygen requirements to be an order of magnitude larger than those reported for the thermotolerant yeast Kluyveromyces marxianus. High oxygen requirements of O. parapolymorpha coincided with a near absence of glycerol, a key NADH/NAD⁺ redox-cofactor-balancing product in many other yeasts, in oxygen-limited cultures. Genome analysis indicated absence of orthologs of the Saccharomyces cerevisiae glycerol-3-phosphate-phosphatase genes GPP1 and GPP2. Co-feeding of acetoin, whose conversion to 2,3-butanediol enables reoxidation of cytosolic NADH, supported a 2.5-fold increase of the biomass concentration in oxygen-limited cultures. An O. parapolymorpha strain in which key genes involved in mitochondrial reoxidation of NADH were inactivated did produce glycerol, but transcriptome analysis did not reveal a clear candidate for a responsible phosphatase. Expression of S. cerevisiae GPD2, which encodes NAD⁺-dependent glycerol-3-phosphate dehydrogenase, and GPP1 supported increased glycerol production by oxygen-limited chemostat cultures of O. parapolymorpha. These results identify dependence on respiration for NADH reoxidation as a key contributor to unexpectedly high oxygen requirements of O. parapolymorpha.

The evolution of the GALactose utilization pathway in budding yeasts

The Leloir galactose utilization or GAL pathway of budding yeasts, including that of the baker's yeast Saccharomyces cerevisiae and the opportunistic human pathogen Candida albicans, breaks down the sugar galactose for energy and biomass production. The GAL pathway has long served as a model system for understanding how eukaryotic metabolic pathways, including their modes of regulation, evolve. More recently, the physical linkage of the structural genes GAL1, GAL7, and GAL10 in diverse budding yeast genomes has been used as a model for understanding the evolution of gene clustering. In this review, we summarize exciting recent work on three different aspects of this iconic pathway's evolution: gene cluster organization, GAL gene regulation, and the population genetics of the GAL pathway.

OUP accepted manuscript

Motivation

In recent years, full-genome sequences have become increasingly available and as a result many modern phylogenetic analyses are based on very long sequences, often with over 100 000 sites. Phylogenetic reconstructions of large-scale alignments are challenging for likelihood-based phylogenetic inference programs and usually require using a powerful computer cluster. Current tools for alignment trimming prior to phylogenetic analysis do not promise a significant reduction in the alignment size and are claimed to have a negative effect on the accuracy of the obtained tree.

Results

Here, we propose an artificial-intelligence-based approach, which provides means to select the optimal subset of sites and a formula by which one can compute the log-likelihood of the entire data based on this subset. Our approach is based on training a regularized Lasso-regression model that optimizes the log-likelihood prediction accuracy while putting a constraint on the number of sites used for the approximation. We show that computing the likelihood based on 5% of the sites already provides accurate approximation of the tree likelihood based on the entire data. Furthermore, we show that using this Lasso-based approximation during a tree search decreased running-time substantially while retaining the same tree-search performance.

Availability and implementation

The code was implemented in Python version 3.8 and is available through GitHub (https://github.com/noaeker/lasso_positions_sampling). The datasets used in this paper were retrieved from Zhou et al. (2018) as described in section 3.

Supplementary information

Supplementary data are available at Bioinformatics online.

Should we hail the Red King? Evolutionary consequences of a mutualistic lifestyle in genomes of lichenized ascomycetes

The Red Queen dynamic is often brought into play for antagonistic relationships. However, the coevolutionary effects of mutualistic interactions, which predict slower evolution for interacting organisms (Red King), have been investigated to a lesser extent. Lichens are a stable, mutualistic relationship of fungi and cyanobacteria and/or algae, which originated several times independently during the evolution of fungi. Therefore, they represent a suitable system to investigate the coevolutionary effect of mutualism on the fungal genome. We measured substitution rates and selective pressure of about 2000 protein-coding genes (plus the rDNA region) in two different classes of Ascomycota, each consisting of closely related lineages of lichenized and non-lichenized fungi. Our results show that independent lichenized clades are characterized by significantly slower rates for both synonymous and non-synonymous substitutions. We hypothesize that this evolutionary pattern is connected to the lichen life cycle (longer generation time of lichenized fungi) rather than a result of different selection strengths, which is described as the main driver for the Red Kind dynamic. This first empirical evidence of slower evolution in lichens provides an important insight on how biotic cooperative interactions are able to shape the evolution of symbiotic organisms.