The RIPper, a web-based tool for genome-wide quantification of Repeat-Induced Point (RIP) mutations

Complexity of the lichen symbiosis revealed by metagenome and transcriptome analysis of Xanthoria parietina

Lichens are composite, symbiotic associations of fungi, algae, and bacteria that result in large, anatomically complex organisms adapted to many of the world's most challenging environments. How such intricate, self-replicating lichen architectures develop from simple microbial components remains unknown because of their recalcitrance to experimental manipulation. Here, we report a metagenomic and metatranscriptomic analysis of the lichen Xanthoria parietina at different developmental stages. We identified 168 genomes of symbionts and lichen-associated microbes across the sampled thalli, including representatives of green algae, three different classes of fungi, and 14 bacterial phyla. By analyzing the occurrence of individual species across lichen thalli from diverse environments, we defined both substrate-specific and core microbial components of the lichen. Metatranscriptomic analysis of the principal fungal symbiont from three different developmental stages of a lichen, compared with axenically grown fungus, revealed differential gene expression profiles indicative of lichen-specific transporter functions, specific cell signaling, transcriptional regulation, and secondary metabolic capacity. Putative immunity-related proteins and lichen-specific structurally conserved secreted proteins resembling fungal pathogen effectors were also identified, consistent with a role for immunity modulation in lichen morphogenesis.

Repeat-induced point mutations driving Parastagonospora nodorum genomic diversity are balanced by selection against non-synonymous mutations

Parastagonospora nodorum is necrotrophic fungal pathogen of wheat with significant genomic resources. Population-level pangenome data for 173 isolates, of which 156 were from Western Australia (WA) and 17 were international, were examined for overall genomic diversity and effector gene content. A heterothallic core population occurred across all regions of WA, with asexually-reproducing clonal clusters in dryer northern regions. High potential for SNP diversity in the form of repeat-induced point mutation (RIP)-like transitions, was observed across the genome, suggesting widespread 'RIP-leakage' from transposon-rich repetitive sequences into non-repetitive regions. The strong potential for RIP-like mutations was balanced by negative selection against non-synonymous SNPs, that was observed within protein-coding regions. Protein isoform profiles of known effector loci (SnToxA, SnTox1, SnTox3, SnTox267, and SnTox5) indicated low-levels of non-synonymous and high-levels of silent RIP-like mutations. Effector predictions identified 186 candidate secreted predicted effector proteins (CSEPs), 69 of which had functional annotations and included confirmed effectors. Pangenome-based effector isoform profiles across WA were distinct from global isolates and were conserved relative to population structure, and may enable new approaches for monitoring crop disease pathotypes.

Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa

The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.

Comparative genomics of the extremophile Cryomyces antarcticus and other psychrophilic Dothideomycetes

Over a billion years of fungal evolution has enabled representatives of this kingdom to populate almost all parts of planet Earth and to adapt to some of its most uninhabitable environments including extremes of temperature, salinity, pH, water, light, or other sources of radiation. Cryomyces antarcticus is an endolithic fungus that inhabits rock outcrops in Antarctica. It survives extremes of cold, humidity and solar radiation in one of the least habitable environments on Earth. This fungus is unusual because it produces heavily melanized, meristematic growth and is thought to be haploid and asexual. Due to its growth in the most extreme environment, it has been suggested as an organism that could survive on Mars. However, the mechanisms it uses to achieve its extremophilic nature are not known. Comparative genomics can provide clues to the processes underlying biological diversity, evolution, and adaptation. This effort has been greatly facilitated by the 1000 Fungal Genomes project and the JGI MycoCosm portal where sequenced genomes have been assembled into phylogenetic and ecological groups representing different projects, lifestyles, ecologies, and evolutionary histories. Comparative genomics within and between these groups provides insights into fungal adaptations, for example to extreme environmental conditions. Here, we analyze two Cryomyces genomes in the context of additional psychrophilic fungi, as well as non-psychrophilic fungi with diverse lifestyles selected from the MycoCosm database. This analysis identifies families of genes that are expanded and contracted in Cryomyces and other psychrophiles and may explain their extremophilic lifestyle. Higher GC contents of genes and of bases in the third positions of codons may help to stabilize DNA under extreme conditions. Numerous smaller contigs in C. antarcticus suggest the presence of an alternative haplotype that could indicate the sequenced isolate is diploid or dikaryotic. These analyses provide a first step to unraveling the secrets of the extreme lifestyle of C. antarcticus.

Update on the Basic Understanding of Fusarium graminearum Virulence Factors in Common Wheat Research

