Genome evolution in filamentous plant pathogens: why bigger can be better

Tandem DNA repeats contain cis-regulatory sequences that activate biotrophy-specific expression of Magnaporthe effector gene PWL2

During plant infection, fungi secrete effector proteins in coordination with distinct infection stages. Thus, the success of plant infection is determined by precise control of effector gene expression. We analysed the PWL2 effector gene of the rice blast fungus Magnaporthe oryzae to understand how effector genes are activated specifically during the early biotrophic stages of rice infection. Here, we used confocal live-cell imaging of M. oryzae transformants with various PWL2 promoter fragments fused to sensitive green fluorescent protein reporter genes to determine the expression patterns of PWL2 at the cellular level, together with quantitative reverse transcription PCR analyses at the tissue level. We found PWL2 expression was coupled with sequential biotrophic invasion of rice cells. PWL2 expression was induced in the appressorium upon penetration into a living rice cell but greatly declined in the highly branched hyphae when the first-invaded rice cell was dead. PWL2 expression then increased again as the hyphae penetrate into living adjacent cells. The expression of PWL2 required fungal penetration into living plant cells of either host rice or nonhost onion. Deletion and mutagenesis experiments further revealed that the tandem repeats in the PWL2 promoter contain 12-base pair sequences required for expression. We conclude that PWL2 expression is (a) activated by an unknown signal commonly present in living plant cells, (b) specific to biotrophic stages of fungal infection, and (c) requires 12-base pair cis-regulatory sequences in the promoter.

Genome of the destructive oomycete Phytophthora cinnamomi provides insights into its pathogenicity and adaptive potential

BACKGROUND

Phytophthora cinnamomi is an oomycete pathogen of global relevance. It is considered as one of the most invasive species, which has caused irreversible damage to natural ecosystems and horticultural crops. There is currently a lack of a high-quality reference genome for this species despite several attempts that have been made towards sequencing its genome. The lack of a good quality genome sequence has been a setback for various genetic and genomic research to be done on this species. As a consequence, little is known regarding its genome characteristics and how these contribute to its pathogenicity and invasiveness.

RESULTS

In this work we generated a high-quality genome sequence and annotation for P. cinnamomi using a combination of Oxford Nanopore and Illumina sequencing technologies. The annotation was done using RNA-Seq data as supporting gene evidence. The final assembly consisted of 133 scaffolds, with an estimated genome size of 109.7 Mb, N50 of 1.18 Mb, and BUSCO completeness score of 97.5%. Genome partitioning analysis revealed that P. cinnamomi has a two-speed genome characteristic, similar to that of other oomycetes and fungal plant pathogens. In planta gene expression analysis revealed up-regulation of pathogenicity-related genes, suggesting their important roles during infection and host degradation.

CONCLUSION

This study has provided a high-quality reference genome and annotation for P. cinnamomi. This is among the best assembled genomes for any Phytophthora species assembled to date and thus resulted in improved identification and characterization of pathogenicity-related genes, some of which were undetected in previous versions of genome assemblies. Phytophthora cinnamomi harbours a large number of effector genes which are located in the gene-poor regions of the genome. This unique genomic partitioning provides P. cinnamomi with a high level of adaptability and could contribute to its success as a highly invasive species. Finally, the genome sequence, its annotation and the pathogenicity effectors identified in this study will serve as an important resource that will enable future studies to better understand and mitigate the impact of this important pathogen.

Secreted in Xylem Genes: Drivers of Host Adaptation in Fusarium oxysporum

Fusarium oxysporum (Fo) is a notorious pathogen that significantly contributes to yield losses in crops of high economic status. It is responsible for vascular wilt characterized by the browning of conductive tissue, wilting, and plant death. Individual strains of Fo are host specific (formae speciales), and approximately, 150 forms have been documented so far. The pathogen secretes small effector proteins in the xylem, termed as Secreted in Xylem (Six), that contribute to its virulence. Most of these proteins contain cysteine residues in even numbers. These proteins are encoded by SIX genes that reside on mobile pathogenicity chromosomes. So far, 14 proteins have been reported. However, formae speciales vary in SIX protein profile and their respective gene sequence. Thus, SIX genes have been employed as ideal markers for pathogen identification. Acquisition of SIX-encoding mobile pathogenicity chromosomes by non-pathogenic lines, through horizontal transfer, results in the evolution of new virulent lines. Recently, some SIX genes present on these pathogenicity chromosomes have been shown to be involved in defining variation in host specificity among formae speciales. Along these lines, the review entails the variability (formae speciales, races, and vegetative compatibility groups) and evolutionary relationships among members of F. oxysporum species complex (FOSC). It provides updated information on the diversity, structure, regulation, and (a)virulence functions of SIX genes. The improved understanding of roles of SIX in variability and virulence of Fo has significant implication in establishment of molecular framework and techniques for disease management. Finally, the review identifies the gaps in current knowledge and provides insights into potential research landscapes that can be explored to strengthen the understanding of functions of SIX genes.

Comparative analyses of saprotrophy in Salisapilia sapeloensis and diverse plant pathogenic oomycetes reveal lifestyle-specific gene expression

Oomycetes include many devastating plant pathogens. Across oomycete diversity, plant-infecting lineages are interspersed by non-pathogenic ones. Unfortunately, our understanding of the evolution of lifestyle switches is hampered by a scarcity of data on the molecular biology of saprotrophic oomycetes, ecologically important primary colonizers of dead tissue that can serve as informative reference points for understanding the evolution of pathogens. Here, we established Salisapilia sapeloensis as a tractable system for the study of saprotrophic oomycetes. We generated multiple transcriptomes from S. sapeloensis and compared them with (i) 22 oomycete genomes and (ii) the transcriptomes of eight pathogenic oomycetes grown under 13 conditions. We obtained a global perspective on gene expression signatures of oomycete lifestyles. Our data reveal that oomycete saprotrophs and pathogens use similar molecular mechanisms for colonization but exhibit distinct expression patterns. We identify a S. sapeloensis-specific array and expression of carbohydrate-active enzymes and putative regulatory differences, highlighted by distinct expression levels of transcription factors. Salisapilia sapeloensis expresses only a small repertoire of candidates for virulence-associated genes. Our analyses suggest lifestyle-specific gene regulatory signatures and that, in addition to variation in gene content, shifts in gene regulatory networks underpin the evolution of oomycete lifestyles.

Two new azaphilone pigments from Talaromyces albobiverticillius and their anti-inflammatory activity

Two new azaphilone pigments, talaralbols A and B (3 and 7), along with five known azaphilone metabolites (1, 2, and 4-6), were isolated from the culture of Talaromyces albobiverticillius associated with the isopod Armadillidium vulgare. Their structures were elucidated by a combination of 1 D and 2 D NMR data, ECD calculations, chemical transformations, and NMR data analogy with model compounds. Talaralbol A (3) showed a moderate inhibition on the lipopolysaccharide (LPS)-induced nitric oxide (NO) production in RAW264.7 cells with the inhibitory rate being 31.0% at the concentration of 10 μM.[Formula: see text].

Soybean anthracnose caused by Colletotrichum species: Current status and future prospects

Soybean (Glycine max) is one of the most important cultivated plants worldwide as a source of protein-rich foods and animal feeds. Anthracnose, caused by different lineages of the hemibiotrophic fungus Colletotrichum, is one of the main limiting factors to soybean production. Losses due to anthracnose have been neglected, but their impact may threaten up to 50% of the grain production.

TAXONOMY

While C. truncatum is considered the main species associated with soybean anthracnose, recently other species have been reported as pathogenic on this host. Until now, it has not been clear whether the association of new Colletotrichum species with the disease is related to emerging species or whether it is due to the undergoing changes in the taxonomy of the genus.

DISEASE SYMPTOMS

