Comparative genomics of two inbred lines of the potato cyst nematode Globodera rostochiensis reveals disparate effector family-specific diversification patterns.

BACKGROUND: Potato cyst nematodes belong to the most harmful pathogens in potato, and durable management of these parasites largely depends on host-plant resistances. These resistances are pathotype specific. The current Globodera rostochiensis pathotype scheme that defines five pathotypes (Ro1 - Ro5) is both fundamentally and practically of limited value. Hence, resistant potato varieties are used worldwide in a poorly informed manner. RESULTS: We generated two novel reference genomes of G. rostochiensis inbred lines derived from a Ro1 and a Ro5 population. These genome sequences comprise 173 and 189 scaffolds respectively, marking a ≈ 24-fold reduction in fragmentation as compared to the current reference genome. We provide copy number variations for 19 effector families. Four dorsal gland effector families were investigated in more detail. SPRYSECs, known to be implicated in plant defence suppression, constitute by far the most diversified family studied herein with 60 and 99 variants in Ro1 and Ro5 distributed over 18 and 26 scaffolds. In contrast, CLEs, effectors involved in feeding site induction, show strong physical clustering. The 10 and 16 variants cluster on respectively 2 and 1 scaffolds. Given that pathotypes are defined by their effectoromes, we pinpoint the disparate nature of the contributing effector families in terms of sequence diversification and loss and gain of variants. CONCLUSIONS: Two novel reference genomes allow for nearly complete inventories of effector diversification and physical organisation within and between pathotypes. Combined with insights we provide on effector family-specific diversification patterns, this constitutes a basis for an effectorome-based virulence scheme for this notorious pathogen.

The effector GpRbp-1 of Globodera pallida targets a nuclear HECT E3 ubiquitin ligase to modulate gene expression in the host.

Plant-parasitic nematodes secrete effectors that manipulate plant cell morphology and physiology to achieve host invasion and establish permanent feeding sites. Effectors from the highly expanded SPRYSEC (SPRY domain with a signal peptide for secretion) family in potato cyst nematodes have been implicated in activation and suppression of plant immunity, but the mechanisms underlying these activities remain largely unexplored. To study the host mechanisms used by SPRYSEC effectors, we identified plant targets of GpRbp-1 from the potato cyst nematode Globodera pallida. Here, we show that GpRbp-1 interacts in yeast and in planta with a functional potato homologue of the Homology to E6-AP C-Terminus (HECT)-type ubiquitin E3 ligase UPL3, which is located in the nucleus. Potato lines lacking StUPL3 are not available, but the Arabidopsis mutant upl3-5 displaying a reduced UPL3 expression showed a consistently small but not significant decrease in susceptibility to cyst nematodes. We observed a major impact on the root transcriptome by the lower levels of AtUPL3 in the upl3-5 mutant, but surprisingly only in association with infections by cyst nematodes. To our knowledge, this is the first example that a HECT-type ubiquitin E3 ligase is targeted by a pathogen effector and that a member of this class of proteins specifically regulates gene expression under biotic stress conditions. Together, our data suggest that GpRbp-1 targets a specific component of the plant ubiquitination machinery to manipulate the stress response in host cells.

Distinct Roles of Non-Overlapping Surface Regions of the Coiled-Coil Domain in the Potato Immune Receptor Rx1.

The intracellular immune receptor Rx1 of potato (Solanum tuberosum), which confers effector-triggered immunity to Potato virus X, consists of a central nucleotide-binding domain (NB-ARC) flanked by a carboxyl-terminal leucine-rich repeat (LRR) domain and an amino-terminal coiled-coil (CC) domain. Rx1 activity is strictly regulated by interdomain interactions between the NB-ARC and LRR, but the contribution of the CC domain in regulating Rx1 activity or immune signaling is not fully understood. Therefore, we used a structure-informed approach to investigate the role of the CC domain in Rx1 functionality. Targeted mutagenesis of CC surface residues revealed separate regions required for the intramolecular and intermolecular interaction of the CC with the NB-ARC-LRR and the cofactor Ran GTPase-activating protein2 (RanGAP2), respectively. None of the mutant Rx1 proteins was constitutively active, indicating that the CC does not contribute to the autoinhibition of Rx1 activity. Instead, the CC domain acted as a modulator of downstream responses involved in effector-triggered immunity. Systematic disruption of the hydrophobic interface between the four helices of the CC enabled the uncoupling of cell death and disease resistance responses. Moreover, a strong dominant negative effect on Rx1-mediated resistance and cell death was observed upon coexpression of the CC alone with full-length Rx1 protein, which depended on the RanGAP2-binding surface of the CC. Surprisingly, coexpression of the N-terminal half of the CC enhanced Rx1-mediated resistance, which further indicated that the CC functions as a scaffold for downstream components involved in the modulation of disease resistance or cell death signaling.

