A New Reference Genome Shows the One-Speed Genome Structure of the Barley Pathogen Ramularia collo-cygni

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 soil-borne white root rot pathogen Rosellinia necatrix expresses antimicrobial proteins during host colonization

Rosellinia necatrix is a prevalent soil-borne plant-pathogenic fungus that is the causal agent of white root rot disease in a broad range of host plants. The limited availability of genomic resources for R. necatrix has complicated a thorough understanding of its infection biology. Here, we sequenced nine R. necatrix strains with Oxford Nanopore sequencing technology, and with DNA proximity ligation we generated a gapless assembly of one of the genomes into ten chromosomes. Whereas many filamentous pathogens display a so-called two-speed genome with more dynamic and more conserved compartments, the R. necatrix genome does not display such genome compartmentalization. It has recently been proposed that fungal plant pathogens may employ effectors with antimicrobial activity to manipulate the host microbiota to promote infection. In the predicted secretome of R. necatrix, 26 putative antimicrobial effector proteins were identified, nine of which are expressed during plant colonization. Two of the candidates were tested, both of which were found to possess selective antimicrobial activity. Intriguingly, some of the inhibited bacteria are antagonists of R. necatrix growth in vitro and can alleviate R. necatrix infection on cotton plants. Collectively, our data show that R. necatrix encodes antimicrobials that are expressed during host colonization and that may contribute to modulation of host-associated microbiota to stimulate disease development.

Fungal Effectoromics: A World in Constant Evolution

Effectors are small, secreted molecules that mediate the establishment of interactions in nature. While some concepts of effector biology have stood the test of time, this area of study is ever-evolving as new effectors and associated characteristics are being revealed. In the present review, the different characteristics that underly effector classifications are discussed, contrasting past and present knowledge regarding these molecules to foster a more comprehensive understanding of effectors for the reader. Research gaps in effector identification and perspectives for effector application in plant disease management are also presented, with a focus on fungal effectors in the plant-microbe interaction and interactions beyond the plant host. In summary, the review provides an amenable yet thorough introduction to fungal effector biology, presenting noteworthy examples of effectors and effector studies that have shaped our present understanding of the field.

High Nucleotide Substitution Rates Associated with Retrotransposon Proliferation Drive Dynamic Secretome Evolution in Smut Pathogens

Transposable elements (TEs) play a pivotal role in shaping diversity in eukaryotic genomes. The covered smut pathogen on barley, Ustilago hordei, encountered a recent genome expansion. Using long reads, we assembled genomes of 6 U. hordei strains and 3 sister species, to study this genome expansion. We found that larger genome sizes can mainly be attributed to a higher genome fraction of long terminal repeat retrotransposons (LTR-RTs). In the studied smut genomes, LTR-RTs fractions are the largest in U. hordei and are positively correlated with the mating-type locus sizes, which is up to ~560 kb in U. hordei. Furthermore, LTR-RTs were found to be associated with higher nucleotide substitution levels, as these occur in specific genome regions of smut species with a recent LTR-RT proliferation. Moreover, genes in genome regions with higher nucleotide substitution levels generally reside closer to LTR-RTs than other genome regions. Genome regions with many nucleotide substitutions encountered an especially high fraction of CG substitutions, which is not observed for LTR-RT sequences. The high nucleotide substitution levels particularly accelerate the evolution of secretome genes, as their more accessory nature results in substitutions that often lead to amino acid alterations. IMPORTANCE Genomic alteration can be generated through various means, in which transposable elements (TEs) can play a pivotal role. Their mobility causes mutagenesis in itself and can disrupt the function of the sequences they insert into. They also impact genome evolution as their repetitive nature facilitates nonhomologous recombination. Furthermore, TEs have been linked to specific epigenetic genome organizations. We report a recent TE proliferation in the genome of the barley covered smut fungus, Ustilago hordei. This proliferation is associated with a distinct nucleotide substitution regime that has a higher rate and a higher fraction of CG substitutions. This different regime shapes the evolution of genes in subjected genome regions. We hypothesize that TEs may influence the error-rate of DNA polymerase in a hitherto unknown fashion.