Typical anthracnose symptoms are pre- and postemergence damping-off; dark, depressed, and irregular spots on cotyledons, stems, petioles, and pods; and necrotic laminar veins on leaves that can result in premature defoliation. Symptoms may evolve to pod rot, immature opening of pods, and premature germination of grains.

CHALLENGES

As accurate species identification of the causal agent is decisive for disease control and prevention, in this work we review the taxonomic designation of Colletotrichum isolated from soybean to understand which lineages are pathogenic on this host. We also present a comprehensive literature review of soybean anthracnose, focusing on distribution, symptomatology, epidemiology, disease management, identification, and diagnosis. We consider the knowledge emerging from population studies and comparative genomics of Colletotrichum spp. associated with soybean providing future perspectives in the identification of molecular factors involved in the pathogenicity process.

USEFUL WEBSITE

Updates on Colletotrichum can be found at http://www.colletotrichum.org/. All available Colletotrichum genomes on GenBank can be viewed at http://www.colletotrichum.org/genomics/.

What Do We Know about Botryosphaeriaceae? An Overview of a Worldwide Cured Dataset

Botryosphaeriaceae-related diseases occur worldwide in a wide variety of plant hosts. The number of studies targeting the distribution, diversity, ecology, and pathogenicity of Botryosphaeriaceae species are consistently increasing. However, with the lack of consistency in species delimitation, the name of hosts, and the locations of studies, it is almost impossible to quantify the presence of these species worldwide, or the number of different host–fungus interactions that occur. In this review, we collected and organized Botryosphaeriaceae occurrences in a single cured dataset, allowing us to obtain for the first time a complete perspective on species’ global diversity, dispersion, host association, ecological niches, pathogenicity, communication efficiency of new occurrences, and new host–fungus associations. This dataset is freely available through an interactive and online application. The current release (version 1.0) contains 14,405 cured isolates and 2989 literature references of 12,121 different host–fungus interactions with 1692 different plant species from 149 countries.

Recent loss of the Dim2 DNA methyltransferase decreases mutation rate in repeats and changes evolutionary trajectory in a fungal pathogen

DNA methylation is found throughout all domains of life, yet the extent and function of DNA methylation differ among eukaryotes. Strains of the plant pathogenic fungus Zymoseptoria tritici appeared to lack cytosine DNA methylation (5mC) because gene amplification followed by Repeat-Induced Point mutation (RIP) resulted in the inactivation of the dim2 DNA methyltransferase gene. 5mC is, however, present in closely related sister species. We demonstrate that inactivation of dim2 occurred recently as some Z. tritici isolates carry a functional dim2 gene. Moreover, we show that dim2 inactivation occurred by a different path than previously hypothesized. We mapped the genome-wide distribution of 5mC in strains with or without functional dim2 alleles. Presence of functional dim2 correlates with high levels of 5mC in transposable elements (TEs), suggesting a role in genome defense. We identified low levels of 5mC in strains carrying non-functional dim2 alleles, suggesting that 5mC is maintained over time, presumably by an active Dnmt5 DNA methyltransferase. Integration of a functional dim2 allele in strains with mutated dim2 restored normal 5mC levels, demonstrating de novo cytosine methylation activity of Dim2. To assess the importance of 5mC for genome evolution, we performed an evolution experiment, comparing genomes of strains with high levels of 5mC to genomes of strains lacking functional dim2. We found that presence of a functional dim2 allele alters nucleotide composition by promoting C to T transitions (C→T) specifically at CpA (CA) sites during mitosis, likely contributing to TE inactivation. Our results show that 5mC density at TEs is a polymorphic trait in Z. tritici populations that can impact genome evolution.

Genomic rearrangements generate hypervariable mini-chromosomes in host-specific isolates of the blast fungus

Supernumerary mini-chromosomes–a unique type of genomic structural variation–have been implicated in the emergence of virulence traits in plant pathogenic fungi. However, the mechanisms that facilitate the emergence and maintenance of mini-chromosomes across fungi remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), mini-chromosomes have been first described in the early 1990s but, until very recently, have been overlooked in genomic studies. Here we investigated structural variation in four isolates of the blast fungus M. oryzae from different grass hosts and analyzed the sequences of mini-chromosomes in the rice, foxtail millet and goosegrass isolates. The mini-chromosomes of these isolates turned out to be highly diverse with distinct sequence composition. They are enriched in repetitive elements and have lower gene density than core-chromosomes. We identified several virulence-related genes in the mini-chromosome of the rice isolate, including the virulence-related polyketide synthase Ace1 and two variants of the effector gene AVR-Pik. Macrosynteny analyses around these loci revealed structural rearrangements, including inter-chromosomal translocations between core- and mini-chromosomes. Our findings provide evidence that mini-chromosomes emerge from structural rearrangements and segmental duplication of core-chromosomes and might contribute to adaptive evolution of the blast fungus.

The genomes of plant pathogens often exhibit an architecture that facilitates high rates of dynamic rearrangements and genetic diversification in virulence associated regions. These regions, which tend to be gene sparse and repeat rich, are thought to serve as a cradle for adaptive evolution. Supernumerary chromosomes, i.e. chromosomes that are only present in some but not all individuals of a species, are a special type of structural variation that have been observed in plants, animals, and fungi. Here we identified and studied supernumerary mini-chromosomes in the blast fungus Magnaporthe oryzae, a pathogen that causes some of the most destructive plant diseases. We found that rice, foxtail millet and goosegrass isolates of this pathogen contain mini-chromosomes with distinct sequence composition. All mini-chromosomes are rich in repetitive genetic elements and have lower gene densities than core-chromosomes. Further, we identified virulence-related genes on the mini-chromosome of the rice isolate. We observed large-scale genomic rearrangements around these loci, indicative of a role of mini-chromosomes in facilitating genome dynamics. Taken together, our results indicate that mini-chromosomes contribute to genome rearrangements and possibly adaptive evolution of the blast fungus.

Near-Chromosome-Level Genome Assembly of the Dark Septate Endophyte Laburnicola rhizohalophila: A Model for Investigating Root-Fungus Symbiosis

The novel DSE Laburnicola rhizohalophila (Pleosporales, Ascomycota) is frequently found in the halophytic seepweed (Suaeda salsa). In this article, we report a near-chromosome-level hybrid assembly of this fungus using a combination of short-read Illumina data to polish assemblies generated from long-read Nanopore data. The reference genome for L. rhizohalophila was assembled into 26 scaffolds with a total length of 64.0 Mb and a N50 length of 3.15 Mb. Of them, 17 scaffolds approached the length of intact chromosomes, and 5 had telomeres at one end only. A total of 10,891 gene models were predicted. Intriguingly, 27.5 Mb of repeat sequences that accounted for 42.97% of the genome was identified, and long terminal repeat retrotransposons were the most frequent known transposable elements, indicating that transposable element proliferation contributes to its increased genome size. BUSCO analyses using the Fungi_odb10 data set showed that 95.0% of genes were complete. In addition, 292 carbohydrate active enzymes, 33 secondary metabolite clusters, and 84 putative effectors were identified in silico. The resulting high-quality assembly and genome features are not only an important resource for further research on understanding the mechanism of root-fungi symbiotic interactions but will also contribute to comparative analyses of genome biology and evolution within Pleosporalean species.

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

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

Lineage-Specific Genes and Cryptic Sex: Parallels and Differences between Arbuscular Mycorrhizal Fungi and Fungal Pathogens

Arbuscular mycorrhizal fungi (AMF) live as obligate root symbionts on almost all land plants. They have long been regarded as ancient asexuals that have propagated clonally for millions of years. However, genomic studies in Rhizophagus irregularis and other AMF revealed many features indicative of sex. Surprisingly, comparative genomics of conspecific isolates of R. irregularis revealed an unexpected interstrain diversity, suggesting that AMF carry a high number of lineage-specific (LS) genes. Intriguingly, cryptic sex and LS genomic regions have previously been reported in a number of fungal pathogens of plants and humans. Here, we discuss these genomic similarities and highlight their potential relevance for AMF adaptation to the environment and for symbiotic functioning.

