OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes

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

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

IMA GENOME - F20 A draft genome assembly of Agroatheliarolfsii, Ceratobasidiumpapillatum, Pyrenopezizabrassicae, Neopestalotiopsismacadamiae, Sphaerellopsisfilum and genomic resources for Colletotrichumspaethianum and Colletotrichumfructicola

This is a genome announacment there is no abstract

The plastomes of Cypripedium (Orchidaceae: Cypripedioideae) exhibit atypical GC content and genome size based on different sequencing strategies

Most of the sequenced plastomes of photosynthetic angiosperms exhibit conservation in size, gene content, gene order, and GC content. In contrast, the sequenced plastomes of Cypripedium are distinguished by genome size expansion, AT-biased base composition, structural variation, and a low substitution rate. Additionally, the impact of sequencing methods is seldom addressed in prior studies, and the species represented in these studies are underrepresented. These atypical plastome features render the genus an ideal candidate for investigating plastome evolution. Besides, the backbone relationships within the genus remain poorly resolved. In this study, we sequenced twelve Cypripedium plastomes using three distinct sequencing strategies and obtained an additional 27 sequences from GenBank for comparative analysis. We classified the plastomes of the genus into two types: one resembling those of most other angiosperms, and the other characterized by inverted repeat (IR) expansion and small single copy (SSC) contraction. The plastomes within this genus exhibit significant size variations (∼72 kb), variations in GC content, and structural differences at the genus level. Furthermore, our comparative analysis revealed that the choice of sequencing strategy significantly impacts the assembly results. The uncovered regions in samples sequenced with short-read technology are predominantly AT-rich, suggesting that short-read sequencing may lead to assembly errors in plastomes with AT-rich regions and long repeats. Additionally, we have reconstructed the phylogeny of the genus using plastome-level data. However, the phylogenetic relationships within the genus remain partially solved. This study provides new insights into the evolution of plastomes, particularly those with AT-rich base compositions and genomes containing long repeat regions.

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.

Population-level whole-genome sequencing of Ascochyta rabiei identifies genomic loci associated with isolate aggressiveness

Ascochyta blight caused by the ascomycete Ascochyta rabiei poses a major biotic threat to chickpea (Cicer arietinum) industries worldwide and incurs substantial costs to the Australian multimillion-dollar chickpea industry in both disease control and yield loss. The fungus was introduced to Australia in the 1970s from an unknown source population and, within a few decades, successfully established in all Australian agroecological chickpea-growing regions. Although genetically highly clonal, a broad range of phenotypic variation in terms of aggressiveness exists among the Australian A. rabiei isolates. More recently, highly aggressive isolates capable of causing severe disease symptoms on moderate to highly resistant chickpea cultivars have increased in frequency. To identify genetic loci potentially associated with A. rabiei aggressiveness on Australian chickpea cultivars, we performed deep genome sequencing of 230 isolates collected from a range of agroecological chickpea-growing regions between 2013 and 2020. Population genetic analyses using genome-wide SNP data identified three main clusters of genetically closely related isolates in Australia. Phylogenetic analyses showed that highly aggressive phenotypes developed multiple times independently throughout the phylogeny. The results point to a minor contribution of multiple genetic regions and most likely epigenomic variations to aggressiveness of A. rabiei isolates on Australian chickpea cultivars.

Genome-wide survey of the bipartite structure and pathogenesis-related genes of Neostagonosporella sichuanensis, a causal agent of Fishscale bamboo rhombic-spot disease

Bamboo resources have garnered significant global attention due to their excellent capacity for regeneration and high yield. Rhombic-spot disease, a substantial threat to fishscale bamboo (Phyllostachys heteroclada), is primarily caused by Neostagonosporella sichuanensis. This study first reported the genome assemblies and characteristics of two N. sichuanensis isolates using PacBio and Illumina sequencing platforms. The genomes of N. sichuanensis strain SICAUCC 16-0001 and strain SICAUCC 23-0140, with sizes of 48.0 Mb and 48.4 Mb, respectively, revealed 10,289 and 10,313 protein-coding genes. Additionally, they contained 34.99 and 34.46% repetitive sequences within AT-rich regions, with notable repeat-induced point mutation activity. Comparative genome analysis identified 1,049 contracted and 45 expanded gene families in the genome of N. sichuanensis, including several related to pathogenicity. Several gene families involved in mycotoxin metabolism, secondary metabolism, sterol biosynthesis and transport, and cell wall degradation were contracted. Compared to most analyzed necrotrophic, hemibiotrophic, and phaeosphaeriacous pathogens, the genomes of two N. sichuanensis isolates exhibited fewer secondary metabolite enzymes, carbohydrate-active enzymes, plant cell wall degrading enzymes, secreted proteins, and effectors. Comparative genomics analysis suggested that N. sichuanensis shares more similar characteristics with hemibiotrophic pathogens. Based on single carbon source tests, N. sichuanensis strains demonstrated a higher potential for xylan decomposition than pectin and cellulose. The proportion of cell wall-degrading enzyme effectors occupied a high proportion of the total effectors of the N. sichuanensis genomes. These findings provide valuable insights into uncovering the pathogenesis of N. sichuanensis toward the efficient management of rhombic-spot disease of fishscale bamboo.

Highly Repetitive Genome of Coniella granati (syn. Pilidiella granati), the Causal Agent of Pomegranate Fruit Rot, Encodes a Minimalistic Proteome with a Streamlined Arsenal of Effector Proteins

