Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae

The jet-like chromatin structure defines active secondary metabolism in fungi

Eukaryotic genomes are spatially organized within the nucleus in a nonrandom manner. However, fungal genome arrangement and its function in development and adaptation remain largely unexplored. Here, we show that the high-order chromosome structure of Fusarium graminearum is sculpted by both H3K27me3 modification and ancient genome rearrangements. Active secondary metabolic gene clusters form a structure resembling chromatin jets. We demonstrate that these jet-like domains, which can propagate symmetrically for 54 kb, are prevalent in the genome and correlate with active gene transcription and histone acetylation. Deletion of GCN5, which encodes a core and functionally conserved histone acetyltransferase, blocks the formation of the domains. Insertion of an exogenous gene within the jet-like domain significantly augments its transcription. These findings uncover an interesting link between alterations in chromatin structure and the activation of fungal secondary metabolism, which could be a general mechanism for fungi to rapidly respond to environmental cues, and highlight the utility of leveraging three-dimensional genome organization in improving gene transcription in eukaryotes.

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

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

Meiosis in an asymmetric dikaryotic genome of Tremella fuciformis Tr01 facilitates new chromosome formation

BACKGROUND

The dikaryotic stage dominates most of the life cycle in basidiomycetes, and each cell carries two different haploid nuclei. Accurate phasing of these two nuclear genomes and their interactions have long been of interest.

RESULTS

We combine PacBio HiFi reads, Nanopore ultra-long reads, and Hi-C data to generate a complete, high-quality asymmetric dikaryotic genome of Tremella fuciformis Tr01, including Haplotypes A and B genomes. We assemble a meiotic haploid DBZ04 genome and detect three recombination events in these two haplotypes. We identify several chromosomal rearrangements that lead to differences in chromosome number, length, content, and sequence arrangement between these two haplotypes. Each nucleus contains a two-speed genome, harboring three accessory chromosomes and two accessory compartments that affect horizontal chromatin transfer between nuclei. We find few basidiospores are ejected from fruiting bodies of Tr01. Most monospore isolates sequenced belong to Tr01-Haplotype A genome architecture. More than one-third of monospore isolates carry one or two extra chromosomes including Chr12B and two new chromosomes ChrN1 and ChrN2. We hypothesize that homologous regions of seven sister chromatids pair into a large complex during meiosis, followed by inter-chromosomal recombination at physical contact sites and formation of new chromosomes.

CONCLUSION

We assemble two haplotype genomes of T. fuciformis Tr01 and provide the first overview of basidiomycetous genomes with discrete genomic architecture. Meiotic activities of asymmetric dikaryotic genomes result in formation of new chromosomes, aneuploidy of some daughter cells, and inviability of most other daughter cells. We propose a new approach for breeding of sporeless mushroom.

Conserved satellite DNA motif and lack of interstitial telomeric sites in highly rearranged African Nothobranchius killifish karyotypes

Using African annual killifishes of the genus Nothobranchius from temporary savannah pools with rapid karyotype and sex chromosome evolution, we analysed the chromosomal distribution of telomeric (TTAGGG)n repeat and Nfu-SatC satellite DNA (satDNA; isolated from Nothobranchius furzeri) in 15 species across the Nothobranchius killifish phylogeny, and with Fundulosoma thierryi as an out-group. Our fluorescence in situ hybridization experiments revealed that all analysed taxa share the presence of Nfu-SatC repeat but with diverse organization and distribution on chromosomes. Nfu-SatC landscape was similar in conspecific populations of Nothobranchius guentheri and Nothobranchius melanospilus but slightly-to-moderately differed between populations of Nothobranchius pienaari, and between closely related Nothobranchius kuhntae and Nothobranchius orthonotus. Inter-individual variability in Nfu-SatC patterns was found in N. orthonotus and Nothobranchius krysanovi. We revealed mostly no sex-linked patterns of studied repetitive DNA distribution. Only in Nothobranchius brieni, possessing multiple sex chromosomes, Nfu-SatC repeat occupied a substantial portion of the neo-Y chromosome, similarly as formerly found in the XY sex chromosome system of turquoise killifish N. furzeri and its sister species Nothobranchius kadleci-representatives not closely related to N. brieni. All studied species further shared patterns of expected telomeric repeats at the ends of all chromosomes and no additional interstitial telomeric sites. In summary, we revealed (i) the presence of conserved satDNA class in Nothobranchius clades (a rare pattern among ray-finned fishes); (ii) independent trajectories of Nothobranchius sex chromosome differentiation, with recurrent and convergent accumulation of Nfu-SatC on the Y chromosome in some species; and (iii) genus-wide shared tendency to loss of telomeric repeats during interchromosomal rearrangements. Collectively, our findings advance our understanding of genome structure, mechanisms of karyotype reshuffling, and sex chromosome differentiation in Nothobranchius killifishes from the genus-wide perspective.

Genome-wide expression QTL mapping reveals the highly dynamic regulatory landscape of a major wheat pathogen

BACKGROUND

In agricultural ecosystems, outbreaks of diseases are frequent and pose a significant threat to food security. A successful pathogen undergoes a complex and well-timed sequence of regulatory changes to avoid detection by the host immune system; hence, well-tuned gene regulation is essential for survival. However, the extent to which the regulatory polymorphisms in a pathogen population provide an adaptive advantage is poorly understood.

RESULTS

We used Zymoseptoria tritici, one of the most important pathogens of wheat, to generate a genome-wide map of regulatory polymorphism governing gene expression. We investigated genome-wide transcription levels of 146 strains grown under nutrient starvation and performed expression quantitative trait loci (eQTL) mapping. We identified cis-eQTLs for 65.3% of all genes and the majority of all eQTL loci are within 2kb upstream and downstream of the transcription start site (TSS). We also show that polymorphism in different gene elements contributes disproportionally to gene expression variation. Investigating regulatory polymorphism in gene categories, we found an enrichment of regulatory variants for genes predicted to be important for fungal pathogenesis but with comparatively low effect size, suggesting a separate layer of gene regulation involving epigenetics. We also show that previously reported trait-associated SNPs in pathogen populations are frequently cis-regulatory variants of neighboring genes with implications for the trait architecture.

CONCLUSIONS

Overall, our study provides extensive evidence that single populations segregate large-scale regulatory variation and are likely to fuel rapid adaptation to resistant hosts and environmental change.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain (TAD)-like structures formed by heterochromatic region aggregation. However, insufficient data exists on how histone post-translational modifications, including acetylation, affect genome organization. In Neurospora, the HCHC complex (comprised of the proteins HDA-1, CDP-2, HP1, and CHAP) deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation and alters the methylation of cytosines in DNA. Here, we assess if the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and inter-chromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone post-translational modifications genome-wide: without CDP-2, heterochromatic H4K16 acetylation is increased, yet smaller heterochromatic regions lose H3K9 trimethylation and gain inter-heterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder, as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

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.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