Silent control: microbial plant pathogens evade host immunity without coding sequence changes

Both animals and plants have evolved a robust immune system to surveil and defeat invading pathogenic microbes. Evasion of host immune surveillance is the key for pathogens to initiate successful infection. To evade the host immunity, plant pathogens evolved a variety of strategies such as masking themselves from host immune recognitions, blocking immune signaling transductions, reprogramming immune responses and adapting to immune microenvironmental changes. Gain of new virulence genes, sequence and structural variations enables plant pathogens to evade host immunity through changes in the genetic code. However, recent discoveries demonstrated that variations at the transcriptional, post-transcriptional, post-translational and glycome level enable pathogens to cope with the host immune system without coding sequence changes. The biochemical modification of pathogen associated molecular patterns and silencing of effector genes emerged as potent ways for pathogens to hide from host recognition. Altered processing in mRNA activities provide pathogens with resilience to microenvironment changes. Importantly, these hiding variants are directly or indirectly modulated by catalytic enzymes or enzymatic complexes and cannot be revealed by classical genomics alone. Unveiling these novel host evasion mechanisms in plant pathogens enables us to better understand the nature of plant disease and pinpoints strategies for rational diseases management in global food protection.

The Impact of Blast Disease: Past, Present, and Future

Rice blast disease is both the most explosive and potentially damaging disease of the world's rice (Oryza sativa) crop and a model system for research on the molecular mechanisms that fungi use to cause plant disease. The blast fungus, Magnaporthe oryzae, is highly evolved to sense when it is on a leaf surface; to develop a pressurized cell, the appressorium, to punch through the leaf cuticle; and then to hijack living rice cells to assist it in causing disease. Host specificity, determining which plants particular fungal strains can infect, is also an important topic for research. The blast fungus is a moving target, quickly overcoming rice resistance genes we deploy to control it, and recently emerging to cause devastating disease on an entirely new cereal crop, wheat. M. oryzae is highly adaptable, with multiple examples of genetic instability at certain gene loci and in certain genomic regions. Understanding the biology of the fungus in the field, and its potential for genetic and genome variability, is key to keep it from adapting to life in the research laboratory and losing relevance to the significant impact it has on global food security.

A unique chromatin profile defines adaptive genomic regions in a fungal plant pathogen

Genomes store information at scales beyond the linear nucleotide sequence, which impacts genome function at the level of an individual, while influences on populations and long-term genome function remains unclear. Here, we addressed how physical and chemical DNA characteristics influence genome evolution in the plant pathogenic fungus Verticillium dahliae. We identified incomplete DNA methylation of repetitive elements, associated with specific genomic compartments originally defined as Lineage-Specific (LS) regions that contain genes involved in host adaptation. Further chromatin characterization revealed associations with features such as H3 Lys-27 methylated histones (H3K27me3) and accessible DNA. Machine learning trained on chromatin data identified twice as much LS DNA as previously recognized, which was validated through orthogonal analysis, and we propose to refer to this DNA as adaptive genomic regions. Our results provide evidence that specific chromatin profiles define adaptive genomic regions, and highlight how different epigenetic factors contribute to the organization of these regions.

Partial pathogenicity chromosomes in Fusarium oxysporum are sufficient to cause disease and can be horizontally transferred

In Fusarium oxysporum f.sp. lycopersici, all effector genes reported so far - also called SIX genes - are located on a single accessory chromosome which is required for pathogenicity and can also be horizontally transferred to another strain. To narrow down the minimal region required for virulence, we selected partial pathogenicity chromosome deletion strains by fluorescence-assisted cell sorting of a strain in which the two arms of the pathogenicity chromosome were labelled with GFP and RFP respectively. By testing the virulence of these deletion mutants, we show that the complete long arm and part of the short arm of the pathogenicity chromosome are not required for virulence. In addition, we demonstrate that smaller versions of the pathogenicity chromosome can also be transferred to a non-pathogenic strain and they are sufficient to turn the non-pathogen into a pathogen. Surprisingly, originally non-pathogenic strains that had received a smaller version of the pathogenicity chromosome were much more aggressive than recipients with a complete pathogenicity chromosome. Whole genome sequencing analysis revealed that partial deletions of the pathogenicity chromosome occurred mainly close to repeats, and that spontaneous duplication of sequences in accessory regions is frequent both in chromosome deletion strains and in horizontal transfer strains.

The Effect of Virulence and Resistance Mechanisms on the Interactions between Parasitic Plants and Their Hosts

Parasitic plants have a unique heterotrophic lifestyle based on the extraction of water and nutrients from host plants. Some parasitic plant species, particularly those of the family Orobanchaceae, attack crops and cause substantial yield losses. The breeding of resistant crop varieties is an inexpensive way to control parasitic weeds, but often does not provide a long-lasting solution because the parasites rapidly evolve to overcome resistance. Understanding mechanisms underlying naturally occurring parasitic plant resistance is of great interest and could help to develop methods to control parasitic plants. In this review, we describe the virulence mechanisms of parasitic plants and resistance mechanisms in their hosts, focusing on obligate root parasites of the genera Orobanche and Striga. We noticed that the resistance (R) genes in the host genome often encode proteins with nucleotide-binding and leucine-rich repeat domains (NLR proteins), hence we proposed a mechanism by which host plants use NLR proteins to activate downstream resistance gene expression. We speculated how parasitic plants and their hosts co-evolved and discussed what drives the evolution of virulence effectors in parasitic plants by considering concepts from similar studies of plant-microbe interaction. Most previous studies have focused on the host rather than the parasite, so we also provided an updated summary of genomic resources for parasitic plants and parasitic genes for further research to test our hypotheses. Finally, we discussed new approaches such as CRISPR/Cas9-mediated genome editing and RNAi silencing that can provide deeper insight into the intriguing life cycle of parasitic plants and could potentially contribute to the development of novel strategies for controlling parasitic weeds, thereby enhancing crop productivity and food security globally.

Whole genome analysis of the koa wilt pathogen (Fusarium oxysporum f. sp. koae) and the development of molecular tools for early detection and monitoring

BACKGROUND

Development and application of DNA-based methods to distinguish highly virulent isolates of Fusarium oxysporum f. sp. koae [Fo koae; cause of koa wilt disease on Acacia koa (koa)] will help disease management through early detection, enhanced monitoring, and improved disease resistance-breeding programs.

RESULTS

This study presents whole genome analyses of one highly virulent Fo koae isolate and one non-pathogenic F. oxysporum (Fo) isolate. These analyses allowed for the identification of putative lineage-specific DNA and predicted genes necessary for disease development on koa. Using putative chromosomes and predicted gene comparisons, Fo koae-exclusive, virulence genes were identified. The putative lineage-specific DNA included identified genes encoding products secreted in xylem (e. g., SIX1 and SIX6) that may be necessary for disease development on koa. Unique genes from Fo koae were used to develop pathogen-specific PCR primers. These diagnostic primers allowed target amplification in the characterized highly virulent Fo koae isolates but did not allow product amplification in low-virulence or non-pathogenic isolates of Fo. Thus, primers developed in this study will be useful for early detection and monitoring of highly virulent strains of Fo koae. Isolate verification is also important for disease resistance-breeding programs that require a diverse set of highly virulent Fo koae isolates for their disease-screening assays to develop disease-resistant koa.

CONCLUSIONS

These results provide the framework for understanding the pathogen genes necessary for koa wilt disease and the genetic variation of Fo koae populations across the Hawaiian Islands.

The cytosine methylation landscape of spring barley revealed by a new reduced representation bisulfite sequencing pipeline, WellMeth

