Two-dimensional data binning for the analysis of genome architecture in filamentous plant pathogens and other eukaryotes

Genomic and transcriptomic analyses of Phytophthora cinnamomi reveal complex genome architecture, expansion of pathogenicity factors, and host-dependent gene expression profiles

Phytophthora cinnamomi is a hemibiotrophic oomycete causing Phytophthora root rot in over 5,000 plant species, threatening natural ecosystems, forestry, and agriculture. Genomic studies of P. cinnamomi are limited compared to other Phytophthora spp. despite the importance of this destructive and highly invasive pathogen. The genome of two genetically and phenotypically distinct P. cinnamomi isolates collected from avocado orchards in California were sequenced using PacBio and Illumina sequencing. Genome sizes were estimated by flow cytometry and assembled de novo to 140-141 Mb genomes with 21,111-21,402 gene models. Genome analyses revealed that both isolates exhibited complex heterozygous genomes fitting the two-speed genome model. The more virulent isolate encodes a larger secretome and more RXLR effectors when compared to the less virulent isolate. Transcriptome analysis after P. cinnamomi infection in Arabidopsis thaliana, Nicotiana benthamiana, and Persea americana de Mill (avocado) showed that this pathogen deploys common gene repertoires in all hosts and host-specific subsets, especially among effectors. Overall, our results suggested that clonal P. cinnamomi isolates employ similar strategies as other Phytophthora spp. to increase phenotypic diversity (e.g., polyploidization, gene duplications, and a bipartite genome architecture) to cope with environmental changes. Our study also provides insights into common and host-specific P. cinnamomi infection strategies and may serve as a method for narrowing and selecting key candidate effectors for functional studies to determine their contributions to plant resistance or susceptibility.

Effector Repertoire of the Sweetpotato Black Rot Fungal Pathogen Ceratocystis fimbriata

In 2015, sweetpotato producers in the United States experienced one of the worst outbreaks of black rot recorded in history, with up to 60% losses reported in the field and packing houses and at shipping ports. Host resistance remains the ideal management tool to decrease crop losses. Lack of knowledge of Ceratocystis fimbriata biology represents a critical barrier for the deployment of resistance to black rot in sweetpotato. In this study, we scanned the recent near chromosomal-level assembly for putative secreted effectors in the sweetpotato C. fimbriata isolate AS236 using a custom fungal effector annotation pipeline. We identified a set of 188 putative effectors on the basis of secretion signal and in silico prediction in EffectorP. We conducted a deep RNA time-course sequencing experiment to determine whether C. fimbriata modulates effectors in planta and to define a candidate list of effectors expressed during infection. We examined the expression profile of two C. fimbriata isolates, a pre-epidemic (1990s) isolate and a post-epidemic (2015) isolate. Our in planta expression profiling revealed clusters of co-expressed secreted effector candidates. Based on fold-change differences of putative effectors in both isolates and over the course of infection, we suggested prioritization of 31 effectors for functional characterization. Among this set, we identified several effectors that provide evidence for a marked biotrophic phase in C. fimbriata during infection of sweetpotato storage roots. Our study revealed a catalog of effector proteins that provide insight into C. fimbriata infection mechanisms and represent a core catalog to implement effector-assisted breeding in sweetpotato. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The genomes of Scedosporium between environmental challenges and opportunism