Wheat is one of the most important food crops, both in China and worldwide. Wheat production is facing extreme stresses posed by different diseases, including Fusarium head blight (FHB), which has recently become an increasingly serious concerns. FHB is one of the most significant and destructive diseases affecting wheat crops all over the world. Recent advancements in genomic tools provide a new avenue for the study of virulence factors in relation to the host plants. The current review focuses on recent progress in the study of different strains of Fusarium infection. The presence of genome-wide repeat-induced point (RIP) mutations causes genomic mutations, eventually leading to host plant susceptibility against Fusarium invasion. Furthermore, effector proteins disrupt the host plant resistance mechanism. In this study, we proposed systematic modification of the host genome using modern biological tools to facilitate plant resistance against foreign invasion. We also suggested a number of scientific strategies, such as gene cloning, developing more powerful functional markers, and using haplotype marker-assisted selection, to further improve FHB resistance and associated breeding methods.

The massive 340 megabase genome of Anisogramma anomala, a biotrophic ascomycete that causes eastern filbert blight of hazelnut

BACKGROUND

The ascomycete fungus Anisogramma anomala causes Eastern Filbert Blight (EFB) on hazelnut (Corylus spp.) trees. It is a minor disease on its native host, the American hazelnut (C. americana), but is highly destructive on the commercially important European hazelnut (C. avellana). In North America, EFB has historically limited commercial production of hazelnut to west of the Rocky Mountains. A. anomala is an obligately biotrophic fungus that has not been grown in continuous culture, rendering its study challenging. There is a 15-month latency before symptoms appear on infected hazelnut trees, and only a sexual reproductive stage has been observed. Here we report the sequencing, annotation, and characterization of its genome.

RESULTS

The genome of A. anomala was assembled into 108 scaffolds totaling 342,498,352 nt with a GC content of 34.46%. Scaffold N50 was 33.3 Mb and L50 was 5. Nineteen scaffolds with lengths over 1 Mb constituted 99% of the assembly. Telomere sequences were identified on both ends of two scaffolds and on one end of another 10 scaffolds. Flow cytometry estimated the genome size of A. anomala at 370 Mb. The genome exhibits two-speed evolution, with 93% of the assembly as AT-rich regions (32.9% GC) and the other 7% as GC-rich (57.1% GC). The AT-rich regions consist predominantly of repeats with low gene content, while 90% of predicted protein coding genes were identified in GC-rich regions. Copia-like retrotransposons accounted for more than half of the genome. Evidence of repeat-induced point mutation (RIP) was identified throughout the AT-rich regions, and two copies of the rid gene and one of dim-2, the key genes in the RIP mutation pathway, were identified in the genome. Consistent with its homothallic sexual reproduction cycle, both MAT1-1 and MAT1-2 idiomorphs were found. We identified a large suite of genes likely involved in pathogenicity, including 614 carbohydrate active enzymes, 762 secreted proteins and 165 effectors.

CONCLUSIONS

This study reveals the genomic structure, composition, and putative gene function of the important pathogen A. anomala. It provides insight into the molecular basis of the pathogen's life cycle and a solid foundation for studying EFB.

The genome of Salmacisia buchloëana, the parasitic puppet master pulling strings of sexual phenotypic monstrosities in buffalograss

To complete its parasitic lifecycle, Salmacisia buchloëana, a biotrophic fungus, manipulates reproductive organ development, meristem determinacy, and resource allocation in its dioecious plant host, buffalograss (Bouteloua dactyloides; Poaceae). To gain insight into S. buchloëana's ability to manipulate its host, we sequenced and assembled the 20.1 Mb genome of S. buchloëana into 22 chromosome-level pseudomolecules. Phylogenetic analysis suggests that S. buchloëana is nested within the genus Tilletia and diverged from Tilletia caries and Tilletia walkeri ∼40 MYA. We find that S. buchloëana contains a novel chromosome arm with no syntenic relationship to other publicly available Tilletia genomes, and that genes on the novel arm are upregulated upon infection, suggesting that this unique chromosomal segment may have played a critical role in S. buchloëana's evolution and host specificity. Salmacisia buchloëana has one of the largest fractions of serine peptidases (1.53% of the proteome) and one of the highest GC contents (62.3%) in all classified fungi. Analysis of codon base composition indicated that GC content is controlled more by selective constraints than directional mutation, and that S. buchloëana has a unique bias for the serine codon UCG. Finally, we identify 3 inteins within the S. buchloëana genome, 2 of which are located in a gene often used in fungal taxonomy. The genomic and transcriptomic resources generated here will aid plant pathologists and breeders by providing insight into the extracellular components contributing to sex determination in dioecious grasses.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Chromosome-Level Assemblies for the Pine Pitch Canker Pathogen Fusarium circinatum