Patterns and level of cytosine methylation vary widely among plant species and are associated with genome size as well as the proportion of transposons and other repetitive elements in the genome. We explored epigenetic patterns and diversity in a representative proportion of the spring barley (Hordeum vulgare L.) genome across several commercial and historical cultivars. This study adapted a genotyping-by-sequencing (GBS) approach for the detection of methylated cytosines in genomic DNA. To analyze the data, we developed WellMeth, a complete pipeline for analysis of reduced representation bisulfite sequencing. WellMeth enabled quantification of context-specific DNA methylation at the single-base resolution as well as identification of differentially methylated sites (DMCs) and regions (DMRs). On average, DNA methylation levels were significantly higher than what is commonly observed in many plants species, reaching over 10-fold higher levels than those in Arabidopsis thaliana (L.) Heynh. in the CHH methylation. Preferential methylation was observed within and at the edges of long-terminal repeats (LTR) retrotransposons Gypsy and Copia. From a pairwise comparison of cultivars, numerous DMRs could be identified of which more than 5,000 were conserved within the analyzed set of barley cultivars. The subset of regions overlapping with genes showed enrichment in gene ontology (GO) categories associated with chromatin and cellular structure and organization. A significant correlation between genetic and epigenetic distances suggests that a considerable portion of methylated regions is under strict genetic control in barley. The data presented herein represents the first step in efforts toward a better understanding of genome-level structural and functional aspects of methylation in barley.

Horizontal Gene Transfer and Tandem Duplication Shape the Unique CAZyme Complement of the Mycoparasitic Oomycetes Pythium oligandrum and Pythium periplocum

Crop protection strategies that are effective but that reduce our reliance on chemical pesticides are urgently needed to meet the UN sustainable development goals for global food security. Mycoparasitic oomycetes such as Pythium oligandrum and Pythium periplocum, have potential for the biological control of plant diseases that threaten crops and have attracted much attention due to their abilities to antagonize plant pathogens and modulate plant immunity. Studies of the molecular and genetic determinants of mycoparasitism in these species have been less well developed than those of their fungal counterparts. Carbohydrate-active enzymes (CAZymes) from P. oligandrum and P. periplocum are predicted to be important components of mycoparasitism, being involved in the degradation of the cell wall of their oomycete and fungal prey species. To explore the evolution of CAZymes of these species we performed an in silico identification and comparison of the full CAZyme complement (CAZyome) of the two mycoparasitic Pythium species (P. oligandrum and P. periplocum), with seven other Pythium species, and four Phytophthora species. Twenty CAZy gene families involved in the degradation of cellulose, hemicellulose, glucan, and chitin were expanded in, or unique to, mycoparasitic Pythium species and several of these genes were expressed during mycoparasitic interactions with either oomycete or fungal prey, as revealed by RNA sequencing and quantitative qRT-PCR. Genes from three of the cellulose and chitin degrading CAZy families (namely AA9, GH5_14, and GH19) were expanded via tandem duplication and predominantly located in gene sparse regions of the genome, suggesting these enzymes are putative pathogenicity factors able to undergo rapid evolution. In addition, five of the CAZy gene families were likely to have been obtained from other microbes by horizontal gene transfer events. The mycoparasitic species are able to utilize complex carbohydrates present in fungal cell walls, namely chitin and N-acetylglucosamine for growth, in contrast to their phytopathogenic counterparts. Nonetheless, a preference for the utilization of simple sugars for growth appears to be a common trait within the oomycete lineage.

Comparative Genomic Analysis of Dactylonectria torresensis Strains from Grapevine, Soil and Weed Highlights Potential Mechanisms in Pathogenicity and Endophytic Lifestyle

The soil-borne fungus Dactylonectria torresensis is the most common causal agent of black-foot disease in Europe. However, there is a lack of understanding on how this fungus can provoke plant symptoms. In this study, we sequenced, annotated and analyzed the genomes of three isolates of D. torresensis collected from asymptomatic vine, weed and soil. Sequenced genomes were further compared to those of 27 fungal species including root and aerial pathogens, white rot degraders, indoor biodeterioration agents, saprotrophs, dark septate endophytes and mycorrhiza. Strains of D. torresensis present genomes with between 64 and 65 Mbp and with up to 18,548 predicted genes for each strain. Average Nucleotide Identity (ANI) shows that strains are different according to genome contents. Clusters of orthologous groups were compared, and clusters of genes related to necroses were particularly detected in all strains of D. torresensis (necrosis inducing peptides and proteins, and ethylene inducing peptides) as well as several genes involved in resistance against fungicides frequently used in viticulture such as copper. Interestingly, an expanded high number of genes related to carbohydrate-active enzymes were detected in each Dactylonectria strain, especially those related to glycoside hydrolases that could be involved in penetration of plant tissues or pathogenicity. An increased number of candidate genes for CAZyme classes AA9 and AA3-1 supports the ability of strains to efficiently degrade plant material. High numbers of genes of D. torresensis related to secretome and small secreted proteins were further characterized. Moreover, the presence of several gene clusters such as fujikurin-like genes was detected and were normally found in Fusariumfujikuroi, that have been linked to fungal pathogenicity. The phenotypes of the three strains investigated showed further difference in light response. We found that Dactylonectria strains have an increased number of photoreceptor encoding genes and we showed sequence alterations. Altogether, the results highlight several gene clusters present in D. torresensis strains that could be linked to endophytic lifestyle, pathogenicity, plant maceration and degradation of plant tissues as well as adaptation to soil contaminated with metals and metalloids and light response.

Threat to Asian wild apple trees posed by gene flow from domesticated apple trees and their "pestified" pathogens

Secondary contact between crops and their wild relatives poses a threat to wild species, not only through gene flow between plants, but also through the dispersal of crop pathogens and genetic exchanges involving these pathogens, particularly those that have become more virulent by indirect selection on resistant crops, a phenomenon known as "pestification." Joint analyses of wild and domesticated hosts and their pathogens are essential to address this issue, but such analyses remain rare. We used population genetics approaches, demographic inference and pathogenicity tests on host-pathogen pairs of wild or domesticated apple trees from Central Asia and their main fungal pathogen, Venturia inaequalis, which itself has differentiated agricultural and wild-type populations. We confirmed the occurrence of gene flow from cultivated (Malus domestica) to wild (Malus sieversii) apple trees in Asian forests, potentially threatening the persistence of Asian wild apple trees. Pathogenicity tests demonstrated the pestification of V. inaequalis, the agricultural-type population being more virulent on both wild and domesticated trees. Single nucleotide polymorphism (SNP) markers and the demographic modelling of pathogen populations revealed hybridization following secondary contact between agricultural and wild-type fungal populations, and dispersal of the agricultural-type pathogen population in wild forests, increasing the threat of disease in the wild apple species. We detected an SNP potentially involved in pathogen pestification, generating an early stop codon in a gene encoding a small secreted protein in the agricultural-type fungal population. Our findings, based on joint analyses of paired host and pathogen data sets, highlight the threat posed by cultivating a crop near its centre of origin, in terms of pestified pathogen invasions in wild plant populations and introgression in the wild-type pathogen population.

Comparative Genomics Analyses of Lifestyle Transitions at the Origin of an Invasive Fungal Pathogen in the Genus Cryphonectria

Forest and agroecosystems, as well as animal and human health, are threatened by emerging pathogens. Following decimation of chestnuts in the United States, the fungal pathogen Cryphonectria parasitica colonized Europe. After establishment, the pathogen population gave rise to a highly successful lineage that spread rapidly across the continent. Core to our understanding of what makes a successful pathogen is the genetic repertoire enabling the colonization and exploitation of host species. Here, we have assembled >100 genomes across two related genera to identify key genomic determinants leading to the emergence of chestnut blight. We found subtle yet highly specific changes in the transition from saprotrophy to latent pathogenicity mostly determined by enzymes involved in carbohydrate metabolism. Large-scale genomic analyses of genes underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