Unveiling the Core Effector Proteins of Oil Palm Pathogen Ganoderma boninense via Pan-Secretome Analysis

Ganoderma boninense is the major causal agent of basal stem rot (BSR) disease in oil palm, causing the progressive rot of the basal part of the stem. Despite its prominence, the key pathogenicity determinants for the aggressive nature of hemibiotrophic infection remain unknown. In this study, genome sequencing and the annotation of G. boninense T10 were carried out using the Illumina sequencing platform, and comparative genome analysis was performed with previously reported G. boninense strains (NJ3 and G3). The pan-secretome of G. boninense was constructed and comprised 937 core orthogroups, 243 accessory orthogroups, and 84 strain-specific orthogroups. In total, 320 core orthogroups were enriched with candidate effector proteins (CEPs) that could be classified as carbohydrate-active enzymes, hydrolases, and non-catalytic proteins. Differential expression analysis revealed an upregulation of five CEP genes that was linked to the suppression of PTI signaling cascade, while the downregulation of four CEP genes was linked to the inhibition of PTI by preventing host defense elicitation. Genome architecture analysis revealed the one-speed architecture of the G. boninense genome and the lack of preferential association of CEP genes to transposable elements. The findings obtained from this study aid in the characterization of pathogenicity determinants and molecular biomarkers of BSR disease.

Fungal endophytes in plants and their relationship to plant disease

The enigmatic endophytic fungi are beginning to reveal their secrets. Like pathogens, they can manipulate the host for their own benefit to create their own optimal habitat. Some endophytic manipulations induce resistance or otherwise outcompete pathogens and can thus be exploited for biological control. Like pathogens and other symbionts, endophytes produce effector proteins and other molecules, ranging from specialised metabolites, phytohormones and microRNAs, to manipulate their hosts and other microorganisms they meet. There is a continuum from endophyte to pathogen: some organisms can infest or cause disease in some hosts, but not in others. Molecular genetics approaches coupled with functional characterisation have demonstrated their worth for understanding the biological phenomena underlying endophytic fungal interactions.

The contribution of DNA repair pathways to genome editing and evolution in filamentous pathogens

DNA double-strand breaks require repair or risk corrupting the language of life. To ensure genome integrity and viability, multiple DNA double-strand break repair pathways function in eukaryotes. Two such repair pathways, canonical non-homologous end joining and homologous recombination, have been extensively studied, while other pathways such as microhomology-mediated end joint and single-strand annealing, once thought to serve as back-ups, now appear to play a fundamental role in DNA repair. Here, we review the molecular details and hierarchy of these four DNA repair pathways, and where possible, a comparison for what is known between animal and fungal models. We address the factors contributing to break repair pathway choice, and aim to explore our understanding and knowledge gaps regarding mechanisms and regulation in filamentous pathogens. We additionally discuss how DNA double-strand break repair pathways influence genome engineering results, including unexpected mutation outcomes. Finally, we review the concept of biased genome evolution in filamentous pathogens, and provide a model, termed Biased Variation, that links DNA double-strand break repair pathways with properties of genome evolution. Despite our extensive knowledge for this universal process, there remain many unanswered questions, for which the answers may improve genome engineering and our understanding of genome evolution.

Detection of Ramularia collo-cygni from barley in Australia using triplex quantitative and droplet digital PCR

BACKGROUND

Ramularia leaf spot (RLS), caused by Ramularia collo-cygni, is an emerging threat to barley (Hordeum vulgare L.) production. RLS has been reported in Australia, however only minimal information is available regarding its detection and distribution. Due to initial asymptomatic growth in planta, slow growth in vitro and symptomatic similarities to net blotch and physiological leaf spots, detection of this pathogen can be challenging. Quantitative polymerase chain reaction (PCR)-based methods for R. collo-cygni-specific identification and detection have been described, however these assays have been demonstrated to lack specificity. False-positive detections may have serious implications, thus we aimed to design a robust R. collo-cygni-specific PCR method.

RESULTS