This study describes the first genome sequence and analysis of Coniella granati, a fungal pathogen with a broad host range, which is responsible for postharvest crown rot, shoot blight, and canker diseases in pomegranates. C. granati is a geographically widespread pathogen which has been reported across Europe, Asia, the Americas, and Africa. Our analysis revealed a 46.8 Mb genome with features characteristic of hemibiotrophic fungi. Approximately one third of its genome was compartmentalised within 'AT-rich' regions exhibiting a low GC content (30 to 45%). These regions primarily comprised transposable elements that are repeated at a high frequency and interspersed throughout the genome. Transcriptome-supported gene annotation of the C. granati genome revealed a streamlined proteome, mirroring similar observations in other pathogens with a latent phase. The genome encoded a relatively compact set of 9568 protein-coding genes with a remarkable 95% having assigned functional annotations. Despite this streamlined nature, a set of 40 cysteine-rich candidate secreted effector-like proteins (CSEPs) was predicted as well as a gene cluster involved in the synthesis of a pomegranate-associated toxin. These potential virulence factors were predominantly located near repeat-rich and AT-rich regions, suggesting that the pathogen evades host defences through Repeat-Induced Point mutation (RIP)-mediated pseudogenisation. Furthermore, 23 of these CSEPs exhibited homology to known effector and pathogenicity genes found in other hemibiotrophic pathogens. The study establishes a foundational resource for the study of the genetic makeup of C. granati, paving the way for future research on its pathogenicity mechanisms and the development of targeted control strategies to safeguard pomegranate production.

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.

Genomic and Pathogenic Characterization of Akanthomyces muscarius Isolated from Living Mite Infesting Hazelnut Big Buds

The capability of entomopathogenic fungi to live as plant endophytes is well established. However, their presence in undiscovered environmental niches represents the beginning of a new challenging research journey. Recently, Akanthomyces muscarius (Ascomycota, Cordycipitaceae) (Petch) Spatafora, Kepler & B. Shrestha was isolated from hazelnut buds infested by the big bud mite pest Phytoptus avellanae Nalepa, which makes the buds swollen, reddish, and unable to further develop. Gall formation is known to be regulated by a consortium of microbes and mites, and to better understand the possible role of A. muscarius within the infested gall, its whole genome sequence was obtained using a hybrid approach of Illumina and Nanopore reads. The functional and comparative genomics analysis provided within this study may help answer questions related to the ecology and the entomopathogenicity of this fungus.

Pan-genome survey of Septoria pistaciarum, causal agent of Septoria leaf spot of pistachios, across three Aegean sub-regions of Greece

Septoria pistaciarum, a causal agent of Septoria leaf spot disease of pistachio, is a fungal pathogen that causes substantial losses in the cultivation, worldwide. This study describes the first pan-genome-based survey of this phytopathogen-comprising a total of 27 isolates, with 9 isolates each from 3 regional units of Greece (Pieria, Larissa and Fthiotida). The reference isolate (SPF8) assembled into a total of 43.1 Mb, with 38.6% contained within AT-rich regions of approximately 37.5% G:C. The genomes of the 27 isolates exhibited on average 42% gene-coding and 20% repetitive regions. The genomes of isolates from the southern Fthiotida region appeared to more diverged from each other than the other regions based on SNP-derived trees, and also contained isolates similar to both the Pieria and Larissa regions. In contrast, isolates of the Pieria and Larissa were less diverse and distinct from one another. Asexual reproduction appeared to be typical, with no MAT1-2 locus detected in any isolate. Genome-based prediction of infection mode indicated hemibiotrophic and saprotrophic adaptations, consistent with its long latent phase. Gene prediction and orthology clustering generated a pan-genome-wide gene set of 21,174 loci. A total of 59 ortholog groups were predicted to contain candidate effector proteins, with 36 (61%) of these either having homologs to known effectors from other species or could be assigned predicted functions from matches to conserved domains. Overall, effector prediction suggests that S. pistaciarum employs a combination of defensive effectors with roles in suppression of host defenses, and offensive effectors with a range of cytotoxic activities. Some effector-like ortholog groups presented as divergent versions of the same protein, suggesting region-specific adaptations may have occurred. These findings provide insights and future research directions in uncovering the pathogenesis and population dynamics of S. pistaciarum toward the efficient management of Septoria leaf spot of pistachio.

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.

Template-Based Modelling of the Structure of Fungal Effector Proteins

The discovery of new fungal effector proteins is necessary to enable the screening of cultivars for disease resistance. Sequence-based bioinformatics methods have been used for this purpose, but only a limited number of functional effector proteins have been successfully predicted and subsequently validated experimentally. A significant obstacle is that many fungal effector proteins discovered so far lack sequence similarity or conserved sequence motifs. The availability of experimentally determined three-dimensional (3D) structures of a number of effector proteins has recently highlighted structural similarities amongst groups of sequence-dissimilar fungal effectors, enabling the search for similar structural folds amongst effector sequence candidates. We have applied template-based modelling to predict the 3D structures of candidate effector sequences obtained from bioinformatics predictions and the PHI-BASE database. Structural matches were found not only with ToxA- and MAX-like effector candidates but also with non-fungal effector-like proteins-including plant defensins and animal venoms-suggesting the broad conservation of ancestral structural folds amongst cytotoxic peptides from a diverse range of distant species. Accurate modelling of fungal effectors were achieved using RaptorX. The utility of predicted structures of effector proteins lies in the prediction of their interactions with plant receptors through molecular docking, which will improve the understanding of effector-plant interactions.

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.