Emerging fungal pathogens are a threat to forest and agroecosystems, as well as animal and human health. How pathogens evolve from nonpathogenic ancestors is still poorly understood, making the prediction of future outbreaks challenging. Most pathogens have evolved lifestyle adaptations, which were enabled by specific changes in the gene content of the species. Hence, understanding transitions in the functions encoded by genomes gives valuable insight into the evolution of pathogenicity. Here, we studied lifestyle evolution in the genus Cryphonectria, including the prominent invasive pathogen Cryphonectria parasitica, the causal agent of chestnut blight on Castanea species. We assembled and compared the genomes of pathogenic and putatively nonpathogenic Cryphonectria species, as well as sister group pathogens in the family Cryphonectriaceae (Diaporthales, Ascomycetes), to investigate the evolution of genome size and gene content. We found a striking loss of genes associated with carbohydrate metabolism (CAZymes) in C. parasitica compared to other Cryphonectriaceae. Despite substantial CAZyme gene loss, experimental data suggest that C. parasitica has retained wood colonization abilities shared with other Cryphonectria species. Putative effectors substantially varied in number, cysteine content, and protein length among species. In contrast, secondary metabolite gene clusters show a high degree of conservation within the genus. Overall, our results underpin the recent lifestyle transition of C. parasitica toward a more pathogenic lifestyle. Our findings suggest that a CAZyme loss may have promoted pathogenicity of C. parasitica on Castanea species. Analyzing gene complements underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

IMPORTANCE Forest and agroecosystems, as well as animal and human health, are threatened by emerging pathogens. Following decimation of chestnuts in the United States, the fungal pathogen Cryphonectria parasitica colonized Europe. After establishment, the pathogen population gave rise to a highly successful lineage that spread rapidly across the continent. Core to our understanding of what makes a successful pathogen is the genetic repertoire enabling the colonization and exploitation of host species. Here, we have assembled >100 genomes across two related genera to identify key genomic determinants leading to the emergence of chestnut blight. We found subtle yet highly specific changes in the transition from saprotrophy to latent pathogenicity mostly determined by enzymes involved in carbohydrate metabolism. Large-scale genomic analyses of genes underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.IMPORTANCE Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

Proteinaceous effector discovery and characterization in filamentous plant pathogens

The complicated interplay of plant-pathogen interactions occurs on multiple levels as pathogens evolve to constantly evade the immune responses of their hosts. Many economically important crops fall victim to filamentous pathogens that produce small proteins called effectors to manipulate the host and aid infection/colonization. Understanding the effector repertoires of pathogens is facilitating an increased understanding of the molecular mechanisms underlying virulence as well as guiding the development of disease control strategies. The purpose of this review is to give a chronological perspective on the evolution of the methodologies used in effector discovery from physical isolation and in silico predictions, to functional characterization of the effectors of filamentous plant pathogens and identification of their host targets.

De novo sequencing, assembly and functional annotation of Armillaria borealis genome

BACKGROUND

Massive forest decline has been observed almost everywhere as a result of negative anthropogenic and climatic effects, which can interact with pests, fungi and other phytopathogens and aggravate their effects. Climatic changes can weaken trees and make fungi, such as Armillaria more destructive. Armillaria borealis (Marxm. & Korhonen) is a fungus from the Physalacriaceae family (Basidiomycota) widely distributed in Eurasia, including Siberia and the Far East. Species from this genus cause the root white rot disease that weakens and often kills woody plants. However, little is known about ecological behavior and genetics of A. borealis. According to field research data, A. borealis is less pathogenic than A. ostoyae, and its aggressive behavior is quite rare. Mainly A. borealis behaves as a secondary pathogen killing trees already weakened by other factors. However, changing environment might cause unpredictable effects in fungus behavior.

RESULTS

The de novo genome assembly and annotation were performed for the A. borealis species for the first time and presented in this study. The A. borealis genome assembly contained ~ 68 Mbp and was comparable with ~ 60 and ~ 79.5 Mbp for the A. ostoyae and A. mellea genomes, respectively. The N50 for contigs equaled 50,544 bp. Functional annotation analysis revealed 21,969 protein coding genes and provided data for further comparative analysis. Repetitive sequences were also identified. The main focus for further study and comparative analysis will be on the enzymes and regulatory factors associated with pathogenicity.

CONCLUSIONS

Pathogenic fungi such as Armillaria are currently one of the main problems in forest conservation. A comprehensive study of these species and their pathogenicity is of great importance and needs good genomic resources. The assembled genome of A. borealis presented in this study is of sufficiently good quality for further detailed comparative study on the composition of enzymes in other Armillaria species. There is also a fundamental problem with the identification and classification of species of the Armillaria genus, where the study of repetitive sequences in the genomes of basidiomycetes and their comparative analysis will help us identify more accurately taxonomy of these species and reveal their evolutionary relationships.

Transcriptome Profiling Data of Botrytis cinerea Infection on Whole Plant Solanum lycopersicum

Botrytis cinerea is a foliar necrotrophic fungal-pathogen capable of infecting >580 genera of plants, is often used as model organism for studying fungal-host interactions. We used RNAseq to study transcriptome of B. cinerea infection on a major (worldwide) vegetable crop, tomato (Solanum lycopersicum). Most previous works explored only few infection stages, using RNA extracted from entire leaf-organ diluting the expression of studied infected region. Many studied B. cinerea infection, on detached organs assuming that similar defense/physiological reactions occurs in the intact plant. We analyzed transcriptome of the pathogen and host in 5 infection stages of whole-plant leaves at the infection site. We supply high quality, pathogen-enriched gene count that facilitates future research of the molecular processes regulating the infection process.

Novel Variation and Evolution of AvrPiz-t of Magnaporthe oryzae in Field Isolates

The product of the avirulence (Avr) gene of Magnaporthe oryzae can be detected by the product of the corresponding resistance (R) gene of rice and activates immunity to rice mediated by the R gene. The high degree of variability of M. oryzae isolates in pathogenicity makes the control of rice blast difficult. That resistance of the R gene in rice has been lost has been ascribed to the instability of the Avr gene in M. oryzae. Further study on the variation of the Avr genes in M. oryze field isolates may yield valuable information on the durable and effective deployment of R genes in rice production areas. AvrPiz-t and Piz-t are a pair of valuable genes in the Rice-Magnaporthe pathosystem. AvrPiz-t is detectable by Piz-t and determines the effectiveness of Piz-t. To effectively deploy the R gene Piz-t, the distribution, variation, and evolution of the corresponding Avr gene AvrPiz-t were found among 312 M. oryzae isolates collected from Yunnan rice production areas of China. PCR amplification and pathogenicity assays of AvrPiz-t showed that 202 isolates (64.7%) held AvrPiz-t alleles and were avirulent to IRBLzt-T (holding Piz-t). There were 42.3–83.3% avirulent isolates containing AvrPiz-t among seven regions in Yunnan Province. Meanwhile, 11 haplotypes of AvrPiz-t encoding three novel AvrPiz-t variants were identified among 100 isolates. A 198 bps insertion homologous to solo-LTR of the retrotransposon inago2 in the promoter region of AvrPiz-t in one isolate and a frameshift mutation of CDS in another isolate were identified among 100 isolates, and those two isolates had evolved to virulent from avirulent. Synonymous mutation and non-AUG-initiated N-terminal extensions keeps the AvrPiz-t gene avirulence function in M. oryzae field isolates in Yunnan. A haplotype network showed that H3 was an ancestral haplotype. Structure variance for absence (28.2%) or partial fragment loss (71.8%) of AvrPiz-t was found among 39 virulent isolates and may cause the AvrPiz-t avirulence function to be lost. Overall, AvrPiz-t evolved to virulent from avirulent forms via point mutation, retrotransposon, shift mutation, and structure variance under field conditions.

Host-specificity factors in plant pathogenic fungi