Emerging fungal pathogens are a global challenge for humankind. Many efforts have been made to understand the mechanisms underlying pathogenicity in bacteria, and OMICs techniques are largely responsible for those advancements. By contrast, our limited understanding of opportunism and antifungal resistance is preventing us from identifying, limiting and interpreting the emergence of fungal pathogens. The genus Scedosporium (Microascaceae) includes fungi with high tolerance to environmental pollution, whilst some species can be considered major human pathogens, such as Scedosporium apiospermum and Scedosporium boydii. However, unlike other fungal pathogens, little is known about the genome evolution of these organisms. We sequenced two novel genomes of Scedosporium aurantiacum and Scedosporium minutisporum isolated from extreme, strongly anthropized environments. We compared all the available Scedosporium and Microascaceae genomes, that we systematically annotated and characterized ex novo in most cases. The genomes in this family were integrated in a Phylum-level comparison to infer the presence of putative, shared genomic traits in filamentous ascomycetes with pathogenic potential. The analysis included the genomes of 100 environmental and clinical fungi, revealing poor evolutionary convergence of putative pathogenicity traits. By contrast, several features in Microascaceae and Scedosporium were detected that might have a dual role in responding to environmental challenges and allowing colonization of the human body, including chitin, melanin and other cell wall related genes, proteases, glutaredoxins and magnesium transporters. We found these gene families to be impacted by expansions, orthologous transposon insertions, and point mutations. With RNA-seq, we demonstrated that most of these anciently impacted genomic features responded to the stress imposed by an antifungal compound (voriconazole) in the two environmental strains S. aurantiacum MUT6114 and S. minutisporum MUT6113. Therefore, the present genomics and transcriptomics investigation stands on the edge between stress resistance and pathogenic potential, to elucidate whether fungi were pre-adapted to infect humans. We highlight the strengths and limitations of genomics applied to opportunistic human pathogens, the multifactoriality of pathogenicity and resistance to drugs, and suggest a scenario where pressures other than anthropic contributed to forge filamentous human pathogens.

Short Linear Motifs (SLiMs) in “Core” RxLR Effectors of Phytophthora parasitica var. nicotianae : a Case of PpRxLR1 Effector

Oomycetes of the genus Phytophthora encompass several of the most successful plant pathogens described to date. The success of infection by Phytophthora species is attributed to the pathogens’ ability to secrete effector proteins that alter the host’s physiological processes. Structural analyses of effector proteins mainly from bacterial and viral pathogens have revealed the presence of intrinsically disordered regions that host short linear motifs (SLiMs). These motifs play important biological roles by facilitating protein-protein interactions as well as protein translocation. Nonetheless, SLiMs in Phytophthora species RxLR effectors have not been investigated previously and their roles remain unknown. Using a bioinformatics pipeline, we identified 333 candidate RxLR effectors in the strain INRA 310 of Phytophthora parasitica. Of these, 71 (21%) were also found to be present in 10 other genomes of P. parasitica, and hence, these were designated core RxLR effectors (CREs). Within the CRE sequences, the N terminus exhibited enrichment in intrinsically disordered regions compared to the C terminus, suggesting a potential role of disorder in effector translocation. Although the disorder content was reduced in the C-terminal regions, it is important to mention that most SLiMs were in this terminus. PpRxLR1 is one of the 71 CREs identified in this study, and its genes encode a 6-amino acid (aa)-long SLiM at the C terminus. We showed that PpRxLR1 interacts with several host proteins that are implicated in defense. Structural analysis of this effector using homology modeling revealed the presence of potential ligand-binding sites. Among key residues that were predicted to be crucial for ligand binding, L¹⁰² and Y¹⁰⁶ were of interest since they form part of the 6-aa-long PpRxLR1 SLiM. In silico substitution of these two residues to alanine was predicted to have a significant effect on both the function and the structure of PpRxLR1 effector. Molecular docking simulations revealed possible interactions between PpRxLR1 effector and ubiquitin-associated proteins. The ubiquitin-like SLiM carried in this effector was shown to be a potential mediator of these interactions. Further studies are required to validate and elucidate the underlying molecular mechanism of action.

IMPORTANCE The continuous gain and loss of RxLR effectors makes the control of Phytophthora spp. difficult. Therefore, in this study, we endeavored to identify RxLR effectors that are highly conserved among species, also known as “core” RxLR effectors (CREs). We reason that these highly conserved effectors target conserved proteins or processes; thus, they can be harnessed in breeding for durable resistance in plants. To further understand the mechanisms of action of CREs, structural dissection of these proteins is crucial. Intrinsically disordered regions (IDRs) that do not adopt a fixed, three-dimensional fold carry short linear motifs (SLiMs) that mediate biological functions of proteins. The presence and potential role of these SLiMs in CREs of Phytophthora spp. have been overlooked. To our knowledge, we have effectively identified CREs as well as SLiMs with the potential of promoting effector virulence. Together, this work has advanced our comprehension of Phytophthora RxLR effector function and may facilitate the development of innovative and effective control strategies.

