CLOCI: unveiling cryptic fungal gene clusters with generalized detection.

Gene clusters are genomic loci that contain multiple genes that are functionally and genetically linked. Gene clusters collectively encode diverse functions, including small molecule biosynthesis, nutrient assimilation, metabolite degradation, and production of proteins essential for growth and development. Identifying gene clusters is a powerful tool for small molecule discovery and provides insight into the ecology and evolution of organisms. Current detection algorithms focus on canonical 'core' biosynthetic functions many gene clusters encode, while overlooking uncommon or unknown cluster classes. These overlooked clusters are a potential source of novel natural products and comprise an untold portion of overall gene cluster repertoires. Unbiased, function-agnostic detection algorithms therefore provide an opportunity to reveal novel classes of gene clusters and more precisely define genome organization. We present CLOCI (Co-occurrence Locus and Orthologous Cluster Identifier), an algorithm that identifies gene clusters using multiple proxies of selection for coordinated gene evolution. Our approach generalizes gene cluster detection and gene cluster family circumscription, improves detection of multiple known functional classes, and unveils non-canonical gene clusters. CLOCI is suitable for genome-enabled small molecule mining, and presents an easily tunable approach for delineating gene cluster families and homologous loci.

A global phylogenomic analysis of the shiitake genus Lentinula.

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia-Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.

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.

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.

First Report of Colletotrichum sansevieriae Causing Anthracnose of Snake Plant (Dracaena trifasciata) in Ohio and its Draft Genome.