The pine pitch canker pathogen, Fusarium circinatum, is globally regarded as one of the most important threats to commercial pine-based forestry. Although genome sequences of this fungus are available, these remain highly fragmented or structurally ill-defined. Our overall goal was to provide high-quality assemblies for two notable strains of F. circinatum, and to characterize these in terms of coding content, repetitiveness and the position of telomeres and centromeres. For this purpose, we used Oxford Nanopore Technologies MinION long-read sequences, as well as Illumina short sequence reads. By leveraging the genomic synteny inherent to F. circinatum and its close relatives, these sequence reads were assembled to chromosome level, where contiguous sequences mostly spanned from telomere to telomere. Comparative analyses unveiled remarkable variability in the twelfth and smallest chromosome, which is known to be dispensable. It presented a striking length polymorphism, with one strain lacking substantial portions from the chromosome's distal and proximal regions. These regions, characterized by a lower gene density, G+C content and an increased prevalence of repetitive elements, contrast starkly with the syntenic segments of the chromosome, as well as with the core chromosomes. We propose that these unusual regions might have arisen or expanded due to the presence of transposable elements. A comparison of the overall chromosome structure revealed that centromeric elements often underpin intrachromosomal differences between F. circinatum strains, especially at chromosomal breakpoints. This suggests a potential role for centromeres in shaping the chromosomal architecture of F. circinatum and its relatives. The publicly available genome data generated here, together with the detailed metadata provided, represent essential resources for future studies of this important plant pathogen.

GWAS reveals a rapidly evolving candidate avirulence effector in the Cercospora leaf spot pathogen

The major resistance gene BvCR4 recently bred into sugar beet hybrids provides a high level of resistance to Cercospora leaf spot caused by the fungal pathogen Cercospora beticola. The occurrence of pathogen strains that overcome BvCR4 was studied using field trials in Switzerland conducted under natural disease pressure. Virulence of a subset of these strains was evaluated in a field trial conducted under elevated artificial disease pressure. We created a new C. beticola reference genome and mapped whole genome sequences of 256 isolates collected in Switzerland and Germany. These were combined with virulence phenotypes to conduct three separate genome-wide association studies (GWAS) to identify candidate avirulence genes. We identified a locus associated with avirulence containing a putative avirulence effector gene named AvrCR4. All virulent isolates either lacked AvrCR4 or had nonsynonymous mutations within the gene. AvrCR4 was present in all 74 isolates from non-BvCR4 hybrids, whereas 33 of 89 isolates from BvCR4 hybrids carried a deletion. We also mapped genomic data from 190 publicly available US isolates to our new reference genome. The AvrCR4 deletion was found in only one of 95 unique isolates from non-BvCR4 hybrids in the United States. AvrCR4 presents a unique example of an avirulence effector in which virulent alleles have only recently emerged. Most likely these were selected out of standing genetic variation after deployment of BvCR4. Identification of AvrCR4 will enable real-time screening of C. beticola populations for the emergence and spread of virulent isolates.

Identification and specificity validation of unique and antimicrobial resistance genes to trace suspected pathogenic AMR bacteria and to monitor the development of AMR in non-AMR strains in the environment and clinical settings

The detection of developing antimicrobial resistance (AMR) has become a global issue. The detection of developing antimicrobial resistance has become a global issue. The growing number of AMR bacteria poses a new threat to public health. Therefore, a less laborious and quick confirmatory test becomes important for further investigations into developing AMR in the environment and in clinical settings. This study aims to present a comprehensive analysis and validation of unique and antimicrobial-resistant strains from the WHO priority list of antimicrobial-resistant bacteria and previously reported AMR strains such as Acinetobacter baumannii, Aeromonas spp., Anaeromonas frigoriresistens, Anaeromonas gelatinfytica, Bacillus spp., Campylobacter jejuni subsp. jejuni, Enterococcus faecalis, Escherichia coli, Haemophilus influenzae, Helicobacter pylori, Klebsiella pneumonia subsp. pneumoniae, Pseudomonas aeruginosa, Salmonella enterica subsp. enterica serovar Typhimurium, Thermanaeromonas toyohensis, and Vibrio proteolyticus. Using in-house designed gene-specific primers, 18 different antibiotic resistance genes (algJ, alpB, AQU-1, CEPH-A3, ciaB, CMY-1-MOX-7, CMY-1-MOX-9, CMY-1/MOX, cphA2, cphA5, cphA7, ebpA, ECP_4655, fliC, OXA-51, RfbU, ThiU2, and tolB) from 46 strains were selected and validated. Hence, this study provides insight into the identification of strain-specific, unique antimicrobial resistance genes. Targeted amplification and verification using selected unique marker genes have been reported. Thus, the present detection and validation use a robust method for the entire experiment. Results also highlight the presence of another set of 18 antibiotic-resistant and unique genes (Aqu1, cphA2, cphA3, cphA5, cphA7, cmy1/mox7, cmy1/mox9, asaI, ascV, asoB, oxa-12, acr-2, pepA, uo65, pliI, dr0274, tapY2, and cpeT). Of these sets of genes, 15 were found to be suitable for the detection of pathogenic strains belonging to the genera Aeromonas, Pseudomonas, Helicobacter, Campylobacter, Enterococcus, Klebsiella, Acinetobacter, Salmonella, Haemophilus, and Bacillus. Thus, we have detected and verified sets of unique and antimicrobial resistance genes in bacteria on the WHO Priority List and from published reports on AMR bacteria. This study offers advantages for confirming antimicrobial resistance in all suspected AMR bacteria and monitoring the development of AMR in non-AMR bacteria, in the environment, and in clinical settings.

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

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

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