Pleoardoris graminearum, gen. et sp. nov., a new member of Pleosporales from North American Plains, its biogeography and effects on a foundation grass species

Diverse fungi colonize plant roots worldwide and include species from many orders of the phylum Ascomycota. These fungi include taxa with dark septate hyphae that colonize grass roots and may modulate plant responses to stress. We describe a novel group of fungal isolates and evaluate their effects on the grass Bouteloua gracilis in vitro. We isolated fungi from roots of six native grasses from 24 sites spanning replicated latitudinal gradients in the south-central US grasslands and characterized isolates phylogenetically using a genome analysis. We analyzed 14 isolates representing a novel clade within the family Montagnulaceae (order Pleosporales), here typified as Pleoardoris graminearum, closely related to the genera Didymocrea and Bimuria. This novel species produces asexual, light brown pycnidium-like conidioma, hyaline hyphae, and chlamydospores when cultured on quinoa and kiwicha agar. To evaluate its effects on B. gracilis, seeds were inoculated with one of three isolates (DS304, DS334, and DS1613) and incubated at 25 C for 20 d. We also tested the effect of volatile organic compounds (VOCs) produced by the same isolates on B. gracilis root and stem lengths. Isolates had variable effects on plant growth. One isolate increased B. gracilis root length up to 34% compared with uninoculated controls. VOCs produced by two isolates increased root and stem lengths (P < 0.05) compared with controls. Internal transcribed spacer ITS2 metabarcode data revealed that P. graminearum is distributed across a wide range of sites in North America (22 of 24 sites sampled), and its relative abundance is influenced by host species identity and latitude. Host species identity and site were the most important factors determining P. graminearum relative abundance in drought experiments at the Extreme Drought in the Grasslands Experiment (EDGE) sites. Variable responses of B. gracilis to inoculation highlight the potential importance of nonmycorrhizal root-associated fungi on plant survival in arid ecosystems.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining (MAT) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By generating additional chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere-flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight toward the identification of a sexual cycle in Malassezia, with possible implications for pathogenicity.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining ( MAT ) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By incorporating newly generated chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere- flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight towards the discovery of a sexual cycle in Malassezia , with possible implications for pathogenicity.

Significance Statement

Malassezia , the dominant fungal group of the mammalian skin microbiome, is associated with numerous skin disorders. Sexual development and yeast-to-hyphae transitions, governed by genes at two mating-type ( MAT ) loci, are thought to be important for fungal pathogenicity. However, Malassezia sexual reproduction has never been observed. Here, we used chromosome-level assemblies and comparative genomics to uncover unforeseen transitions in MAT loci organization within Malassezia , possibly related with fragility of centromeric-associated regions. Additionally, by expressing different MAT alleles in the same cell, we show that Malassezia can undergo hyphal development and this phenotype is associated with increased expression of key mating genes along with other genes known to be virulence factors, providing a possible connection between hyphal development, sexual reproduction, and pathogenicity.

Nuclear genome organization in fungi: from gene folding to Rabl chromosomes

Comparative genomics has recently provided unprecedented insights into the biology and evolution of the fungal lineage. In the postgenomics era, a major research interest focuses now on detailing the functions of fungal genomes, i.e. how genomic information manifests into complex phenotypes. Emerging evidence across diverse eukaryotes has revealed that the organization of DNA within the nucleus is critically important. Here, we discuss the current knowledge on the fungal genome organization, from the association of chromosomes within the nucleus to topological structures at individual genes and the genetic factors required for this hierarchical organization. Chromosome conformation capture followed by high-throughput sequencing (Hi-C) has elucidated how fungal genomes are globally organized in Rabl configuration, in which centromere or telomere bundles are associated with opposite faces of the nuclear envelope. Further, fungal genomes are regionally organized into topologically associated domain-like (TAD-like) chromatin structures. We discuss how chromatin organization impacts the proper function of DNA-templated processes across the fungal genome. Nevertheless, this view is limited to a few fungal taxa given the paucity of fungal Hi-C experiments. We advocate for exploring genome organization across diverse fungal lineages to ensure the future understanding of the impact of nuclear organization on fungal genome function.

Hybrid de novo genome assembly and comparative genomics of three different isolates of Gnomoniopsis castaneae

The first genome assemblies of Gnomoniopsis castaneae (syn. G. smithogilvyi), the causal agent of chestnut brown rot of kernels, shoot blight and cankers, are provided here. Specifically, the complete genome of the Italian ex-type MUT401 isolate was compared to the draft genome of a second Italian isolate (GN01) and to the ICMP 14040 isolate from New Zealand. The three genome sequences were obtained through a hybrid assembly using both short Illumina reads and long Nanopore reads, their coding sequences were annotated and compared with each other and with other Diaporthales. The information offered by the genome assembly of the three isolates represents the base of data for further application related to -omics strategies of the fungus and to develop markers for population studies at a local and global scale.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

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.

Comparative genomics of five Valsa species gives insights on their pathogenicity evolution