Investigating the origin of subtelomeric and centromeric AT-rich elements in Aspergillus flavus

An in silico study of Aspergillus flavus genome stability uncovered significant variations in both coding and non-coding regions. The non-coding insertions uniformly consisted of AT-rich sequences that are evolutionarily maintained, albeit distributed at widely different sites in an array of A. flavus strains. A survey of ≥ 2kb AT-rich elements (AT ≥ 70%; ATEs) in non-centromeric regions uncovered two major categories of ATEs. The first category is composed of homologous insertions at ectopic, non-allelic sites that contain homology to transposable elements (TEs; Classes B, C, D, and E). Strains differed significantly in frequency, position, and TE type, but displayed a common enrichment in subtelomeric regions. The TEs were heavily mutated, with patterns consistent with the ancestral activity of repeat-induced point mutations (RIP). The second category consists of a conserved set of novel subtelomeric ATE repeats (Classes A, G, G, H, I and J) which lack discernible TEs and, unlike TEs, display a constant polarity relative to the telomere. Members of one of these classes are derivatives of a progenitor ATE that is predicted to have undergone extensive homologous recombination during evolution. A third category of ATEs consists of ~100 kb regions at each centromere. Centromeric ATEs and TE clusters within these centromeres display a high level of sequence identity between strains. These studies suggest that transposition and RIP are forces in the evolution of subtelomeric and centromeric structure and function.

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.

Folding Features and Dynamics of 3D Genome Architecture in Plant Fungal Pathogens

The folding and dynamics of three-dimensional (3D) genome organization are fundamental for eukaryotes executing genome functions but have been largely unexplored in nonmodel fungi. Using high-throughput sequencing coupled with chromosome conformation capture (Hi-C) data, we generated two chromosome-level assemblies for Puccinia striiformis f. sp. tritici, a fungus causing stripe rust disease on wheat, for studying 3D genome architectures of plant pathogenic fungi. The chromatin organization of the fungus followed a combination of the fractal globule model and the equilibrium globule model. Surprisingly, chromosome compartmentalization was not detected. Dynamics of 3D genome organization during two developmental stages of P. striiformis f. sp. tritici indicated that regulation of gene activities might be independent of the changes of genome organization. In addition, chromatin conformation conservation was found to be independent of genome sequence synteny conservation among different fungi. These results highlighted the distinct folding principles of fungal 3D genomes. Our findings should be an important step toward a holistic understanding of the principles and functions of genome architecture across different eukaryotic kingdoms. IMPORTANCE Previously, our understanding of 3D genome architecture has mainly come from model mammals, insects, and plants. However, the organization and regulatory functions of 3D genomes in fungi are largely unknown. In this study, we comprehensively investigated P. striiformis f. sp. tritici, a plant fungal pathogen, and revealed distinct features of the 3D genome, comparing it with the universal folding feature of 3D genomes in higher eukaryotic organisms. We further suggested that there might be distinct regulatory mechanisms of gene expression that are independent of chromatin organization changes during the developmental stages of this rust fungus. Moreover, we showed that the evolutionary pattern of 3D genomes in this fungus is also different from the cases in mammalian genomes. In addition, the genome assembly pipeline and the generated two chromosome-level genomes will be valuable resources. These results highlighted the unexplored distinct features of 3D genome organization in fungi. Therefore, our study provided complementary knowledge to holistically understand the organization and functions of 3D genomes across different eukaryotes.

Transcriptomic study reveals lncRNA-mediated downregulation of innate immune and inflammatory response in the SARS-CoV-2 vaccination breakthrough infections

INTRODUCTION

The emergence of multiple variants of concerns (VOCs) with higher number of Spike mutations have led to enhanced immune escape by the SARS-CoV-2. With the increasing number of vaccination breakthrough (VBT) infections, it is important to understand the possible reason/s of the breakthrough infections.

METHODS

We performed transcriptome sequencing of 57 VBT and unvaccinated COVID-19 patients, followed by differential expression and co-expression analysis of the lncRNAs and the mRNAs. The regulatory mechanism was highlighted by analysis towards repeat element distribution within the co-expressed lncRNAs, followed by repeats driven homologous interaction between the lncRNAs and the promoter regions of genes from the same topologically associated domains (TAD).

RESULTS

We identified 727 differentially expressed lncRNAs (153 upregulated and 574 downregulated) and 338 mRNAs (34 up- and 334 downregulated) in the VBT patients. This includes LUCAT1, MALAT1, ROR1-AS1, UGDH-AS1 and LINC00273 mediated modulation of immune response, whereas MALAT1, NEAT1 and GAS5 regulated inflammatory response in the VBT. LncRNA-mRNA co-expression analysis highlighted 34 lncRNAs interacting with 267 mRNAs. We also observed a higher abundance of Alu, LINE1 and LTRs within the interacting lncRNAs of the VBT patients. These interacting lncRNAs have higher interaction with the promoter region of the genes from the same TAD, compared to the non-interacting lncRNAs with the enrichment of Alu and LINE1 in the gene promoter.

DISCUSSION

Significant downregulation and GSEA of the TAD gene suggest Alu and LINE1 driven homologous interaction between the lncRNAs and the TAD genes as a possible mechanism of lncRNA-mediated suppression of innate immune/inflammatory responses and activation of adaptive immune response. The lncRNA-mediated suppression of innate immune/inflammatory responses and activation of adaptive immune response might explain the SARS-CoV-2 breakthrough infections with milder symptoms in the VBT. Besides, the study also highlights repeat element mediated regulation of genes in 3D as another possible way of lncRNA-mediated immune-regulation modulating vaccination breakthroughs milder disease phenotype and shorter hospital stay.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

A chromosome-level, fully phased genome assembly of the oat crown rust fungus Puccinia coronata f. sp. avenae: a resource to enable comparative genomics in the cereal rusts

Advances in sequencing technologies as well as development of algorithms and workflows have made it possible to generate fully phased genome references for organisms with nonhaploid genomes such as dikaryotic rust fungi. To enable discovery of pathogen effectors and further our understanding of virulence evolution, we generated a chromosome-scale assembly for each of the 2 nuclear genomes of the oat crown rust pathogen, Puccinia coronata f. sp. avenae (Pca). This resource complements 2 previously released partially phased genome references of Pca, which display virulence traits absent in the isolate of historic race 203 (isolate Pca203) which was selected for this genome project. A fully phased, chromosome-level reference for Pca203 was generated using PacBio reads and Hi-C data and a recently developed pipeline named NuclearPhaser for phase assignment of contigs and phase switch correction. With 18 chromosomes in each haplotype and a total size of 208.10 Mbp, Pca203 has the same number of chromosomes as other cereal rust fungi such as Puccinia graminis f. sp. tritici and Puccinia triticina, the causal agents of wheat stem rust and wheat leaf rust, respectively. The Pca203 reference marks the third fully phased chromosome-level assembly of a cereal rust to date. Here, we demonstrate that the chromosomes of these 3 Puccinia species are syntenous and that chromosomal size variations are primarily due to differences in repeat element content.