CRISPR-Cas9 Gene Editing and Secondary Metabolite Screening Confirm Fusarium graminearum C16 Biosynthetic Gene Cluster Products as Decalin-Containing Diterpenoid Pyrones

Fusarium graminearum is a causal organism of Fusarium head blight in cereals and maize. Although a few secondary metabolites produced by F. graminearum are considered disease virulence factors, many molecular products of biosynthetic gene clusters expressed by F. graminearum during infection and their associated role in the disease are unknown. In particular, the predicted meroterpenoid products of the biosynthetic gene cluster historically designated as "C16" are likely associated with pathogenicity. Presented here are the results of CRISPR-Cas9 gene-editing experiments disrupting the polyketide synthase and terpene synthase genes associated with the C16 biosynthetic gene cluster in F. graminearum. Culture medium screening experiments using transformant strains were profiled by UHPLC-HRMS and targeted MS² experiments to confirm the associated secondary metabolite products of the C16 biosynthetic gene cluster as the decalin-containing diterpenoid pyrones, FDDP-D and FDDP-E. Both decalin-containing diterpenoid pyrones were confirmed to be produced in wheat heads challenged with F. graminearum in growth chamber trials. The extent to which the F. graminearum C16 biosynthetic gene cluster is dispersed within the genus Fusarium is discussed along with a proposed role of the FDDPs as pathogen virulence factors.

Chromosome-level analysis of the Colletotrichum graminicola genome reveals the unique characteristics of core and minichromosomes

The fungal pathogen Colletotrichum graminicola causes the anthracnose of maize (Zea mays) and is responsible for significant yield losses worldwide. The genome of C. graminicola was sequenced in 2012 using Sanger sequencing, 454 pyrosequencing, and an optical map to obtain an assembly of 13 pseudochromosomes. We re-sequenced the genome using a combination of short-read (Illumina) and long-read (PacBio) technologies to obtain a chromosome-level assembly. The new version of the genome sequence has 13 chromosomes with a total length of 57.43 Mb. We detected 66 (23.62 Mb) structural rearrangements in the new assembly with respect to the previous version, consisting of 61 (21.98 Mb) translocations, 1 (1.41 Mb) inversion, and 4 (221 Kb) duplications. We annotated the genome and obtained 15,118 predicted genes and 3,614 new gene models compared to the previous version of the assembly. We show that 25.88% of the new assembly is composed of repetitive DNA elements (13.68% more than the previous assembly version), which are mostly found in gene-sparse regions. We describe genomic compartmentalization consisting of repeat-rich and gene-poor regions vs. repeat-poor and gene-rich regions. A total of 1,140 secreted proteins were found mainly in repeat-rich regions. We also found that ~75% of the three smallest chromosomes (minichromosomes, between 730 and 551 Kb) are strongly affected by repeat-induced point mutation (RIP) compared with 28% of the larger chromosomes. The gene content of the minichromosomes (MCs) comprises 121 genes, of which 83.6% are hypothetical proteins with no predicted function, while the mean percentage of Chr1-Chr10 is 36.5%. No predicted secreted proteins are present in the MCs. Interestingly, only 2% of the genes in Chr11 have homologs in other strains of C. graminicola, while Chr12 and 13 have 58 and 57%, respectively, raising the question as to whether Chrs12 and 13 are dispensable. The core chromosomes (Chr1-Chr10) are very different with respect to the MCs (Chr11-Chr13) in terms of the content and sequence features. We hypothesize that the higher density of repetitive elements and RIPs in the MCs may be linked to the adaptation and/or host co-evolution of this pathogenic fungus.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

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.

Stepwise recombination suppression around the mating-type locus in an ascomycete fungus with self-fertile spores