Valsa is a genus of ascomycetes within the Valsaceae family. This family includes many wood destructive pathogens such as the well known Valsa mali and Valsa pyri which cause canker diseases in fruit trees and threaten the global fruit production. Lack of genomic information of this family is impeding our understandings about their evolution and genetic basis of their pathogenicity divergence. Here, we report genome assemblies of Valsa malicola, Valsa persoonii, and Valsa sordida which represent close relatives of Valsa mali and Valsa pyri with different host preferences. Comparative genomics analysis revealed that segmental rearrangements, inversions, and translocations frequently occurred among Valsa spp. genomes. Gene families that exhibited gene copy expansions tended to be associated with secondary metabolism, transmembrane transport, and pyrophosphatase activities. Orthologous genes in regions lost synteny exhibited significantly higher rate of synonymous substitution (KS) than those in regions retained synteny. Moreover, among these genes, membrane transporter families associated with antidrug (MFS, DHA) activities and nutrient transportation (SP and APCs) activities were significantly over-represented. Lineage specific synonymous substitution (KS) and nonsynonymous substitution (KA) analysis based on the phylogeny constructed from 11 fungal species identified a set of genes with selection signatures in Valsa clade and these genes were significantly enriched in functions associated with fatty acid beta-oxidation, DNA helicase activity, and ATPase activity. Furthermore, unique genes that possessed or retained by each of the five Valsa species are more likely part of the secondary metabolic (SM) gene clusters. SM gene clusters conserved across five Valsa species showed various degrees of diversification in both identity and completeness. All 11 syntenically conserved SM clusters showed differential expression during the infection of apple branch with Valsa mali suggesting involvements of secondary metabolism in the pathogenicity of Valsa species.

Metagenomics Shines Light on the Evolution of "Sunscreen" Pigment Metabolism in the Teloschistales (Lichen-Forming Ascomycota)

Fungi produce a vast number of secondary metabolites that shape their interactions with other organisms and the environment. Characterizing the genes underpinning metabolite synthesis is therefore key to understanding fungal evolution and adaptation. Lichenized fungi represent almost one-third of Ascomycota diversity and boast impressive secondary metabolites repertoires. However, most lichen biosynthetic genes have not been linked to their metabolite products. Here we used metagenomic sequencing to survey gene families associated with production of anthraquinones, UV-protectant secondary metabolites present in various fungi, but especially abundant in a diverse order of lichens, the Teloschistales (class Lecanoromycetes, phylum Ascomycota). We successfully assembled 24 new, high-quality lichenized-fungal genomes de novo and combined them with publicly available Lecanoromycetes genomes from taxa with diverse secondary chemistry to produce a whole-genome tree. Secondary metabolite biosynthetic gene cluster (BGC) analysis showed that whilst lichen BGCs are numerous and highly dissimilar, core enzyme genes are generally conserved across taxa. This suggests metabolite diversification occurs via re-shuffling existing enzyme genes with novel accessory genes rather than BGC gains/losses or de novo gene evolution. We identified putative anthraquinone BGCs in our lichen dataset that appear homologous to anthraquinone clusters from non-lichenized fungi, suggesting these genes were present in the common ancestor of the subphylum Pezizomycotina. Finally, we identified unique transporter genes in Teloschistales anthraquinone BGCs that may explain why these metabolites are so abundant and ubiquitous in these lichens. Our results support the importance of metagenomics for understanding the secondary metabolism of non-model fungi such as lichens.

Isolation of strains and their genome sequencing to analyze the mating system of Ophiocordyceps robertsii

The fungal genus Ophiocordyceps contains a number of insect pathogens. One of the best known of these is Ophiocordyceps sinensis, which is used in Chinese medicine and its overharvesting threatens sustainability; hence, alternative species are being sought. Ophiocordyceps robertsii, found in Australia and New Zealand, has been proposed to be a close relative to O. sinensis, but little is known about this species despite being also of historical significance. Here, O. robertsii strains were isolated into culture and high coverage draft genome sequences obtained and analyzed. This species has a large genome expansion, as also occurred in O. sinensis. The mating type locus was characterized, indicating a heterothallic arrangement whereby each strain has an idiomorphic region of two (MAT1-2-1, MAT1-2-2) or three (MAT1-1-1, MAT1-1-2, MAT1-1-3) genes flanked by the conserved APN2 and SLA2 genes. These resources provide a new opportunity for understanding the evolution of the expanded genome in the homothallic species O. sinensis, as well as capabilities to explore the pharmaceutical potential in a species endemic to Australia and New Zealand.

A global pangenome for the wheat fungal pathogen Pyrenophora tritici-repentis and prediction of effector protein structural homology

The adaptive potential of plant fungal pathogens is largely governed by the gene content of a species, consisting of core and accessory genes across the pathogen isolate repertoire. To approximate the complete gene repertoire of a globally significant crop fungal pathogen, a pan genomic analysis was undertaken for Pyrenophora tritici-repentis (Ptr), the causal agent of tan (or yellow) spot disease in wheat. In this study, 15 new Ptr genomes were sequenced, assembled and annotated, including isolates from three races not previously sequenced. Together with 11 previously published Ptr genomes, a pangenome for 26 Ptr isolates from Australia, Europe, North Africa and America, representing nearly all known races, revealed a conserved core-gene content of 57 % and presents a new Ptr resource for searching natural homologues (orthologues not acquired by horizontal transfer from another species) using remote protein structural homology. Here, we identify for the first time a non-synonymous mutation in the Ptr necrotrophic effector gene ToxB, multiple copies of the inactive toxb within an isolate, a distant natural Pyrenophora homologue of a known Parastagonopora nodorum necrotrophic effector (SnTox3), and clear genomic break points for the ToxA effector horizontal transfer region. This comprehensive genomic analysis of Ptr races includes nine isolates sequenced via long read technologies. Accordingly, these resources provide a more complete representation of the species, and serve as a resource to monitor variations potentially involved in pathogenicity.

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.