Fortunately, no fungus can cause disease on all plant species, and although some plant-pathogenic fungi have quite a broad host range, most are highly limited in the range of plant species or even cultivars that they cause disease in. The mechanisms of host specificity have been extensively studied in many plant-pathogenic fungi, especially in fungal pathogens causing disease on economically important crops. Specifically, genes involved in host specificity have been identified during the last few decades. In this overview, we describe and discuss these host-specificity genes. These genes encode avirulence (Avr) proteins, proteinaceous host-specific toxins or secondary metabolites. We discuss the genomic context of these genes, their expression, polymorphism, horizontal transfer and involvement in pathogenesis.

Genomic and transcriptomic insights into Raffaelea lauricola pathogenesis

BACKGROUND

Laurel wilt caused by Raffaelea lauricola is a lethal vascular disease of North American members of the Lauraceae plant family. This fungus and its primary ambrosia beetle vector Xyleborus glabratus originated from Asia; however, there is no report of laurel wilt causing widespread mortality on native Lauraceae trees in Asia. To gain insight into why R. lauricola is a tree-killing plant pathogen in North America, we generated and compared high quality draft genome assemblies of R. lauricola and its closely related non-pathogenic species R. aguacate.

RESULTS

Relative to R. aguacate, the R. lauricola genome uniquely encodes several small-secreted proteins that are associated with virulence in other pathogens and is enriched in secondary metabolite biosynthetic clusters, particularly polyketide synthase (PKS), non-ribosomal peptide synthetase (NRPS) and PKS-NRPS anchored gene clusters. The two species also exhibit significant differences in secreted proteins including CAZymes that are associated with polysaccharide binding including the chitin binding CBM50 (LysM) domain. Transcriptomic comparisons of inoculated redbay trees and in vitro-grown fungal cultures further revealed a number of secreted protein genes, secondary metabolite clusters and alternative sulfur uptake and assimilation pathways that are coordinately up-regulated during infection.

CONCLUSIONS

Through these comparative analyses we have identified potential adaptations of R. lauricola that may enable it to colonize and cause disease on susceptible hosts. How these adaptations have interacted with co-evolved hosts in Asia, where little to no disease occurs, and non-co-evolved hosts in North America, where lethal wilt occurs, requires additional functional analysis of genes and pathways.

The Effector Repertoire of the Hop Downy Mildew Pathogen Pseudoperonospora humuli

Pseudoperonospora humuli is an obligate biotrophic oomycete that causes downy mildew (DM), one of the most destructive diseases of cultivated hop that can lead to 100% crop loss in susceptible cultivars. We used the published genome of P. humuli to predict the secretome and effectorome and analyze the transcriptome variation among diverse isolates and during infection of hop leaves. Mining the predicted coding genes of the sequenced isolate OR502AA of P. humuli revealed a secretome of 1,250 genes. We identified 296 RXLR and RXLR-like effector-encoding genes in the secretome. Among the predicted RXLRs, there were several WY-motif-containing effectors that lacked canonical RXLR domains. Transcriptome analysis of sporangia from 12 different isolates collected from various hop cultivars revealed 754 secreted proteins and 201 RXLR effectors that showed transcript evidence across all isolates with reads per kilobase million (RPKM) values > 0. RNA-seq analysis of OR502AA-infected hop leaf samples at different time points after infection revealed highly expressed effectors that may play a relevant role in pathogenicity. Quantitative RT-PCR analysis confirmed the differential expression of selected effectors. We identified a set of P. humuli core effectors that showed transcript evidence in all tested isolates and elevated expression during infection. These effectors are ideal candidates for functional analysis and effector-assisted breeding to develop DM resistant hop cultivars.

Composition and biosynthetic machinery of the Blumeria graminis f. sp. hordei conidia cell wall

Infection of barley with the powdery mildew causal agent, Blumeria graminis f. sp. hordei (Bgh), can lead to devastating damage to barley crops. The recent emergence of fungicide resistance imposes a need to develop new antifungal strategies. The enzymes involved in cell wall biosynthesis are ideal targets for the development of fungicides. However, in order to narrow down any target proteins involved in cell wall formation, a greater understanding of the cell wall structure and composition is required. Here, we present a detailed carbohydrate analysis of the Bgh conidial cell wall, a full annotation of Carbohydrate Active enZymes (CAZy) in the Bgh genome, and a comprehensive expression profile of the genes involved in cell wall metabolism. Glycosidic linkage analysis has revealed that the cell wall polysaccharide fraction of Bgh conidia predominantly consists of glucosyl residues (63.1%) and has a greater proportion of galactopyranosyl residues compared to other species (8.5%). Trace amounts of xylosyl residues were also detected, which is unusual in ascomycetes. Transcripts of the genes involved in cell wall metabolism show high expression of chitin deacetylases, which assist fungi in evading the host defence system by deacetylating chitin to chitosan. The data presented suggest that the cell wall components of the conidia and the corresponding obligate biotrophic CAZy gene profile play a key role in the infection process.

The rise and fall of genes: origins and functions of plant pathogen pangenomes

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

Prediction and Characterization of RXLR Effectors in Pythium Species

RXLR effectors, a class of secreted proteins that are transferred into host cells to manipulate host immunity, have been reported to widely exist in oomycetes, including those from genera Phytophthora, Hyaloperonospora, Albugo, and Saprolegnia. However, in Pythium species, no RXLR effector has yet been characterized, and the origin and evolution of such virulent effectors are still unknown. Here, we developed a modified regular expression method for de novo identification of RXLRs and characterized 359 putative RXLR effectors in nine Pythium species. Phylogenetic analysis revealed that all oomycetous RXLRs formed a single superfamily, suggesting that they might have a common ancestor. RXLR effectors from Pythium and Phytophthora species exhibited similar sequence features, protein structures, and genome locations. In particular, there were significantly more RXLR proteins in the mosquito biological control agent P. guiyangense than in the other eight Pythium species, and P. guiyangense RXLRs might be the result of gene duplication and genome rearrangement events, as indicated by synteny analysis. Expression pattern analysis of RXLR-encoding genes in the plant pathogen P. ultimum detected transcripts of the majority of the predicted RXLR genes, with some RXLR effectors induced in infection stages and one RXLR showing necrosis-inducing activity. Furthermore, all predicted RXLR genes were cloned from two biocontrol agents, P. oligandrum and P. periplocum, and three of the RXLR genes were found to induce a defense response in Nicotiana benthamiana. Taken together, our findings represent the first evidence of RXLR effectors in Pythium species, providing valuable information on their evolutionary patterns and the mechanisms of their interactions with diverse hosts.

Differential loss of effector genes in three recently expanded pandemic clonal lineages of the rice blast fungus

BACKGROUND

Understanding the mechanisms and timescales of plant pathogen outbreaks requires a detailed genome-scale analysis of their population history. The fungus Magnaporthe (Syn. Pyricularia) oryzae-the causal agent of blast disease of cereals- is among the most destructive plant pathogens to world agriculture and a major threat to the production of rice, wheat, and other cereals. Although M. oryzae is a multihost pathogen that infects more than 50 species of cereals and grasses, all rice-infecting isolates belong to a single genetically defined lineage. Here, we combined the two largest genomic datasets to reconstruct the genetic history of the rice-infecting lineage of M. oryzae based on 131 isolates from 21 countries.

RESULTS

The global population of the rice blast fungus consists mainly of three well-defined genetic groups and a diverse set of individuals. Multiple population genetic tests revealed that the rice-infecting lineage of the blast fungus probably originated from a recombining diverse group in Southeast Asia followed by three independent clonal expansions that took place over the last ~ 200 years. Patterns of allele sharing identified a subpopulation from the recombining diverse group that introgressed with one of the clonal lineages before its global expansion. Remarkably, the four genetic lineages of the rice blast fungus vary in the number and patterns of presence and absence of candidate effector genes. These genes encode secreted proteins that modulate plant defense and allow pathogen colonization. In particular, clonal lineages carry a reduced repertoire of effector genes compared with the diverse group, and specific combinations of presence and absence of effector genes define each of the pandemic clonal lineages.