Using the phylogenetically informative RNA polymerase II second largest subunit (rpb2) and translation elongation factor 1-alpha (tef1-α) genes, along with the tef1-α gene of H. vulgare, a triplex assay was developed for both quantitative and droplet digital PCR. The triplex assay detected R. collo-cygni DNA in barley leaves from New South Wales, South Australia, Tasmania, Victoria and Western Australia. No R. collo-cygni DNA was detected in barley seed grown in Western Australia.

CONCLUSION

The presence of R. collo-cygni DNA has been confirmed in Australian barley crops, suggesting a distribution ranging across the southern barley growing regions of Australia. The R. collo-cygni-specific assay will be a valuable tool to assist with monitoring the distribution and impact of R. collo-cygni in Australia and other regions. © 2021 Society of Chemical Industry.

Positive Selection of Transcription Factors Is a Prominent Feature of the Evolution of a Plant Pathogenic Genus Originating in the Miocene

Tests based on the d_N/d_S statistic are used to identify positive selection of nonsynonymous polymorphisms. Using these tests on alignments of all orthologs from related species can provide insights into which gene categories have been most frequently positively selected. However, longer alignments have more power to detect positive selection, creating a detection bias that could create misleading results from functional enrichment tests. Most studies of positive selection in plant pathogens focus on genes with specific virulence functions, with little emphasis on broader molecular processes. Furthermore, no studies in plant pathogens have accounted for detection bias due to alignment length when performing functional enrichment tests. To address these research gaps, we analyze 12 genomes of the phytopathogenic fungal genus Botrytis, including two sequenced in this study. To establish a temporal context, we estimated fossil-calibrated divergence times for the genus. We find that Botrytis likely originated 16–18 Ma in the Miocene and underwent continuous radiation ending in the Pliocene. An untargeted scan of Botrytis single-copy orthologs for positive selection with three different statistical tests uncovered evidence for positive selection among proteases, signaling proteins, CAZymes, and secreted proteins. There was also a strong overrepresentation of transcription factors among positively selected genes. This overrepresentation was still apparent after two complementary controls for detection bias due to sequence length. Positively selected sites were depleted within DNA-binding domains, suggesting changes in transcriptional responses to internal and external cues or protein–protein interactions have undergone positive selection more frequently than changes in promoter fidelity.

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.

Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Stripe rust of wheat, caused by the obligate biotrophic fungus Puccinia striiformis f.sp. tritici, is a major threat to wheat production worldwide with an estimated yearly loss of US $1 billion. The recent advances in long-read sequencing technologies and tailored-assembly algorithms enabled us to disentangle the two haploid genomes of Pst. This provides us with haplotype-specific information at a whole-genome level. Exploiting this novel information, we perform whole-genome comparative genomics of two P. striiformis f.sp. tritici isolates with contrasting life histories. We compare one isolate of the old European lineage (PstS0), which has been asexual for over 50 years, and a Warrior isolate (PstS7 lineage) from a novel incursion into Europe in 2011 from a sexual population in the Himalayan region. This comparison provides evidence that long-term asexual evolution leads to genome expansion, accumulation of transposable elements, and increased heterozygosity at the single nucleotide, structural, and allele levels. At the whole-genome level, candidate effectors are not compartmentalized and do not exhibit reduced levels of synteny. Yet we were able to identify two subsets of candidate effector populations. About 70% of candidate effectors are invariant between the two isolates, whereas 30% are hypervariable. The latter might be involved in host adaptation on wheat and explain the different phenotypes of the two isolates. Overall, this detailed comparative analysis of two haplotype-aware assemblies of P. striiformis f.sp. tritici is the first step in understanding the evolution of dikaryotic rust fungi at a whole-genome level.

Distinct Life Histories Impact Dikaryotic Genome Evolution in the Rust Fungus Puccinia striiformis Causing Stripe Rust in Wheat