Two independent approaches converge to the cloning of a new Leptosphaeria maculans avirulence effector gene, AvrLmS-Lep2

Brassica napus (oilseed rape, canola) seedling resistance to Leptosphaeria maculans, the causal agent of blackleg (stem canker) disease, follows a gene-for-gene relationship. The avirulence genes AvrLmS and AvrLep2 were described to be perceived by the resistance genes RlmS and LepR2, respectively, present in B. napus 'Surpass 400'. Here we report cloning of AvrLmS and AvrLep2 using two independent methods. AvrLmS was cloned using combined in vitro crossing between avirulent and virulent isolates with sequencing of DNA bulks from avirulent or virulent progeny (bulked segregant sequencing). AvrLep2 was cloned using a biparental cross of avirulent and virulent L. maculans isolates and a classical map-based cloning approach. Taking these two approaches independently, we found that AvrLmS and AvrLep2 are the same gene. Complementation of virulent isolates with this gene confirmed its role in inducing resistance on Surpass 400, Topas-LepR2, and an RlmS-line. The gene, renamed AvrLmS-Lep2, encodes a small cysteine-rich protein of unknown function with an N-terminal secretory signal peptide, which is a common feature of the majority of effectors from extracellular fungal plant pathogens. The AvrLmS-Lep2/LepR2 interaction phenotype was found to vary from a typical hypersensitive response through intermediate resistance sometimes towards susceptibility, depending on the inoculation conditions. AvrLmS-Lep2 was nevertheless sufficient to significantly slow the systemic growth of the pathogen and reduce the stem lesion size on plant genotypes with LepR2, indicating the potential efficiency of this resistance to control the disease in the field.

Genomic Variations and Mutational Events Associated with Plant-Pathogen Interactions

Phytopathologists are actively researching the molecular basis of plant-pathogen interactions. The mechanisms of responses to pathogens have been studied extensively in model crop plant species and natural populations. Today, with the rapid expansion of genomic technologies such as DNA sequencing, transcriptomics, proteomics, and metabolomics, as well as the development of new methods and protocols, data analysis, and bioinformatics, it is now possible to assess the role of genetic variation in plant-microbe interactions and to understand the underlying molecular mechanisms of plant defense and microbe pathogenicity with ever-greater resolution and accuracy. Genetic variation is an important force in evolution that enables organisms to survive in stressful environments. Moreover, understanding the role of genetic variation and mutational events is essential for crop breeders to produce improved cultivars. This review focuses on genetic variations and mutational events associated with plant-pathogen interactions and discusses how these genome compartments enhance plants' and pathogens' evolutionary processes.

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

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

Hybrid de novo genome-reassembly reveals new insights on pathways and pathogenicity determinants in rice blast pathogen Magnaporthe oryzae RMg_Dl

Blast disease incited by Magnaporthe oryzae is a major threat to sustain rice production in all rice growing nations. The pathogen is widely distributed in all rice paddies and displays rapid aerial transmissions, and seed-borne latent infection. In order to understand the genetic variability, host specificity, and molecular basis of the pathogenicity-associated traits, the whole genome of rice infecting Magnaporthe oryzae (Strain RMg_Dl) was sequenced using the Illumina and PacBio (RSII compatible) platforms. The high-throughput hybrid assembly of short and long reads resulted in a total of 375 scaffolds with a genome size of 42.43 Mb. Furthermore, comparative genome analysis revealed 99% average nucleotide identity (ANI) with other oryzae genomes and 83% against M. grisea, and 73% against M. poe genomes. The gene calling identified 10,553 genes with 10,539 protein-coding sequences. Among the detected transposable elements, the LTR/Gypsy and Type LINE showed high occurrence. The InterProScan of predicted protein sequences revealed that 97% protein family (PFAM), 98% superfamily, and 95% CDD were shared among RMg_Dl and reference 70-15 genome, respectively. Additionally, 550 CAZymes with high GH family content/distribution and cell wall degrading enzymes (CWDE) such endoglucanase, beta-glucosidase, and pectate lyase were also deciphered in RMg_Dl. The prevalence of virulence factors determination revealed that 51 different VFs were found in the genome. The biochemical pathway such as starch and sucrose metabolism, mTOR signaling, cAMP signaling, MAPK signaling pathways related genes were identified in the genome. The 49,065 SNPs, 3267 insertions and 3611 deletions were detected, and majority of these varinats were located on downstream and upstream region. Taken together, the generated information will be useful to develop a specific marker for diagnosis, pathogen surveillance and tracking, molecular taxonomy, and species delineation which ultimately leads to device improved management strategies for blast disease.

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.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

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.

The "Bipartite" Structure of the First Genome of Ampelomyces quisqualis, a Common Hyperparasite and Biocontrol Agent of Powdery Mildews, May Point to Its Evolutionary Origin from Plant Pathogenic Fungi