Recombination is often suppressed at sex-determining loci in plants and animals, and at self-incompatibility or mating-type loci in plants and fungi. In fungal ascomycetes, recombination suppression around the mating-type locus is associated with pseudo-homothallism, i.e. the production of self-fertile dikaryotic sexual spores carrying the two opposite mating types. This has been well studied in two species complexes from different families of Sordariales: Podospora anserina and Neurospora tetrasperma. However, it is unclear whether this intriguing association holds in other species. We show here that Schizothecium tetrasporum, a fungus from a third family in the order Sordariales, also produces mostly self-fertile dikaryotic spores carrying the two opposite mating types. This was due to a high frequency of second meiotic division segregation at the mating-type locus, indicating the occurrence of a single and systematic crossing-over event between the mating-type locus and the centromere, as in P. anserina. The mating-type locus has the typical Sordariales organization, plus a MAT1-1-1 pseudogene in the MAT1-2 haplotype. High-quality genome assemblies of opposite mating types and segregation analyses revealed a suppression of recombination in a region of 1.47 Mb around the mating-type locus. We detected three evolutionary strata, indicating a stepwise extension of recombination suppression. The three strata displayed no rearrangement or transposable element accumulation but gene losses and gene disruptions were present, and precisely at the strata margins. Our findings indicate a convergent evolution of self-fertile dikaryotic sexual spores across multiple ascomycete fungi. The particular pattern of meiotic segregation at the mating-type locus was associated with recombination suppression around this locus, that had extended stepwise. This association between pseudo-homothallism and recombination suppression across lineages and the presence of gene disruption at the strata limits are consistent with a recently proposed mechanism of sheltering deleterious alleles to explain stepwise recombination suppression.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

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

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

Publicly Available and Validated DNA Reference Sequences Are Critical to Fungal Identification and Global Plant Protection Efforts: A Use-Case in Colletotrichum

Publicly available and validated DNA reference sequences useful for phylogeny estimation and identification of fungal pathogens are an increasingly important resource in the efforts of plant protection organizations to facilitate safe international trade of agricultural commodities. Colletotrichum species are among the most frequently encountered and regulated plant pathogens at U.S. ports-of-entry. The RefSeq Targeted Loci (RTL) project at NCBI (BioProject no. PRJNA177353) contains a database of curated fungal internal transcribed spacer (ITS) sequences that interact extensively with NCBI Taxonomy, resulting in verified name-strain-sequence type associations for >12,000 species. We present a publicly available dataset of verified and curated name-type strain-sequence associations for all available Colletotrichum species. This includes an updated GenBank Taxonomy for 238 species associated with up to 11 protein coding loci and an updated RTL ITS dataset for 226 species. We demonstrate that several marker loci are well suited for phylogenetic inference and identification. We improve understanding of phylogenetic relationships among verified species, verify or improve phylogenetic circumscriptions of 14 species complexes, and reveal that determining relationships among these major clades will require additional data. We present detailed comparisons between phylogenetic and similarity-based approaches to species identification, revealing complex patterns among single marker loci that often lead to misidentification when based on single-locus similarity approaches. We also demonstrate that species-level identification is elusive for a subset of samples regardless of analytical approach, which may be explained by novel species diversity in our dataset and incomplete lineage sorting and lack of accumulated synapomorphies at these loci.

Comparative Genomics of Fusarium circinatum Isolates Used to Screen Southern Pines for Pitch Canker Resistance

Pitch canker, caused by the fungal pathogen Fusarium circinatum, is a global disease affecting many Pinus spp. Often fatal, this disease causes significant mortality in both commercially grown and natural pine forests and is an issue of current and growing concern. F. circinatum isolates collected from three locations in the U.S. state of Florida were shown to be virulent on both slash and loblolly pine, with two of the isolates causing equivalent and significantly larger lesions than those caused by the third isolate during pathogenicity trials. In addition, significant genetic variation in lesion length in the pedigreed slash pine population was evident and rankings of parents for lesion length were similar across isolates. Experimental data demonstrate that both host and pathogen genetics contribute to disease severity. High-quality genomic assemblies of all three isolates were created and compared for structural differences and gene content. No major structural differences were observed among the isolates; however, missing or altered genes do contribute to genomic variation in the pathogen population. This work evaluates in planta virulence among three isolates of F. circinatum, provides genomic resources to facilitate study of this organism, and details comparative genomic methods that may be used to explore the pathogen's contribution to disease development.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

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

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

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

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

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.

Histone H3K27 methylation perturbs transcriptional robustness and underpins dispensability of highly conserved genes in fungi

Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 underpins the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation, we used a framework of 23 species across the Fusarium genus including three species covered with histone methylation data. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. In contrast, genes marked by H3K27me3 do not show such associations. Furthermore, highly conserved genes marked by H3K27me3 encode smaller proteins, exhibit weaker codon usage bias, higher levels of hydrophobicity, show lower intrinsically disordered regions, and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival during host infection. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related Fusarium fungi are associated with transcriptional perturbation at the microevolutionary scale. In conclusion, we show how repressive histone marks are entangled in the evolutionary fate of highly conserved genes across evolutionary timescales.

