Transposable Elements versus the Fungal Genome: Impact on Whole-Genome Architecture and Transcriptional Profiles

The dynamics of fungal genome organization and its impact on host adaptation and antifungal resistance

Fungi are a diverse kingdom characterized by remarkable genomic plasticity that facilitates pathogenicity and adaptation to adverse environmental conditions. In this review, we delve into the dynamic organization of fungal genomes and its implications for host adaptation and antifungal resistance. We examine key features and the heterogeneity of genomes across different fungal species, including but not limited to their chromosome content, DNA composition, distribution and arrangement of their content across chromosomes, and other major traits. We further highlight how this variability in genomic traits influences their virulence and adaptation to adverse conditions. Fungal genomes exhibit large variations in size, gene content, and structural features, such as the abundance of transposable elements (TEs), compartmentalization into gene-rich and TE-rich regions, and the presence or absence of dispensable chromosomes. Genomic structural variations are equally diverse in fungi, ranging from whole-chromosome duplications that may enhance tolerance to antifungal compounds, to targeted deletion of effector encoding genes that may promote virulence. Finally, the often-overlooked fungal mitochondrial genomes can also affect virulence and resistance to fungicides. Such and other features of fungal genome organization are reviewed and discussed in the context of host-microbe interactions and antifungal resistance.

How genomics can help unravel the evolution of endophytic fungi

Endophytic fungi (EFs) form intimate associations with plants, residing within their tissues without causing apparent harm. Understanding the evolution of endophytic fungal genomes is essential for uncovering the mechanisms that drive their symbiotic relationships with host plants. This review explores the dynamic interactions between EFs and host plants, focusing on the evolutionary processes that shape their genomes. We highlighted key genomic adaptations promoting their endophytic lifestyle, including genes involved in plant cell wall degradation, secondary metabolite production, and stress tolerance. By combining genomic data with ecological and physiological information, this review provides a comprehensive understanding of the coevolutionary dynamics between EFs and host plants. Moreover, it provides insights that help elucidate the complex interdependencies governing their symbiotic interactions.

Analyses of Transposable Elements in Arbuscular Mycorrhizal Fungi Support Evolutionary Parallels With Filamentous Plant Pathogens

Transposable elements are repetitive DNA sequences that excise or create copies that are inserted elsewhere in the genome. Their expansion shapes genome variability and evolution by impacting gene expression and rearrangement rates. Arbuscular mycorrhizal fungi are beneficial plant symbionts with large, transposable element-rich genomes, and recent findings showed these elements vary significantly in abundance, evolution, and regulation among model strains. Here, we aimed to obtain a more comprehensive understanding of transposable element function and evolution in arbuscular mycorrhizal fungi by investigating assembled genomes from representatives of all known families. We uncovered multiple, family-specific bursts of insertions in different species, indicating variable past and ongoing transposable element activity contributing to the diversification of arbuscular mycorrhizal fungi lineages. We also found that transposable elements are preferentially located within and around candidate effectors/secreted proteins, as well as in proximity to promoters. Altogether, these findings support the role of transposable elements in promoting the diversity in proteins involved in molecular dialogs with hosts and, more generally, in driving gene regulation. The mechanisms of transposable element evolution we observed in these prominent plant symbionts bear striking similarities to those of many filamentous plant pathogens.

Genomic features of lichen-associated black fungi

Lichens are mutualistic associations consisting of a primary fungal host, and one to few primary phototrophic symbiont(s), usually a green alga and/or a cyanobacterium. They form complex thallus structures, which provide unique and stable habitats for many other microorganisms. Frequently isolated from lichens are the so-called black fungi, or black yeasts, which are mainly characterized by melanized cell walls and extremophilic lifestyles. It is presently unclear in which ways these fungi interact with other members of the lichen symbiosis. Genomic resources of lichen-associated black fungi are needed to better understand the physiological potential of these fungi and shed light on the complexity of the lichen consortium. Here, we present high-quality genomes of 14 black fungal lineages, isolated from lichens of the rock-dwelling genus Umbilicaria. Nine of the lineages belong to the Eurotiomycetes (Chaetothyriales), four to the Dothideomycetes, and one to the Arthoniomycetes, representing the first genome of a black fungus in this class. The PacBio-based assemblies are highly contiguous (5-42 contigs per genome, mean coverage of 79-502, N50 of 1.0-7.3 mega-base-pair (Mb), Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness generally ≥95.4%). Most contigs are flanked by a telomere sequence, suggesting we achieved near chromosome-level assemblies. Genome sizes range between 26 and 44 Mb. Transcriptome-based annotations yielded ~11,000-18,000 genes per genome. We analyzed genome content with respect to repetitive elements, biosynthetic genes, and effector genes. Each genome contained a polyketide synthase gene related to the dihydroxynaphthalene-melanin pathway. This research provides insights into genome content and metabolic potential of these relatively unknown, but frequently encountered lichen associates.

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.

Comparative genome analysis and the genome-shaping role of long terminal repeat retrotransposons in the evolutionary divergence of fungal pathogens Blastomyces dermatitidis and Blastomyces gilchristii

Blastomyces dermatitidis and Blastomyces gilchristii are cryptic species of fungi that cause blastomycosis, an often severe disease involving pulmonary infection capable of systemic dissemination. While these species appear morphologically identical, differences exist in the genetic makeup, geographical range, and possibly the clinical presentation of infection. Here, we show genetic divergence between the cryptic species through both a Blastomyces species tree constructed from orthologous protein sequences and whole genome single-nucleotide variant phylogenomic analysis. Following linked-read sequencing and de novo genome assembly, we characterized and compared the genomes of 3 B. dermatitidis and 3 B. gilchristii isolates. The B. gilchristii genomes (73.25-75.4 Mb) were ∼8 Mb larger than the B. dermatitidis genomes (64.88-66.61 Mb). Average nucleotide identity was lower between genomes of different species than genomes of the same species, yet functional classification of genes suggested similar proteomes. The most striking difference involved long terminal repeat retrotransposons. Although the same retrotransposon elements were detected in the genomes, the quantity of elements differed between the 2 species. Gypsy retrotransposon content was significantly higher in B. gilchristii (38.04-39.26 Mb) than in B. dermatitidis (30.85-32.40 Mb), accounting for the majority of genome size difference between species. Age estimation and phylogenetic analysis of the reverse transcriptase domains suggested that these retrotransposons are relatively ancient, with genome insertion predating the speciation of B. dermatitidis and B. gilchristii. We postulate that different trajectories of genome contraction led to genetic incompatibility, reproductive isolation, and speciation, highlighting the role of transposable elements in fungal evolution.

Genome analysis of the esca-associated Basidiomycetes Fomitiporia mediterranea, Fomitiporia polymorpha, Inonotus vitis, and Tropicoporus texanus reveals virulence factor repertoires characteristic of white-rot fungi

Some Basidiomycete fungi are important plant pathogens, and certain species have been associated with the grapevine trunk disease esca. We present the genomes of 4 species associated with esca: Fomitiporia mediterranea, Fomitiporia polymorpha, Tropicoporus texanus, and Inonotus vitis. We generated high-quality phased genome assemblies using long-read sequencing. The genomic and functional comparisons identified potential virulence factors, suggesting their roles in disease development. Similar to other white-rot fungi known for their ability to degrade lignocellulosic substrates, these 4 genomes encoded a variety of lignin peroxidases and carbohydrate-active enzymes (CAZymes) such as CBM1, AA9, and AA2. The analysis of gene family expansion and contraction revealed dynamic evolutionary patterns, particularly in genes related to secondary metabolite production, plant cell wall decomposition, and xenobiotic degradation. The availability of these genomes will serve as a reference for further studies of diversity and evolution of virulence factors and their roles in esca symptoms and host resistance.