Powdery mildews are among the most important plant pathogens worldwide, which are often attacked in the field by mycoparasitic fungi belonging to the genus Ampelomyces. The taxonomy of the genus Ampelomyces is unresolved, but well-supported molecular operational taxonomic units were repeatedly defined suggesting that the genus may include at least four to seven species. Some Ampelomyces strains were commercialized as biocontrol agents of crop pathogenic powdery mildews. However, the genomic mechanisms underlying their mycoparasitism are still poorly understood. To date, the draft genome of a single Ampelomyces strain, designated as HMLAC 05119, has been released. We report a high-quality, annotated hybrid draft genome assembly of A. quisqualis strain BRIP 72107, which, based on phylogenetic analyses, is not conspecific with HMLAC 05119. The constructed genome is 40.38 Mb in size, consisting of 24 scaffolds with an N50 of 2.99 Mb and 96.2% completeness. Our analyses revealed "bipartite" structure of Ampelomyces genomes, where GC-balanced genomic regions are interspersed by longer or shorter stretches of AT-rich regions. This is also a hallmark of many plant pathogenic fungi and provides further evidence for evolutionary affinity of Ampelomyces species to plant pathogenic fungi. The high-quality genome and annotation produced here provide an important resource for future genomic studies of mycoparasitisim to decipher molecular mechanisms underlying biocontrol processes and natural tritrophic interactions.

Gene Validation and Remodelling Using Proteogenomics of Phytophthora cinnamomi, the Causal Agent of Dieback

Phytophthora cinnamomi is a pathogenic oomycete that causes plant dieback disease across a range of natural ecosystems and in many agriculturally important crops on a global scale. An annotated draft genome sequence is publicly available (JGI Mycocosm) and suggests 26,131 gene models. In this study, soluble mycelial, extracellular (secretome), and zoospore proteins of P. cinnamomi were exploited to refine the genome by correcting gene annotations and discovering novel genes. By implementing the diverse set of sub-proteomes into a generated proteogenomics pipeline, we were able to improve the P. cinnamomi genome annotation. Liquid chromatography mass spectrometry was used to obtain high confidence peptides with spectral matching to both the annotated genome and a generated 6-frame translation. Two thousand seven hundred sixty-four annotations from the draft genome were confirmed by spectral matching. Using a proteogenomic pipeline, mass spectra were used to edit the P. cinnamomi genome and allowed identification of 23 new gene models and 60 edited gene features using high confidence peptides obtained by mass spectrometry, suggesting a rate of incorrect annotations of 3% of the detectable proteome. The novel features were further validated by total peptide support, alongside functional analysis including the use of Gene Ontology and functional domain identification. We demonstrated the use of spectral data in combination with our proteogenomics pipeline can be used to improve the genome annotation of important plant diseases and identify missed genes. This study presents the first use of spectral data to edit and manually annotate an oomycete pathogen.

Reference genome assembly for Australian Ascochyta lentis isolate Al4

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use, and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

Genome sequence of Monilinia vaccinii-corymbosi sheds light on mummy berry disease infection of blueberry and mating type

Mummy berry disease, caused by the fungal pathogen Monilinia vaccinii-corymbosi (Mvc), is one of the most economically important diseases of blueberries in North America. Mvc is capable of inducing two separate blighting stages during its life cycle. Infected fruits are rendered mummified and unmarketable. Genomic data for this pathogen is lacking, but could be useful in understanding the reproductive biology of Mvc and the mechanisms it deploys to facilitate host infection. In this study, PacBio sequencing and Hi-C interaction data were utilized to create a chromosome-scale reference genome for Mvc. The genome comprises nine chromosomes with a total length of 30 Mb, an N50 length of 4.06 Mb, and an average 413X sequence coverage. A total of 9399 gene models were predicted and annotated, and BUSCO analysis revealed that 98% of 1,438 searched conserved eukaryotic genes were present in the predicted gene set. Potential effectors were identified, and the mating-type (MAT) locus was characterized. Biotrophic effectors allow the pathogen to avoid recognition by the host plant and evade or mitigate host defense responses during the early stages of fruit infection. Following locule colonization, necrotizing effectors promote the mummification of host tissues. Potential biotrophic effectors utilized by Mvc include chorismate mutase for reducing host salicylate and necrotrophic effectors include necrosis-inducing proteins and hydrolytic enzymes for macerating host tissue. The MAT locus sequences indicate the potential for homothallism in the reference genome, but a deletion allele of the MAT locus, characterized in a second isolate, indicates heterothallism. Further research is needed to verify the roles of individual effectors in virulence and to determine the role of the MAT locus in outcrossing and population genotypic diversity.

Chromosome-level genome assembly and manually-curated proteome of model necrotroph Parastagonospora nodorum Sn15 reveals a genome-wide trove of candidate effector homologs, and redundancy of virulence-related functions within an accessory chromosome

Background

The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3.

Results

A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions.

Conclusions

We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07699-8.

Austropuccinia psidii, causing myrtle rust, has a gigabase-sized genome shaped by transposable elements

Austropuccinia psidii, originating in South America, is a globally invasive fungal plant pathogen that causes rust disease on Myrtaceae. Several biotypes are recognized, with the most widely distributed pandemic biotype spreading throughout the Asia-Pacific and Oceania regions over the last decade. Austropuccinia psidii has a broad host range with more than 480 myrtaceous species. Since first detected in Australia in 2010, the pathogen has caused the near extinction of at least three species and negatively affected commercial production of several Myrtaceae. To enable molecular and evolutionary studies into A. psidii pathogenicity, we assembled a highly contiguous genome for the pandemic biotype. With an estimated haploid genome size of just over 1 Gb (gigabases), it is the largest assembled fungal genome to date. The genome has undergone massive expansion via distinct transposable element (TE) bursts. Over 90% of the genome is covered by TEs predominantly belonging to the Gypsy superfamily. These TE bursts have likely been followed by deamination events of methylated cytosines to silence the repetitive elements. This in turn led to the depletion of CpG sites in TEs and a very low overall GC content of 33.8%. Compared to other Pucciniales, the intergenic distances are increased by an order of magnitude indicating a general insertion of TEs between genes. Overall, we show how TEs shaped the genome evolution of A. psidii and provide a greatly needed resource for strategic approaches to combat disease spread.