An Improved Assembly of the Albugo candida Ac2V Genome Reveals the Expansion of the "CCG" Class of Effectors

Albugo candida is an obligate oomycete pathogen that infects many plants in the Brassicaceae family. We resequenced the genome of isolate Ac2V using PacBio long reads and constructed an assembly augmented by Illumina reads. The Ac2VPB genome assembly is 10% larger and more contiguous compared with a previous version. Our annotation of the new assembly, aided by RNA-sequencing information, revealed a 175% expansion (40 to 110) in the CHxC effector class, which we redefined as "CCG" based on motif analysis. This class of effectors consist of arrays of phylogenetically related paralogs residing in gene sparse regions, and shows signatures of positive selection and presence/absence polymorphism. This work provides a resource that allows the dissection of the genomic components underlying A. candida adaptation and, particularly, the role of CCG effectors in virulence and avirulence on different hosts.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

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.

Emergence and diversification of a highly invasive chestnut pathogen lineage across southeastern Europe

Invasive microbial species constitute a major threat to biodiversity, agricultural production and human health. Invasions are often dominated by one or a small number of genotypes, yet the underlying factors driving invasions are poorly understood. The chestnut blight fungus Cryphonectria parasitica first decimated the North American chestnut, and a more recent outbreak threatens European chestnut stands. To unravel the chestnut blight invasion of southeastern Europe, we sequenced 230 genomes of predominantly European strains. Genotypes outside of the invasion zone showed high levels of diversity with evidence for frequent and ongoing recombination. The invasive lineage emerged from the highly diverse European genotype pool rather than a secondary introduction from Asia or North America. The expansion across southeastern Europe was mostly clonal and is dominated by a single mating type, suggesting a fitness advantage of asexual reproduction. Our findings show how an intermediary, highly diverse bridgehead population gave rise to an invasive, largely clonally expanding pathogen.

A High-Quality Grapevine Downy Mildew Genome Assembly Reveals Rapidly Evolving and Lineage-Specific Putative Host Adaptation Genes

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high contiguity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant–pathogen interactions. Plasmopara viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.

Genomic Plasticity Mediated by Transposable Elements in the Plant Pathogenic Fungus Colletotrichum higginsianum

Phytopathogen genomes are under constant pressure to change, as pathogens are locked in an evolutionary arms race with their hosts, where pathogens evolve effector genes to manipulate their hosts, whereas the hosts evolve immune components to recognize the products of these genes. Colletotrichum higginsianum (Ch), a fungal pathogen with no known sexual morph, infects Brassicaceae plants including Arabidopsis thaliana. Previous studies revealed that Ch differs in its virulence toward various Arabidopsis thaliana ecotypes, indicating the existence of coevolutionary selective pressures. However, between-strain genomic variations in Ch have not been studied. Here, we sequenced and assembled the genome of a Ch strain, resulting in a highly contiguous genome assembly, which was compared with the chromosome-level genome assembly of another strain to identify genomic variations between strains. We found that the two closely related strains vary in terms of large-scale rearrangements, the existence of strain-specific regions, and effector candidate gene sets and that these variations are frequently associated with transposable elements (TEs). Ch has a compartmentalized genome consisting of gene-sparse, TE-dense regions with more effector candidate genes and gene-dense, TE-sparse regions harboring conserved genes. Additionally, analysis of the conservation patterns and syntenic regions of effector candidate genes indicated that the two strains vary in their effector candidate gene sets because of de novo evolution, horizontal gene transfer, or gene loss after divergence. Our results reveal mechanisms for generating genomic diversity in this asexual pathogen, which are important for understanding its adaption to hosts.

Phytophthora methylomes are modulated by 6mA methyltransferases and associated with adaptive genome regions

BACKGROUND

Filamentous plant pathogen genomes often display a bipartite architecture with gene-sparse, repeat-rich compartments serving as a cradle for adaptive evolution. The extent to which this two-speed genome architecture is associated with genome-wide DNA modifications is unknown.

RESULTS