Discovering the hidden function in fungal genomes

New molecular technologies have helped unveil previously unexplored facets of the genome beyond the canonical proteome, including microproteins and short ORFs, products of alternative splicing, regulatory non-coding RNAs, as well as transposable elements, cis-regulatory DNA, and other highly repetitive regions of DNA. In this Review, we highlight what is known about this 'hidden genome' within the fungal kingdom. Using well-established model systems as a contextual framework, we describe key elements of this hidden genome in diverse fungal species, and explore how these factors perform critical functions in regulating fungal metabolism, stress tolerance, and pathogenesis. Finally, we discuss new technologies that may be adapted to further characterize the hidden genome in fungi.

Comparative Genomics of Fungi in Nectriaceae Reveals Their Environmental Adaptation and Conservation Strategies

This study presents the first genome assembly of the freshwater saprobe fungus Neonectria lugdunensis and a comprehensive phylogenomics analysis of the Nectriaceae family, examining genomic traits according to fungal lifestyles. The Nectriaceae family, one of the largest in Hypocreales, includes fungi with significant ecological roles and economic importance as plant pathogens, endophytes, and saprobes. The phylogenomics analysis identified 2684 single-copy orthologs, providing a robust evolutionary framework for the Nectriaceae family. We analyzed the genomic characteristics of 17 Nectriaceae genomes, focusing on their carbohydrate-active enzymes (CAZymes), biosynthetic gene clusters (BGCs), and adaptations to environmental temperatures. Our results highlight the adaptation mechanisms of N. lugdunensis, emphasizing its capabilities for plant litter degradation and enzyme activity in varying temperatures. The comparative genomics of different Nectriaceae lifestyles revealed significant differences in genome size, gene content, repetitive elements, and secondary metabolite production. Endophytes exhibited larger genomes, more effector proteins, and BGCs, while plant pathogens had higher thermo-adapted protein counts, suggesting greater resilience to global warming. In contrast, the freshwater saprobe shows less adaptation to warmer temperatures and is important for conservation goals. This study underscores the importance of understanding fungal genomic adaptations to predict ecosystem impacts and conservation targets in the face of climate change.

Comparative Analysis of the Mitochondrial Genome Sequences of Diaporthe longicolla (syn. Phomopsis longicolla) Isolates Causing Phomopsis Seed Decay in Soybean

Diaporthe longicolla (syn. Phomopsis longicolla) is an important seed-borne fungal pathogen and the primary cause of Phomopsis seed decay (PSD) in soybean. PSD is one of the most devastating seed diseases, reducing soybean seed quality and yield worldwide. As part of a genome sequencing project on the fungal Diaporthe-Phomopsis complex, draft genomes of eight D. longicolla isolates were sequenced and assembled. Sequences of mitochondrial genomes were extracted and analyzed. The circular mitochondrial genomes ranged from 52,534 bp to 58,280 bp long, with a mean GC content of 34%. A total of 14 core protein-coding genes, 23 tRNA, and 2 rRNA genes were identified. Introns were detected in the genes of atp6, cob, cox1, cox2, cox3, nad1, nad2, nad5, and rnl. Three isolates (PL7, PL10, and PL185E) had more introns than other isolates. Approximately 6.4% of the mitochondrial genomes consist of repetitive elements. Moreover, 48 single-nucleotide polymorphisms (SNPs) and were identified. The mitochondrial genome sequences of D. longicolla will be useful to further study the molecular basis of seed-borne pathogens causing seed diseases, investigate genetic variation among isolates, and develop improved control strategies for Phomopsis seed decay of soybean.

Whole-genome sequencing of marine water-derived Curvularia verruculosa KHW-7: a pioneering study

Marine microorganisms are renowned for being a rich source of new secondary metabolites that are significant to humans. The fungi strain KHW-7 was isolated from the seawater collected from the Gulf of Khambhat, India, and identified as Curvularia verruculosa KHW-7. On a next-generation sequencing platform, C. verruculosa KHW-7's whole-genome sequencing (WGS) and gene annotation were carried out using several bioinformatic methods. The 31.59 MB genome size, 52.3% GC, and 158 bp mean read length were discovered using WGS. This genome also contained 9,745 protein-coding genes, including 852 secreted proteins and 2048 transmembrane proteins. The antiSMASH algorithm used to analyze genomes found 25 secondary metabolite biosynthetic gene clusters (BGCs) that are abundant in terpene, non-ribosomal peptide synthetase (NRPS), and polyketides type 1 (T1PKS). To our knowledge, this is the first whole-genome sequence report of C. verruculosa. The WGS analysis of C. verruculosa KHW-7 indicated that this marine-derived fungus could be an efficient generator of bioactive secondary metabolites and an important industrial enzyme, both of which demand further investigation and development.

Transposable elements impact the population divergence of rice blast fungus Magnaporthe oryzae

Dynamic transposition of transposable elements (TEs) in fungal pathogens has significant impact on genome stability, gene expression, and virulence to the host. In Magnaporthe oryzae, genome plasticity resulting from TE insertion is a major driving force leading to the rapid evolution and diversification of this fungus. Despite their importance in M. oryzae population evolution and divergence, our understanding of TEs in this context remains limited. Here, we conducted a genome-wide analysis of TE transposition dynamics in the 11 most abundant TE families in M. oryzae populations. Our results show that these TEs have specifically expanded in recently isolated M. oryzae rice populations, with the presence/absence polymorphism of TE insertions highly concordant with population divergence on Geng/Japonica and Xian/Indica rice cultivars. Notably, the genes targeted by clade-specific TEs showed clade-specific expression patterns and are involved in the pathogenic process, suggesting a transcriptional regulation of TEs on targeted genes. Our study provides a comprehensive analysis of TEs in M. oryzae populations and demonstrates a crucial role of recent TE bursts in adaptive evolution and diversification of the M. oryzae rice-infecting lineage.

IMPORTANCE

Magnaporthe oryzae is the causal agent of the destructive blast disease, which caused massive loss of yield annually worldwide. The fungus diverged into distinct clades during adaptation toward the two rice subspecies, Xian/Indica and Geng/Japonica. Although the role of TEs in the adaptive evolution was well established, mechanisms underlying how TEs promote the population divergence of M. oryzae remain largely unknown. In this study, we reported that TEs shape the population divergence of M. oryzae by differentially regulating gene expression between Xian/Indica-infecting and Geng/Japonica-infecting populations. Our results revealed a TE insertion-mediated gene expression adaption that led to the divergence of M. oryzae population infecting different rice subspecies.

Deciphering the Genomic Landscape and Virulence Mechanisms of the Wheat Powdery Mildew Pathogen Blumeria graminis f. sp. tritici Wtn1: Insights from Integrated Genome Assembly and Conidial Transcriptomics

