Transcriptional activity and epigenetic regulation of transposable elements in the symbiotic fungus Rhizophagus irregularis

Strain-specific evolution and host-specific regulation of transposable elements in the model plant symbiont Rhizophagus irregularis

Transposable elements (TEs) are repetitive DNA that can create genome structure and regulation variability. The genome of Rhizophagus irregularis, a widely studied arbuscular mycorrhizal fungus (AMF), comprises ∼50% repetitive sequences that include TEs. Despite their abundance, two-thirds of TEs remain unclassified, and their regulation among AMF life stages remains unknown. Here, we aimed to improve our understanding of TE diversity and regulation in this model species by curating repeat datasets obtained from chromosome-level assemblies and by investigating their expression across multiple conditions. Our analyses uncovered new TE superfamilies and families in this model symbiont and revealed significant differences in how these sequences evolve both within and between R. irregularis strains. With this curated TE annotation, we also found that the number of upregulated TE families in colonized roots is 4 times higher than in the extraradical mycelium, and their overall expression differs depending on the plant host. This work provides a fine-scale view of TE diversity and evolution in model plant symbionts and highlights their transcriptional dynamism and specificity during host-microbe interactions. We also provide Hidden Markov Model profiles of TE domains for future manual curation of uncharacterized sequences (https://github.com/jordana-olive/TE-manual-curation/tree/main).

A journey into the world of small RNAs in the arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal (AM) symbiosis is a mutualistic interaction between fungi and most land plants that is underpinned by a bidirectional exchange of nutrients. AM development is a tightly regulated process that encompasses molecular communication for reciprocal recognition, fungal accommodation in root tissues and activation of symbiotic function. As such, a complex network of transcriptional regulation and molecular signaling underlies the cellular and metabolic reprogramming of host cells upon AM fungal colonization. In addition to transcription factors, small RNAs (sRNAs) are emerging as important regulators embedded in the gene network that orchestrates AM development. In addition to controlling cell-autonomous processes, plant sRNAs also function as mobile signals capable of moving to different organs and even to different plants or organisms that interact with plants. AM fungi also produce sRNAs; however, their function in the AM symbiosis remains largely unknown. Here, we discuss the contribution of host sRNAs in the development of AM symbiosis by considering their role in the transcriptional reprogramming of AM fungal colonized cells. We also describe the characteristics of AM fungal-derived sRNAs and emerging evidence for the bidirectional transfer of functional sRNAs between the two partners to mutually modulate gene expression and control the symbiosis.

Epigenetic regulation influenced by soil microbiota and nutrients: Paving road to epigenome editing in plants

Soil is a complex ecosystem that houses microbes and nutrients that are necessary for plant development. Edaphic properties of the soil and environmental conditions influence microbial growth and nutrient accessibility. Various environmental stimuli largely affect the soil microbes and ionic balance, in turn influencing plants. Soil microflora helps decompose organic matter and is involved in mineral uptake. The combination of soil microbes and mineral nutrients notably affects plant growth. Recent advancements have enabled a deeper understanding of plant genetic/molecular regulators. Deficiencies/sufficiencies of soil minerals and microbes also alter plant gene regulation. Gene regulation mediated by epigenetic mechanisms comprises conformational alterations in chromatin structure, DNA/histone modifications, or involvement of small RNAs. Epigenetic regulation is unique due to its potential to inherit without involving alteration of the DNA sequence. Thus, the compilation study of heritable epigenetic changes driven by nutrient imbalances and soil microbes would facilitate understanding this molecular phenomenon in plants. This information can aid in epigenome editing, which has recently emerged as a promising technology for plant non-transgenic/non-mutagenic modification. Potential epigenetic marks induced by biotic and abiotic stresses in plants could be explored as target sites for epigenome editing. This review discusses novel ways of epigenome editing to create epigenome edited plants with desirable and heritable phenotypes. As plants are sessile and in constant exposure to the soil microbiome and nutrients, epigenetic changes induced by these factors could provide more effective, stable and a sustainable molecular solution for crop improvement.

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

Exploring the potential of endophyte-plant interactions for improving crop sustainable yields in a changing climate

Climate change poses a major threat to global food security, significantly reducing crop yields as cause of abiotic stresses, and for boosting the spread of new and old pathogens and pests. Sustainable crop management as a route to mitigation poses the challenge of recruiting an array of solutions and tools for the new aims. Among these, the deployment of positive interactions between the micro-biotic components of agroecosystems and plants can play a highly significant role, as part of the agro-ecological revolution. Endophytic microorganisms have emerged as a promising solution to tackle this challenge. Among these, Arbuscular Mycorrhizal Fungi (AMF) and endophytic bacteria and fungi have demonstrated their potential to alleviate abiotic stresses such as drought and heat stress, as well as the impacts of biotic stresses. They can enhance crop yields in a sustainable way also by other mechanisms, such as improving the nutrient uptake, or by direct effects on plant physiology. In this review we summarize and update on the main types of endophytes, we highlight several studies that demonstrate their efficacy in improving sustainable yields and explore possible avenues for implementing crop-microbiota interactions. The mechanisms underlying these interactions are highly complex and require a comprehensive understanding. For this reason, omic technologies such as genomics, transcriptomics, proteomics, and metabolomics have been employed to unravel, by a higher level of information, the complex network of interactions between plants and microorganisms. Therefore, we also discuss the various omic approaches and techniques that have been used so far to study plant-endophyte interactions.

Fungi of the order Mucorales express a "sealing-only" tRNA ligase

Some eukaryotic pre-tRNAs contain an intron that is removed by a dedicated set of enzymes. Intron-containing pre-tRNAs are cleaved by tRNA splicing endonuclease, followed by ligation of the two exons and release of the intron. Fungi use a "heal and seal" pathway that requires three distinct catalytic domains of the tRNA ligase enzyme, Trl1. In contrast, humans use a "direct ligation" pathway carried out by RTCB, an enzyme completely unrelated to Trl1. Because of these mechanistic differences, Trl1 has been proposed as a promising drug target for fungal infections. To validate Trl1 as a broad-spectrum drug target, we show that fungi from three different phyla contain Trl1 orthologs with all three domains. This includes the major invasive human fungal pathogens, and these proteins can each functionally replace yeast Trl1. In contrast, species from the order Mucorales, including the pathogens Rhizopus arrhizus and Mucor circinelloides, have an atypical Trl1 that contains the sealing domain but lacks both healing domains. Although these species contain fewer tRNA introns than other pathogenic fungi, they still require splicing to decode three of the 61 sense codons. These sealing-only Trl1 orthologs can functionally complement defects in the corresponding domain of yeast Trl1 and use a conserved catalytic lysine residue. We conclude that Mucorales use a sealing-only enzyme together with unidentified nonorthologous healing enzymes for their heal and seal pathway. This implies that drugs that target the sealing activity are more likely to be broader-spectrum antifungals than drugs that target the healing domains.

Fungi of the order Mucorales express a "sealing-only" tRNA ligase

Some eukaryotic pre-tRNAs contain an intron that is removed by a dedicated set of enzymes. Intron-containing pre-tRNAs are cleaved by tRNA splicing endonuclease (TSEN), followed by ligation of the two exons and release of the intron. Fungi use a "heal and seal" pathway that requires three distinct catalytic domains of the tRNA ligase enzyme, Trl1. In contrast, humans use a "direct ligation" pathway carried out by RTCB, an enzyme completely unrelated to Trl1. Because of these mechanistic differences, Trl1 has been proposed as a promising drug target for fungal infections. To validate Trl1 as a broad-spectrum drug target, we show that fungi from three different phyla contain Trl1 orthologs with all three domains. This includes the major invasive human fungal pathogens, and these proteins each can functionally replace yeast Trl1. In contrast, species from the order Mucorales, including the pathogens Rhizopus arrhizus and Mucor circinelloides, contain an atypical Trl1 that contains the sealing domain, but lack both healing domains. Although these species contain fewer tRNA introns than other pathogenic fungi, they still require splicing to decode three of the 61 sense codons. These sealing-only Trl1 orthologs can functionally complement defects in the corresponding domain of yeast Trl1 and use a conserved catalytic lysine residue. We conclude that Mucorales use a sealing-only enzyme together with unidentified non-orthologous healing enzymes for their heal and seal pathway. This implies that drugs that target the sealing activity are more likely to be broader-spectrum antifungals than drugs that target the healing domains.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

Arbuscular mycorrhizal fungi heterokaryons have two nuclear populations with distinct roles in host-plant interactions

Arbuscular mycorrhizal fungi (AMF) are prominent root symbionts that can carry thousands of nuclei deriving from two parental strains in a large syncytium. These co-existing genomes can also vary in abundance with changing environmental conditions. Here we assemble the nuclear genomes of all four publicly available AMF heterokaryons using PacBio high-fidelity and Hi-C sequencing. We find that the two co-existing genomes of these strains are phylogenetically related but differ in structure, content and epigenetics. We confirm that AMF heterokaryon genomes vary in relative abundance across conditions and show this can lead to nucleus-specific differences in expression during interactions with plants. Population analyses also reveal signatures of genetic exchange indicative of past events of sexual reproduction in these strains. This work uncovers the origin and contribution of two nuclear genomes in AMF heterokaryons and opens avenues for the improvement and environmental application of these strains.

A perspective on cross-kingdom RNA interference in mutualistic symbioses

RNA interference (RNAi) is arguably one of the more versatile mechanisms in cell biology, facilitating the fine regulation of gene expression and protection against mobile genomic elements, whilst also constituting a key aspect of induced plant immunity. More recently, the use of this mechanism to regulate gene expression in heterospecific partners - cross-kingdom RNAi (ckRNAi) - has been shown to form a critical part of bidirectional interactions between hosts and endosymbionts, regulating the interplay between microbial infection mechanisms and host immunity. Here, we review the current understanding of ckRNAi as it relates to interactions between plants and their pathogenic and mutualistic endosymbionts, with particular emphasis on evidence in support of ckRNAi in the arbuscular mycorrhizal symbiosis.

Lessons from arbuscular mycorrhizal fungal genomes

Arbuscular mycorrhizal fungi (AMF) have accompanied the majority of land plants since their evolution in the Devonian period with a symbiotic alliance centered on nutrient exchanges. The exploration of AMF genomes is providing clues to explain major questions about their biology, evolution, and ecology. The dynamics of nuclei across the fungal life cycle, the abundance of transposable elements, and the epigenome landscape are emerging as sources of intraspecific variability, which can be especially important in organisms with no or rare sexual reproduction such as AMF. These features have been hypothesized to support AMF adaptability to a wide host range and to environmental changes. New insights on plant-fungus communication and on the iconic function of phosphate transport were also recently obtained that overall contribute to a better understanding of this ancient and fascinating symbiosis.

A highly contiguous genome assembly reveals sources of genomic novelty in the symbiotic fungus Rhizophagus irregularis

The root systems of most plant species are aided by the soil-foraging capacities of symbiotic arbuscular mycorrhizal (AM) fungi of the Glomeromycotina subphylum. Despite recent advances in our knowledge of the ecology and molecular biology of this mutualistic symbiosis, our understanding of the AM fungi genome biology is just emerging. Presented here is a close to T2T genome assembly of the model AM fungus Rhizophagus irregularis DAOM197198, achieved through Nanopore long-read DNA sequencing and Hi-C data. This haploid genome assembly of R. irregularis, alongside short- and long-read RNA-Sequencing data, was used to produce a comprehensive annotation catalog of gene models, repetitive elements, small RNA loci, and DNA cytosine methylome. A phylostratigraphic gene age inference framework revealed that the birth of genes associated with nutrient transporter activity and transmembrane ion transport systems predates the emergence of Glomeromycotina. While nutrient cycling in AM fungi relies on genes that existed in ancestor lineages, a burst of Glomeromycotina-restricted genetic innovation is also detected. Analysis of the chromosomal distribution of genetic and epigenetic features highlights evolutionarily young genomic regions that produce abundant small RNAs, suggesting active RNA-based monitoring of genetic sequences surrounding recently evolved genes. This chromosome-scale view of the genome of an AM fungus genome reveals previously unexplored sources of genomic novelty in an organism evolving under an obligate symbiotic life cycle.

A pathogen effector FOLD diversified in symbiotic fungi

Pathogenic fungi use secreted effector proteins to suppress immunity and support their infection, but effectors have also been reported from fungi that engage in nutritional symbioses with plants. Sequence-based effector comparisons between pathogens and symbiotic arbuscular mycorrhizal (AM) fungi are hampered by the huge diversity of effector sequences even within closely related microbes. To find sequence-divergent but structurally similar effectors shared between symbiotic and pathogenic fungi, we compared secreted protein structure models of the AM fungus Rhizophagus irregularis to known pathogen effectors. We identified proteins with structural similarity to known Fusarium oxysporum f. sp. lycopersici dual domain (FOLD) effectors, which occur in low numbers in several fungal pathogens. Contrastingly, FOLD genes from AM fungi (MycFOLDs) are found in enlarged and diversified gene families with higher levels of positive selection in their C-terminal domains. Our structure model comparison suggests that MycFOLDs are similar to carbohydrate-binding motifs. Different MycFOLD genes are expressed during colonisation of different hosts and MycFOLD-17 transcripts accumulate in plant intracellular arbuscules. The exclusive presence of MycFOLDs across unrelated plant-colonising fungi, their inducible expression, lineage-specific sequence diversification and transcripts in arbuscules suggest that FOLD proteins act as effectors during plant colonisation of symbiotic and pathogenic fungi.

Whole-genome DNA methylome analysis of different developmental stages of the entomopathogenic fungus Beauveria bassiana NCHU-157 by nanopore sequencing

The entomopathogenic fungus (EPF), Beauveria bassiana, is an important and commonly used EPF for microbial control. However, the role of DNA methylation has not been thoroughly studied. Therefore, the whole genomic DNA methylome of one promising EPF isolate, B. bassiana NCHU-157 (Bb-NCHU-157), was investigated by Oxford Nanopore Technologies (ONT). First, the whole genome of Bb-NCHU-157 was sequenced by next-generation sequencing (NGS) and ONT. The genome of Bb-NCHU-157 contains 16 contigs with 34.19 Mb and 50% GC content, which are composed of 10,848 putative protein-coding genes. Two putative DNA methyltransferases (DNMTs) were found, including Dim-2 and C-5 cytosine-specific DNA methylases. Both DNMTs showed higher expression levels in the mycelium stage than in the conidia stage, indicating that development of DNA methylation in Bb-NCHU-157 might occur in the mycelium stage. The global methylation level of the mycelium stage (5 mC = 4.56%, CG = 3.33%, CHG = 0.74%, CHH = 0.49%) was higher than that of the conidial stage (5 mC = 2.99%, CG = 1.99%, CHG = 0.63%, CHH = 0.37%) in both the gene and transposable element (TE) regions. Furthermore, the TE regions showed higher methylation frequencies than the gene regions, especially for CHH site methylation, suggesting regulation of genomic stabilization during mycelium development. In the gene regions, high methylation frequencies were found around the transcription start site (TSS) and transcription end site (TES). Moreover, CG and CHG methylation mainly occur in the promoter and intergenic regions, while CHH methylation occurs in the TE region. Among the methylated regions, 371, 661, and 756 differentially DNA methylated regions (DMRs) were hypermethylated in the mycelium in CG, CHG, and CHH, while only 13 and 7 DMRs were hypomethylated in the mycelium in CHG, and CHH, respectively. Genes located in the DMR shared the GO terms, DNA binding (GO: 0003677), and sequence-specific DNA binding (GO: 0043565) for hypermethylation in the mycelium, suggesting that methylation might regulate gene expression from the initial process. Evaluation of the DNA methylome in Bb-NCHU-157 by ONT provided new insight into this field. These data will be further validated, and epigenetic regulation during the development of B. bassiana will be explored.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

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.