Systematic identification of cargo-mobilizing genetic elements reveals new dimensions of eukaryotic diversity.

Cargo-mobilizing mobile elements (CMEs) are genetic entities that faithfully transpose diverse protein coding sequences. Although common in bacteria, we know little about eukaryotic CMEs because no appropriate tools exist for their annotation. For example, Starships are giant fungal CMEs whose functions are largely unknown because they require time-intensive manual curation. To address this knowledge gap, we developed starfish, a computational workflow for high-throughput eukaryotic CME annotation. We applied starfish to 2 899 genomes of 1 649 fungal species and found that starfish recovers known Starships with 95% combined precision and recall while expanding the number of annotated elements ten-fold. Extant Starship diversity is partitioned into 11 families that differ in their enrichment patterns across fungal classes. Starship cargo changes rapidly such that elements from the same family differ substantially in their functional repertoires, which are predicted to contribute to diverse biological processes such as metabolism. Many elements have convergently evolved to insert into 5S rDNA and AT-rich sequence while others integrate into random locations, revealing both specialist and generalist strategies for persistence. Our work establishes a framework for advancing mobile element biology and provides the means to investigate an emerging dimension of eukaryotic genetic diversity, that of genomes within genomes.

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.

The Architecture of Metabolism Maximizes Biosynthetic Diversity in the Largest Class of Fungi.

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or "specialized" metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

The Canadian Fungal Research Network: current challenges and future opportunities.