Unique patterns of mating pheromone presence and absence could result in the ambiguous sexual behaviors of Colletotrichum species

Colletotrichum species are known to engage in unique sexual behaviors that differ significantly from the mating strategies of other filamentous ascomycete species. For example, most ascomycete fungi require the expression of both the MAT1-1-1 and MAT1-2-1 genes to induce sexual reproduction. In contrast, all isolates of Colletotrichum harbor only the MAT1-2-1 gene and yet, are capable of recognizing suitable mating partners and producing sexual progeny. The molecular mechanisms contributing to mating types and behaviors in Colletotrichum are, however, unknown. A comparative genomics approach analyzing 35 genomes, representing 31 Colletotrichum species and two Verticillium species, was used to elucidate a putative molecular mechanism underlying the unique sexual behaviors observed in Colletotrichum species. The existence of only the MAT1-2 idiomorph was confirmed across all species included in this study. Comparisons of the loci harboring the two mating pheromones and their cognate receptors revealed interesting patterns of gene presence and absence. The results showed that these genes have been lost multiple, independent times over the evolutionary history of this genus. These losses indicate that the pheromone pathway no longer plays an active role in mating type determination, suggesting an undiscovered mechanism by which mating partner recognition is controlled in these species. This further suggests that there has been a redirection of the underlying genetic mechanisms that regulate sexual development in Colletotrichum species. This research thus provides a foundation from which further interrogation of this topic can take place.

Remarkable recent changes in the genetic diversity of the avirulence gene AvrStb6 in global populations of the wheat pathogen Zymoseptoria tritici

Septoria tritici blotch (STB), caused by the fungus Zymoseptoria tritici, is one of the most economically important diseases of wheat. Recently, both factors of a gene-for-gene interaction between Z. tritici and wheat, the wheat receptor-like kinase Stb6 and the Z. tritici secreted effector protein AvrStb6, have been identified. Previous analyses revealed a high diversity of AvrStb6 haplotypes present in earlier Z. tritici isolate collections, with up to c.18% of analysed isolates possessing the avirulence isoform of AvrStb6 identical to that originally identified in the reference isolate IPO323. With Stb6 present in many commercial wheat cultivars globally, we aimed to assess potential changes in AvrStb6 genetic diversity and the incidence of haplotypes allowing evasion of Stb6-mediated resistance in more recent Z. tritici populations. Here we show, using targeted resequencing of AvrStb6, that this gene is universally present in field isolates sampled from major wheat-growing regions of the world in 2013-2017. However, in contrast to the data from previous AvrStb6 population studies, we report a complete absence of the originally described avirulence isoform of AvrStb6 amongst modern Z. tritici isolates. Moreover, a remarkably small number of haplotypes, each encoding AvrStb6 protein isoforms conditioning virulence on Stb6-containing wheat, were found to predominate among modern Z. tritici isolates. A single virulence isoform of AvrStb6 was found to be particularly abundant throughout the global population. These findings indicate that, despite the ability of Z. tritici to sexually reproduce on resistant hosts, AvrStb6 avirulence haplotypes tend to be eliminated in subsequent populations.

Increased abundance of secreted hydrolytic enzymes and secondary metabolite gene clusters define the genomes of latent plant pathogens in the Botryosphaeriaceae

BACKGROUND

The Botryosphaeriaceae are important plant pathogens, but also have the ability to establish asymptomatic infections that persist for extended periods in a latent state. In this study, we used comparative genome analyses to shed light on the genetic basis of the interactions of these fungi with their plant hosts. For this purpose, we characterised secreted hydrolytic enzymes, secondary metabolite biosynthetic gene clusters and general trends in genomic architecture using all available Botryosphaeriaceae genomes, and selected Dothideomycetes genomes.

RESULTS

The Botryosphaeriaceae genomes were rich in carbohydrate-active enzymes (CAZymes), proteases, lipases and secondary metabolic biosynthetic gene clusters (BGCs) compared to other Dothideomycete genomes. The genomes of Botryosphaeria, Macrophomina, Lasiodiplodia and Neofusicoccum, in particular, had gene expansions of the major constituents of the secretome, notably CAZymes involved in plant cell wall degradation. The Botryosphaeriaceae genomes were shown to have moderate to high GC contents and most had low levels of repetitive DNA. The genomes were not compartmentalized based on gene and repeat densities, but genes of secreted enzymes were slightly more abundant in gene-sparse regions.

CONCLUSION