CONCLUSIONS

Our analyses reconstruct the genetic history of the rice-infecting lineage of M. oryzae revealing three clonal lineages associated with rice blast pandemics. Each of these lineages displays a specific pattern of presence and absence of effector genes that may have shaped their adaptation to the rice host and their evolutionary history.

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

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

In vitro evidence of root colonization suggests ecological versatility in the genus Mycena

The root-associated habit has evolved on numerous occasions in different fungal lineages, suggesting a strong evolutionary pressure for saprotrophic fungi to switch to symbiotic associations with plants. Species within the ubiquitous, saprotrophic genus Mycena are frequently major components in molecular studies of root-associated fungal communities, suggesting that an evaluation of their trophic status is warranted. Here, we report on interactions between a range of Mycena species and the plant Betula pendula. In all, 17 Mycena species were inoculated onto B. pendula seedlings. Physical interactions between hyphae and fine roots were examined using differential staining and fluorescence microscopy. Physiological interactions were investigated using ¹⁴ C and ³² P to show potential transfer between symbionts. All Mycena species associated closely with fine roots, showing hyphal penetration into the roots, which in some cases were intracellular. Seven species formed mantle-like structures around root tips, but none formed a Hartig net. Mycena pura and Mycena galopus both enhanced seedling growth, with M. pura showing significant transfer of ³² P to the seedlings. Our results support the view that several Mycena species can associate closely with plant roots and some may potentially occupy a transitional state between saprotrophy and biotrophy.

Genetic Variation Bias toward Noncoding Regions and Secreted Proteins in the Rice Blast Fungus Magnaporthe oryzae

The genomes of plant pathogens are highly variable and plastic. Pathogen gene repertoires change quickly with the plant environment, which results in a rapid loss of plant resistance shortly after the pathogen emerges in the field. Extensive studies have evaluated natural pathogen populations to understand their evolutionary effects; however, the number of studies that have examined the dynamic processes of the mutation and adaptation of plant pathogens to host plants remains limited. Here, we applied experimental evolution and high-throughput pool sequencing to Magnaporthe oryzae, a fungal pathogen that causes massive losses in rice production, to observe the evolution of genome variation. We found that mutations, including single-nucleotide variants (SNVs), insertions and deletions (indels), and transposable element (TE) insertions, accumulated very rapidly throughout the genome of M. oryzae during sequential plant inoculation and preferentially in noncoding regions, while such mutations were not frequently found in coding regions. However, we also observed that new TE insertions accumulated with time and preferentially accumulated at the proximal region of secreted protein (SP) coding genes in M. oryzae populations. Taken together, these results revealed a bias in genetic variation toward noncoding regions and SP genes in M. oryzae and may contribute to the rapid adaptive evolution of the blast fungal effectors under host selection.IMPORTANCE Plants "lose" resistance toward pathogens shortly after their widespread emergence in the field because plant pathogens mutate and adapt rapidly under resistance selection. Thus, the rapid evolution of pathogens is a serious threat to plant health. Extensive studies have evaluated natural pathogen populations to understand their evolutionary effects; however, the study of the dynamic processes of the mutation and adaptation of plant pathogens to host plants remains limited. Here, by performing an experimental evolution study, we found a bias in genetic variation toward noncoding regions and SPs in the rice blast fungus Magnaporthe oryzae, which explains the ability of the rice blast fungus to maintain high virulence variation to overcome rice resistance in the field.

Characterization of Two Satellite DNA Families in the Genome of the Oomycete Plant Pathogen Phytophthora parasitica

Satellite DNA is a class of repetitive sequences that are organized in long arrays of tandemly repeated units in most eukaryotes. Long considered as selfish DNA, satellite sequences are now proposed to contribute to genome integrity. Despite their potential impact on the architecture and evolution of the genome, satellite DNAs have not been investigated in oomycetes due to the paucity of genomic data and the difficulty of assembling highly conserved satellite arrays. Yet gaining knowledge on the structure and evolution of genomes of oomycete pathogens is crucial to understanding the mechanisms underlying adaptation to their environment and to proposing efficient disease control strategies. A de novo assembly of the genome of Phytophthora parasitica, an important oomycete plant pathogen, led to the identification of several families of tandemly repeated sequences varying in size, copy number, and sequence conservation. Among them, two abundant families, designated as PpSat1 and PpSat2, displayed typical features of satellite DNA and were collectively designated as PpSat. These two satellite families differ by their length, sequence, organization, genomic environment, and evolutionary dynamics. PpSat1, but not PpSat2, presented homologs among oomycetes. This observation, as well as the characterization of transcripts of PpSat families, suggested that these satellite DNA families likely play a conserved role within this important group of pathogens.

Genome-Wide Increased Copy Number is Associated with Emergence of Dominant Clones of the Irish Potato Famine Pathogen Phytophthora infestans

The plant pathogen that caused the Irish potato famine, Phytophthora infestans, continues to reemerge globally. These modern epidemics are caused by clonally reproducing lineages. In contrast, a sexual mode of reproduction is observed at its center of origin in Mexico. We conducted a comparative genomic analysis of 47 high-coverage genomes to infer changes in genic copy number. We included samples from sexual populations at the center of origin as well as several dominant clonal lineages sampled worldwide. We conclude that sexual populations at the center of origin are diploid, as was the lineage that caused the famine, while modern clonal lineages showed increased copy number (3×). Copy number variation (CNV) was found genome-wide and did not to adhere to the two-speed genome hypothesis. Although previously reported, tetraploidy was not found in any of the genomes evaluated. We propose a model of dominant clone emergence supported by the epidemiological record (e.g., EU_13_A2, US-11, US-23) whereby a higher copy number provides fitness, leading to replacement of prior clonal lineages.IMPORTANCE The plant pathogen implicated in the Irish potato famine, Phytophthora infestans, continues to reemerge globally. Understanding changes in the genome during emergence can provide insights useful for managing this pathogen. Previous work has relied on studying individuals from the United States, South America, Europe, and China reporting that these can occur as diploids, triploids, or tetraploids and are clonal. We studied variation in sexual populations at the pathogen's center of origin, in Mexico, where it has been reported to reproduce sexually as well as within clonally reproducing, dominant clones from the United States and Europe. Our results newly show that sexual populations at the center of origin are diploid, whereas populations elsewhere are more variable and show genome-wide variation in gene copy number. We propose a model of evolution whereby new pathogen clones emerge predominantly by increasing the gene copy number genome-wide.

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.

Mycosis is a Disease State Encountered Rarely in Shore Crabs, Carcinus maenas

There is a paucity of knowledge regarding the diversity and impact(s) of disease-causing fungi in marine animals, especially shellfish. In efforts to address this knowledge gap for the shore crab Carcinus maenas, a year-long disease screen was carried out across two sites in Swansea Bay (Wales, UK) with a view to characterising putative fungal infections. Crabs were sampled between November 2017 and October 2018, and screened systematically for disease signatures using haemolymph (blood) preparations, targeted PCR and tissue histopathology. Strikingly, mycosis was confirmed in ~0.4% of total crabs tested (n = 1191) and restricted to one location only (Mumbles Pier). Clinical infections were observed in four out of four infected crabs. In these animals, the gills and hepatopancreas were congested with fungal morphotypes. In addition, some evidence indicates haemocyte (immune cell) reactivity toward the fungi. Phylogenetic placement of the partial internal transcribed spacer (ITS1) gene regions amplified from three mycotic crabs revealed the causative agent to be related to hypocrealean fungi, thereby representing a novel species.

Effector Repertoire of Phytophthora betacei: In Search of Possible Virulence Factors Responsible for Its Host Specificity