Stripe rust of wheat, caused by the obligate biotrophic fungus Puccinia striiformis f.sp. tritici, is a major threat to wheat production worldwide with an estimated yearly loss of US $1 billion. The recent advances in long-read sequencing technologies and tailored-assembly algorithms enabled us to disentangle the two haploid genomes of Pst. This provides us with haplotype-specific information at a whole-genome level. Exploiting this novel information, we perform whole-genome comparative genomics of two P. striiformis f.sp. tritici isolates with contrasting life histories. We compare one isolate of the old European lineage (PstS0), which has been asexual for over 50 years, and a Warrior isolate (PstS7 lineage) from a novel incursion into Europe in 2011 from a sexual population in the Himalayan region. This comparison provides evidence that long-term asexual evolution leads to genome expansion, accumulation of transposable elements, and increased heterozygosity at the single nucleotide, structural, and allele levels. At the whole-genome level, candidate effectors are not compartmentalized and do not exhibit reduced levels of synteny. Yet we were able to identify two subsets of candidate effector populations. About 70% of candidate effectors are invariant between the two isolates, whereas 30% are hypervariable. The latter might be involved in host adaptation on wheat and explain the different phenotypes of the two isolates. Overall, this detailed comparative analysis of two haplotype-aware assemblies of P. striiformis f.sp. tritici is the first step in understanding the evolution of dikaryotic rust fungi at a whole-genome level.

"CATAStrophy," a Genome-Informed Trophic Classification of Filamentous Plant Pathogens - How Many Different Types of Filamentous Plant Pathogens Are There?

The traditional classification of fungal and oomycete phytopathogens into three classes – biotrophs, hemibiotrophs, or necrotrophs – is unsustainable. This study highlights multiple phytopathogen species for which these labels have been inappropriately applied. We propose a novel and reproducible classification based solely on genome-derived analysis of carbohydrate-active enzyme (CAZyme) gene content called CAZyme-Assisted Training And Sorting of -trophy (CATAStrophy). CATAStrophy defines four major divisions for species associated with living plants. These are monomertrophs (Mo) (corresponding to biotrophs), polymertrophs (P) (corresponding to necrotrophs), mesotrophs (Me) (corresponding to hemibiotrophs), and vasculartrophs (including species commonly described as wilts, rots, or anthracnoses). The Mo class encompasses symbiont, haustorial, and non-haustorial species. Me are divided into the subclasses intracellular and extracellular Me, and the P into broad and narrow host sub-classes. This gives a total of seven discrete plant-pathogenic classes. The classification provides insight into the properties of these species and offers a facile route to develop control measures for newly recognized diseases. Software for CATAStrophy is available online at https://github.com/ccdmb/catastrophy. We present the CATAStrophy method for the prediction of trophic phenotypes based on CAZyme gene content, as a complementary method to the traditional tripartite “biotroph–hemibiotroph–necrotroph” classifications that may encourage renewed investigation and revision within the fungal biology community.

The Current Epidemic of the Barley Pathogen Ramularia collo-cygni Derives from a Population Expansion and Shows Global Admixture

Ramularia leaf spot is becoming an ever-increasing problem in main barley-growing regions since the 1980s, causing up to 70% yield loss in extreme cases. Yet, the causal agent Ramularia collo-cygni remains poorly studied. The diversity of the pathogen in the field thus far remains unknown. Furthermore, it is unknown to what extent the pathogen has a sexual reproductive cycle. The teleomorph of R. collo-cygni has not been observed. To study the genetic diversity of R. collo-cygni and get more insights in its biology, we sequenced the genomes of 19 R. collo-cygni isolates from multiple geographic locations and diverse hosts. Nucleotide polymorphism analyses of all isolates shows that R. collo-cygni is genetically diverse worldwide, with little geographic or host specific differentiation. Next, we used two different methods to detect signals of recombination in our sample set. Both methods find putative recombination events, which indicate that sexual reproduction happens or has happened in the global R. collo-cygni population. Lastly, we used these data on recombination to perform historic population size analyses. These suggest that the effective population size of R. collo-cygni decreased during the domestication of barley and subsequently grew with the rise of agriculture. Our findings deepen our understanding of R. collo-cygni biology and can help us to understand the current epidemic. We discuss how our findings support possible global spread through seed transfer, and we highlight how recombination, clonal spreading, and lack of host specificity could amplify global epidemics of this increasingly important disease and suggest specific approaches to combat the pathogen.