The abundance of secreted hydrolytic enzymes and secondary metabolite BGCs in the genomes of Botryosphaeria, Macrophomina, Lasiodiplodia, and Neofusicoccum were similar to those in necrotrophic plant pathogens and some endophytes of woody plants. The results provide a foundation for comparative genomic analyses and hypotheses to explore the mechanisms underlying Botryosphaeriaceae host-plant interactions.

Apicidin biosynthesis is linked to accessory chromosomes in Fusarium poae isolates

BACKGROUND

Fusarium head blight is a disease of global concern that reduces crop yields and renders grains unfit for consumption due to mycotoxin contamination. Fusarium poae is frequently associated with cereal crops showing symptoms of Fusarium head blight. While previous studies have shown F. poae isolates produce a range of known mycotoxins, including type A and B trichothecenes, fusarins and beauvericin, genomic analysis suggests that this species may have lineage-specific accessory chromosomes with secondary metabolite biosynthetic gene clusters awaiting description.

METHODS

We examined the biosynthetic potential of 38 F. poae isolates from Eastern Canada using a combination of long-read and short-read genome sequencing and untargeted, high resolution mass spectrometry metabolome analysis of extracts from isolates cultured in multiple media conditions.

RESULTS

A high-quality assembly of isolate DAOMC 252244 (Fp157) contained four core chromosomes as well as seven additional contigs with traits associated with accessory chromosomes. One of the predicted accessory contigs harbours a functional biosynthetic gene cluster containing homologs of all genes associated with the production of apicidins. Metabolomic and genomic analyses confirm apicidins are produced in 4 of the 38 isolates investigated and genomic PCR screening detected the apicidin synthetase gene APS1 in approximately 7% of Eastern Canadian isolates surveyed.

CONCLUSIONS

Apicidin biosynthesis is linked to isolate-specific putative accessory chromosomes in F. poae. The data produced here are an important resource for furthering our understanding of accessory chromosome evolution and the biosynthetic potential of F. poae.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

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.

Historical genomics reveals the evolutionary mechanisms behind multiple outbreaks of the host-specific coffee wilt pathogen Fusarium xylarioides

BACKGROUND

Nearly 50% of crop yields are lost to pests and disease, with plants and pathogens locked in an amplified co-evolutionary process of disease outbreaks. Coffee wilt disease, caused by Fusarium xylarioides, decimated coffee production in west and central Africa following its initial outbreak in the 1920s. After successful management, it later re-emerged and by the 2000s comprised two separate epidemics on arabica coffee in Ethiopia and robusta coffee in east and central Africa.

RESULTS

Here, we use genome sequencing of six historical culture collection strains spanning 52 years to identify the evolutionary processes behind these repeated outbreaks. Phylogenomic reconstruction using 13,782 single copy orthologs shows that the robusta population arose from the initial outbreak, whilst the arabica population is a divergent sister clade to the other strains. A screen for putative effector genes involved in pathogenesis shows that the populations have diverged in gene content and sequence mainly by vertical processes within lineages. However, 15 putative effector genes show evidence of horizontal acquisition, with close homology to genes from F. oxysporum. Most occupy small regions of homology within wider scaffolds, whereas a cluster of four genes occupy a 20Kb scaffold with strong homology to a region on a mobile pathogenicity chromosome in F. oxysporum that houses known effector genes. Lacking a match to the whole mobile chromosome, we nonetheless found close associations with DNA transposons, especially the miniature impala type previously proposed to facilitate horizontal transfer of pathogenicity genes in F. oxysporum. These findings support a working hypothesis that the arabica and robusta populations partly acquired distinct effector genes via transposition-mediated horizontal transfer from F. oxysporum, which shares coffee as a host and lives on other plants intercropped with coffee.

CONCLUSION

Our results show how historical genomics can help reveal mechanisms that allow fungal pathogens to keep pace with our efforts to resist them. Our list of putative effector genes identifies possible future targets for fungal control. In turn, knowledge of horizontal transfer mechanisms and putative donor taxa might help to design future intercropping strategies that minimize the risk of transfer of effector genes between closely-related Fusarium taxa.

Rapid sequence evolution driven by transposable elements at a virulence locus in a fungal wheat pathogen

BACKGROUND

Plant pathogens cause substantial crop losses in agriculture production and threaten food security. Plants evolved the ability to recognize virulence factors and pathogens have repeatedly escaped recognition due rapid evolutionary change at pathogen virulence loci (i.e. effector genes). The presence of transposable elements (TEs) in close physical proximity of effector genes can have important consequences for gene regulation and sequence evolution. Species-wide investigations of effector gene loci remain rare hindering our ability to predict pathogen evolvability.

RESULTS

