Whole-genome sequencing of Ganoderma boninense, the causal agent of basal stem rot disease in oil palm, via combined short- and long-read sequencing.

The hemibiotrophic Basidiomycete pathogen Ganoderma boninense (Gb) is the dominant causal agent of oil palm basal stem rot disease. Here, we report a complete chromosomal genome map of Gb using a combination of short-read Illumina and long-read Pacific Biosciences (PacBio) sequencing platforms combined with chromatin conformation capture data from the Chicago and Hi-C platforms. The genome was 55.87 Mb in length and assembled to a high contiguity (N50: 304.34 kb) of 12 chromosomes built from 112 scaffolds, with a total of only 4.34 Mb (~ 7.77%) remaining unplaced. The final assemblies were evaluated for completeness of the genome by using Benchmarking Universal Single Copy Orthologs (BUSCO) v4.1.4, and based on 4464 total BUSCO polyporales group searches, the assemblies yielded 4264 (95.52%) of the conserved orthologs as complete and only a few fragmented BUSCO of 42 (0.94%) as well as a missing BUSCO of 158 (3.53%). Genome annotation predicted a total of 21,074 coding genes, with a GC content ratio of 59.2%. The genome features were analyzed with different databases, which revealed 2471 Gene Ontology/GO (11.72%), 5418 KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthologous/KO (25.71%), 13,913 Cluster of Orthologous Groups of proteins/COG (66.02%), 60 ABC transporter (0.28%), 1049 Carbohydrate-Active Enzymes/CAZy (4.98%), 4005 pathogen-host interactions/PHI (19%), and 515 fungal transcription factor/FTFD (2.44%) genes. The results obtained in this study provide deep insight for further studies in the future.

High intragenomic, intergenomic, and phenotypic diversity in pulcherrimin-producing Metschnikowia yeasts indicates a special mode of genome evolution.

In molecular systematics, the delimitation of yeast species is based on the notion that the barcode differences are smaller within species than between them. The most widely used barcodes are segments of the chromosomal repeats coding for ribosomal RNAs that are homogenised in yeasts. The analysis of these segments of the type strains of ten species recently merged in Metschnikowia pulcherrima and 37 new isolates demonstrated that this is not the case in this species. The intragenomic diversity significantly exceeded the threshold gaps used to differentiate related yeast species. Large segments of the D1/D2 domains were not diverse within the genomes and could therefore be used to determine the taxonomic affiliation of the isolates. The genome structures of the isolates were compared by RAPD and the RFLP of the mitochondrial DNA. Both patterns were highly heterogeneous. The sequence analysis of the PUL4 gene (a member of the PUL gene cluster involved in pulcherrimin production) revealed very high intragenomic differences, suggesting that the genomes may be chimerised. Three phenotypic traits related to the antimicrobial antagonism characteristic of the species were also highly diverse and prone to reversible segregation resembling epigenetic processes (silencing and reactivation of regulators) rather than mutations and back-mutations. These features make M. pulcherrima unique among yeasts and indicate that it evolves in a non-standard way.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Complete genome sequencing of Hortaea werneckii M-3 for identifying polyester polyurethane degrading enzymes.

Hortaea werneckii M-3, a black yeast isolated from the marine sediment of the West Pacific, can utilize polyester polyurethane (PU, Impranil DLN) as a sole carbon source. Here, we present the complete genome of Hortaea werneckii M-3 with the focus on PU degradation enzymes. The total genome size is 38,167,921 bp, consisting of 186 contigs with a N50 length of 651,266 bp and a GC content of 53.06%. Genome annotation analysis predicts a total of 13,462 coding genes, which include 99 tRNAs and 105 rRNAs. Some genes encoding PU degrading enzymes including cutinase and urease are identified in this genome. The genome analysis of Hortaea werneckii M-3 will be helpful for further understanding the degradation mechanism of polyester PU by marine yeasts.

Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components.

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae.

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

An Optimized Genotyping Workflow for Identifying Highly SCRaMbLEd Synthetic Yeasts.