Telomere-to-Telomere Genome Sequences across a Single Genus Reveal Highly Variable Chromosome Rearrangement Rates but Absolute Stasis of Chromosome Number

Genome rearrangements in filamentous fungi are prevalent but little is known about the modalities of their evolution, in part because few complete genomes are available within a single genus. To address this, we have generated and compared 15 complete telomere-to-telomere genomes across the phylogeny of a single genus of filamentous fungi, Epichloë. We find that the striking distinction between gene-rich and repeat-rich regions previously reported for isolated species is ubiquitous across the Epichloë genus. We built a species phylogeny from single-copy gene orthologs to provide a comparative framing to study chromosome composition and structural change through evolutionary time. All Epichloë genomes have exactly seven nuclear chromosomes, but despite this conserved ploidy, analyses reveal low synteny and substantial rearrangement of gene content across the genus. These rearrangements are highly lineage-dependent, with most occurring over short evolutionary distances, with long periods of structural stasis. Quantification of chromosomal rearrangements shows they are uncorrelated with numbers of substitutions and evolutionary distances, suggesting that different modes of evolution are acting to create nucleotide and chromosome-scale changes.

The genome organization of Neurospora crassa at high resolution uncovers principles of fungal chromosome topology

The eukaryotic genome must be precisely organized for its proper function, as genome topology impacts transcriptional regulation, cell division, replication, and repair, among other essential processes. Disruptions to human genome topology can lead to diseases, including cancer. The advent of chromosome conformation capture with high-throughput sequencing (Hi-C) to assess genome organization has revolutionized the study of nuclear genome topology; Hi-C has elucidated numerous genomic structures, including chromosomal territories, active/silent chromatin compartments, Topologically Associated Domains, and chromatin loops. While low-resolution heatmaps can provide important insights into chromosomal level contacts, high-resolution Hi-C datasets are required to reveal folding principles of individual genes. Of particular interest are high-resolution chromosome conformation datasets of organisms modeling the human genome. Here, we report the genome topology of the fungal model organism Neurospora crassa at a high resolution. Our composite Hi-C dataset, which merges 2 independent datasets generated with restriction enzymes that monitor euchromatin (DpnII) and heterochromatin (MseI), along with our DpnII/MseI double digest dataset, provide exquisite detail for both the conformation of entire chromosomes and the folding of chromatin at the resolution of individual genes. Within constitutive heterochromatin, we observe strong yet stochastic internal contacts, while euchromatin enriched with either activating or repressive histone post-translational modifications associates with constitutive heterochromatic regions, suggesting intercompartment contacts form to regulate transcription. Consistent with this, a strain with compromised heterochromatin experiences numerous changes in gene expression. Our high-resolution Neurospora Hi-C datasets are outstanding resources to the fungal community and provide valuable insights into higher organism genome topology.

Cross-species transcriptomics identifies core regulatory changes differentiating the asymptomatic asexual and virulent sexual life cycles of grass-symbiotic Epichloë fungi

Fungi from the genus Epichloë form systemic endobiotic infections of cool season grasses, producing a range of host-protective natural products in return for access to nutrients. These infections are asymptomatic during vegetative host growth, with associations between asexual Epichloë spp. and their hosts considered mutualistic. However, the sexual cycle of Epichloë spp. involves virulent growth, characterized by the envelopment and sterilization of a developing host inflorescence by a dense sheath of mycelia known as a stroma. Microscopic analysis of stromata revealed a dramatic increase in hyphal propagation and host degradation compared with asymptomatic tissues. RNAseq was used to identify differentially expressed genes in asymptomatic vs stromatized tissues from 3 diverse Epichloë-host associations. Comparative analysis identified a core set of 135 differentially expressed genes that exhibited conserved transcriptional changes across all 3 associations. The core differentially expressed genes more strongly expressed during virulent growth encode proteins associated with host suppression, digestion, adaptation to the external environment, a biosynthetic gene cluster, and 5 transcription factors that may regulate Epichloë stroma formation. An additional 5 transcription factor encoding differentially expressed genes were suppressed during virulent growth, suggesting they regulate mutualistic processes. Expression of biosynthetic gene clusters for natural products that suppress herbivory was universally suppressed during virulent growth, and additional biosynthetic gene clusters that may encode production of novel host-protective natural products were identified. A comparative analysis of 26 Epichloë genomes found a general decrease in core differentially expressed gene conservation among asexual species, and a specific decrease in conservation for the biosynthetic gene cluster expressed during virulent growth and an unusual uncharacterized gene.

Nucleic acids movement and its relation to genome dynamics of repetitive DNA: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?

There is growing evidence of evolutionary genome plasticity. The evolution of repetitive DNA elements, the major components of most eukaryotic genomes, involves the amplification of various classes of mobile genetic elements, the expansion of satellite DNA, the transfer of fragments or entire organellar genomes and may have connections with viruses. In addition to various repetitive DNA elements, a plethora of large and small RNAs migrate within and between cells during individual development as well as during evolution and contribute to changes of genome structure and function. Such migration of DNA and RNA molecules often results in horizontal gene transfer, thus shaping the whole genomic network of interconnected species. Here, we propose that a high evolutionary dynamism of repetitive genome components is often related to the migration/movement of DNA or RNA molecules. We speculate that the cytoplasm is probably an ideal compartment for such evolutionary experiments.

Epichloë scottii sp. nov., a new endophyte isolated from Melica uniflora is the missing ancestor of Epichloë disjuncta

Here we describe a new, haploid and stroma forming species within the genus Epichloë, as Epichloë scottii sp. nov. The fungus was isolated from Melica uniflora growing in Bad Harzburg, Germany. Phylogenetic reconstruction using a combined dataset of the tubB and tefA genes strongly support that E. scottii is a distinct species and the so far unknown ancestor species of the hybrid E. disjuncta. A distribution analysis showed a high infection rate in close vicinity of the initial sampling site and only two more spots with low infection rates. Genetic variations in key genes required for alkaloid production suggested that E. scottii sp. nov. might not be capable of producing any of the major alkaloids including ergot alkaloid, loline, indole-diterpene and peramine. All isolates and individuals found in the distribution analysis were identified as mating-type B explaining the lack of mature stromata during this study. We further release a telomere-to-telomere de novo assembly of all seven chromosomes and the mitogenome of E. scottii sp. nov.