Fungi critically impact the health and function of global ecosystems and economies. In Canada, fungal researchers often work within silos defined by subdiscipline and institutional type, complicating the collaborations necessary to understand the impacts fungi have on the environment, economy, and plant and animal health. Here, we announce the establishment of the Canadian Fungal Research Network (CanFunNet, https://fungalresearch.ca), whose mission is to strengthen and promote fungal research in Canada by facilitating dialogue among scientists. We summarize the challenges and opportunities for Canadian fungal research that were discussed at CanFunNet's inaugural meeting in 2019, and identify 4 priorities for our community: (i) increasing collaboration among scientists, (ii) studying diversity in the context of ecological disturbance, (iii) preserving culture collections in the absence of sustained funding, and (iv) leveraging diverse expertise to attract trainees. We have gathered additional information to support our recommendations, including a survey identifying underrepresentation of fungal-related courses at Canadian universities, a list of Canadian fungaria and culture collections, and a case study of a human fungal pathogen outbreak. We anticipate that these discussions will help prioritize fungal research in Canada, and we welcome all researchers to join this nationwide effort to enhance knowledge dissemination and funding advocacy.

Fungal adaptation to plant defences through convergent assembly of metabolic modules.

The ongoing diversification of plant defence compounds exerts dynamic selection pressures on the microorganisms that colonize plant tissues. Evolutionary processes that generate resistance towards these compounds increase microbial fitness by giving access to plant resources and increasing pathogen virulence. These processes entail sequence-based mechanisms that result in adaptive gene functions, and combinatorial mechanisms that result in novel syntheses of existing gene functions. However, the priority and interactions among these processes in adaptive resistance remain poorly understood. Using a combination of molecular genetic and computational approaches, we investigated the contributions of sequence-based and combinatorial processes to the evolution of fungal metabolic gene clusters encoding stilbene cleavage oxygenases (SCOs), which catalyse the degradation of biphenolic plant defence compounds known as stilbenes into monophenolic molecules. We present phylogenetic evidence of convergent assembly among three distinct types of SCO gene clusters containing alternate combinations of phenolic catabolism. Multiple evolutionary transitions between different cluster types suggest recurrent selection for distinct gene assemblages. By comparison, we found that the substrate specificities of heterologously expressed SCO enzymes encoded in different clusters types were all limited to stilbenes and related molecules with a 4'-OH group, and differed modestly in substrate range and activity under the experimental conditions. Together, this work suggests a primary role for genome structural rearrangement, and the importance of enzyme modularity, in promoting fungal metabolic adaptation to plant defence chemistry.

Endophyte genomes support greater metabolic gene cluster diversity compared with non-endophytes in Trichoderma.

Trichoderma is a cosmopolitan genus with diverse lifestyles and nutritional modes, including mycotrophy, saprophytism, and endophytism. Previous research has reported greater metabolic gene repertoires in endophytic fungal species compared to closely-related non-endophytes. However, the extent of this ecological trend and its underlying mechanisms are unclear. Some endophytic fungi may also be mycotrophs and have one or more mycoparasitism mechanisms. Mycotrophic endophytes are prominent in certain genera like Trichoderma, therefore, the mechanisms that enable these fungi to colonize both living plants and fungi may be the result of expanded metabolic gene repertoires. Our objective was to determine what, if any, genomic features are overrepresented in endophytic fungi genomes in order to undercover the genomic underpinning of the fungal endophytic lifestyle. Here we compared metabolic gene cluster and mycoparasitism gene diversity across a dataset of thirty-eight Trichoderma genomes representing the full breadth of environmental Trichoderma's diverse lifestyles and nutritional modes. We generated four new Trichoderma endophyticum genomes to improve the sampling of endophytic isolates from this genus. As predicted, endophytic Trichoderma genomes contained, on average, more total biosynthetic and degradative gene clusters than non-endophytic isolates, suggesting that the ability to create/modify a diversity of metabolites potential is beneficial or necessary to the endophytic fungi. Still, once the phylogenetic signal was taken in consideration, no particular class of metabolic gene cluster was independently associated with the Trichoderma endophytic lifestyle. Several mycoparasitism genes, but no chitinase genes, were associated with endophytic Trichoderma genomes. Most genomic differences between Trichoderma lifestyles and nutritional modes are difficult to disentangle from phylogenetic divergences among species, suggesting that Trichoderma genomes maybe particularly well-equipped for lifestyle plasticity. We also consider the role of endophytism in diversifying secondary metabolism after identifying the horizontal transfer of the ergot alkaloid gene cluster to Trichoderma.

Giant mobile elements: Agents of multivariate phenotypic evolution in fungi.

Fungal mobile genetic elements are typically small, vertically inherited, and are not known to encode adaptive traits. A new study documents HEPHAESTUS - a large, horizontally transferred, cargo-carrying mobile element that confers tolerance to several metals in fungi.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom-forming fungi (Agaricomycetes).

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a 'developmental hourglass,' act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Repeated gain and loss of a single gene modulates the evolution of vascular plant pathogen lifestyles.

Vascular plant pathogens travel long distances through host veins, leading to life-threatening, systemic infections. In contrast, nonvascular pathogens remain restricted to infection sites, triggering localized symptom development. The contrasting features of vascular and nonvascular diseases suggest distinct etiologies, but the basis for each remains unclear. Here, we show that the hydrolase CbsA acts as a phenotypic switch between vascular and nonvascular plant pathogenesis. cbsA was enriched in genomes of vascular phytopathogenic bacteria in the family Xanthomonadaceae and absent in most nonvascular species. CbsA expression allowed nonvascular Xanthomonas to cause vascular blight, while cbsA mutagenesis resulted in reduction of vascular or enhanced nonvascular symptom development. Phylogenetic hypothesis testing further revealed that cbsA was lost in multiple nonvascular lineages and more recently gained by some vascular subgroups, suggesting that vascular pathogenesis is ancestral. Our results overall demonstrate how the gain and loss of single loci can facilitate the evolution of complex ecological traits.

Metabolic gene clusters, fungal diversity, and the generation of accessory functions.

Ecological interactions are largely determined by adaptive traits deemed 'accessory'. In plants, fungi, and bacteria, such traits mainly comprise metabolic pathways that produce or transform diverse molecules. While accessory metabolic pathways are pervasive, it is often difficult to identify their genetic bases. Recently, in-depth descriptions of metabolic gene clusters (MGCs), which encode discrete metabolic pathways, have greatly simplified the characterization of genotype-phenotype maps, yet questions of how this genome architecture relates to the evolution of accessory functions remain. Fungi are uniquely positioned to spearhead investigations into these dynamics because they display gradients in clustering across pathways and taxa. This review will focus on the role of MGCs as both agents and consequences of the accessory function evolution that underpins fungal diversification.

Whole-Genome Sequence of the Phlox Powdery Mildew Pathogen Golovinomyces magnicellulatus Strain FPH2017-1.

Powdery mildew (PM) fungi are obligate biotrophs capable of infecting diverse plant hosts, ranging from monocotyledonous agricultural crops to dicotyledonous ornamental crops. The PM lifestyle poses significant challenges for studying these pathogens in isolation from their host. We present a draft genome of Golovinomyces magnicellulatus, a host-specific PM on Phlox species.

Psychoactive plant- and mushroom-associated alkaloids from two behavior modifying cicada pathogens.

Entomopathogenic fungi routinely kill their hosts before releasing infectious spores, but a few species keep insects alive while sporulating, which enhances dispersal. Transcriptomics- and metabolomics-based studies of entomopathogens with post-mortem dissemination from their parasitized hosts have unraveled infection processes and host responses. However, the mechanisms underlying active spore transmission by Entomophthoralean fungi in living insects remain elusive. Here we report the discovery, through metabolomics, of the plant-associated amphetamine, cathinone, in four Massospora cicadina-infected periodical cicada populations, and the mushroom-associated tryptamine, psilocybin, in annual cicadas infected with Massospora platypediae or Massospora levispora, which likely represent a single fungal species. The absence of some fungal enzymes necessary for cathinone and psilocybin biosynthesis along with the inability to detect intermediate metabolites or gene orthologs are consistent with possibly novel biosynthesis pathways in Massospora. The neurogenic activities of these compounds suggest the extended phenotype of Massospora that modifies cicada behavior to maximize dissemination is chemically-induced.

Specialized plant biochemistry drives gene clustering in fungi.

The fitness and evolution of prokaryotes and eukaryotes are affected by the organization of their genomes. In particular, the physical clustering of genes can coordinate gene expression and can prevent the breakup of co-adapted alleles. Although clustering may thus result from selection for phenotype optimization and persistence, the impact of environmental selection pressures on eukaryotic genome organization has rarely been systematically explored. Here, we investigated the organization of fungal genes involved in the degradation of phenylpropanoids, a class of plant-produced secondary metabolites that mediate many ecological interactions between plants and fungi. Using a novel gene cluster detection method, we identified 1110 gene clusters and many conserved combinations of clusters in a diverse set of fungi. We demonstrate that congruence in genome organization over small spatial scales is often associated with similarities in ecological lifestyle. Additionally, we find that while clusters are often structured as independent modules with little overlap in content, certain gene families merge multiple modules into a common network, suggesting they are important components of phenylpropanoid degradation strategies. Together, our results suggest that phenylpropanoids have repeatedly selected for gene clustering in fungi, and highlight the interplay between genome organization and ecological evolution in this ancient eukaryotic lineage.

Horizontal gene cluster transfer increased hallucinogenic mushroom diversity.

Secondary metabolites are a heterogeneous class of chemicals that often mediate interactions between species. The tryptophan-derived secondary metabolite, psilocin, is a serotonin receptor agonist that induces altered states of consciousness. A phylogenetically disjunct group of mushroom-forming fungi in the Agaricales produce the psilocin prodrug, psilocybin. Spotty phylogenetic distributions of fungal compounds are sometimes explained by horizontal transfer of metabolic gene clusters among unrelated fungi with overlapping niches. We report the discovery of a psilocybin gene cluster in three hallucinogenic mushroom genomes, and evidence for its horizontal transfer between fungal lineages. Patterns of gene distribution and transmission suggest that synthesis of psilocybin may have provided a fitness advantage in the dung and late wood-decay fungal niches, which may serve as reservoirs of fungal indole-based metabolites that alter behavior of mycophagous and wood-eating invertebrates. These hallucinogenic mushroom genomes will serve as models in neurochemical ecology, advancing the (bio)prospecting and synthetic biology of novel neuropharmaceuticals.

Genes of the de novo and Salvage Biosynthesis Pathways of Vitamin B6 are Regulated under Oxidative Stress in the Plant Pathogen Rhizoctonia solani.

Vitamin B6 is recognized as an important cofactor required for numerous metabolic enzymes, and has been shown to act as an antioxidant and play a role in stress responses. It can be synthesized through two different routes: salvage and de novo pathways. However, little is known about the possible function of the vitamin B6 pathways in the fungal plant pathogen Rhizoctonia solani. Using genome walking, the de novo biosynthetic pathway genes; RsolPDX1 and RsolPDX2 and the salvage biosynthetic pathway gene, RsolPLR were sequenced. The predicted amino acid sequences of the three genes had high degrees of similarity to other fungal PDX1, PDX2, and PLR proteins and are closely related to other R. solani anastomosis groups. We also examined their regulation when subjected to reactive oxygen species (ROS) stress inducers, the superoxide generator paraquat, or H2O2, and compared it to the well-known antioxidant genes, catalase and glutathione-S-transferase (GST). The genes were differentially regulated with transcript levels as high as 33 fold depending on the gene and type of stress reflecting differences in the type of damage induced by ROS. Exogenous addition of the vitamers PN or PLP in culture medium significantly induced the transcription of the vitamin B6 de novo encoding genes as early as 0.5 hour post treatment (HPT). On the other hand, transcription of RsolPLR was vitamer-specific; a down regulation upon supplementation of PN and upregulation with PLP. Our results suggest that accumulation of ROS in R. solani mycelia is linked to transcriptional regulation of the three genes and implicate the vitamin B6 biosynthesis machinery in R. solani, similar to catalases and GST, as an antioxidant stress protector against oxidative stress.