Synthetic Sc2.0 yeast strains contain hundreds to thousands of loxPsym recombination sites that allow restructuring of the Saccharomyces cerevisiae genome by SCRaMbLE. Thus, a highly diverse yeast population can arise from a single genotype. The selection of genetically diverse candidates with rearranged synthetic chromosomes for downstream analysis requires an efficient and straightforward workflow. Here we present loxTags, a set of qPCR primers for genotyping across loxPsym sites to detect not only deletions but also inversions and translocations after SCRaMbLE. To cope with the large number of amplicons, we generated qTagGer, a qPCR genotyping primer prediction tool. Using loxTag-based genotyping and long-read sequencing, we show that light-inducible Cre recombinase L-SCRaMbLE can efficiently generate diverse recombination events when applied to Sc2.0 strains containing a linear or a circular version of synthetic chromosome III.

Genome and secretome insights: unravelling the lignocellulolytic potential of Myceliophthora verrucosa for enhanced hydrolysis of lignocellulosic biomass.

Lignocellulolytic enzymes from a novel Myceliophthora verrucosa (5DR) strain was found to potentiate the efficacy of benchmark cellulase during saccharification of acid/alkali treated bagasse by ~ 2.24 fold, indicating it to be an important source of auxiliary enzymes. The De-novo sequencing and analysis of M. verrucosa genome (31.7 Mb) revealed to encode for 7989 putative genes, representing a wide array of CAZymes (366) with a high proportions of auxiliary activity (AA) genes (76). The LC/MS QTOF based secretome analysis of M. verrucosa showed high abundance of glycosyl hydrolases and AA proteins with cellobiose dehydrogenase (CDH) (AA8), being the most prominent auxiliary protein. A gene coding for lytic polysaccharide monooxygenase (LPMO) was expressed in Pichia pastoris and CDH produced by M. verrucosa culture on rice straw based solidified medium were purified and characterized. The mass spectrometry of LPMO catalyzed hydrolytic products of avicel showed the release of both C1/C4 oxidized products, indicating it to be type-3. The lignocellulolytic cocktail comprising of in-house cellulase produced by Aspergillus allahabadii strain spiked with LPMO & CDH exhibited enhanced and better hydrolysis of mild alkali deacetylated (MAD) and unwashed acid pretreated rice straw slurry (UWAP), when compared to Cellic CTec3 at high substrate loading rate.

Methodological advances enabled by the construction of a synthetic yeast genome.

The international Synthetic Yeast Project (Sc2.0) aims to construct the first synthetic designer eukaryote genome. Over the past few years, the Sc2.0 consortium has achieved several significant milestones by synthesizing and characterizing all 16 nuclear chromosomes of the yeast Saccharomyces cerevisiae, as well as a 17thde novo neochromosome containing all nuclear tRNA genes. In this commentary, we discuss the recent technological advances achieved in this project and provide a perspective on how they will impact the emerging field of synthetic genomics in the future.

Draft Genome Sequence of Candida saopaulonensis from a Very Premature Infant with Sepsis.

The rare fungus Candida saopaulonensis has never been reported to be associated with human infection. We report the draft genome sequence of the first clinical isolate of C. saopaulonensis, which was isolated from a very premature infant with sepsis. This is the first genome assembly reaching the near-complete chromosomal level with structural annotation for this species, opening up avenues for exploring evolutionary patterns and genetic mechanisms of pathogenesis.

Machine learning enables identification of an alternative yeast galactose utilization pathway.