Sequence Exchange between Homologous NB-LRR Genes Converts Virus Resistance into Nematode Resistance, and Vice Versa.

Plants have evolved a limited repertoire of NB-LRR disease resistance (R) genes to protect themselves against myriad pathogens. This limitation is thought to be counterbalanced by the rapid evolution of NB-LRR proteins, as only a few sequence changes have been shown to be sufficient to alter resistance specificities toward novel strains of a pathogen. However, little is known about the flexibility of NB-LRR R genes to switch resistance specificities between phylogenetically unrelated pathogens. To investigate this, we created domain swaps between the close homologs Gpa2 and Rx1, which confer resistance in potato (Solanum tuberosum) to the cyst nematode Globodera pallida and Potato virus X, respectively. The genetic fusion of the CC-NB-ARC of Gpa2 with the LRR of Rx1 (Gpa2CN/Rx1L) results in autoactivity, but lowering the protein levels restored its specific activation response, including extreme resistance to Potato virus X in potato shoots. The reciprocal chimera (Rx1CN/Gpa2L) shows a loss-of-function phenotype, but exchange of the first three LRRs of Gpa2 by the corresponding region of Rx1 was sufficient to regain a wild-type resistance response to G. pallida in the roots. These data demonstrate that exchanging the recognition moiety in the LRR is sufficient to convert extreme virus resistance in the leaves into mild nematode resistance in the roots, and vice versa. In addition, we show that the CC-NB-ARC can operate independently of the recognition specificities defined by the LRR domain, either aboveground or belowground. These data show the versatility of NB-LRR genes to generate resistance to unrelated pathogens with completely different lifestyles and routes of invasion.

The genome of the yellow potato cyst nematode, Globodera rostochiensis, reveals insights into the basis of parasitism and virulence.

BACKGROUND: The yellow potato cyst nematode, Globodera rostochiensis, is a devastating plant pathogen of global economic importance. This biotrophic parasite secretes effectors from pharyngeal glands, some of which were acquired by horizontal gene transfer, to manipulate host processes and promote parasitism. G. rostochiensis is classified into pathotypes with different plant resistance-breaking phenotypes. RESULTS: We generate a high quality genome assembly for G. rostochiensis pathotype Ro1, identify putative effectors and horizontal gene transfer events, map gene expression through the life cycle focusing on key parasitic transitions and sequence the genomes of eight populations including four additional pathotypes to identify variation. Horizontal gene transfer contributes 3.5 % of the predicted genes, of which approximately 8.5 % are deployed as effectors. Over one-third of all effector genes are clustered in 21 putative 'effector islands' in the genome. We identify a dorsal gland promoter element motif (termed DOG Box) present upstream in representatives from 26 out of 28 dorsal gland effector families, and predict a putative effector superset associated with this motif. We validate gland cell expression in two novel genes by in situ hybridisation and catalogue dorsal gland promoter element-containing effectors from available cyst nematode genomes. Comparison of effector diversity between pathotypes highlights correlation with plant resistance-breaking. CONCLUSIONS: These G. rostochiensis genome resources will facilitate major advances in understanding nematode plant-parasitism. Dorsal gland promoter element-containing effectors are at the front line of the evolutionary arms race between plant and parasite and the ability to predict gland cell expression a priori promises rapid advances in understanding their roles and mechanisms of action.

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.

Structural determinants at the interface of the ARC2 and leucine-rich repeat domains control the activation of the plant immune receptors Rx1 and Gpa2.