Long reads and Hi-C sequencing illuminate the two-compartment genome of the model arbuscular mycorrhizal symbiont Rhizophagus irregularis

Chromosome folding links genome structure with gene function by generating distinct nuclear compartments and topologically associating domains. In mammals, these undergo preferential interactions and regulate gene expression. However, their role in fungal genome biology is unclear. Here, we combine Nanopore (ONT) sequencing with chromatin conformation capture sequencing (Hi-C) to reveal chromosome and epigenetic diversity in a group of obligate plant symbionts: the arbuscular mycorrhizal fungi (AMF). We find that five phylogenetically distinct strains of the model AMF Rhizophagus irregularis carry 33 chromosomes with substantial within-species variability in size, as well as in gene and repeat content. Strain-specific Hi-C contact maps reveal a 'checkerboard' pattern that underline two dominant euchromatin (A) and heterochromatin (B) compartments. Each compartment differs in the level of gene transcription, regulation of candidate effectors and methylation frequencies. The A-compartment is more gene-dense and contains most core genes, while the B-compartment is more repeat-rich and has higher rates of chromosomal rearrangement. While the B-compartment is transcriptionally repressed, it has significantly more secreted proteins and in planta upregulated candidate effectors, suggesting a possible host-induced change in chromosome conformation. Overall, this study provides a fine-scale view into the genome biology and evolution of model plant symbionts, and opens avenues to study the epigenetic mechanisms that modify chromosome folding during host-microbe interactions.

Metabolic Potential of Epichloë Endophytes for Host Grass Fungal Disease Resistance

Asexual species of the genus Epichloë (Clavicipitaceae, Ascomycota) form endosymbiotic associations with Pooidae grasses. This association is important both ecologically and to the pasture and turf industries, as the endophytic fungi confer a multitude of benefits to their host plant that improve competitive ability and performance such as growth promotion, abiotic stress tolerance, pest deterrence and increased host disease resistance. Biotic stress tolerance conferred by the production of bioprotective metabolites has a critical role in an industry context. While the known antimammalian and insecticidal toxins are well characterized due to their impact on livestock welfare, antimicrobial metabolites are less studied. Both pasture and turf grasses are challenged by many phytopathogenic diseases that result in significant economic losses and impact livestock health. Further investigations of Epichloë endophytes as natural biocontrol agents can be conducted on strains that are safe for animals. With the additional benefits of possessing host disease resistance, these strains would increase their commercial importance. Field reports have indicated that pasture grasses associated with Epichloë endophytes are superior in resisting fungal pathogens. However, only a few antifungal compounds have been identified and chemically characterized, and these from sexual (pathogenic) Epichloë species, rather than those utilized to enhance performance in turf and pasture industries. This review provides insight into the various strategies reported in identifying antifungal activity from Epichloë endophytes and, where described, the associated antifungal metabolites responsible for the activity.

Genetic Diversity of the Symbiotic Fungus Epichloë festucae in Naturally Occurring Host Grass Populations

Epichloë festucae is a common symbiont of the perennial and widely distributed cool season grass, Festuca rubra. The symbiosis is highly integrated involving systemic growth of the fungus throughout above-ground host parts and vertical transmission from plant to its offspring via host seeds. However, the nature of symbiosis is labile ranging from antagonistic to mutualistic depending on prevailing selection pressures. Both the loss of fungus in the maternal host lineage and horizontal transmission through sexual spores within the host population may partly explain the detected variation in symbiosis in wild grass populations. Epichloë species are commonly considered as pathogens when they produce sexual spores and partly castrate their host plant. This is the pathogenic end of the continuum from antagonistic to mutualistic interactions. Here we examined the population genetic structure of E. festucae to reveal the gene flow, importance of reproduction modes, and alkaloid potential of the symbiotic fungus in Europe. Epichloë-species are highly dependent on the host in survival and reproduction whilst benefits to the host are largely linked to defensive mutualism attributable to fungal-origin bioactive alkaloids that negatively affect vertebrate and/or invertebrate herbivores. We detected decreased genetic diversity in previously glaciated areas compared to non-glaciated regions during the last glacial maximum period and found three major genetic clusters in E. festucae populations: southern, northeastern and northwestern Europe. Sexual reproduction may have a higher role than expected in Spanish E. festucae populations due to the predominance of unique genotypes and presence of both mating types in the region. In contrast, asexual reproduction via host seeds predominates in the Faroe Island and Finland in northern Europe due to the presence of biased mating-type ratios and large dominant genotypes in the E. festucae populations within the region. A substantially larger variation of alkaloid genotypes was observed in the fungal populations than expected, although the variability of the alkaloid genotypes within populations is considerably lower in northern than Spanish populations in southern Europe. E. festucae populations consist of different combinations of alkaloid classes from the gene clusters of ergot alkaloid and indole-terpenes, and from pyrrolopyrazine alkaloid gene. We suggest that the postglacial distribution history of the host grass, prevailing reproduction strategies of E. festucae, and local selection pressures likely explain a large part of the genetic variation observed in fungal populations among geographic regions. The identified alkaloid genotypes can be used by turfgrass breeders to improve resistance against herbivores in red fescue varieties and to develop new sustainable cultivars in Europe.

Neuron-specific chromosomal megadomain organization is adaptive to recent retrotransposon expansions

Regulatory mechanisms associated with repeat-rich sequences and chromosomal conformations in mature neurons remain unexplored. Here, we map cell-type specific chromatin domain organization in adult mouse cerebral cortex and report strong enrichment of Endogenous Retrovirus 2 (ERV2) repeat sequences in the neuron-specific heterochromatic B₂^NeuN+ megabase-scaling subcompartment. Single molecule long-read sequencing and comparative Hi-C chromosomal contact mapping in wild-derived SPRET/EiJ (Mus spretus) and laboratory inbred C57BL/6J (Mus musculus) reveal neuronal reconfigurations tracking recent ERV2 expansions in the murine germline, with significantly higher B₂^NeuN+ contact frequencies at sites with ongoing insertions in Mus musculus. Neuronal ablation of the retrotransposon silencer Kmt1e/Setdb1 triggers B₂^NeuN+ disintegration and rewiring with open chromatin domains enriched for cellular stress response genes, along with severe neuroinflammation and proviral assembly with infiltration of dendrites . We conclude that neuronal megabase-scale chromosomal architectures include an evolutionarily adaptive heterochromatic organization which, upon perturbation, results in transcriptional dysregulation and unleashes ERV2 proviruses with strong neuronal tropism.

Chromosome-level genomes provide insights into genome evolution, organization and size in Epichloe fungi