A high-quality genome sequence from an Indian isolate of Blumeria graminis f. sp. tritici Wtn1, a persistent threat in wheat farming, was obtained using a hybrid method. The assembly of over 9.24 million DNA-sequence reads resulted in 93 contigs, totaling a 140.61 Mb genome size, potentially encoding 8480 genes. Notably, more than 73.80% of the genome, spanning approximately 102.14 Mb, comprises retro-elements, LTR elements, and P elements, influencing evolution and adaptation significantly. The phylogenomic analysis placed B. graminis f. sp. tritici Wtn1 in a distinct monocot-infecting clade. A total of 583 tRNA anticodon sequences were identified from the whole genome of the native virulent strain B. graminis f. sp. tritici, which comprises distinct genome features with high counts of tRNA anticodons for leucine (70), cysteine (61), alanine (58), and arginine (45), with only two stop codons (Opal and Ochre) present and the absence of the Amber stop codon. Comparative InterProScan analysis unveiled "shared and unique" proteins in B. graminis f. sp. tritici Wtn1. Identified were 7707 protein-encoding genes, annotated to different categories such as 805 effectors, 156 CAZymes, 6102 orthologous proteins, and 3180 distinct protein families (PFAMs). Among the effectors, genes like Avra10, Avrk1, Bcg-7, BEC1005, CSEP0105, CSEP0162, BEC1016, BEC1040, and HopI1 closely linked to pathogenesis and virulence were recognized. Transcriptome analysis highlighted abundant proteins associated with RNA processing and modification, post-translational modification, protein turnover, chaperones, and signal transduction. Examining the Environmental Information Processing Pathways in B. graminis f. sp. tritici Wtn1 revealed 393 genes across 33 signal transduction pathways. The key pathways included yeast MAPK signaling (53 genes), mTOR signaling (38 genes), PI3K-Akt signaling (23 genes), and AMPK signaling (21 genes). Additionally, pathways like FoxO, Phosphatidylinositol, the two-component system, and Ras signaling showed significant gene representation, each with 15-16 genes, key SNPs, and Indels in specific chromosomes highlighting their relevance to environmental responses and pathotype evolution. The SNP and InDel analysis resulted in about 3.56 million variants, including 3.45 million SNPs, 5050 insertions, and 5651 deletions within the whole genome of B. graminis f. sp. tritici Wtn1. These comprehensive genome and transcriptome datasets serve as crucial resources for understanding the pathogenicity, virulence effectors, retro-elements, and evolutionary origins of B. graminis f. sp. tritici Wtn1, aiding in developing robust strategies for the effective management of wheat powdery mildew.

Orthoptera-TElib: a library of Orthoptera transposable elements for TE annotation

Transposable elements (TEs) are a major component of eukaryotic genomes and are present in almost all eukaryotic organisms. TEs are highly dynamic between and within species, which significantly affects the general applicability of the TE databases. Orthoptera is the only known group in the class Insecta with a significantly enlarged genome (0.93-21.48 Gb). When analyzing the large genome using the existing TE public database, the efficiency of TE annotation is not satisfactory. To address this limitation, it becomes imperative to continually update the available TE resource library and the need for an Orthoptera-specific library as more insect genomes are publicly available. Here, we used the complete genome data of 12 Orthoptera species to de novo annotate TEs, then manually re-annotate the unclassified TEs to construct a non-redundant Orthoptera-specific TE library: Orthoptera-TElib. Orthoptera-TElib contains 24,021 TE entries including the re-annotated results of 13,964 unknown TEs. The naming of TE entries in Orthoptera-TElib adopts the same naming as RepeatMasker and Dfam and is encoded as the three-level form of "level1/level2-level3". Orthoptera-TElib can be directly used as an input reference database and is compatible with mainstream repetitive sequence analysis software such as RepeatMasker and dnaPipeTE. When analyzing TEs of Orthoptera species, Orthoptera-TElib performs better TE annotation as compared to Dfam and Repbase regardless of using low-coverage sequencing or genome assembly data. The most improved TE annotation result is Angaracris rhodopa, which has increased from 7.89% of the genome to 53.28%. Finally, Orthoptera-TElib is stored in Sqlite3 for the convenience of data updates and user access.

PiggyBac Transposon Mining in the Small Genomes of Animals

TEs, including DNA transposons, are major contributors of genome expansions, and have played a very significant role in shaping the evolution of animal genomes, due to their capacity to jump from one genomic position to the other. In this study, we investigated the evolution landscapes of PB transposons, including their distribution, diversity, activity and structure organization in 79 species of small (compact) genomes of animals comprising both vertebrate and invertebrates. Overall, 212 PB transposon types were detected from almost half (37) of the total number of the small genome species (79) investigated. The detected PB transposon types, which were unevenly distributed in various genera and phyla, have been classified into seven distinct clades or families with good bootstrap support (>80%). The PB transposon types that were identified have a length ranging from 1.23 kb to 9.51 kb. They encode transposases of approximately ≥500 amino acids in length, and possess terminal inverted repeats (TIRs) ranging from 4 bp to 24 bp. Though some of the transposon types have long TIRs (528 bp), they still maintain the consistent and reliable 4 bp target site duplication (TSD) of TTAA. However, PiggyBac-2_Cvir transposon originating from the Crassostrea virginica species exhibits a unique TSD of TATG. The TIRs of the transposons in all the seven families display high divergence, with a highly conserved 5' end motif. The core transposase domains (DDD) were better conserved among the seven different families compared to the other protein domains, which were less prevalent in the vertebrate genome. The divergent evolution dynamics analysis also indicated that the majority of the PB transposon types identified in this study are either relatively young or old, with some being active. Additionally, numerous invasions of PB transposons were found in the genomes of both vertebrate and invertebrate animals. The data reveals that the PB superfamily is widely distributed in these species. PB transposons exhibit high diversity and activity in the small genomes of animals, and might play a crucial role in shaping the evolution of these small genomes of animals.

The genomes of Scedosporium between environmental challenges and opportunism

Emerging fungal pathogens are a global challenge for humankind. Many efforts have been made to understand the mechanisms underlying pathogenicity in bacteria, and OMICs techniques are largely responsible for those advancements. By contrast, our limited understanding of opportunism and antifungal resistance is preventing us from identifying, limiting and interpreting the emergence of fungal pathogens. The genus Scedosporium (Microascaceae) includes fungi with high tolerance to environmental pollution, whilst some species can be considered major human pathogens, such as Scedosporium apiospermum and Scedosporium boydii. However, unlike other fungal pathogens, little is known about the genome evolution of these organisms. We sequenced two novel genomes of Scedosporium aurantiacum and Scedosporium minutisporum isolated from extreme, strongly anthropized environments. We compared all the available Scedosporium and Microascaceae genomes, that we systematically annotated and characterized ex novo in most cases. The genomes in this family were integrated in a Phylum-level comparison to infer the presence of putative, shared genomic traits in filamentous ascomycetes with pathogenic potential. The analysis included the genomes of 100 environmental and clinical fungi, revealing poor evolutionary convergence of putative pathogenicity traits. By contrast, several features in Microascaceae and Scedosporium were detected that might have a dual role in responding to environmental challenges and allowing colonization of the human body, including chitin, melanin and other cell wall related genes, proteases, glutaredoxins and magnesium transporters. We found these gene families to be impacted by expansions, orthologous transposon insertions, and point mutations. With RNA-seq, we demonstrated that most of these anciently impacted genomic features responded to the stress imposed by an antifungal compound (voriconazole) in the two environmental strains S. aurantiacum MUT6114 and S. minutisporum MUT6113. Therefore, the present genomics and transcriptomics investigation stands on the edge between stress resistance and pathogenic potential, to elucidate whether fungi were pre-adapted to infect humans. We highlight the strengths and limitations of genomics applied to opportunistic human pathogens, the multifactoriality of pathogenicity and resistance to drugs, and suggest a scenario where pressures other than anthropic contributed to forge filamentous human pathogens.