How genomic differences contribute to phenotypic differences is a major question in biology. The recently characterized genomes, isolation environments, and qualitative patterns of growth on 122 sources and conditions of 1,154 strains from 1,049 fungal species (nearly all known) in the yeast subphylum Saccharomycotina provide a powerful, yet complex, dataset for addressing this question. We used a random forest algorithm trained on these genomic, metabolic, and environmental data to predict growth on several carbon sources with high accuracy. Known structural genes involved in assimilation of these sources and presence/absence patterns of growth in other sources were important features contributing to prediction accuracy. By further examining growth on galactose, we found that it can be predicted with high accuracy from either genomic (92.2%) or growth data (82.6%) but not from isolation environment data (65.6%). Prediction accuracy was even higher (93.3%) when we combined genomic and growth data. After the GALactose utilization genes, the most important feature for predicting growth on galactose was growth on galactitol, raising the hypothesis that several species in two orders, Serinales and Pichiales (containing the emerging pathogen Candida auris and the genus Ogataea, respectively), have an alternative galactose utilization pathway because they lack the GAL genes. Growth and biochemical assays confirmed that several of these species utilize galactose through an alternative oxidoreductive D-galactose pathway, rather than the canonical GAL pathway. Machine learning approaches are powerful for investigating the evolution of the yeast genotype-phenotype map, and their application will uncover novel biology, even in well-studied traits.

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

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

Phylogenomics reveals extensive misidentification of fungal strains from the genus Aspergillus.

Modern taxonomic classification is often based on phylogenetic analyses of a few molecular markers, although single-gene studies are still common. Here, we leverage genome-scale molecular phylogenetics (phylogenomics) of species and populations to reconstruct evolutionary relationships in a dense data set of 710 fungal genomes from the biomedically and technologically important genus Aspergillus. To do so, we generated a novel set of 1,362 high-quality molecular markers specific for Aspergillus and provided profile Hidden Markov Models for each, facilitating their use by others. Examining the resulting phylogeny helped resolve ongoing taxonomic controversies, identified new ones, and revealed extensive strain misidentification (7.59% of strains were previously misidentified), underscoring the importance of population-level sampling in species classification. These findings were corroborated using the current standard, taxonomically informative loci. These findings suggest that phylogenomics of species and populations can facilitate accurate taxonomic classifications and reconstructions of the Tree of Life.IMPORTANCEIdentification of fungal species relies on the use of molecular markers. Advances in genomic technologies have made it possible to sequence the genome of any fungal strain, making it possible to use genomic data for the accurate assignment of strains to fungal species (and for the discovery of new ones). We examined the usefulness and current limitations of genomic data using a large data set of 710 publicly available genomes from multiple strains and species of the biomedically, agriculturally, and industrially important genus Aspergillus. Our evolutionary genomic analyses revealed that nearly 8% of publicly available Aspergillus genomes are misidentified. Our work highlights the usefulness of genomic data for fungal systematic biology and suggests that systematic genome sequencing of multiple strains, including reference strains (e.g., type strains), of fungal species will be required to reduce misidentification errors in public databases.

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.

FGCD: a database of fungal gene clusters related to secondary metabolism.

Fungal secondary metabolites are not necessary for growth, but they are important for fungal metabolism and ecology because they provide selective advantages for competition, survival and interactions with the environment. These various metabolites are widely used as medicinal precursors and insecticides. Secondary metabolism genes are commonly arranged in clusters along chromosomes, which allow for the coordinate control of complete pathways. In this study, we created the Fungal Gene Cluster Database to store, retrieve, and visualize secondary metabolite gene cluster information across fungal species. The database was created by merging data from RNA sequencing, Basic Local Alignment Search Tool, genome browser, enrichment analysis and the R Shiny web framework to visualize and query putative gene clusters. This database facilitated the rapid and thorough examination of significant gene clusters across fungal species by detecting, defining and graphically displaying the architecture, organization and expression patterns of secondary metabolite gene clusters. In general, this genomic resource makes use of the tremendous chemical variety of the products of these ecologically and biotechnologically significant gene clusters to our further understanding of fungal secondary metabolism. Database URL: https://www.hebaubioinformatics.cn/FungalGeneCluster/.

An evolutionary view of the Fusarium core genome.