Epichloe fungi are endophytes of cool season grasses, both wild species and commercial cultivars, where they may exhibit mutualistic or pathogenic lifestyles. The Epichloe-grass symbiosis is of great interest to agricultural research for the fungal bioprotective properties conferred to host grasses but also serves as an ideal system to study the evolution of fungal plant-pathogens in natural environments. Here, we assembled and annotated gapless chromosome-level genomes of two pathogenic Epichloe sibling species. Both genomes have a bipartite genome organization, with blocks of highly syntenic gene-rich regions separated by blocks of AT-rich DNA. The AT-rich regions show an extensive signature of RIP (repeat-induced point mutation) and the expansion of this compartment accounts for the large difference in genome size between the two species. This study reveals how the rapid evolution of repeat structure can drive divergence between closely related taxa and highlights the evolutionary role of dynamic compartments in fungal genomes.

Three-dimensional genome organization via triplex-forming RNAs

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF.

Emerging Roles of Repetitive and Repeat-Containing RNA in Nuclear and Chromatin Organization and Gene Expression

Genomic repeats have been intensely studied as regulatory elements controlling gene transcription, splicing and genome architecture. Our understanding of the role of the repetitive RNA such as the RNA coming from genomic repeats, or repetitive sequences embedded in mRNA/lncRNAs, in nuclear and cellular functions is instead still limited. In this review we discuss evidence supporting the multifaceted roles of repetitive RNA and RNA binding proteins in nuclear organization, gene regulation, and in the formation of dynamic membrane-less aggregates. We hope that our review will further stimulate research in the consolidating field of repetitive RNA biology.

Chitin Deacetylases Are Required for Epichloë festucae Endophytic Cell Wall Remodeling During Establishment of a Mutualistic Symbiotic Interaction with Lolium perenne

Epichloë festucae forms a mutualistic symbiotic association with Lolium perenne. This biotrophic fungus systemically colonizes the intercellular spaces of aerial tissues to form an endophytic hyphal network and also grows as an epiphyte. However, little is known about the cell wall-remodeling mechanisms required to avoid host defense and maintain intercalary growth within the host. Here, we use a suite of molecular probes to show that the E. festucae cell wall is remodeled by conversion of chitin to chitosan during infection of L. perenne seedlings, as the hyphae switch from free-living to endophytic growth. When hyphae transition from endophytic to epiphytic growth, the cell wall is remodeled from predominantly chitosan to chitin. This conversion from chitin to chitosan is catalyzed by chitin deacetylase. The genome of E. festucae encodes three putative chitin deacetylases, two of which (cdaA and cdaB) are expressed in planta. Deletion of either of these genes results in disruption of fungal intercalary growth in the intercellular spaces of plants infected with these mutants. These results establish that these two genes are required for maintenance of the mutualistic symbiotic interaction between E. festucae and L. perenne.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Dynamics of the compartmentalized Streptomyces chromosome during metabolic differentiation

Bacteria of the genus Streptomyces are prolific producers of specialized metabolites, including antibiotics. The linear chromosome includes a central region harboring core genes, as well as extremities enriched in specialized metabolite biosynthetic gene clusters. Here, we show that chromosome structure in Streptomyces ambofaciens correlates with genetic compartmentalization during exponential phase. Conserved, large and highly transcribed genes form boundaries that segment the central part of the chromosome into domains, whereas the terminal ends tend to be transcriptionally quiescent compartments with different structural features. The onset of metabolic differentiation is accompanied by a rearrangement of chromosome architecture, from a rather 'open' to a 'closed' conformation, in which highly expressed specialized metabolite biosynthetic genes form new boundaries. Thus, our results indicate that the linear chromosome of S. ambofaciens is partitioned into structurally distinct entities, suggesting a link between chromosome folding, gene expression and genome evolution.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

Global Changes in Asexual Epichloë Transcriptomes during the Early Stages, from Seed to Seedling, of Symbiotum Establishment

Asexual Epichloë fungi are strictly seed-transmitted endophytic symbionts of cool-season grasses and spend their entire life cycle within the host plant. Endophyte infection can confer protective benefits to its host through the production of bioprotective compounds. Inversely, plants provide nourishment and shelter to the resident endophyte in return. Current understanding of the changes in global gene expression of asexual Epichloë endophytes during the early stages of host-endophyte symbiotum is limited. A time-course study using a deep RNA-sequencing approach was performed at six stages of germination, using seeds infected with one of three endophyte strains belonging to different representative taxa. Analysis of the most abundantly expressed endophyte genes identified that most were predicted to have a role in stress and defence responses. The number of differentially expressed genes observed at early time points was greater than those detected at later time points, suggesting an active transcriptional reprogramming of endophytes at the onset of seed germination. Gene ontology enrichment analysis revealed dynamic changes in global gene expression consistent with the developmental processes of symbiotic relationships. Expression of pathway genes for biosynthesis of key secondary metabolites was studied comprehensively and fuzzy clustering identified some unique expression patterns. Furthermore, comparisons of the transcriptomes from three endophyte strains in planta identified genes unique to each strain, including genes predicted to be associated with secondary metabolism. Findings from this study highlight the importance of better understanding the unique properties of individual endophyte strains and will serve as an excellent resource for future studies of host-endophyte interactions.

Toward incorporating epigenetics into regulation of gene expression in the parasite Trichomonas vaginalis

Trichomonas vaginalis is an extracellular parasite that colonizes the human urogenital tract, causing a highly prevalent sexually transmitted infection. The parasite must change its transcriptional profile in order to establish and maintain infection. However, few core regulatory elements and transcription factors have been identified to date and little is known about other mechanisms that may control these rapid changes in gene expression during parasite infection. In the last years, epigenetic mechanisms involved in the regulation of gene expression have been gaining major attention in this parasite. In this review, we summarize and discuss the major advances of the last few years with regard to epigenetics (DNA methylation, post-translational histone modifications, and histone variants) in the parasite T. vaginalis. These studies can shed light into our current understanding of this parasite's biology with far-reaching implications for the prognosis and treatment of trichomoniasis.

Trichoderma reesei Rad51 tolerates mismatches in hybrid meiosis with diverse genome sequences

Sexual eukaryotes fall into two groups with respect to their RecA-like recombinases. The first group possesses Rad51 (ubiquitous) and Dmc1 (meiosis-specific), which cooperate to conduct interhomolog recombination in zygotes with high sequence heterogeneity. Interestingly, Dmc1 was lost from the second group of eukaryotic organisms. Here we used the industrial workhorse fungus Trichoderma reesei to address if and how Rad51-only eukaryotes carry out hybrid meiosis. We show that T. reesei Rad51 (TrRad51) is indispensable for interhomolog recombination during meiosis and that TrRad51, like Saccharomyces cerevisiae Dmc1, possesses a better mismatch tolerability than S. cerevisiae Rad51. Our results indicate that the ancestral TrRad51 evolved to acquire Dmc1-like properties by adopting multiple structural variations in the L1 and L2 DNA-binding loops.