Intraspecific Variation of Transposable Elements Reveals Differences in the Evolutionary History of Fungal Phytopathogen Pathotypes

Transposable elements (TEs) contribute to intraspecific variation and play important roles in the evolution of fungal genomes. However, our understanding of the processes that shape TE landscapes is limited, as is our understanding of the relationship between TE content, population structure, and evolutionary history of fungal species. Fungal plant pathogens, which often have host-specific populations, are useful systems in which to study intraspecific TE content diversity. Here, we describe TE dynamics in five lineages of Magnaporthe oryzae, the fungus that causes blast disease of rice, wheat, and many other grasses. We identified differences in TE content across these lineages and showed that recent lineage-specific expansions of certain TEs have contributed to overall greater TE content in rice-infecting and Setaria-infecting lineages. We reconstructed the evolutionary histories of long terminal repeat-retrotransposon expansions and found that in some cases they were caused by complex proliferation dynamics of one element and in others by multiple elements from an older population of TEs multiplying in parallel. Additionally, we found evidence suggesting the recent transfer of a DNA transposon between rice- and wheat-infecting M. oryzae lineages and a region showing evidence of homologous recombination between those lineages, which could have facilitated such a transfer. By investigating intraspecific TE content variation, we uncovered key differences in the proliferation dynamics of TEs in various pathotypes of a fungal plant pathogen, giving us a better understanding of the evolutionary history of the pathogen itself.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

Chromosome-scale assembly uncovers genomic compartmentation of Fusarium oxysporum f. sp. albedinis, the causal agent of Bayoud disease in date palm

Date palm (Phoenixdactylifera) is the most significant crop across North Africa and the Middle East. However, the crop faces a severe threat from Bayoud disease caused by the fungal pathogen Fusarium oxysporum f. sp. albedinis (FOA). FOA is a soil-borne fungus that infects the roots and vascular system of date palms, leading to widespread destruction of date palm plantations in North Africa over the last century. This is considered the most devastating pathogen of oasis agriculture in North Africa and responsible for loss of 13 million trees in Algeria and Morocco alone. In this study, we present a chromosome-scale high-quality genome assembly of the virulent isolate Foa 44, which provides valuable insights into understanding the genetic basis of Bayoud disease. The genome assembly consists of 11 chromosomes and 40 unplaced contigs, totalling 65,971,825 base pairs in size. It exhibits a GC ratio of 47.77% and a TE (transposable element) content of 17.30%. Through prediction and annotation, we identified 20,416 protein-coding genes. By combining gene and repeat densities analysis with alignment to Fusarium oxysporum f. sp. lycopersici (FOL) 4287 isolate genome sequence, we determined the core and lineage-specific compartments in Foa 44, shedding light on the genome structure of this pathogen. Furthermore, a phylogenomic analysis based on the 3,292 BUSCOs core genome revealed a distinct clade of FOA isolates within the Fusarium oxysporum species complex (FOSC). Notably, the genealogies of the five identified Secreted In Xylem (SIX) genes (1, 6, 9, 11 and 14) in FOA displayed a polyphyletic pattern, suggesting a horizontal inheritance of these effectors. These findings provide a valuable genomics toolbox for further research aimed at combatting the serious biotic constraints posed by FOA to date palm. This will pave the way for a deeper understanding of Bayoud disease and facilitate the development of effective diagnostic tools and control measures.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

A chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Genome characterization of two novel deep-sea sediment fungi, Penicillium pacificagyrus sp. nov. and Penicillium pacificasedimenti sp. nov., from South Pacific Gyre subseafloor sediments, highlights survivability

BACKGROUND

Marine deep subsurface sediments were once thought to be devoid of eukaryotic life, but advances in molecular technology have unlocked the presence and activity of well-known closely related terrestrial and marine fungi. Commonly detected fungi in deep marine sediment environments includes Penicillium, Aspergillus, Cladosporium, Fusarium, and Schizophyllum, which could have important implications in carbon and nitrogen cycling in this isolated environment. In order to determine the diversity and unknown metabolic capabilities of fungi in deep-sea sediments, their genomes need to be fully analyzed. In this study, two Penicillium species were isolated from South Pacific Gyre sediment enrichments during Integrated Ocean Drilling Program Expedition 329. The inner gyre has very limited productivity, organic carbon, and nutrients.

RESULTS

Here, we present high-quality genomes of two proposed novel Penicillium species using Illumina HiSeq and PacBio sequencing technologies. Single-copy homologues within the genomes were compared to other closely related genomes using OrthoMCL and maximum-likelihood estimation, which showed that these genomes were novel species within the genus Penicillium. We propose to name isolate SPG-F1 as Penicillium pacificasedimenti sp. nov. and SPG-F15 as Penicillium pacificagyrus sp. nov. The resulting genome sizes were 32.6 Mbp and 36.4 Mbp, respectively, and both genomes were greater than 98% complete as determined by the presence of complete single-copy orthologs. The transposable elements for each genome were 4.87% for P. pacificasedimenti and 10.68% for P. pacificagyrus. A total of 12,271 genes were predicted in the P. pacificasedimenti genome and 12,568 genes in P. pacificagyrus. Both isolates contained genes known to be involved in the degradation of recalcitrant carbon, amino acids, and lignin-derived carbon.

CONCLUSIONS

Our results provide the first constructed genomes of novel Penicillium isolates from deep marine sediments, which will be useful for future studies of marine subsurface fungal diversity and function. Furthermore, these genomes shed light on the potential impact fungi in marine sediments and the subseafloor could have on global carbon and nitrogen biogeochemical cycles and how they may be persisting in the most energy-limited sedimentary biosphere.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Haplotype-phased and chromosome-level genome assembly of Puccinia polysora, a giga-scale fungal pathogen causing southern corn rust

Rust fungi are characterized by large genomes with high repeat content and have two haploid nuclei in most life stages, which makes achieving high-quality genome assemblies challenging. Here, we described a pipeline using HiFi reads and Hi-C data to assemble a gigabase-sized fungal pathogen, Puccinia polysora f.sp. zeae, to haplotype-phased and chromosome-scale. The final assembled genome is 1.71 Gbp, with ~850 Mbp and 18 chromosomes in each haplotype, being currently one of the two giga-scale fungi assembled to chromosome level. Transcript-based annotation identified 47,512 genes for the dikaryotic genome with a similar number for each haplotype. A high level of interhaplotype variation was found with 10% haplotype-specific BUSCO genes, 5.8 SNPs/kbp, and structural variation accounting for 3% of the genome size. The P. polysora genome displayed over 85% repeat contents, with genome-size expansion and copy number increasing of species-specific orthogroups. Interestingly, these features did not affect overall synteny with other Puccinia species having smaller genomes. Fine-time-point transcriptomics revealed seven clusters of coexpressed secreted proteins that are conserved between two haplotypes. The fact that candidate effectors interspersed with all genes indicated the absence of a "two-speed genome" evolution in P. polysora. Genome resequencing of 79 additional isolates revealed a clonal population structure of P. polysora in China with low geographic differentiation. Nevertheless, a minor population differentiated from the major population by having mutations on secreted proteins including AvrRppC, indicating the ongoing virulence to evade recognition by RppC, a major resistance gene in Chinese corn cultivars. The high-quality assembly provides valuable genomic resources for future studies on disease management and the evolution of P. polysora.