Many plant and animal immune receptors have a modular nucleotide-binding-leucine-rich repeat (NB-LRR) architecture in which a nucleotide-binding switch domain, NB-ARC, is tethered to a LRR sensor domain. The cooperation between the switch and sensor domains, which regulates the activation of these proteins, is poorly understood. Here, we report structural determinants governing the interaction between the NB-ARC and LRR in the highly homologous plant immune receptors Gpa2 and Rx1, which recognize the potato cyst nematode Globodera pallida and Potato virus X, respectively. Systematic shuffling of polymorphic sites between Gpa2 and Rx1 showed that a minimal region in the ARC2 and N-terminal repeats of the LRR domain coordinate the activation state of the protein. We identified two closely spaced amino acid residues in this region of the ARC2 (positions 401 and 403) that distinguish between autoactivation and effector-triggered activation. Furthermore, a highly acidic loop region in the ARC2 domain and basic patches in the N-terminal end of the LRR domain were demonstrated to be required for the physical interaction between the ARC2 and LRR. The NB-ARC and LRR domains dissociate upon effector-dependent activation, and the complementary-charged regions are predicted to mediate a fast reassociation, enabling multiple rounds of activation. Finally, we present a mechanistic model showing how the ARC2, NB, and N-terminal half of the LRR form a clamp, which regulates the dissociation and reassociation of the switch and sensor domains in NB-LRR proteins.

The effector SPRYSEC-19 of Globodera rostochiensis suppresses CC-NB-LRR-mediated disease resistance in plants.

The potato cyst nematode Globodera rostochiensis invades roots of host plants where it transforms cells near the vascular cylinder into a permanent feeding site. The host cell modifications are most likely induced by a complex mixture of proteins in the stylet secretions of the nematodes. Resistance to nematodes conferred by nucleotide-binding-leucine-rich repeat (NB-LRR) proteins usually results in a programmed cell death in and around the feeding site, and is most likely triggered by the recognition of effectors in stylet secretions. However, the actual role of these secretions in the activation and suppression of effector-triggered immunity is largely unknown. Here we demonstrate that the effector SPRYSEC-19 of G. rostochiensis physically associates in planta with the LRR domain of a member of the SW5 resistance gene cluster in tomato (Lycopersicon esculentum). Unexpectedly, this interaction did not trigger defense-related programmed cell death and resistance to G. rostochiensis. By contrast, agroinfiltration assays showed that the coexpression of SPRYSEC-19 in leaves of Nicotiana benthamiana suppresses programmed cell death mediated by several coiled-coil (CC)-NB-LRR immune receptors. Furthermore, SPRYSEC-19 abrogated resistance to Potato virus X mediated by the CC-NB-LRR resistance protein Rx1, and resistance to Verticillium dahliae mediated by an unidentified resistance in potato (Solanum tuberosum). The suppression of cell death and disease resistance did not require a physical association of SPRYSEC-19 and the LRR domains of the CC-NB-LRR resistance proteins. Altogether, our data demonstrated that potato cyst nematodes secrete effectors that enable the suppression of programmed cell death and disease resistance mediated by several CC-NB-LRR proteins in plants.

A genome-wide genetic map of NB-LRR disease resistance loci in potato.

Like all plants, potato has evolved a surveillance system consisting of a large array of genes encoding for immune receptors that confer resistance to pathogens and pests. The majority of these so-called resistance or R proteins belong to the super-family that harbour a nucleotide binding and a leucine-rich-repeat domain (NB-LRR). Here, sequence information of the conserved NB domain was used to investigate the genome-wide genetic distribution of the NB-LRR resistance gene loci in potato. We analysed the sequences of 288 unique BAC clones selected using filter hybridisation screening of a BAC library of the diploid potato clone RH89-039-16 (S. tuberosum ssp. tuberosum) and a physical map of this BAC library. This resulted in the identification of 738 partial and full-length NB-LRR sequences. Based on homology of these sequences with known resistance genes, 280 and 448 sequences were classified as TIR-NB-LRR (TNL) and CC-NB-LRR (CNL) sequences, respectively. Genetic mapping revealed the presence of 15 TNL and 32 CNL loci. Thirty-six are novel, while three TNL loci and eight CNL loci are syntenic with previously identified functional resistance genes. The genetic map was complemented with 68 universal CAPS markers and 82 disease resistance trait loci described in literature, providing an excellent template for genetic studies and applied research in potato.

Comparative sequence analysis of the potato cyst nematode resistance locus H1 reveals a major lack of co-linearity between three haplotypes in potato (Solanum tuberosum ssp.).