Fusarium, a member of the Ascomycota fungi, encompasses several pathogenic species significant to plants and animals. Some phytopathogenic species have received special attention due to their negative economic impact on the agricultural industry around the world. Traditionally, identification and taxonomic analysis of Fusarium have relied on morphological and phenotypic features, including the fungal host, leading to taxonomic conflicts that have been solved using molecular systematic technologies. In this work, we applied a phylogenomic approach that allowed us to resolve the evolutionary history of the species complexes of the genus and present evidence that supports the F. ventricosum species complex as the most basal lineage of the genus. Additionally, we present evidence that proposes modifications to the previous hypothesis of the evolutionary history of the F. staphyleae, F. newnesense, F. nisikadoi, F. oxysporum, and F. fujikuroi species complexes. Evolutionary analysis showed that the genome GC content tends to be lower in more modern lineages, in both, the whole-genome and core-genome coding DNA sequences. In contrast, genome size gain and losses are present during the evolution of the genus. Interestingly, core genome duplication events positively correlate with genome size. Evolutionary and genome conservation analysis supports the F3 hypothesis of Fusarium as a more compact and conserved group in terms of genome conservation. By contrast, outside of the F3 hypothesis, the most basal clades only share 8.8% of its genomic sequences with the F3 clade.

Genome biology and evolution of mating-type loci in four cereal rust fungi.

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Annotation of 2,507 Saccharomyces cerevisiae genomes.

Saccharomyces cerevisiae (baker's yeast, budding yeast) is one of the most important model organisms for biological research and is a crucial microorganism in industry. Currently, a huge number of Saccharomyces cerevisiae genome sequences are available at the public domain. However, these genomes are distributed at different websites and a large number of them are released without annotation information. To provide one complete annotated genome data resource, we collected 2,507 Saccharomyces cerevisiae genome assemblies and re-annotated 2,506 assemblies using a custom annotation pipeline, producing a total of 15,407,164 protein-coding gene models. With a custom pipeline, all these gene sequences were clustered into families. A total of 1,506 single-copy genes were selected as marker genes, which were then used to evaluate the genome completeness and base qualities of all assemblies. Pangenomic analyses were performed based on a selected subset of 847 medium-high-quality genomes. Statistical comparisons revealed a number of gene families showing copy number variations among different organism sources. To the authors' knowledge, this study represents the largest genome annotation project of S. cerevisiae so far, providing rich genomic resources for the future studies of the model organism S. cerevisiae and its relatives.IMPORTANCESaccharomyces cerevisiae (baker's yeast, budding yeast) is one of the most important model organisms for biological research and is a crucial microorganism in industry. Though a huge number of Saccharomyces cerevisiae genome sequences are available at the public domain, these genomes are distributed at different websites and most are released without annotation, hindering the efficient reuse of these genome resources. Here, we collected 2,507 genomes for Saccharomyces cerevisiae, performed genome annotation, and evaluated the genome qualities. All the obtained data have been deposited at public repositories and are freely accessible to the community. This study represents the largest genome annotation project of S. cerevisiae so far, providing one complete annotated genome data set for S. cerevisiae, an important workhorse for fundamental biology, biotechnology, and industry.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi.

Fungi biosynthesize diverse secondary metabolites, small organic bioactive molecules with key roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce a cadre of secondary metabolites, some of which are useful (e.g. the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and others harmful (e.g. the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. Fungal genomes often also encode resistance genes that confer protection against toxic secondary metabolites. Some Penicillium species, such as Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity. To investigate the evolutionary conservation of homologs of the gliotoxin BGC and of genes involved in gliotoxin resistance in Penicillium, we analyzed 35 Penicillium genomes from 23 species. Homologous, lesser fragmented gliotoxin BGCs were found in 12 genomes, mostly fragmented remnants of the gliotoxin BGC were found in 21 genomes, whereas the remaining 2 Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes was observed. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes, suggestive of constraint and potential functional significance or more recent decay. Gene tree-species tree reconciliation analyses suggested that the history of homologs in the gliotoxin BGC across the genus Penicillium likely involved multiple duplications, losses, and horizontal gene transfers. Our analyses suggest that genes encoded in BGCs can have complex evolutionary histories and be retained in genomes long after the loss of secondary metabolite biosynthesis.

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).