Most eukaryotes possess two RecA-like recombinases (ubiquitous Rad51 and meiosis-specific Dmc1) to promote interhomolog recombination during meiosis. However, some eukaryotes have lost Dmc1. Given that mammalian and yeast Saccharomyces cerevisiae (Sc) Dmc1 have been shown to stabilize recombination intermediates containing mismatches better than Rad51, we used the Pezizomycotina filamentous fungus Trichoderma reesei to address if and how Rad51-only eukaryotes conduct interhomolog recombination in zygotes with high sequence heterogeneity. We applied multidisciplinary approaches (next- and third-generation sequencing technology, genetics, cytology, bioinformatics, biochemistry, and single-molecule biophysics) to show that T. reesei Rad51 (TrRad51) is indispensable for interhomolog recombination during meiosis and, like ScDmc1, TrRad51 possesses better mismatch tolerance than ScRad51 during homologous recombination. Our results also indicate that the ancestral TrRad51 evolved to acquire ScDmc1-like properties by creating multiple structural variations, including via amino acid residues in the L1 and L2 DNA-binding loops.

Telomere length de novo assembly of all 7 chromosomes and mitogenome sequencing of the model entomopathogenic fungus, Metarhizium brunneum, by means of a novel assembly pipeline

Background

More accurate and complete reference genomes have improved understanding of gene function, biology, and evolutionary mechanisms. Hybrid genome assembly approaches leverage benefits of both long, relatively error-prone reads from third-generation sequencing technologies and short, accurate reads from second-generation sequencing technologies, to produce more accurate and contiguous de novo genome assemblies in comparison to using either technology independently. In this study, we present a novel hybrid assembly pipeline that allowed for both mitogenome de novo assembly and telomere length de novo assembly of all 7 chromosomes of the model entomopathogenic fungus, Metarhizium brunneum.

Results

The improved assembly allowed for better ab initio gene prediction and a more BUSCO complete proteome set has been generated in comparison to the eight current NCBI reference Metarhizium spp. genomes. Remarkably, we note that including the mitogenome in ab initio gene prediction training improved overall gene prediction. The assembly was further validated by comparing contig assembly agreement across various assemblers, assessing the assembly performance of each tool. Genomic synteny and orthologous protein clusters were compared between Metarhizium brunneum and three other Hypocreales species with complete genomes, identifying core proteins, and listing orthologous protein clusters shared uniquely between the two entomopathogenic fungal species, so as to further facilitate the understanding of molecular mechanisms underpinning fungal-insect pathogenesis.

Conclusions

The novel assembly pipeline may be used for other haploid fungal species, facilitating the need to produce high-quality reference fungal genomes, leading to better understanding of fungal genomic evolution, chromosome structuring and gene regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07390-y.

The Epichloë festucae Antifungal Protein Efe-AfpA Is also a Possible Effector Protein Required for the Interaction of the Fungus with Its Host Grass Festuca rubra subsp. rubra

Strong creeping red fescue (Festuca rubra subsp. rubra) is a commercially important low-maintenance turfgrass and is often naturally infected with the fungal endophyte Epichloë festucae. Epichloë spp. are endophytes of several cool-season grass species, often conferring insect resistance to the grass hosts due to the production of toxic alkaloids. In addition to insect resistance, a unique feature of the strong creeping red fescue/E. festucae symbiosis is the endophyte-mediated disease resistance to the fungal pathogen Clarireedia jacksonii, the causal agent of dollar spot disease. Such disease resistance is not a general feature of other grass/ Epichloë interactions. E. festucae isolates infecting red fescue have an antifungal protein gene Efe-afpA, whereas most other Epichloë spp. do not have a similar gene. The uniqueness of this gene suggests it may, therefore, be a component of the unique disease resistance seen in endophyte-infected red fescue. Here, we report the generation of CRISPR-Cas9 Efe-afpA gene knockouts with the goal of determining if absence of the protein in endophyte-infected Festuca rubra leads to disease susceptibility. However, it was not possible to infect plants with the knockout isolates, although infection was possible with the wild type E. festucae and with complemented isolates. This raises the interesting possibility that, in addition to having antifungal activity, the protein is required for the symbiotic interaction. The antifungal protein is a small secreted protein with high expression in planta relative to its expression in culture, all characteristics consistent with effector proteins. If Efe-AfpA is an effector protein it must be specific to certain interactions, since most Epichloë spp. do not have such a gene in their genomes.

Regulation of host-infection ability in the grass-symbiotic fungus Epichloë festucae by histone H3K9 and H3K36 methyltransferases

Recent studies have identified key genes that control the symbiotic interaction between Epichloë festucae and Lolium perenne. Here we report on the identification of specific E. festucae genes that control host infection. Deletion of setB, which encodes a homologue of the H3K36 histone methyltransferase Set2/KMT3, reduced histone H3K36 trimethylation and led to severe defects in colony growth and hyphal development. The E. festucae ΔclrD mutant, which lacks the gene encoding the homologue of the H3K9 methyltransferase KMT1, displays similar developmental defects. Both mutants are completely defective in their ability to infect L. perenne. Alleles that complement the culture and plant phenotypes of both mutants also complement the histone methylation defects. Co-inoculation of either ΔsetB or ΔclrD with the wild-type strain enables these mutants to colonize the host. However, successful colonization by the mutants resulted in death or stunting of the host plant. Transcriptome analysis at the early infection stage identified four fungal candidate genes, three of which encode small-secreted proteins, that are differentially regulated in these mutants compared to wild type. Deletion of crbA, which encodes a putative carbohydrate binding protein, resulted in significantly reduced host infection rates by E. festucae.

Repetitive Elements Contribute to the Diversity and Evolution of Centromeres in the Fungal Genus Verticillium