Dracaena trifasciata (Prain) Mabb. is a popular houseplant in the United States. In September 2021, two diseased samples from two Ohio homeowners were received by the Ornamental Pathology Laboratory at The Ohio State University. Each sample included one or two detached leaves displaying circular gray water-soaked lesions scattered throughout the lamina and blighted areas with concentric rings bearing brown to black acervuli. Lesions covered between 25 and 50% of the leaf surface. Isolations were made by excising small portions of leaf tissue from the margin of the lesions, surface-disinfesting in 10% bleach for 45 s, rinsing in sterile water, and plating on potato dextrose agar (PDA). Plates were incubated at 23°C for one week. Two representative isolates, one per sample (FPH2021-5 and -6), were obtained by transferring hyphal tips to fresh PDA plates. Mycelia of both isolates were aerial, cottony, grayish-white, producing spores in a gelatinous orange matrix, and appeared gray to olivaceous-gray on the plate underside. Conidia produced by both isolates were cylindrical, single-celled, hyaline, measuring 12.02 to 18.11 (15.51) × 5.03 to 7.29 (6.14) µm (FPH2021-5; n=50) and 15.58 to 20.90 (18.39) × 5.63 to 8.27 (7.05) µm (FPH2021-6; n=50). Appressoria were globose to subglobose, single-celled, dark brown to sepia, measuring 6.62 to 13.98 (8.97) × 5.05 to 6.58 (6.58) µm (FPH2021-5; n=50), and 6.54 to 11.32 (8.63) × 4.54 to 8.94 (7.09) µm (FPH2021-6; n=50). Genomic DNA (gDNA) samples were extracted from both isolates and the internal transcribed spacer (ITS) region was amplified using primers ITS1F/ITS4 (Gardes and Bruns, 1993; White et al. 1990). GenBank BLAST sequence analysis resulted in 99.83% (FPH2021-5; GenBank Acc. No. OP410918.1) and 100% (FPH2021-6; OP410917.1) identity with 100% query coverage to the type strain of Colletotrichum sansevieriae Miho Nakam. & Ohzono MAFF239721 or Sa-1-2 (NR_152313.1; Nakamura et al. 2006). Whole genome sequencing was conducted for FPH2021-6 and the assembly was deposited in GenBank (JAOQIF000000000.1). The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and ß-tubulin (ß-tub) regions were either extracted from the genome of FPH2021-6 (OP414603.1 and OP414601.1, respectively) or amplified from FPH2021-5 gDNA using primers GDF/GDR (OP414604.1) and Bt-2b/T1 (OP414602.1), respectively (Templeton et al. 1992; Glass and Donaldson 1995; O'Donnell and Cigelnik 1997). A multilocus partitioned analysis (Chernomor et al. 2016) based on concatenated sequences of ITS, GAPDH, and ß-tub using ModelFinder (Kalyaanamoorthy et al. 2017) was performed to build a maximum likelihood tree (IQ-TREE v2.0.3; Nguyen et al. 2015), suggesting that these two isolates are phylogenetically closer to the type strain from Japan than to a previously reported isolate 1047 from Florida (Palmateer et al. 2012). To fulfill Koch's postulates, two parallel leaf sections from one 10-inch D. trifasciata 'Laurentii' plant maintained in a 1.3-liter container were selected. Three wounds were made in each section using a sterile syringe needle. A 10-µl drop of either a 1×106 conidia/ml suspension of isolate FPH2021-6 or sterile water was placed on each wound. The plant was covered with a plastic bag for two days post-inoculation (DPI) and maintained in a greenhouse at 25°C with a 12- h photoperiod. The experiment was conducted twice. Grayish water-soaked lesions, acervuli, and leaf blight were observed on the inoculated sections 3, 10, and 14 DPI, respectively, while no symptoms appeared on the sections treated with sterile water. C. sansevieriae was re-isolated from the lesions and confirmed to be identical to the original isolate based on ITS sequencing and morphological examinations. To the best of our knowledge, this is the first report of C. sansevieriae on D. trifasciata in Ohio and the first genome draft of an isolate from the United States. Availability of whole-genome sequence data is paramount for resolving species identification in this highly diverse fungal genus, and a powerful tool to conduct comparative genomic analyses in the future.

A metabolomics-guided approach to discover Fusarium graminearum metabolites after removal of a repressive histone modification.

Many secondary metabolites are produced by biosynthetic gene clusters (BGCs) that are repressed during standard growth conditions, which complicates the discovery of novel bioactive compounds. In the genus Fusarium, many BGCs reside in chromatin enriched for trimethylated histone 3 lysine 27 (H3K27me3), a modification correlated with transcriptional gene silencing. Here we report on our progress in assigning metabolites to genes by using a strain lacking the H3K27 methyltransferase, Kmt6. To guide isolation efforts, we coupled genetics to multivariate analysis of liquid chromatography-mass spectrometry (LCMS) data from both wild type and kmt6, which allowed identification of compounds previously unknown from F. graminearum. We found low molecular weight, amino acid-derived metabolites (N-ethyl anthranilic acid, N-phenethylacetamide, N-acetyltryptamine). We identified one new compound, protofusarin, as derived from fusarin biosynthesis. Similarly, we isolated large amounts of fusaristatin A, gibepyrone A, and fusarpyrones A and B, simply by using the kmt6 mutant, instead of having to optimize growth media. To increase the abundance of metabolites underrepresented in wild type, we generated kmt6 fus1 double mutants and discovered tricinolone and tricinolonoic acid, two new sesquiterpenes belonging to the tricindiol class. Our approach allows rapid visualization and analyses of the genetically induced changes in metabolite production, and discovery of new molecules by a combination of chemical and genetic dereplication. Of 22 fungal metabolites identified here, 10 compounds had not been reported from F. graminearum before. We show that activating silent metabolic pathways by mutation of a repressive chromatin modification enzyme can result in the discovery of new chemistry even in a well-studied organism, and helps to connect new or known small molecules to the BGCs responsible for their production.

Live tracking of a plant pathogen outbreak reveals rapid and successive, multidecade plasmid reduction.

Quickly understanding the genomic changes that lead to pathogen emergence is necessary to launch mitigation efforts and reduce harm. In this study, we tracked in real time a 2022 bacterial plant disease outbreak in U.S. geraniums (Pelargonium × hortorum) caused by Xhp2022, a novel lineage of Xanthomonas hortorum. Genomes from 31 Xhp2022 isolates from seven states showed limited chromosomal variation and all contained a single plasmid (p93). Time tree and single nucleotide polymorphism whole-genome analysis estimated that Xhp2022 emerged within the last decade. The phylogenomic analysis determined that p93 resulted from the cointegration of three plasmids (p31, p45, and p66) found sporadically across isolates from previous outbreaks. Although p93 had a 49 kb nucleotide reduction, it retained putative fitness genes, which became predominant in the 2022 outbreak. Overall, we demonstrated, through rapid whole-genome sequencing and analysis, a recent, traceable event of genome reduction for niche adaptation typically observed over millennia in obligate and fastidious pathogens.IMPORTANCEThe geranium industry, valued at $4 million annually, faces an ongoing Xanthomonas hortorum pv. pelargonii (Xhp) pathogen outbreak. To track and describe the outbreak, we compared the genome structure across historical and globally distributed isolates. Our research revealed Xhp population has not had chromosome rearrangements since 1974 and has three distinct plasmids. In 2012, we found all three plasmids in individual Xhp isolates. However, in 2022, the three plasmids co-integrated into one plasmid named p93. p93 retained putative fitness genes but lost extraneous genomic material. Our findings show that the 2022 strain group of the bacterial plant pathogen Xanthomonas hortorum underwent a plasmid reduction. We also observed several Xanthomonas species from different years, hosts, and continents have similar plasmids to p93, possibly due to shared agricultural settings. We noticed parallels between genome efficiency and reduction that we see across millennia with obligate parasites with increased niche specificity.

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.

Whole-Genome Sequence Data for the Holotype Strain of Diaporthe ilicicola, a Fungus Associated with Latent Fruit Rot in Deciduous Holly.

Diaporthe ilicicola is a newly described fungal species that is associated with latent fruit rot in deciduous holly. This announcement provides a whole-genome assembly and annotation for this plant pathogen, which will inform research on its parasitism and identification of gene clusters involved in the production of bioactive metabolites.

Draft Genome Sequence of the Termite-Associated "Cuckoo Fungus," Athelia (Fibularhizoctonia) sp. TMB Strain TB5.

Atheliales is a diverse order of crust-forming Basidiomycota fungi. Here, we report the draft genome of the "cuckoo fungus," Athelia (Fibularhizoctonia) sp. TMB strain TB5 (Atheliales), which forms termite-egg-mimicking sclerotia that termites tend. We further compare its repertoire of psilocybin gene homologs to homologs previously reported for Fibularhizoctonia psychrophila.