We show that the oomycetes Phytophthora infestans and Phytophthora sojae possess functional adenine N6-methylation (6mA) methyltransferases that modulate patterns of 6mA marks across the genome. In contrast, 5-methylcytosine could not be detected in these species. Methylated DNA IP sequencing (MeDIP-seq) of each species reveals 6mA is depleted around the transcription start sites (TSSs) and is associated with lowly expressed genes, particularly transposable elements. Genes occupying the gene-sparse regions have higher levels of 6mA in both genomes, possibly implicating the methylome in adaptive evolution. All six putative adenine methyltransferases from P. infestans and P. sojae, except PsDAMT2, display robust enzymatic activities. Surprisingly, single knockouts in P. sojae significantly reduce in vivo 6mA levels, indicating that the three enzymes are not fully redundant. MeDIP-seq of the psdamt3 mutant reveals uneven 6mA methylation reduction across genes, suggesting that PsDAMT3 may have a preference for gene body methylation after the TSS. Furthermore, transposable elements such as DNA elements are more active in the psdamt3 mutant. A large number of genes, particularly those from the adaptive genomic compartment, are differentially expressed.

CONCLUSIONS

Our findings provide evidence that 6mA modification is potentially an epigenetic mark in Phytophthora genomes, and complex patterns of 6mA methylation may be associated with adaptive evolution in these important plant pathogens.

Draft Genome Sequence for the Tree Pathogen Phytophthora plurivora

Species from the genus Phytophthora are well represented among organisms causing serious diseases on trees. Phytophthora plurivora has been implicated in long-term decline of woodland trees across Europe. Here we present a draft genome sequence of P. plurivora, originally isolated from diseased European beech (Fagus sylvatica) in Malmö, Sweden. When compared with other sequenced Phytophthora species, the P. plurivora genome assembly is relatively compact, spanning 41 Mb. This is organized in 1,919 contigs and 1,898 scaffolds, encompassing 11,741 predicted genes, and has a repeat content of approximately 15%. Comparison of allele frequencies revealed evidence for tetraploidy in the sequenced isolate. As in other sequenced Phytophthora species, P. plurivora possesses genes for pathogenicity-associated RXLR and Crinkle and Necrosis effectors, predominantly located in gene-sparse genomic regions. Comparison of the P. plurivora RXLR effectors with orthologs in other sequenced species in the same clade (Phytophthora multivora and Phytophthora capsici) revealed that the orthologs were likely to be under neutral or purifying selection.

De Novo Assembly and Phasing of Dikaryotic Genomes from Two Isolates of Puccinia coronata f. sp. avenae, the Causal Agent of Oat Crown Rust

Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenaeIMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae, which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae, resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae.

Secretome Analysis from the Ectomycorrhizal Ascomycete Cenococcum geophilum

Cenococcum geophilum is an ectomycorrhizal fungus with global distribution in numerous habitats and associates with a large range of host species including gymnosperm and angiosperm trees. Moreover, C. geophilum is the unique ectomycorrhizal species within the clade Dothideomycetes, the largest class of Ascomycetes containing predominantly saprotrophic and many devastating phytopathogenic fungi. Recent studies highlight that mycorrhizal fungi, as pathogenic ones, use effectors in form of Small Secreted Proteins (SSPs) as molecular keys to promote symbiosis. In order to better understand the biotic interaction of C. geophilum with its host plants, the goal of this work was to characterize mycorrhiza-induced small-secreted proteins (MiSSPs) that potentially play a role in the ectomycorrhiza formation and functioning of this ecologically very important species. We combined different approaches such as gene expression profiling, genome localization and conservation of MiSSP genes in different C. geophilum strains and closely related species as well as protein subcellular localization studies of potential targets of MiSSPs in interacting plants using in tobacco leaf cells. Gene expression analyses of C. geophilum interacting with Pinus sylvestris (pine) and Populus tremula × Populus alba (poplar) showed that similar sets of genes coding for secreted proteins were up-regulated and only few were specific to each host. Whereas pine induced more carbohydrate active enzymes (CAZymes), the interaction with poplar induced the expression of specific SSPs. We identified a set of 22 MiSSPs, which are located in both, gene-rich, repeat-poor or gene-sparse, repeat-rich regions of the C. geophilum genome, a genome showing a bipartite architecture as seen for some pathogens but not yet for an ectomycorrhizal fungus. Genome re-sequencing data of 15 C. geophilum strains and two close relatives Glonium stellatum and Lepidopterella palustris were used to study sequence conservation of MiSSP-encoding genes. The 22 MiSSPs showed a high presence-absence polymorphism among the studied C. geophilum strains suggesting an evolution through gene gain/gene loss. Finally, we showed that six CgMiSSPs target four distinct sub-cellular compartments such as endoplasmic reticulum, plasma membrane, cytosol and tonoplast. Overall, this work presents a comprehensive analysis of secreted proteins and MiSSPs in different genetic level of C. geophilum opening a valuable resource to future functional analysis.