Transposable elements in the transcriptome of the velvetbean caterpillar Anticarsia gemmatalis Hübner, 1818 (Lepidoptera: Erebidae)

Transposable elements (TEs) are DNA sequences that possess the ability to move from one genomic location to another. These sequences contribute to a significant fraction of the genomes of most eukaryotes and can impact their architecture and regulation. In this paper, we present the first data related to the identification and characterization of TEs present in the transcriptome of Anticarsia gemmatalis. Approximately, 835 transcripts showed significant similarity to TEs and (or) characteristic domains. Retrotransposons accounted for 71.2% (595 sequences) of the identified elements, while DNA transposons were less abundant, with 240 annotations (28.8%). TEs were classified into 30 superfamilies, with SINE3/5S and Gypsy being the most abundant. Based on the sequences of TEs found in the transcriptome, we were able to locate conserved regions in the chromosomes of this species. The analysis of differential expression of TEs in susceptible and resistant strains, challenged and not challenged with Bacillus thuringiensis ( Bt) from in silico analysis, indicated that exposure to Bt can regulate the transcription of mobile genetic elements in the velvetbean caterpillar. Thus, these data contribute significantly to the knowledge of the structure and composition of these elements in the genome of this species, and suggest the role of stress on their expression.

Insights into the genomic evolution and the alkali tolerance mechanisms of Agaricus sinodeliciosus by comparative genomic and transcriptomic analyses

Agaricus sinodeliciosus is a rare wild edible mushroom from northwest China, and grows naturally in mild saline-alkali soil, which is also unusual in mushrooms. A. sinodeliciosus represents a potential model organism for explaining saline-alkali tolerance mechanisms and revealing related physiological processes in mushrooms. Here, we provide a high-quality genome of A. sinodeliciosus. Comparative genomic analyses reveal A. sinodeliciosus has numerous changes to its genome organization after a solitary evolutionary history under saline-alkali environments, such as gene family contraction, retrotransposon expansion and rapid evolution of adaptative genes. Our saline and alkali tolerance tests show that mycelium growth and fruit body formation of this species are effected by mild alkalinity. Transcriptomic analyses reveal that genes involved in carbon and nitrogen utilization, cell stability and fruit body formation of A. sinodeliciosus could be activated under mildly alkaline conditions. In particular, the 'starch and sucrose metabolism', 'biosynthesis of amino acids' and 'phenylpropanoid biosynthesis' pathways are important for mildly alkaline tolerance of A. sinodeliciosus. Like plants and arbuscular mycorrhizal fungi, in the rot fungus A. sinodeliciosus, the biosynthesis of intracellular small molecules could be enhanced to counter osmotic and oxidative stresses caused by mild alkalinity, and the biosynthesis of monolignol could be suppressed to increase cell wall infiltrates under mildly alkaline conditions. This research provides an understanding of the genomic evolution and mechanisms of A. sinodeliciosus in tolerance to saline-alkali environments. The A. sinodeliciosus genome constitutes a valuable resource for evolutionary and ecological studies of Agaricus.

Biocatalyst collection and heterologous expression of sesquiterpene synthases from basidiomycetous fungi: Discovery of a novel sesquiterpene hydrocarbon

Basidiomycetes produce a wide variety of sesquiterpenoids, which attract significant interest in pharmaceutical and industrial applications. Structural diversification of sesquiterpenoids is performed by sesquiterpene synthases (STSs), which produce a wide array of backbone structures; therefore, functional characterization and increased biocatalyst collection of STSs are important for expanding scientific knowledge and meeting the needs of advanced biotechnology. Gene identification and functional annotation of STSs from the basidiomycetous fungi Agaricus bisporus, Auriscalpium vulgare, Lepista nuda, Pleurotus ostreatus and Trametes versicolor were conducted. Through these investigations, the catalytic functions of 30 STSs were revealed using recombinant enzymes heterologously expressed in Saccharomyces cerevisiae. Furthermore, the unique function of an STS from P. ostreatus, PoSTS-06, was revealed to be the production of a novel sesquiterpene hydrocarbon that we named pleostene. The absolute structure of pleostene was determined by NMR spectroscopy and X-ray crystallography using the crystalline sponge method.

Transcriptome Analysis of Antrodia cinnamomea Mycelia from Different Wood Substrates

Antrodia cinnamomea, an edible and medicinal fungus with significant economic value and application prospects, is rich in terpenoids, benzenoids, lignans, polysaccharides, and benzoquinone, succinic and maleic derivatives. In this study, the transcriptome of A. cinnamomea cultured on the wood substrates of Cinnamomum glanduliferum (YZM), C. camphora (XZM), and C. kanehirae (NZM) was sequenced using the high-throughput sequencing technology Illumina HiSeq 2000, and the data were assembled by de novo strategy to obtain 78,729 Unigenes with an N50 of 4,463 bp. Compared with public databases, about 11,435, 6,947, and 5,994 Unigenes were annotated to the Non-Redundant (NR), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genome (KEGG), respectively. The comprehensive analysis of the mycelium terpene biosynthesis-related genes in A. cinnamomea revealed that the expression of acetyl-CoA acetyltransferase (AACT), acyl-CoA dehydrogenase (MCAD), 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA), mevalonate pyrophosphate decarboxylase (MVD), and isopentenyl diphosphate isomerase (IDI) was significantly higher on NZM compared to the other two wood substrates. Similarly, the expression of geranylgeranyltransferase (GGT) was significantly higher on YZM compared to NZM and XZM, and the expression of farnesyl transferase (FTase) was significantly higher on XZM. Furthermore, the expressions of 2,3-oxidized squalene cyclase (OCS), squalene synthase (SQS), and squalene epoxidase (SE) were significantly higher on NZM. Overall, this study provides a potential approach to explore the molecular regulation mechanism of terpenoid biosynthesis in A. cinnamomea.

Gapless Genome Assembly of Puccinia triticina Provides Insights into Chromosome Evolution in Pucciniales

Chromosome evolution drives species evolution, speciation, and adaptive radiation. Accurate genome assembly is crucial to understanding chromosome evolution of species, such as dikaryotic fungi. Rust fungi (Pucciniales) in dikaryons represent the largest group of plant pathogens, but the evolutionary process of adaptive radiation in Pucciniales remains poorly understood. Here, we report a gapless genome for the wheat leaf rust fungus Puccinia triticina determined using PacBio high-fidelity (HiFi) sequencing. This gapless assembly contains two sets of chromosomes, showing that one contig represents one chromosome. Comparisons of homologous chromosomes between the phased haplotypes revealed that highly frequent small-scale sequence divergence shapes haplotypic variation. Genome analyses of Puccinia triticina along with other rusts revealed that recent transposable element bursts and extensive segmental gene duplications synergistically highlight the evolution of chromosome structures. Comparative analysis of chromosomes indicated that frequent chromosomal rearrangements may act as a major contributor to rapid radiation of Pucciniales. This study presents the first gapless, phased assembly for a dikaryotic rust fungus and provides insights into adaptive evolution and species radiation in Pucciniales. IMPORTANCE Rust fungi (Pucciniales) are the largest group of plant pathogens. Adaptive radiation is a predominant feature in Pucciniales evolution. Chromosome evolution plays an important role in adaptive evolution. Accurate chromosome-scale assembly is required to understand the role of chromosome evolution in Pucciniales. We took advantage of HiFi sequencing to construct a gapless, phased genome for Puccinia triticina. Further analyses revealed that the evolution of chromosome structures in rust lineage is shaped by the combination of transposable element bursts and segmental gene duplications. Chromosome comparisons of Puccinia triticina and other rusts suggested that frequent chromosomal arrangements may make remarkable contributions to high species diversity of rust fungi. Our results present the first gapless genome for Pucciniales and shed light on the feature of chromosome evolution in Pucciniales.