Phytophthora betacei is an oomycete plant pathogen closely related to Phytophthora infestans. It infects tree tomato (Solanum betaceum) in northern South America, but is, under natural conditions, unable to infect potatoes or tomatoes, the main hosts of its sister species P. infestans. We characterized, and compared the effector repertoires of P. betacei and other Phytophthora species. To this end, we used in silico approaches to predict and describe the repertoire of secreted proteins in Phytophthora species and determine unique and core effectors. P. betacei has the largest proteome and secretome of all Phytophthora species evaluated. We identified between 450 and 1933 candidate effector genes in Phytophthora ramorum, Phytophthora sojae, Phytophthora cactorum, Phytophthora parasitica, Phytophthora palmivora, P. infestans, and P. betacei genomes. The P. betacei predicted secretome contains 5653 proteins, 1126 of which are apoplastic effectors and 807cytoplasmic effectors. Genes encoding cytoplasmic effectors include 791 genes with an RxLR domain (the largest number known so far in a Phytophthora species) and 16 with a Crinkler (CRN) domain. We detected homologs of previously described avirulence gene (Avr) present in Phytophthora spp., such as Avr1, Avr3b, Avr4, and Avrblb1, suggesting a high level of effector gene conservation among Phytophthora species. Nonetheless, fewer CRN effectors were obtained in P. betacei compared to all other Phytophthora species analyzed. The comparison between P. infestans and P. betacei effector profiles shows unique features in P. betacei that might be involved in pathogenesis and host preference. Indeed, 402 unique predicted effector genes were detected in P. betacei, corresponding to 197 apoplastic effector genes, 203 RxLR cytoplasmic effector genes, and 2 effector genes with CRN domain. This is the first characterization of the effector profile of P. betacei and the broadest comparison of predicted effector repertoires in the genus Phytophthora following a standardized prediction pipeline. The resultant P. betacei putative effector repertoire provides a reasonable set of proteins whose experimental validation could lead to understand the specific virulence factors responsible for the host specificity of this species.

Threats to global food security from emerging fungal and oomycete crop pathogens

Emerging fungal and oomycete pathogens infect staple calorie crops and economically important commodity crops, thereby posing a significant risk to global food security. Our current agricultural systems - with emphasis on intensive monoculture practices - and globalized markets drive the emergence and spread of new pathogens and problematic traits, such as fungicide resistance. Climate change further promotes the emergence of pathogens on new crops and in new places. Here we review the factors affecting the introduction and spread of pathogens and current disease control strategies, illustrating these with the historic example of the Irish potato famine and contemporary examples of soybean rust, wheat blast and blotch, banana wilt and cassava root rot. Our Review looks to the future, summarizing what we see as the main challenges and knowledge gaps, and highlighting the direction that research must take to face the challenge of emerging crop pathogens.

Phylotranscriptomics of the Pentapetalae Reveals Frequent Regulatory Variation in Plant Local Responses to the Fungal Pathogen Sclerotinia sclerotiorum

Quantitative disease resistance (QDR) is a conserved form of plant immunity that limits infections caused by a broad range of pathogens. QDR has a complex genetic determinism. The extent to which molecular components of the QDR response vary across plant species remains elusive. The fungal pathogen Sclerotinia sclerotiorum, causal agent of white mold diseases on hundreds of plant species, triggers QDR in host populations. To document the diversity of local responses to S. sclerotiorum at the molecular level, we analyzed the complete transcriptomes of six species spanning the Pentapetalae (Phaseolus vulgaris, Ricinus communis, Arabidopsis [Arabidopsis thaliana], Helianthus annuus, Solanum lycopersicum, and Beta vulgaris) inoculated with the same strain of S. sclerotiorum About one-third of plant transcriptomes responded locally to S. sclerotiorum, including a high proportion of broadly conserved genes showing frequent regulatory divergence at the interspecific level. Evolutionary inferences suggested a trend toward the acquisition of gene induction relatively recently in several lineages. Focusing on a group of ABCG transporters, we propose that exaptation by regulatory divergence contributed to the evolution of QDR. This evolutionary scenario has implications for understanding the QDR spectrum and durability. Our work provides resources for functional studies of gene regulation and QDR molecular mechanisms across the Pentapetalae.

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.

The Many Questions about Mini Chromosomes in Colletotrichum spp

Many fungal pathogens carry accessory regions in their genome, which are not required for vegetative fitness. Often, although not always, these regions occur as relatively small chromosomes in different species. Such mini chromosomes appear to be a typical feature of many filamentous plant pathogens. Since these regions often carry genes coding for effectors or toxin-producing enzymes, they may be directly related to virulence of the respective pathogen. In this review, we outline the situation of small accessory chromosomes in the genus Colletotrichum, which accounts for ecologically important plant diseases. We summarize which species carry accessory chromosomes, their gene content, and chromosomal makeup. We discuss the large variation in size and number even between different isolates of the same species, their potential roles in host range, and possible mechanisms for intra- and interspecies exchange of these interesting genetic elements.

Genome reconstruction of the non-culturable spinach downy mildew Peronospora effusa by metagenome filtering

Peronospora effusa (previously known as P. farinosa f. sp. spinaciae, and here referred to as Pfs) is an obligate biotrophic oomycete that causes downy mildew on spinach (Spinacia oleracea). To combat this destructive many disease resistant cultivars have been bred and used. However, new Pfs races rapidly break the employed resistance genes. To get insight into the gene repertoire of Pfs and identify infection-related genes, the genome of the first reference race, Pfs1, was sequenced, assembled, and annotated. Due to the obligate biotrophic nature of this pathogen, material for DNA isolation can only be collected from infected spinach leaves that, however, also contain many other microorganisms. The obtained sequences can, therefore, be considered a metagenome. To filter and obtain Pfs sequences we utilized the CAT tool to taxonomically annotate ORFs residing on long sequences of a genome pre-assembly. This study is the first to show that CAT filtering performs well on eukaryotic contigs. Based on the taxonomy, determined on multiple ORFs, contaminating long sequences and corresponding reads were removed from the metagenome. Filtered reads were re-assembled to provide a clean and improved Pfs genome sequence of 32.4 Mbp consisting of 8,635 scaffolds. Transcript sequencing of a range of infection time points aided the prediction of a total of 13,277 gene models, including 99 RxLR(-like) effector, and 14 putative Crinkler genes. Comparative analysis identified common features in the predicted secretomes of different obligate biotrophic oomycetes, regardless of their phylogenetic distance. Their secretomes are generally smaller, compared to hemi-biotrophic and necrotrophic oomycete species. We observe a reduction in proteins involved in cell wall degradation, in Nep1-like proteins (NLPs), proteins with PAN/apple domains, and host translocated effectors. The genome of Pfs1 will be instrumental in studying downy mildew virulence and for understanding the molecular adaptations by which new isolates break spinach resistance.

Genetic diversity of Alternaria species associated with black point in wheat grains

The genus Alternaria is a widely distributed major plant pathogen that can act as a saprophyte in plant debris. Fungi of this genus frequently infect cereal crops and cause such diseases as black point and wheat leaf blight, which decrease the yield and quality of cereal products. A total of 25 Alternaria sp. isolates were collected from germ grains of various wheat cultivars from different geographic regions in Kazakhstan. We investigated the genetic relationships of the main Alternaria species related to black point disease of wheat in Kazakhstan, using the inter-primer binding site (iPBS) DNA profiling technique. We used 25 retrotransposon-based iPBS primers to identify the differences among and within Alternaria species populations, and analyzed the variation using clustering (UPGMA) and statistical approaches (AMOVA). Isolates of Alternaria species clustered into two main genetic groups, with species of A.alternata and A.tennuissima forming one cluster, and isolates of A. infectoria forming another. The genetic diversity found using retrotransposon profiles was strongly correlated with geographic data. Overall, the iPBS fingerprinting technique is highly informative and useful for the evaluation of genetic diversity and relationships of Alternaria species.