Centromeres are chromosomal regions that are crucial for chromosome segregation during mitosis and meiosis, and failed centromere formation can contribute to chromosomal anomalies. Despite this conserved function, centromeres differ significantly between and even within species. Thus far, systematic studies into the organization and evolution of fungal centromeres remain scarce. In this study, we identified the centromeres in each of the 10 species of the fungal genus Verticillium and characterized their organization and evolution. Chromatin immunoprecipitation of the centromere-specific histone CenH3 (ChIP-seq) and chromatin conformation capture (Hi-C) followed by high-throughput sequencing identified eight conserved, large (∼150-kb), AT-, and repeat-rich regional centromeres that are embedded in heterochromatin in the plant pathogen Verticillium dahliae Using Hi-C, we similarly identified repeat-rich centromeres in the other Verticillium species. Strikingly, a single degenerated long terminal repeat (LTR) retrotransposon is strongly associated with centromeric regions in some but not all Verticillium species. Extensive chromosomal rearrangements occurred during Verticillium evolution, of which some could be linked to centromeres, suggesting that centromeres contributed to chromosomal evolution. The size and organization of centromeres differ considerably between species, and centromere size was found to correlate with the genome-wide repeat content. Overall, our study highlights the contribution of repetitive elements to the diversity and rapid evolution of centromeres within the fungal genus VerticilliumIMPORTANCE The genus Verticillium contains 10 species of plant-associated fungi, some of which are notorious pathogens. Verticillium species evolved by frequent chromosomal rearrangements that contribute to genome plasticity. Centromeres are instrumental for separation of chromosomes during mitosis and meiosis, and failed centromere functionality can lead to chromosomal anomalies. Here, we used a combination of experimental techniques to identify and characterize centromeres in each of the Verticillium species. Intriguingly, we could strongly associate a single repetitive element to the centromeres of some of the Verticillium species. The presence of this element in the centromeres coincides with increased centromere sizes and genome-wide repeat expansions. Collectively, our findings signify a role of repetitive elements in the function, organization, and rapid evolution of centromeres in a set of closely related fungal species.

Third-generation sequencing-based mapping and visualization of single nucleotide polymorphism, meiotic recombination, illegitimate mutation and repeat-induced point mutation

Generation of new genetic diversity by crossover (CO) and non-crossover (NCO) is a fundamental process in eukaryotes. Fungi have played critical roles in studying this process because they permit tetrad analysis, which has been used by geneticists for several decades to determine meiotic recombination products. New genetic variations can also be generated in zygotes via illegitimate mutation (IM) and repeat-induced point mutation (RIP). RIP is a genome defense mechanism for preventing harmful expansion of transposable elements or duplicated sequences in filamentous fungi. Although the exact mechanism of RIP is unknown, the C:G to T:A mutations might result from DNA cytosine methylation. A comprehensive approach for understanding the molecular mechanisms underlying these important processes is to perform high-throughput mapping of CO, NCO, RIP and IM in zygotes bearing large numbers of heterozygous variant markers. To this aim, we developed ‘TSETA’, a versatile and user-friendly pipeline that utilizes high-quality and chromosome-level genome sequences involved in a single meiotic event of the industrial workhorse fungus Trichoderma reesei. TSETA not only can be applied to most sexual eukaryotes for genome-wide tetrad analysis, it also outcompetes most currently used methods for calling out single nucleotide polymorphisms between two or more intraspecies strains or isolates.

90 YEARS OF PROGESTERONE: Molecular mechanisms of progesterone receptor action on the breast cancer genome

Gene regulation by steroid hormones has been at the forefront in elucidating the intricacies of transcriptional regulation in eukaryotes ever since the discovery by Karlson and Clever that the insect steroid hormone ecdysone induces chromatin puffs in giant chromosomes. After the successful cloning of the hormone receptors toward the end of the past century, detailed mechanistic insight emerged in some model systems, in particular the MMTV provirus. With the arrival of next generation DNA sequencing and the omics techniques, we have gained even further insight into the global cellular response to steroid hormones that in the past decades also extended to the function of the 3D genome topology. More recently, advances in high resolution microcopy, single cell genomics and the new vision of liquid-liquid phase transitions in the context of nuclear space bring us closer than ever to unravelling the logic of gene regulation and its complex integration of global cellular signaling networks. Using the function of progesterone and its cellular receptor in breast cancer cells, we will briefly summarize the history and describe the present extent of our knowledge on how regulatory proteins deal with the chromatin structure to gain access to DNA sequences and interpret the genomic instructions that enable cells to respond selectively to external signals by reshaping their gene regulatory networks.

Lolium perenne apoplast metabolomics for identification of novel metabolites produced by the symbiotic fungus Epichloë festucae

Epichloë festucae is an endophytic fungus that forms a symbiotic association with Lolium perenne. Here we analysed how the metabolome of the ryegrass apoplast changed upon infection of this host with sexual and asexual isolates of E. festucae. A metabolite fingerprinting approach was used to analyse the metabolite composition of apoplastic wash fluid from uninfected and infected L. perenne. Metabolites enriched or depleted in one or both of these treatments were identified using a set of interactive tools. A genetic approach in combination with tandem MS was used to identify a novel product of a secondary metabolite gene cluster. Metabolites likely to be present in the apoplast were identified using MarVis in combination with the BioCyc and KEGG databases, and an in-house Epichloë metabolite database. We were able to identify the known endophyte-specific metabolites, peramine and epichloëcyclins, as well as a large number of unknown markers. To determine whether these methods can be applied to the identification of novel Epichloë-derived metabolites, we deleted a gene encoding a NRPS (lgsA) that is highly expressed in planta. Comparative MS analysis of apoplastic wash fluid from wild-type- vs mutant-infected plants identified a novel Leu/Ile glycoside metabolite present in the former.

Molecular mechanisms in grass-Epichloë interactions: towards endophyte driven farming to improve plant fitness and immunity

All plants harbor many microbial species including bacteria and fungi in their tissues. The interactions between the plant and these microbes could be symbiotic, mutualistic, parasitic or commensalistic. Mutualistic microorganisms are endophytic in nature and are known to play a role in plant growth, development and fitness. Endophytes display complex diversity depending upon the agro-climatic conditions and this diversity could be exploited for crop improvement and sustainable agriculture. Plant-endophyte partnerships are highly specific, several genetic and molecular cascades play a key role in colonization of endophytes in host plants leading to rapid changes in host and endophyte metabolism. This results in the accumulation of secondary metabolites, which play an important role in plant defense against biotic and abiotic stress conditions. Alkaloids are one of the important class of metabolites produced by Epichloë genus and other related classes of endophytes and confer protection against insect and mammalian herbivory. In this context, this review discusses the evolutionary aspects of the Epichloë genus along with key molecular mechanisms determining the lifestyle of Epichloë endophytes in host system. Novel hypothesis is proposed to outline the initial cellular signaling events during colonization of Epichloë in cool season grasses. Complex clustering of alkaloid biosynthetic genes and molecular mechanisms involved in the production of alkaloids have been elaborated in detail. The natural defense and advantages of the endophyte derived metabolites have also been extensively discussed. Finally, this review highlights the importance of endophyte-arbitrated plant immunity to develop novel approaches for eco-friendly agriculture.

Polycomb Repression without Bristles: Facultative Heterochromatin and Genome Stability in Fungi