Heterochromatin and RNAi act independently to ensure genome stability in Mucorales human fungal pathogens

Chromatin modifications play a fundamental role in controlling transcription and genome stability and yet despite their importance, are poorly understood in early-diverging fungi. We present a comprehensive study of histone lysine and DNA methyltransferases across the Mucoromycota, emphasizing heterochromatin formation pathways that rely on the Clr4 complex involved in H3K9-methylation, the Polycomb-repressive complex 2 driving H3K27-methylation, or DNMT1-like methyltransferases that catalyze 5mC DNA methylation. Our analysis uncovered H3K9-methylated heterochromatin as the major chromatin modification repressing transcription in these fungi, which lack both Polycomb silencing and cytosine methylation. Although small RNAs generated by RNA interference (RNAi) pathways facilitate the formation of heterochromatin in many eukaryotic organisms, we show that RNAi is not required to maintain either genomic or centromeric heterochromatin in Mucor. H3K9-methylation and RNAi act independently to control centromeric regions, suggesting a functional subspecialization. Whereas the H3K9 methyltransferase Clr4 and heterochromatin formation are essential for cell viability, RNAi is dispensable for viability yet acts as the main epigenetic, regulatory force repressing transposition of centromeric GremLINE1 elements. Mutations inactivating canonical RNAi lead to rampant transposition and insertional inactivation of targets resulting in antimicrobial drug resistance. This fine-tuned, Rdrp2-dependent RNAi activity is critical for genome stability, restricting GremLINE1 retroelements to the centromeres where they occupy long heterochromatic islands. Taken together, our results suggest that RNAi and heterochromatin formation are independent genome defense and regulatory mechanisms in the Mucorales, contributing to a paradigm shift from the cotranscriptional gene silencing observed in fission yeasts to models in which heterochromatin and RNAi operate independently in early-diverging fungi.

Purification and Characterization of Class III Lipase from a White-Rot Fungus Pleurotus ostreatus

Pleurotus ostreatus is an edible white-rot fungus with lignocellulosic biomass degrading enzymes that have been studied extensively. However, until now, lipolytic enzymes from P. ostreatus, which degrade extractives in lignocellulosic biomass, have not been purified and characterized. In this study, P. ostreatus was inoculated into the rapeseed oil containing culture to induce lipase. The lipase in the culture broth was successfully purified to homogeneity by chromatographic methods. The molecular weight of the purified lipase was 27 kDa, and its optimal pH and temperature were 5.0 and 30 °C, respectively. The purified lipase showed high activity with the substrates 4-methylumbelliferyl (4-MU) decanoate (C10:0) and 4-MU oleate (C18:1), and no activity with 4-MU acetate (C2:0) and 4-MU butyrate (C4:0). The amino acid sequences and substrate specificities of the purified lipase suggested that it belonged to class III. Kinetic parameters measurements (Km and Vmax) showed that 4-MU palmitate had a high affinity for the purified lipase, and it was the substrate most efficiently hydrolyzed by the purified lipase.

Two-speed genome evolution drives pathogenicity in fungal pathogens of animals

The origins and evolution of virulence in amphibian-infecting chytrids Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are largely unknown. Here, we use deep nanopore sequencing of Bsal and comparative genomics against 21 high-quality genome assemblies that span the fungal Chytridiomycota. We discover that Bsal has the most repeat-rich genome of the Chytridiomycota, comprising 40.9% repetitive elements; this genome has expanded to more than 3× the length of its conspecific Bd, with autonomous and fully functional LTR/Gypsy elements contributing significantly to the expansion. The M36 metalloprotease virulence factors are highly expanded (n = 177) in Bsal, most of which (53%) are flanked by transposable elements, suggesting they have a repeat-associated expansion. We find enrichment upstream of M36 metalloprotease genes of three novel repeat families belonging to the repeat superfamily of LINEs that are implicated with gene copy number variations. Additionally, Bsal has a highly compartmentalized genome architecture, with virulence factors enriched in gene-sparse/repeat-rich compartments, while core conserved genes are enriched in gene-rich/repeat-poor compartments. Genes upregulated during infection are primarily found in the gene-sparse/repeat-rich compartment in both Bd and Bsal. Furthermore, genes with signatures of positive selection in Bd are enriched in repeat-rich regions, suggesting these regions are a cradle for the evolution of chytrid pathogenicity. These are the hallmarks of two-speed genome evolution, and this study provides evidence of two-speed genomes in an animal pathogen, shedding light on the evolution of fungal pathogens of vertebrates driving global declines and extinctions.

Virulence and pathogenicity determinants in whole genome sequence of Fusarium udum causing wilt of pigeon pea

The present study deals with whole genome analysis of Fusarium udum, a wilt causing pathogen of pigeon pea. The de novo assembly identified a total of 16,179 protein-coding genes, of which 11,892 genes (73.50%) were annotated using BlastP and 8,928 genes (55.18%) from KOG annotation. In addition, 5,134 unique InterPro domains were detected in the annotated genes. Apart from this, we also analyzed genome sequence for key pathogenic genes involved in virulence, and identified 1,060 genes (6.55%) as virulence genes as per the PHI-BASE database. The secretome profiling of these virulence genes indicated the presence of 1,439 secretory proteins. Of those, an annotation of 506 predicted secretory proteins through CAZyme database indicated maximum abundance of Glycosyl hydrolase (GH, 45%) family proteins followed by auxiliary activity (AA) family proteins. Interestingly, the presence of effectors for cell wall degradation, pectin degradation, and host cell death was found. The genome comprised approximately 895,132 bp of repetitive elements, which includes 128 long terminal repeats (LTRs), and 4,921 simple sequence repeats (SSRs) of 80,875 bp length. The comparative mining of effector genes among different Fusarium species revealed five common and two specific effectors in F. udum that are related to host cell death. Furthermore, wet lab experiment validated the presence of effector genes like SIX (for Secreted in Xylem). We conclude that deciphering the whole genome of F. udum would be instrumental in understanding evolution, virulence determinants, host-pathogen interaction, possible control strategies, ecological behavior, and many other complexities of the pathogen.

Draft genome sequencing and secretome profiling of Sclerotinia sclerotiorum revealed effector repertoire diversity and allied broad-host range necrotrophy