Effector Biology in Focus: A Primer for Computational Prediction and Functional Characterization

Plant-pathogen interactions are controlled by a multilayered immune system, which is activated by pathogen recognition in the host. Pathogens secrete effector molecules to interfere with the immune recognition or signaling network and reprogram cell structure or metabolism. Understanding the effector repertoires of diverse pathogens will contribute to unraveling the molecular mechanism of virulence and developing sustainable disease-control strategies for crops and natural ecosystems. Effector functionality has been investigated extensively in only a small number of pathogen species. However, many more pathogen genomes are becoming available, and much can be learned from a broader view of effector biology in diverse pathosystems. The purpose of this review is to summarize methodology for computational prediction of protein effectors, functional characterization of effector proteins and their targets, and the use of effectors as probes to screen for new sources of host resistance. Although these techniques were generally developed in model pathosystems, many of the approaches are directly applicable for exploration and exploitation of effector biology in pathosystems that are less well studied. We hope to facilitate such exploration, which will broaden understanding of the mechanisms that underpin the biological diversity of plant-pathogen interactions, and maximize the impact of new approaches that leverage effector biology for disease control.

Comparative analysis of the predicted secretomes of Rosaceae scab pathogens Venturia inaequalis and V. pirina reveals expanded effector families and putative determinants of host range

BACKGROUND

Fungal plant pathogens belonging to the genus Venturia cause damaging scab diseases of members of the Rosaceae. In terms of economic impact, the most important of these are V. inaequalis, which infects apple, and V. pirina, which is a pathogen of European pear. Given that Venturia fungi colonise the sub-cuticular space without penetrating plant cells, it is assumed that effectors that contribute to virulence and determination of host range will be secreted into this plant-pathogen interface. Thus the predicted secretomes of a range of isolates of Venturia with distinct host-ranges were interrogated to reveal putative proteins involved in virulence and pathogenicity.

RESULTS

Genomes of Venturia pirina (one European pear scab isolate) and Venturia inaequalis (three apple scab, and one loquat scab, isolates) were sequenced and the predicted secretomes of each isolate identified. RNA-Seq was conducted on the apple-specific V. inaequalis isolate Vi1 (in vitro and infected apple leaves) to highlight virulence and pathogenicity components of the secretome. Genes encoding over 600 small secreted proteins (candidate effectors) were identified, most of which are novel to Venturia, with expansion of putative effector families a feature of the genus. Numerous genes with similarity to Leptosphaeria maculans AvrLm6 and the Verticillium spp. Ave1 were identified. Candidates for avirulence effectors with cognate resistance genes involved in race-cultivar specificity were identified, as were putative proteins involved in host-species determination. Candidate effectors were found, on average, to be in regions of relatively low gene-density and in closer proximity to repeats (e.g. transposable elements), compared with core eukaryotic genes.

CONCLUSIONS

Comparative secretomics has revealed candidate effectors from Venturia fungal plant pathogens that attack pome fruit. Effectors that are putative determinants of host range were identified; both those that may be involved in race-cultivar and host-species specificity. Since many of the effector candidates are in close proximity to repetitive sequences this may point to a possible mechanism for the effector gene family expansion observed and a route to diversification via transposition and repeat-induced point mutation.

Genome Sequence and Architecture of the Tobacco Downy Mildew Pathogen Peronospora tabacina