Genome integrity is essential to maintain cellular function and viability. Consequently, genome instability is frequently associated with dysfunction in cells and associated with plant, animal, and human diseases. One consequence of relaxed genome maintenance that may be less appreciated is an increased potential for rapid adaptation to changing environments in all organisms. Here, we discuss evidence for the control and function of facultative heterochromatin, which is delineated by methylation of histone H3 lysine 27 (H3K27me) in many fungi. Aside from its relatively well understood role in transcriptional repression, accumulating evidence suggests that H3K27 methylation has an important role in controlling the balance between maintenance and generation of novelty in fungal genomes. We present a working model for a minimal repressive network mediated by H3K27 methylation in fungi and outline challenges for future research.

Active and repressed biosynthetic gene clusters have spatially distinct chromosome states

While colocalization within a bacterial operon enables coexpression of the constituent genes, the mechanistic logic of clustering of nonhomologous monocistronic genes in eukaryotes is not immediately obvious. Biosynthetic gene clusters that encode pathways for specialized metabolites are an exception to the classical eukaryote rule of random gene location and provide paradigmatic exemplars with which to understand eukaryotic cluster dynamics and regulation. Here, using 3C, Hi-C, and Capture Hi-C (CHi-C) organ-specific chromosome conformation capture techniques along with high-resolution microscopy, we investigate how chromosome topology relates to transcriptional activity of clustered biosynthetic pathway genes in Arabidopsis thaliana Our analyses reveal that biosynthetic gene clusters are embedded in local hot spots of 3D contacts that segregate cluster regions from the surrounding chromosome environment. The spatial conformation of these cluster-associated domains differs between transcriptionally active and silenced clusters. We further show that silenced clusters associate with heterochromatic chromosomal domains toward the periphery of the nucleus, while transcriptionally active clusters relocate away from the nuclear periphery. Examination of chromosome structure at unrelated clusters in maize, rice, and tomato indicates that integration of clustered pathway genes into distinct topological domains is a common feature in plant genomes. Our results shed light on the potential mechanisms that constrain coexpression within clusters of nonhomologous eukaryotic genes and suggest that gene clustering in the one-dimensional chromosome is accompanied by compartmentalization of the 3D chromosome.

Mycelial biomass and concentration of loline alkaloids driven by complex population structure in Epichloë uncinata and meadow fescue (Schedonorus pratensis)

Many efforts have been made to select and isolate naturally occurring animal-friendly Epichloë strains for later reinfection into elite cultivars. Often this process involves large-scale screening of Epichloë-infected wild grass populations where strains are characterized and alkaloids measured. Here, we describe for the first time the use of genotyping-by-sequencing (GBS) on a collection of 217 Epichloë-infected grasses (7 S. arundinaceum, 4 L. perenne, and 206 S. pratensis). This genotyping strategy is cheaper than complete genome sequencing, is suitable for a large number of individuals, and, when applied to endophyte-infected grasses, conveniently genotypes both organisms. In total, 6273 single nucleotide polymorphisms (SNPs) in the endophyte data set and 38 323 SNPs in the host data set were obtained. Our findings reveal a composite structure with three distinct endophyte clusters unrelated to the three main S. pratensis gene pools that have most likely spread from different glacial refugia in Eurasia. All three gene pools can establish symbiosis with E. uncinata. A comparison of the endophyte clusters with microsatellite-based fingerprinting of the same samples allows a quick test to discriminate between these clusters using two simple sequence repeats (SSRs). Concentrations of loline alkaloids and mycelial biomass are correlated and differ significantly among the plant and endophyte subpopulations; one endophyte strain has higher levels of lolines than others, and one specific host genotype is particularly suitable to host E. uncinata. These findings pave the way for targeted artificial inoculations of specific host-endophyte combinations to boost loline production in the symbiota and for genome association studies with the aim of isolating genes involved in the compatibility between meadow fescue and E. uncinata.

Sterol Demethylation Inhibitor Fungicide Resistance in Leptosphaeria maculans is Caused by Modifications in the Regulatory Region of ERG11

Blackleg is a worldwide disease of canola (Brassica napus), caused by a complex of fungal species in the genus Leptosphaeria, that impacts canola production and seed quality. Demethylation inhibitor (DMI) fungicides that target sterol 14α-demethylase are an integral part of disease control. Here, we report six DMI-resistant isolates of Leptosphaeria maculans and two different types of genetic modification related to the resistance. Analysis of the regulatory region of the DMI target gene ERG11 (also known as CYP51) revealed a 275-bp insertion in two of the isolates and three long terminal repeat retrotransposons (5,263, 5,267, and 5,248 bp) inserted in the promoter region of three resistant isolates. Genetic approaches confirmed that these elements are responsible for DMI resistance in L. maculans and crosses show segregation consistent with a single locus. Reverse-transcription quantitative PCR assays demonstrated that the 275-bp insertion increases ERG11 gene expression, conferring DMI fungicide resistance both in vitro and in planta. Moreover, transformation of a susceptible isolate of L. maculans with ERG11 driven by a promoter containing the 275-bp insertion increased resistance to tebuconazole. A minimal shift of the values of concentration whereby 50% of the mycelial growth is inhibited in vitro was observed in resistant isolates containing long terminal repeat retrotransposons; nevertheless, these isolates were able to develop significant lesions on cotyledons from fungicide-treated seedlings. This is the first report of genetic modifications in L. maculans relating to DMI fungicide resistance.

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome.

Generation of Epichloë Strains Expressing Fluorescent Proteins Suitable for Studying Host-Endophyte Interactions and Characterisation of a T-DNA Integration Event

Methods for the identification and localisation of endophytic fungi are required to study the establishment, development, and progression of host-symbiont interactions, as visible reactions or disease symptoms are generally absent from host plants. Fluorescent proteins have proved valuable as reporter gene products, allowing non-invasive detection in living cells. This study reports the introduction of genes for two fluorescent proteins, green fluorescent protein (GFP) and red fluorescent protein, DsRed, into the genomes of two distinct perennial ryegrass (Lolium perenne L.)-associated Epichloë endophyte strains using A. tumefaciens-mediated transformation. Comprehensive characterisation of reporter gene-containing endophyte strains was performed using molecular genetic, phenotypic, and bioinformatic tools. A combination of long read and short read sequencing of a selected transformant identified a single complex T-DNA insert of 35,530 bp containing multiple T-DNAs linked together. This approach allowed for comprehensive characterisation of T-DNA integration to single-base resolution, while revealing the unanticipated nature of T-DNA integration in the transformant analysed. These reporter gene endophyte strains were able to establish and maintain stable symbiotum with the host. In addition, the same endophyte strain labelled with two different fluorescent proteins were able to cohabit the same plant. This knowledge can be used to provide the basis to develop strategies to gain new insights into the host-endophyte interaction through independent and simultaneous monitoring in planta throughout its life cycle in greater detail.