White mold commonly known as Sclerotinia sclerotiorum causes stem rot disease and has emerged as one of the major fungal pathogens of oilseed Brassica across the world. In the present study, consistently virulent S. sclerotiorum isolate "ESR-01" was sequenced and an assembly size of ~ 41 Mb with 328 scaffolds having N50 of 447,128 was obtained. Additionally, 27,450 single nucleotide polymorphisms (SNPs) were identified from 155 scaffolds against S. sclerotiorum 1980 isolate, with an average SNP density of ~ 1.5 per kb genome. 667 repetitive elements were identified and approximately comprised 7% of the total annotated genes. The DDE_1 with 454 in numbers was found to be the most abundant and accounts for 68% of the total predicted repetitive elements. In total, 3844 simple sequence repeats are identified in the 328 scaffolds. A total of 9469 protein-coding genes were predicted from the whole genome assembly with an average gene length of 1587 bp and their distribution as 230.95 genes per Mb in the genome. Out of 9469 predicted protein-coding genes, 529 genes were observed encoding the CAZymes (Carbohydrate-Active enzymes) capable of degradation of the complex polysaccharides. Glycosyltransferase (GT) families were most abundant (49.71%) among the predicted CAZymes and GT2 (23%), GT4 (20%), and glycoside hydrolase (GH) 23% with GH18 (11%) were the prominent cell wall degrading enzyme families in the ESR-01 secretome. Besides this, 156 genes essential for the pathogen-host interactions were also identified. The effector analysis in the whole genome proteomics dataset revealed a total of 57 effector candidates (ECs) and 27 of them were having their analogs whereas the remaining 30 were novel ones. Eleven selected ECs were validated experimentally by analyzing the expression profile of the ESR-01 isolate of S. sclerotiorum. Together, the present investigation offers a better understanding of the S. sclerotiorum genome, secretome, and its effector repertoire which will help in refining the present knowledge on S. sclerotiorum-Brassica interactions and necrotrophic lifestyle of the phytopathogen in general.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

Genome-wide comparison deciphers lifestyle adaptation and glass biodeterioration property of Curvularia eragrostidis C52

Glass biodeterioration by fungi has caused irreversible damage to valuable glass materials such as cultural heritages and optical devices. To date, knowledge about metabolic potential and genomic profile of biodeteriorative fungi is still scarce. Here, we report for the first time the whole genome sequence of Curvularia eragrostidis C52 that strongly degraded silica-based glasses coated with fluorine and hafnium, as expressed by the hyphal surface coverage of 46.16 ± 3.3% and reduced light transmission of 50.93 ± 1.45%. The genome of C. eragrostidis C52 is 36.9 Mb long with a GC content of 52.1% and contains 14,913 protein-coding genes, which is the largest genome ever recorded in the genus Curvularia. Phylogenomic analysis revealed C. eragrostidis C52 formed a distinct cluster with Curvularia sp. IFB-Z10 and was not evolved from compared genomes. Genome-wide comparison showed that strain C52 harbored significantly higher proportion of proteins involved in carbohydrate-active enzymes, peptidases, secreted proteins, and transcriptional factors, which may be potentially attributed to a lifestyle adaptation. Furthermore, 72 genes involved in the biosynthesis of 6 different organic acids were identified and expected to be crucial for the fungal survival in the glass environment. To form biofilm against stress, the fungal strain utilized 32 genes responsible for exopolysaccharide production. These findings will foster a better understanding of the biology of C. eragrostidis and the mechanisms behind fungal biodeterioration in the future.

A Bacterial Type Three Secretion-Based Delivery System for Functional Characterization of Sporisorium scitamineum Plant Immune Suppressing Effector Proteins

The facultative biotrophic basidiomycete Sporisorium scitamineum causes smut disease in sugarcane. This study applied an assay to identify S. scitamineum candidate effectors (CEs) with plant immunity suppression activities by delivering them into Nicotiana benthamiana cells via the type-three secretion system of Pseudomonas fluorescens EtHAn. Six CEs were individually cloned into the pEDV6 vector and expressed by P. fluorescens EtHAn for translocation into the plant cells. Three CEs (g1052, g3890, and g5159) could suppress pattern-triggered immunity (PTI) responses with high reproducibility in different coinfiltration experiments with P. syringae pv. tomato DC3000. In addition, three CEs (g1052, g4549, and g5159) were also found to be AvrB-induced suppressors of effector-triggered immunity (ETI), demonstrating for the first time that S. scitamineum can defeat both PTI and ETI responses. A transcriptomic analysis at different stages of infection by the smut fungus of three sugarcane cultivars with contrasting responses to the pathogen revealed that suppressors g1052, g3890, g4549, and g5159 were induced at the early stage of infection. By contrast, the two CEs (g2666 and g6610) that did not exhibit suppression activities expressed only at the late stage of infection. Moreover, genomic structures of the CEs and searches for orthologs in other smut species suggested duplication events and further divergence in CEs evolution of S. scitamineum. Thus, the transient assay applied here demonstrated the potential of pEDV6 and P. fluorescens EtHAn as biological tools for identifying plant immune suppressors from S. scitamineum.

Exploring a novel genomic safe-haven site in the human pathogenic mould Aspergillus fumigatus

Aspergillus fumigatus is the most important airborne fungal pathogen and allergen of humans causing high morbidity and mortality worldwide. The factors that govern pathogenicity of this organism are multi-factorial and are poorly understood. Molecular tools to dissect the mechanisms of pathogenicity in A. fumigatus have improved significantly over the last 20 years however many procedures have not been standardised for A. fumigatus. Here, we present a new genomic safe-haven locus at the site of an inactivated transposon, named SH-aft4, which can be used to insert DNA sequences in the genome of this fungus without impacting its phenotype. We show that we are able to effectively express a transgene construct from the SH-aft4 and that natural regulation of promoter function is conserved at this site. Furthermore, the SH-aft4 locus is highly conserved in the genome of a wide range of clinical and environmental isolates including the isolates commonly used by many laboratories CEA10, Af293 and ATCC46645, allowing a wide range of isolates to be manipulated. Our results show that the aft4 locus can serve as a site for integration of a wide range of genetic constructs to aid functional genomics studies of this important human fungal pathogen.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

A chromosomal-level reference genome of the widely utilized Coccidioides posadasii laboratory strain "Silveira"

Coccidioidomycosis is a common fungal disease that is endemic to arid and semi-arid regions of both American continents. Coccidioides immitis and Coccidioides posadasii are the etiological agents of the disease, also known as Valley Fever. For several decades, the C. posadasii strain Silveira has been used widely in vaccine studies, is the source strain for production of diagnostic antigens, and is a widely used experimental strain for functional studies. In 2009, the genome was sequenced using Sanger sequencing technology, and a draft assembly and annotation were made available. In this study, the genome of the Silveira strain was sequenced using single molecule real-time sequencing PacBio technology, assembled into chromosomal-level contigs, genotyped, and the genome was reannotated using sophisticated and curated in silico tools. This high-quality genome sequencing effort has improved our understanding of chromosomal structure, gene set annotation, and lays the groundwork for identification of structural variants (e.g. transversions, translocations, and copy number variants), assessment of gene gain and loss, and comparison of transposable elements in future phylogenetic and population genomics studies.

A chromosome-scale genome assembly of the tomato pathogen Cladosporium fulvum reveals a compartmentalized genome architecture and the presence of a dispensable chromosome

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the 'two-speed genome' model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

Distribution and regulatory roles of oxidized 5-methylcytosines in DNA and RNA of the basidiomycete fungi Laccaria bicolor and Coprinopsis cinerea