Large-scale transcriptomics to dissect 2 years of the life of a fungal phytopathogen interacting with its host plant

BACKGROUND

The fungus Leptosphaeria maculans has an exceptionally long and complex relationship with its host plant, Brassica napus, during which it switches between different lifestyles, including asymptomatic, biotrophic, necrotrophic, and saprotrophic stages. The fungus is also exemplary of "two-speed" genome organisms in the genome of which gene-rich and repeat-rich regions alternate. Except for a few stages of plant infection under controlled conditions, nothing is known about the genes mobilized by the fungus throughout its life cycle, which may last several years in the field.

RESULTS

We performed RNA-seq on samples corresponding to all stages of the interaction of L. maculans with its host plant, either alive or dead (stem residues after harvest) in controlled conditions or in field experiments under natural inoculum pressure, over periods of time ranging from a few days to months or years. A total of 102 biological samples corresponding to 37 sets of conditions were analyzed. We show here that about 9% of the genes of this fungus are highly expressed during its interactions with its host plant. These genes are distributed into eight well-defined expression clusters, corresponding to specific infection lifestyles or to tissue-specific genes. All expression clusters are enriched in effector genes, and one cluster is specific to the saprophytic lifestyle on plant residues. One cluster, including genes known to be involved in the first phase of asymptomatic fungal growth in leaves, is re-used at each asymptomatic growth stage, regardless of the type of organ infected. The expression of the genes of this cluster is repeatedly turned on and off during infection. Whatever their expression profile, the genes of these clusters are enriched in heterochromatin regions associated with H3K9me3 or H3K27me3 repressive marks. These findings provide support for the hypothesis that part of the fungal genes involved in niche adaptation is located in heterochromatic regions of the genome, conferring an extreme plasticity of expression.

CONCLUSION

This work opens up new avenues for plant disease control, by identifying stage-specific effectors that could be used as targets for the identification of novel durable disease resistance genes, or for the in-depth analysis of chromatin remodeling during plant infection, which could be manipulated to interfere with the global expression of effector genes at crucial stages of plant infection.

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.

Comparative Genome-Wide Characterization of Microsatellites in Candida albicans and Candida dubliniensis Leading to the Development of Species-Specific Marker

BACKGROUND

Microsatellites or simple sequence repeats (SSR) are related to genomic structure, function, and certain diseases of taxonomically different organisms.

OBJECTIVE

To characterize microsatellites in two closely related Candida species by searching and comparing 1-6 bp nucleotide motifs and utilizing them to develop species-specific markers.

METHODS

Whole-genome sequence was downloaded from the public domain, microsatellites were mined and analyzed, and primers were synthesized.

RESULTS

A total of 15,821 and 7,868 microsatellites, with mono-nucleotides (8,679) and trinucleotides (3,156) as most frequent microsatellites, were mined in Candida dubliniensis and Candida albicans, respectively. Chromosome size was found positively correlated with microsatellite number in both the species, whereas it was negatively correlated with the relative abundance and density of microsatellites. A number of unique motifs were also found in both the species. Overall, microsatellite frequencies of each chromosome in C. dubliniensis were higher than in C. albicans.

CONCLUSION

The features of microsatellite distribution in the two species' genomes revealed that it is probably not conserved in the genus Candida. Data generated in this article could be used for comparative genome mapping and understanding the distribution of microsatellites and genome structure between these closely related and phenotypically misidentified species and may provide a foundation for the development of a new set of species-specific microsatellite markers. Here, we also report a novel microsatellite-based marker for C. dubliniensis-specific identification.

Codon usage signatures in the genus Cryptococcus: A complex interplay of gene expression, translational selection and compositional bias

The fungal genus Cryptococcus comprises of several diverse species. The pathogens forming Cryptococcus neoformans/ Cryptococcus gatti species complex are of immense clinical significance owing to the high frequency of infections and deaths globally. Three closely related non-pathogenic species namely, Cryptococcus amylolentus, Cryptococcus wingfieldii and Cryptococcus depauperatus are the non-pathogenic ancestral species from which pathogenic lineages have diverged. In the current study, a comprehensive analysis of factors influencing the codon and amino acid usage bias in six pathogenic and three non-pathogenic species was performed. Our results revealed that though compositional bias played a crucial role, translational selection and gene expression were the key determinants of codon usage variations. Analysis of relative dinucleotide abundance and codon context signatures revealed strict avoidance of TpA dinucleotide across genomes. Multivariate statistical analysis based on codon usage data resulted in discrete clustering of pathogens and non-pathogens which correlated with previous reports on their phylogenetic distribution.

Substitution scoring matrices for proteins - An overview