The H1 locus confers resistance to the potato cyst nematode Globodera rostochiensis pathotypes 1 and 4. It is positioned at the distal end of chromosome V of the diploid Solanum tuberosum genotype SH83-92-488 (SH) on an introgression segment derived from S. tuberosum ssp. andigena. Markers from a high-resolution genetic map of the H1 locus (Bakker et al. in Theor Appl Genet 109:146-152, 2004) were used to screen a BAC library to construct a physical map covering a 341-kb region of the resistant haplotype coming from SH. For comparison, physical maps were also generated of the two haplotypes from the diploid susceptible genotype RH89-039-16 (S. tuberosum ssp. tuberosum/S. phureja), spanning syntenic regions of 700 and 319 kb. Gene predictions on the genomic segments resulted in the identification of a large cluster consisting of variable numbers of the CC-NB-LRR type of R genes for each haplotype. Furthermore, the regions were interspersed with numerous transposable elements and genes coding for an extensin-like protein and an amino acid transporter. Comparative analysis revealed a major lack of gene order conservation in the sequences of the three closely related haplotypes. Our data provide insight in the evolutionary mechanisms shaping the H1 locus and will facilitate the map-based cloning of the H1 resistance gene.

Nucleocytoplasmic distribution is required for activation of resistance by the potato NB-LRR receptor Rx1 and is balanced by its functional domains.

The Rx1 protein, as many resistance proteins of the nucleotide binding-leucine-rich repeat (NB-LRR) class, is predicted to be cytoplasmic because it lacks discernable nuclear targeting signals. Here, we demonstrate that Rx1, which confers extreme resistance to Potato virus X, is located both in the nucleus and cytoplasm. Manipulating the nucleocytoplasmic distribution of Rx1 or its elicitor revealed that Rx1 is activated in the cytoplasm and cannot be activated in the nucleus. The coiled coil (CC) domain was found to be required for accumulation of Rx1 in the nucleus, whereas the LRR domain promoted the localization in the cytoplasm. Analyses of structural subdomains of the CC domain revealed no autonomous signals responsible for active nuclear import. Fluorescence recovery after photobleaching and nuclear fractionation indicated that the CC domain binds transiently to large complexes in the nucleus. Disruption of the Rx1 resistance function and protein conformation by mutating the ATP binding phosphate binding loop in the NB domain, or by silencing the cochaperone SGT1, impaired the accumulation of Rx1 protein in the nucleus, while Rx1 versions lacking the LRR domain were not affected in this respect. Our results support a model in which interdomain interactions and folding states determine the nucleocytoplasmic distribution of Rx1.

Identification and characterization of the most abundant cellulases in stylet secretions from Globodera rostochiensis.

Plant-parasitic cyst nematodes secrete cell wall modifying proteins during their invasion of host plants. In this study, we used a monoclonal antibody to immunopurify and to sequence the N terminus of the most abundant cellulases in stylet secretions of preparasitic juveniles of Globodera rostochiensis. The N-terminal amino acid sequence perfectly matched the sequence of an expressed sequence tag of two nearly identical genes, named Gr-eng3 and Gr-eng4, which show relatively low similarity with the previously identified Gr-eng1 and Gr-eng2 (i.e., 62% similarity and 42% identity). The recombinantly produced proteins from Gr-eng3 and Gr-eng4 demonstrated specific activity on carboxymethylcellulose, indicating that these genes encode active cellulases. To date, the cellulases in cyst nematodes are comprised of three possible domain structure variants with different types of ancillary domains at the C terminus of the glycosyl hydrolase family 5 (GHF5) domain. We used Bayesian inference to show that the phylogeny of the GHF5 domain based on currently available data suggest that the extant nematode cellulases arose through reshuffling of the GHF5 domain with different types of ancillary domains as relatively independent units. Knocking-down Gr-eng3 and Gr-eng4 using RNA interference resulted in a reduction of nematode infectivity by 57%. Our observations show that the reduced infectivity of the nematodes can be attributed to poor penetration of the host's root system at the onset of parasitism.

Origin, distribution and 3D-modeling of Gr-EXPB1, an expansin from the potato cyst nematode Globodera rostochiensis.

Southern analysis showed that Gr-EXPB1, a functional expansin from the potato cyst nematode Globodera rostochiensis, is member of a multigene family, and EST data suggest expansins to be present in other plant parasitic nematodes as well. Homology modeling predicted that Gr-EXPB1 domain 1 (D1) has a flat beta-barrel structure with surface-exposed aromatic rings, whereas the 3D structure of Gr-EXPB1-D2 was remarkably similar to plant expansins. Gr-EXPB1 shows highest sequence similarity to two extracellular proteins from saprophytic soil-inhabiting Actinobacteria, and includes a bacterial type II carbohydrate-binding module. These results support the hypothesis that a number of pathogenicity factors of cyst nematodes is of procaryotic origin and were acquired by horizontal gene transfer.