The formation of three oxidative DNA 5-methylcytosine (5mC) modifications (oxi-mCs)-5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC)-by the TET/JBP family of dioxygenases prompted intensive studies of their functional roles in mammalian cells. However, the functional interplay of these less abundant modified nucleotides in other eukaryotic lineages remains poorly understood. We carried out a systematic study of the content and distribution of oxi-mCs in the DNA and RNA of the basidiomycetes Laccaria bicolor and Coprinopsis cinerea, which are established models to study DNA methylation and developmental and symbiotic processes. Quantitative liquid chromatography-tandem mass spectrometry revealed persistent but uneven occurrences of 5hmC, 5fC and 5caC in the DNA and RNA of the two organisms, which could be upregulated by vitamin C. 5caC in RNA (5carC) was predominantly found in non-ribosomal RNA, which potentially includes non-coding, messenger and small RNA species. Genome-wide mapping of 5hmC and 5fC using the single CG analysis techniques hmTOP-seq and foTOP-seq pointed at involvement of oxi-mCs in the regulation of gene expression and silencing of transposable elements. The implicated diverse roles of 5mC and oxi-mCs in the two fungi highlight the epigenetic importance of the latter modifications, which are often neglected in standard whole-genome bisulfite analyses.

Genome Analysis of the Broad Host Range Necrotroph Nalanthamala psidii Highlights Genes Associated With Virulence

Guava wilt disease is caused by the fungus Nalanthamala psidii. The wilt disease results in large-scale destruction of orchards in South Africa, Taiwan, and several Southeast Asian countries. De novo assembly, annotation, and in-depth analysis of the N. psidii genome were carried out to facilitate the identification of characteristics associated with pathogenicity and pathogen evolution. The predicted secretome revealed a range of CAZymes, proteases, lipases and peroxidases associated with plant cell wall degradation, nutrient acquisition, and disease development. Further analysis of the N. psidii carbohydrate-active enzyme profile exposed the broad-spectrum necrotrophic lifestyle of the pathogen, which was corroborated by the identification of putative effectors and secondary metabolites with the potential to induce tissue necrosis and cell surface-dependent immune responses. Putative regulatory proteins including transcription factors and kinases were identified in addition to transporters potentially involved in the secretion of secondary metabolites. Transporters identified included important ABC and MFS transporters involved in the efflux of fungicides. Analysis of the repetitive landscape and the detection of mechanisms linked to reproduction such as het and mating genes rendered insights into the biological complexity and evolutionary potential of N. psidii as guava pathogen. Hence, the assembly and annotation of the N. psidii genome provided a valuable platform to explore the pathogenic potential and necrotrophic lifestyle of the guava wilt pathogen.

Sexual Crossing, Chromosome-Level Genome Sequences, and Comparative Genomic Analyses for the Medicinal Mushroom Taiwanofungus Camphoratus (Syn. Antrodia Cinnamomea, Antrodia Camphorata)

Taiwanofungus camphoratus mushrooms are a complementary and alternative medicine for hangovers, cancer, hypertension, obesity, diabetes, and inflammation. Though Taiwanofungus camphoratus has attracted considerable biotechnological and pharmacological attention, neither classical genetic nor genomic approaches have been properly established for it. We isolated four sexually competent monokaryons from two T. camphoratus dikaryons used for the commercial cultivation of orange-red (HC1) and milky-white (SN1) mushrooms, respectively. We also sequenced, annotated, and comparatively analyzed high-quality and chromosome-level genome sequences of these four monokaryons. These genomic resources represent a valuable basis for understanding the biology, evolution, and secondary metabolite biosynthesis of this economically important mushrooms. We demonstrate that T. camphoratus has a tetrapolar mating system and that HC1 and SN1 represent two intraspecies isolates displaying karyotypic variation. Compared with several edible mushroom model organisms, T. camphoratus underwent a significant contraction in the gene family and individual gene numbers, most notably for plant, fungal, and bacterial cell-wall-degrading enzymes, explaining why T. camphoratus mushrooms are rare in natural environments, are difficult and time-consuming to artificially cultivate, and are susceptible to fungal and bacterial infections. Our results lay the foundation for an in-depth T. camphoratus study, including precise genetic manipulation, improvements to mushroom fruiting, and synthetic biology applications for producing natural medicinal products. IMPORTANCETaiwanofungus camphoratus (Tc) is a basidiomycete fungus that causes brown heart rot of the aromatic tree Cinnamomum kanehirae. The Tc fruiting bodies have been used to treat hangovers, abdominal pain, diarrhea, hypertension, and other diseases first by aboriginal Taiwanese and later by people in many countries. To establish classical genetic and genomic approaches for this economically important medicinal mushroom, we first isolated and characterized four sexually competent monokaryons from two dikaryons wildly used for commercial production of Tc mushrooms. We applied PacBio single molecule, real-time sequencing technology to determine the near-completed genome sequences of four monokaryons. These telomere-to-telomere and gapless haploid genome sequences reveal all genomic variants needed to be studied and discovered, including centromeres, telomeres, retrotransposons, mating type loci, biosynthetic, and metabolic gene clusters. Substantial interspecies diversities are also discovered between Tc and several other mushroom model organisms, including Agrocybe aegerita, Coprinopsis cinerea, and Schizophyllum commune, and Ganoderma lucidum.

The role of DNA methylation in genome defense in Cnidaria and other invertebrates

Considerable attention has recently been focused on the potential involvement of DNA methylation in regulating gene expression in cnidarians. Much of this work has been centered on corals, in the context of changes in methylation perhaps facilitating adaptation to higher seawater temperatures and other stressful conditions. Although first proposed more than 30 years ago, the possibility that DNA methylation systems function in protecting animal genomes against the harmful effects of transposon activity has largely been ignored since that time. Here, we show that transposons are specifically targeted by the DNA methylation system in cnidarians, and that the youngest transposons (i.e., those most likely to be active) are most highly methylated. Transposons in longer and highly active genes were preferentially methylated and, as transposons aged, methylation levels declined, reducing the potentially harmful side effects of CpG methylation. In Cnidaria and a range of other invertebrates, correlation between the overall extent of methylation and transposon content was strongly supported. Present transposon burden is the dominant factor in determining overall level of genomic methylation in a range of animals that diverged in or before the early Cambrian, suggesting that genome defense represents the ancestral role of CpG methylation.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Fungi as mediators linking organisms and ecosystems

Fungi form a major and diverse component of most ecosystems on Earth. They are both micro and macroorganisms with high and varying functional diversity as well as great variation in dispersal modes. With our growing knowledge of microbial biogeography, it has become increasingly clear that fungal assembly patterns and processes differ from other microorganisms such as bacteria, but also from macroorganisms such as plants. The success of fungi as organisms and their influence on the environment lies in their ability to span multiple dimensions of time, space, and biological interactions, that is not rivalled by other organism groups. There is also growing evidence that fungi mediate links between different organisms and ecosystems, with the potential to affect the macroecology and evolution of those organisms. This suggests that fungal interactions are an ecological driving force, interconnecting different levels of biological and ecological organisation of their hosts, competitors, and antagonists with the environment and ecosystem functioning. Here we review these emerging lines of evidence by focusing on the dynamics of fungal interactions with other organism groups across various ecosystems. We conclude that the mediating role of fungi through their complex and dynamic ecological interactions underlie their importance and ubiquity across Earth's ecosystems.