Sequence analysis is the primary and simplest approach to discover structural, functional and evolutionary details of related proteins. All the alignment based approaches of sequence analysis make use of amino acid substitution matrices, and the accuracy of the results largely depends on the type of scoring matrices used to perform alignment tasks. An amino acid substitution matrix is a 20 × 20 matrix in which the individual elements encapsulate the rates at which each of the 20 amino acid residues in proteins are substituted by other amino acid residues over time. In contrast to most globular/ordered proteins whose amino acids composition is considered as standard, there are several classes of proteins (e.g., transmembrane proteins) in which certain types of amino acid (e.g., hydrophobic residues) are enriched. These compositional differences among various classes of proteins are manifested in their underlying residue substitution frequencies. Therefore, each of the compositionally distinct class of proteins or protein segments should be studied using specific scoring matrices that reflect their distinct residue substitution pattern. In this review, we describe the development and application of various substitution scoring matrices peculiar to proteins with standard and biased compositions. Along with most commonly used standard matrices (PAM, BLOSUM, MD and VTML) that act as default parameters in various homologs search and alignment tools, different substitution scoring matrices specific to compositionally distinct class of proteins are discussed in detail.

Reference Genome Assembly for Australian Ascochyta rabiei Isolate ArME14

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

An Avirulence Gene Cluster in the Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) Identified through Genetic Mapping and Whole-Genome Sequencing of a Sexual Population

Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe (yellow) rust, is an obligate, biotrophic fungus. It was difficult to study the genetics of the pathogen due to the lack of sexual reproduction. The recent discovery of alternate hosts for P. striiformis f. sp. tritici makes it possible to study inheritance and map genes involved in its interaction with plant hosts. To identify avirulence (Avr) genes in P. striiformis f. sp. tritici, we developed a segregating population by selfing isolate 12-368 on barberry (Berberis vulgaris) plants under controlled conditions. The dikaryotic sexual population segregated for avirulent/virulent phenotypes on nine Yr single-gene lines. The parental and progeny isolates were whole-genome sequenced at >30× coverage using Illumina HiSeq PE150 technology. A total of 2,637 high-quality markers were discovered by mapping the whole-genome sequencing (WGS) reads to the reference genome of strain 93-210 and used to construct a genetic map, consisting of 41 linkage groups, spanning 7,715.0 centimorgans (cM) and covering 68 Mb of the reference genome. The recombination rate was estimated to be 1.81 ± 2.32 cM/10 kb. Quantitative trait locus analysis mapped six Avr gene loci to the genetic map, including an Avr cluster harboring four Avr genes, AvYr7, AvYr43, AvYr44, and AvYrExp2 Aligning the genetic map to the reference genome identified Avr candidates and narrowed them to a small genomic region (<200 kb). The discovery of the Avr gene cluster is useful for understanding pathogen evolution, and the identification of candidate genes is an important step toward cloning Avr genes for studying molecular mechanisms of pathogen-host interactions.IMPORTANCE Stripe rust is a destructive disease of wheat worldwide. Growing resistant cultivars is the most effective, easy-to-use, economical, and environmentally friendly strategy for the control of the disease. However, P. striiformis f. sp. tritici can produce new virulent races that may circumvent race-specific resistance. Therefore, understanding the genetic basis of the interactions between wheat genes for resistance and P. striiformis f. sp. tritici genes for avirulence is useful for improving cultivar resistance for more effective control of the disease. This study developed a high-quality map that facilitates genomic and genetic studies of important traits related to pathogen pathogenicity and adaptation to different environments and crop cultivars carrying different resistance genes. The information on avirulence/virulence genes identified in this study can be used for guiding breeding programs to select combinations of genes for developing new cultivars with effective resistance to mitigate this devastating disease.

Genome mining of the citrus pathogen Elsinoë fawcettii; prediction and prioritisation of candidate effectors, cell wall degrading enzymes and secondary metabolite gene clusters

Elsinoë fawcettii, a necrotrophic fungal pathogen, causes citrus scab on numerous citrus varieties around the world. Known pathotypes of E. fawcettii are based on host range; additionally, cryptic pathotypes have been reported and more novel pathotypes are thought to exist. E. fawcettii produces elsinochrome, a non-host selective toxin which contributes to virulence. However, the mechanisms involved in potential pathogen-host interactions occurring prior to the production of elsinochrome are unknown, yet the host-specificity observed among pathotypes suggests a reliance upon such mechanisms. In this study we have generated a whole genome sequencing project for E. fawcettii, producing an annotated draft assembly 26.01 Mb in size, with 10,080 predicted gene models and low (0.37%) coverage of transposable elements. A small proportion of the assembly showed evidence of AT-rich regions, potentially indicating genomic regions with increased plasticity. Using a variety of computational tools, we mined the E. fawcettii genome for potential virulence genes as candidates for future investigation. A total of 1,280 secreted proteins and 276 candidate effectors were predicted and compared to those of other necrotrophic (Botrytis cinerea, Parastagonospora nodorum, Pyrenophora tritici-repentis, Sclerotinia sclerotiorum and Zymoseptoria tritici), hemibiotrophic (Leptosphaeria maculans, Magnaporthe oryzae, Rhynchosporium commune and Verticillium dahliae) and biotrophic (Ustilago maydis) plant pathogens. Genomic and proteomic features of known fungal effectors were analysed and used to guide the prioritisation of 120 candidate effectors of E. fawcettii. Additionally, 378 carbohydrate-active enzymes were predicted and analysed for likely secretion and sequence similarity with known virulence genes. Furthermore, secondary metabolite prediction indicated nine additional genes potentially involved in the elsinochrome biosynthesis gene cluster than previously described. A further 21 secondary metabolite clusters were predicted, some with similarity to known toxin producing gene clusters. The candidate virulence genes predicted in this study provide a comprehensive resource for future experimental investigation into the pathogenesis of E. fawcettii.