Here, we performed genome-wide association studies (GWAS) on a highly polymorphic mapping population of 120 isolates of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe. We identified a major locus underlying significant variation in reproductive success of the pathogen and damage caused on the wheat cultivar Claro. The most strongly associated locus is intergenic and flanked by genes encoding a predicted effector and a serine-type endopeptidase. The center of the locus contained a highly dynamic region consisting of multiple families of TEs. Based on a large global collection of assembled genomes, we show that the virulence locus has undergone substantial recent sequence evolution. Large insertion and deletion events generated length variation between the flanking genes by a factor of seven (5-35 kb). The locus showed also strong signatures of genomic defenses against TEs (i.e. RIP) contributing to the rapid diversification of the locus.

CONCLUSIONS

In conjunction, our work highlights the power of combining GWAS and population-scale genome analyses to investigate major effect loci in pathogens.

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

Genome-Wide Analyses of Repeat-Induced Point Mutations in the Ascomycota

The Repeat-Induced Point (RIP) mutation pathway is a fungus-specific genome defense mechanism that mitigates the deleterious consequences of repeated genomic regions and transposable elements (TEs). RIP mutates targeted sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence and extent of RIP with a sliding-window approach. Using genome-wide RIP data and two sets of control groups, the association between RIP, TEs, and GC content were contrasted in organisms capable and incapable of RIP. Based on these data, we then set out to determine the extent and occurrence of RIP in 58 representatives of the Ascomycota. The findings were summarized by placing each of the fungi investigated in one of six categories based on the extent of genome-wide RIP. In silico RIP analyses, using a sliding-window approach with stringent RIP parameters, implemented simultaneously within the same genetic context, on high quality genome assemblies, yielded superior results in determining the genome-wide RIP among the Ascomycota. Most Ascomycota had RIP and these mutations were particularly widespread among classes of the Pezizomycotina, including the early diverging Orbiliomycetes and the Pezizomycetes. The most extreme cases of RIP were limited to representatives of the Dothideomycetes and Sordariomycetes. By contrast, the genomes of the Taphrinomycotina and Saccharomycotina contained no detectable evidence of RIP. Also, recent losses in RIP combined with controlled TE proliferation in the Pezizomycotina subphyla may promote substantial genome enlargement as well as the formation of sub-genomic compartments. These findings have broadened our understanding of the taxonomic range and extent of RIP in Ascomycota and how this pathway affects the genomes of fungi harboring it.

Genome comparisons suggest an association between Ceratocystis host adaptations and effector clusters in unique transposable element families

Ceratocystis fimbriata is a host specific fungal pathogen of sweet potato (Ipomoea batatas). The closely related species, C. manginecans, is an important pathogen of trees (e.g. Acacia mangium and Mangifera indica) but has never been isolated from tuber crops. The genetic factors that determine the host range and host specificity of these species have not been determined. The aim of this study was to compare the genomes of C. fimbriata and C. manginecans in order to identify species-specific genetic differences that could be associated with host specificity. This included whole-genome alignments as well as comparisons of gene content and transposable elements (TEs). The genomes of the two species were found to be very similar, sharing similar catalogues of CAZymes, peptidases and lipases. However, the genomes of the two species also varied, harbouring species-specific genes (e.g. small secreted effectors, nutrient processing proteins and stress response proteins). A portion of the TEs identified (17%) had a unique distribution in each species. Transposable elements appeared to have played a prominent role in the divergence of the two species because they were strongly associated with chromosomal translocations and inversions as well as with unique genomic regions containing species-specific genes. Two large effector clusters, with unique TEs in each species, were identified. These effectors displayed non-synonymous mutations and deletions, conserved within a species, and could serve as mutational hot-spots for the development of host specificity in the two species.

Repeat-Induced Point Mutations Drive Divergence between Fusarium circinatum and Its Close Relatives

The Repeat-Induced Point (RIP) mutation pathway is a fungal-specific genome defense mechanism that counteracts the deleterious effects of transposable elements. This pathway permanently mutates its target sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence of RIP in the pitch canker pathogen, Fusarium circinatum, and its close relatives in the Fusarium fujikuroi species complex (FFSC). Our results showed that the examined fungi all exhibited hallmarks of RIP, but that they differed in terms of the extent to which their genomes were affected by this pathway. RIP mutations constituted a large proportion of all the FFSC genomes, including both core and dispensable chromosomes, although the latter were generally more extensively affected by RIP. Large RIP-affected genomic regions were also much more gene sparse than the rest of the genome. Our data further showed that RIP-directed sequence diversification increased the variability between homologous regions of related species, and that RIP-affected regions can interfere with homologous recombination during meiosis, thereby contributing to post-mating segregation distortion. Taken together, these findings suggest that RIP can drive the independent divergence of chromosomes, alter chromosome architecture, and contribute to the divergence among F. circinatum and other members of this economically important group of fungi.