Peronospora tabacina is an obligate biotrophic oomycete that causes blue mold or downy mildew on tobacco (Nicotiana tabacum). It is an economically important disease occurring frequently in tobacco-growing regions worldwide. We sequenced and characterized the genomes of two P. tabacina isolates and mined them for pathogenicity-related proteins and effector-encoding genes. De novo assembly of the genomes using Illumina reads resulted in 4,016 (63.1 Mb, N50 = 79 kb) and 3,245 (55.3 Mb, N50 = 61 kb) scaffolds for isolates 968-J2 and 968-S26, respectively, with an estimated genome size of 68 Mb. The mitochondrial genome has a similar size (approximately 43 kb) and structure to those of other oomycetes, plus several minor unique features. Repetitive elements, primarily retrotransposons, make up approximately 24% of the nuclear genome. Approximately 18,000 protein-coding gene models were predicted. Mining the secretome revealed approximately 120 candidate RxLR, six CRN (candidate effectors that elicit crinkling and necrosis), and 61 WY domain-containing proteins. Candidate RxLR effectors were shown to be predominantly undergoing diversifying selection, with approximately 57% located in variable gene-sparse regions of the genome. Aligning the P. tabacina genome to Hyaloperonospora arabidopsidis and Phytophthora spp. revealed a high level of synteny. Blocks of synteny show gene inversions and instances of expansion in intergenic regions. Extensive rearrangements of the gene-rich genomic regions do not appear to have occurred during the evolution of these highly variable pathogens. These assemblies provide the basis for studies of virulence in this and other downy mildew pathogens.

Genomic characterisation of the effector complement of the potato cyst nematode Globodera pallida

Background

The potato cyst nematode Globodera pallida has biotrophic interactions with its host. The nematode induces a feeding structure – the syncytium – which it keeps alive for the duration of the life cycle and on which it depends for all nutrients required to develop to the adult stage. Interactions of G. pallida with the host are mediated by effectors, which are produced in two sets of gland cells. These effectors suppress host defences, facilitate migration and induce the formation of the syncytium.

Results

The recent completion of the G. pallida genome sequence has allowed us to identify the effector complement from this species. We identify 128 orthologues of effectors from other nematodes as well as 117 novel effector candidates. We have used in situ hybridisation to confirm gland cell expression of a subset of these effectors, demonstrating the validity of our effector identification approach. We have examined the expression profiles of all effector candidates using RNAseq; this analysis shows that the majority of effectors fall into one of three clusters of sequences showing conserved expression characteristics (invasive stage nematode only, parasitic stage only or invasive stage and adult male only). We demonstrate that further diversity in the effector pool is generated by alternative splicing. In addition, we show that effectors target a diverse range of structures in plant cells, including the peroxisome. This is the first identification of effectors from any plant pathogen that target this structure.

Conclusion

This is the first genome scale search for effectors, combined to a life-cycle expression analysis, for any plant-parasitic nematode. We show that, like other phylogenetically unrelated plant pathogens, plant parasitic nematodes deploy hundreds of effectors in order to parasitise plants, with different effectors required for different phases of the infection process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-923) contains supplementary material, which is available to authorized users.

Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum

BACKGROUND

The white mold fungus Sclerotinia sclerotiorum is a devastating necrotrophic plant pathogen with a remarkably broad host range. The interaction of necrotrophs with their hosts is more complex than initially thought, and still poorly understood.

RESULTS

We combined bioinformatics approaches to determine the repertoire of S. sclerotiorum effector candidates and conducted detailed sequence and expression analyses on selected candidates. We identified 486 S. sclerotiorum secreted protein genes expressed in planta, many of which have no predicted enzymatic activity and may be involved in the interaction between the fungus and its hosts. We focused on those showing (i) protein domains and motifs found in known fungal effectors, (ii) signatures of positive selection, (iii) recent gene duplication, or (iv) being S. sclerotiorum-specific. We identified 78 effector candidates based on these properties. We analyzed the expression pattern of 16 representative effector candidate genes on four host plants and revealed diverse expression patterns.

CONCLUSIONS

These results reveal diverse predicted functions and expression patterns in the repertoire of S. sclerotiorum effector candidates. They will facilitate the functional analysis of fungal pathogenicity determinants and should prove useful in the search for plant quantitative disease resistance components active against the white mold.