Skip to main content

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



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.


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.


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.


Plant-parasitic nematodes (PPN) cause severe damage to crops throughout the world and are an important constraint on delivering food security. Establishing the value of damage and control costs for these pathogens is difficult but has been calculated as being in excess of US$80 billion each year [1]. The largest economic losses are caused by the sedentary endoparasitic root-knot and cyst nematodes of the genera Meloidogyne and Heterodera/Globodera[2]. These nematodes have complex, biotrophic interactions with their hosts and induce the formation of feeding structures from which they derive all the food required for development. The feeding site needs to be maintained and protected from host defence responses, as well as from potential pathogens, for several weeks while the nematode matures to the adult stage.

The white potato cyst nematode, Globodera pallida, is an important pathogen of potato wherever it is grown. Globodera pallida originated in South America [3] and was introduced into Europe in the 19th century with wild potato material used for resistance breeding against late blight. It is now widely distributed in Europe and in regions that have imported seed potato from Europe [4]. Yield losses in excess of 50% have been reported and the lack of major gene resistance for G. pallida, coupled to the increasing legislative restrictions on the use of nematicides, mean that new control strategies are required.

Globodera pallida hatches as a second stage juvenile (J2) in response to diffusates from roots of suitable host plants. The J2 then locates the root, invades and migrates destructively through root cells until it reaches the inner cortex layers. During this migration the nematode uses the stylet to mechanically disrupt host cells and releases a cocktail of plant cell wall degrading enzymes that soften the cell wall [5]. At the inner cortex the behaviour of the nematode changes and it probes individual cells with the stylet until a cell that does not respond adversely is detected (reviewed in [6]). The nematode then secretes proteins into this initial syncytial cell which is transformed into a large multinucleate syncytium. Cell wall openings are formed, initially by widening of pre-existing plasmodesmata, followed by controlled breakdown of the plant cell wall in these regions. The cytoplasm of the initial syncytial cell proliferates, the central vacuole breaks down and the nucleus becomes enlarged. These changes are also observed in the cells surrounding the initial syncytial cell which implies communication between these cells. Eventually the protoplasts of the initial syncytial cell and its neighbours fuse at the cell wall openings and this process is repeated with further layers of cells until 200–300 cells are incorporated into the syncytium. Each nematode can only create one syncytium and depends on it for all the nutrients required for development to the adult stage, a process that takes 4–6 weeks. This prolonged period of biotrophy is almost unparalleled in plant-pathogen interactions and demonstrates a clear need for the nematode to suppress host defences throughout the life cycle.

The interactions of PPN, including G. pallida, are mediated by effectors – defined broadly here as secreted nematode proteins that manipulate the host to the benefit of the nematode. Effectors of PPNs are mainly secreted from two sets of gland cells, the dorsal and subventral gland cells, through the stylet into the host. These sets of gland cells show distinct developmental profiles. The two subventral gland cells are large and full of secretory granules in invasive stage J2, decrease in size and activity during the sedentary parasitic stages but become active again in adult male nematodes, which leave the root in order to locate females. In contrast, the dorsal gland cell is small in J2 but increases in size and activity in the sedentary stages [7]. It has therefore been widely suggested that proteins produced in the subventral gland cells are important for invasion, migration and processes occurring at the early stages of parasitism whereas those produced in the dorsal gland cell play a role in the later stages of the parasitic process.

Because of their key roles in parasitism, effectors of plant parasitic nematodes have been studied in detail. Effectors from cyst nematodes have been identified through EST sequencing e.g.[8], expression profiling [9] and, in particular, through sequencing of mRNA extracted from aspirated gland cell cytoplasm [10]. More recently, a protocol for purification of oesophageal gland cells and sequencing of purified mRNA has been described [11]. In each case, confirmation that a sequence of interest represents a candidate effector is obtained by in situ hybridisation to demonstrate expression in the subventral or dorsal oesophageal gland cells [12]. Following on from these projects, many effectors have been the subject of further functional characterisation. In addition to the cell wall modifying enzymes described above, effectors that play a wide variety of roles at various stages in the host parasite interaction have been described (reviewed in [13]. Effectors that suppress PAMP triggered immunity (PTI) and Effector triggered immunity (ETI) [14] have recently been identified from both cyst nematodes and root-knot nematodes [15, 16]. Other effectors may have roles in syncytium induction, either by targeting host auxin transporters [17] or by mimicking plant CLAVATA3 peptides [18]. An effector has also been identified that may play a role in the extensive remodelling of the host cell walls that takes place in the syncytium through its interaction with a host pectin methylesterase [19]. In other cases the host target of an effector has been identified along with a demonstration that the presence of the effector is of benefit to the nematode, although the precise function of the effector may remain uncertain (e.g.[20]).

Although enormous progress has clearly been made in terms of unravelling the functions of some nematode effectors, little genome scale analysis of the effectors present in any PPN has yet been reported. The genomes of two root knot nematodes, Meloidogyne incognita and M. hapla have been sequenced and some of the effectors present in each of these genomes were analysed as part of these projects [21, 22]. Further publications have described the large scale analysis of secreted proteins collected from one of these species [23, 24]. In addition, the genome and the secretome of a migratory endoparasitic nematode, Bursaphelenchus xylophilus have been sequenced [25, 26]. In contrast to cyst nematodes, B. xylophilus is not a biotrophic pathogen (in that it does not need to keep its host alive during the feeding process) and the effectors identified in this species were substantially different from those present in other PPN, with the exception of pectate lyases and expansins. This is likely to be a reflection of both the different modes of parasitism between B. xylophilus and cyst/root-knot nematodes and the taxonomic distance between them.

We recently reported the full genome sequence of G. pallida[27] including a fully replicated life cycle transcriptome analysis and a preliminary description of the effectors likely to be present in this species. Here we report the detailed genome scale analysis of the effector complement present in G. pallida, including analysis of expression profiles across the life cycle. We also examine the location of various effector gene families within the genome. In addition, we demonstrate the gland cell expression of effector candidates and in functional studies show that the effectors of G. pallida target a wide range of host subcellular compartments.

Results & discussion

Cell wall modifying proteins

PPN are almost unique among animals in having the capacity to metabolise the plant cell wall using endogenous cell wall degrading and modifying enzymes. The first effectors identified from PPN were cellulases, secreted from the subventral gland cells of cyst nematodes [28]. The similarity of these genes to bacterial cellulases and their absence from other animals led to the suggestion that they had been acquired from bacteria by horizontal gene transfer. Subsequent studies have led to the identification of numerous other horizontally acquired cell wall modifying proteins in a wide range of PPN, including cellulases, pectate lyases, xylanases, polygalacturonases, arabinogalactan galactosidases, arabinanases, carbohydrate binding modules (CBM) and expansin-like proteins (reviewed in [5]). The genome project for Meloidogyne incognita showed that more than 90 genes from seven protein families are present that could be involved in metabolism of the cell wall ([21] – Table 1). The M. hapla genome contains fewer genes albeit encoding a similar range of proteins [22], while non-plant parasites such as Caenorhabditis elegans or Brugia malayi have none (Table 1). The only noted exception to this rule is the presence of sequences similar to GH5 cellulases in Pristionchus pacificus[29] and these are thought to be part of a subclass of these enzymes that is distinct from those of the plant parasitic species. The G. pallida genome also contains a large number of genes that could encode cell wall modifying proteins ([27]; Table 1). The availability of the G. pallida transcriptome sequences allowed us to examine the expression profiles of such genes across the life cycle. The J2 (which needs to invade the host) and adult male (which needs to leave the root to locate females at the root surface) are the stages of the nematode that migrate through host tissues and many of the genes encoding cell wall modifying proteins, including a subset of the expansins and cellulases, were expressed specifically at these two life stages (Figure 1A and B). The pectate lyases that had detectable expression in the life stages tested were largely restricted to the J2 stage (Figure 1C), possibly also reflecting a role in migration. Some of the cellulases and expansins of G. pallida consist of a catalytic domain fused to a CBM. Like other PPN, G. pallida also has genes encoding proteins consisting of the CBM domain alone. All CBM protein domains present in G. pallida and other PPN, regardless of whether they are “stand alone” proteins or fused to cellulases or expansins are from CBM family 2. Some of these CBM proteins (not fused to any catylytic domain) are thought to be involved in migration. However, in H. glycines one CBM protein has been shown to interact with a host pectin methylesterase which is involved in the regulation of cell growth and expansion. This H. glycines CBM may thus be involved in syncytium expansion [19]. The CBMs present in G. pallida showed expression profiles reflecting these two functional roles; one CBM gene (GPLIN_000536400 – blue line in Figure 1D) is upregulated in J2 while another two CBM genes (GPLIN_000707900 and GPLIN_000706300 – green and red lines in Figure 1D) are upregulated at parasitic stages suggesting they could be the functional orthologues of the H. glycines genes involved in syncytium development (Figure 1D). Intriguingly, a similar pattern of expression was observed for the GH53 (arabinogalactan endo 1,4 β-galactosidase) genes; while one gene (GPLIN_000142900 – blue line in Figure 1E) was upregulated in J2 and particularly in males, expression of the other (GPLIN_000143000 – red line in Figure 1E) was restricted to the parasitic stages (Figure 1E). It was previously thought that nematode cell wall degrading enzymes are used solely for migration, with the plant’s own cell wall degrading machinery activated to allow the controlled modifications to the host cell wall that occur during syncytium formation [30, 31]. However, these expression data suggest that some nematode enzymes may also play a role in this process.

Table 1 Genes potentially encoding cell wall modifying proteins in various nematode species
Figure 1

Expression profiles of cell wall modifying proteins across the Globodera pallida life cycle. Life stages indicated on X axis, Y-axis indicates reads per kilobase per million. Each line represents the expression profile for an individual gene. (A) Expansins: Expression profiles of nine sequences are shown. GPLIN_000293700 (orange line), GPLIN_000536200 (darker green line) and GPLIN_000599100 (light blue line) are upregulated at J2 and adult male stages while the other six sequences (GPLIN_000092400, GPLIN_000293400, GPLIN_000590900, GPLIN_000599200, GPLIN_001571600 and GPLIN_001621500) are only weakly expressed at the life stages tested; (B) Cellulases: Expression profiles of fifteen sequences are shown. One sequence GPLIN_000536400 – dark blue) is upregulated at J2 and male, five sequences (GPLIN_001111200 – purple, GPLIN_001111300 – light blue, GPLIN_001215600 – light olive, GPLIN_000313600 – pink and GPLIN_000552400 – brick red) are upregulated at J2, two sequences (GPLIN_000755200 – orange, and GPLIN_000755100 – light blue) are upregulated at male only. The remaining cellulase-like sequences (GPLIN_000616300, GPLIN_000694900, GPLIN_000779000, GPLIN_779200, GPLIN_000827200, GPLIN_001185800 and GPLIN_000304900) are weakly expressed at the life stages shown here; (C) Pectate Lyases: Expression profiles of seven sequences are shown. Three sequences (GPLIN_000467400 – blue line, GPLIN_000142600 and GPLIN_000412300 – overlaying graphs appearing green in the figure) are upregulated at J2 while four others (GPLIN_000673000, GPLIN_000294400, GPLIN_000294500 and GPLIN_000322300) are weakly expressed at the stages shown here; (D) Carbohydrate Binding Modules; (E) GH53 arabinogalactan endo 1,4 ß galactosidases. For Figure 1E the left hand Y axis shows scale of expression for GPLIN_000142900 (blue line) while the right hand Y axis shows scale of expression for GPLIN_000143000 (red line).

Orthologues of other previously characterised nematode effectors

The G. pallida predicted protein sequences were BLAST searched with a list of previously characterised effectors from the root-knot nematode M. incognita and the cyst nematodes G. rostochiensis and H. glycines. This analysis revealed 128 putative G. pallida orthologues of 37 effectors from these species (Additional file 1: Table S3). In 6 cases (G. rostochiensis sequences A42, 66P1 and 747, and H. glycines effectors Hgg6, G8H07 and G28B03), G. pallida orthologues were not present among the gene predictions from the final assembly but were present as uncalled sequences in the assembly and/or identified from transcriptome sequences. Several of these sequences were cloned from cDNA in order to ensure that they were genuinely present in G. pallida and for further functional analysis (below). In addition, 2 gene families similar to effectors IA7 and IVG9, previously identified from G. pallida were also present, containing 7 and 5 members respectively (Additional file 1: Table S3). A BLAST search of the G. rostochiensis assembled transcriptome (Eves van den Akker, unpublished) showed that, as expected, sequences similar to all of those identified in this analysis were present in G. rostochiensis. In addition to these sequences, two substantial gene families were present that could encode homologues of SPRYSEC effectors and orthologues of a sequence annotated as “Heterodera avenae gland cell protein” (Genbank HM147943.1). These families have been described in detail as part of the G. pallida genome analysis and are not considered further here. Sixteen effectors from other cyst nematodes appeared to have no orthologue in the current assembly of the G. pallida genome or transcriptome (Additional file 2: Table S4). As previously described [27] almost no overlap (other than the cell wall modifying proteins and chorismate mutase) was found with effectors from M. incognita.

In order to demonstrate that the cyst nematode orthologues identified represent genuine effector candidates rather than secreted proteins with functions within the nematode body, we performed in situ hybridisation to confirm expression of selected genes in the gland cells of the nematode. This analysis identified genes expressed in the subventral gland cells and in the dorsal gland cell (Figure 2). Staining patterns in G. pallida matched those reported for the H. glycines orthologues, where these had previously been described ([10]; Figure 2A,C). In common with all other SPRYSEC sequences examined to date, the G. pallida SPRYSEC sequences tested were expressed in the dorsal gland cell (two examples shown in Figure 2H and I). In addition to confirming that the products of these genes represent genuine effector candidates, this analysis can provide information about the part of the life cycle in which the effectors are important as the subventral gland cells are thought to produce proteins required at the early stages of the parasitic process while the dorsal gland cells produce proteins that are important at later stages. On a larger scale, we were able to undertake a cluster analysis of the temporal expression profiles of all effector candidates using the RNAseq data generated for G. pallida. This analysis revealed that five different clusters of genes sharing similar expression profiles were present: J2 (30 sequences), J2 and male (5 sequences), parasitic stages (61 sequences), constitutive (20 sequences) and parasitic and male (4 sequences) (Figure 3). Comparing the in situ hybridisation data with this cluster information gave conflicting results. All the genes expressed in the subventral gland cells showed expression peaks at J2 (e.g. GPLIN_000662500 and many of the cell wall degrading enzymes), suggesting that expression in these tissues is usually indicative of a role at the earliest stages of the host-parasite interaction. However, genes expressed in the dorsal gland cell showed a variety of different temporal expression profiles with some showing a clear peak of expression at the J2 stage while others peaked in expression during the parasitic stages as expected. This implies that control of the temporal expression of effectors is complex and is not, as previously suggested, simply a question of expression in one gland cell or another.

Figure 2

Localisation of gene expression by in situ hybridisation of Globodera pallida SPRYSEC genes and orthologues of effectors from other PPN. Sections of G. pallida preparasitic stage 2 juveniles (B, C, E, H and I) or 14dpi parasitic nematodes (A, D, F and G) were incubated with antisense probes based on the DNA coding sequence for the following gene loci: (A) GPLIN_001203000 (ortholog of Heterodera glycines 10C02 sequence); (B) GPLIN_000854400 (ortholog of H. glycines 16H02 sequence); (C) GPLIN_000662500 (ortholog of H. glycines G20E03 sequence); (D) GPLIN_000235400 (ortholog of Globodera rostochiensis 1106 sequence); (E) GPLIN_000243800 (ortholog of H. glycines 4D06 sequence); (F) G. pallida ortholog of G. rostochiensis A42 sequence; (G) GPLIN_000201400 (G. pallida ortholog of G. rostochiensis E9 sequence); (H) SPRYSEC GPLIN_001082900; (I) SPRYSEC GPLIN_000657200. st, stylet; dg, dorsal gland cell; svg, subventral gland cells; m, metacorpal bulb.

Figure 3

Cluster analysis of co-regulated candidate effectors during Globodera pallida life cycle. The y-axis represents fold change in expression values, determined by calculating fold changes over mean expression values across all samples from RNAseq data. Each cluster contains more than one sub group of genes with each sub group showing subtle differences in expression profile.

G. pallidaeffectors target various host subcellular structures

Understanding the subcellular localisation of effectors can provide information about their putative function(s), can be used to prioritise future functional studies and is also of value when analysing the results of yeast two-hybrid screens that aim to identify targets of effectors. Effectors from a range of plant pathogens have been shown to target diverse subcellular structures in plant cells. Several studies have shown in particular that some effectors from either cyst or root-knot nematodes can target the host cell nucleus and/or nucleolus, suggesting that PPN may be hijacking key nuclear functions in their host [8, 3234].

Using eGFP-fusions in transient expression assays in Nicotiana benthamiana, we observed that putative effectors from G. pallida targeted a wide range of plant cell structures. While many were localised throughout the cytoplasm in a manner similar to free eGFP (Figure 4A and B), others showed more restricted localisation. Although many effectors are small enough to allow passive diffusion into the nucleus as eGFP fusions (e.g. Figure 4B), these are excluded from the nucleolus. However, several effectors were identified that were localised in the nucleolus (Figure 4C and D) and these are likely to represent genuine localisations rather than passive diffusion. These two effectors have nuclear localisation signals (NLS) predicted by PSORT. One SPRYSEC effector (GPLIN_001465500) looked to be excluded from the cytoplasm but was localised specifically in the nucleus and was often seen surrounding the nucleolus (Figure 4E and F). No NLS was predicted for this effector sequence. Similar localisations have been observed for some viral proteins [35]. Effectors were also identified that targeted peroxisomes and the peroxisome membrane (Figure 4G-I). This localisation was confirmed using an mRFP-tagged marker for peroxisomes.

Figure 4

Subcellular localisations of Globodera pallida effector fusions with the enhanced green fluorescent protein (eGFP) in leaves of Nicotiana benthamiana . (A) Localisation of free eGFP. (B) Localisation of GPLIN_000015300 (ortholog of Heterodera glycines G7E05 sequence) in the cytoplasm and nucleoplasm but excluded from the nucleolus. (C) Localisation of GpA42 _ (ortholog of Globodera rostochiensis A42 sequence) in the nucleoli. (D) Faint expression of GPLIN_000235400 (ortholog of G. rostochiensis 1106 sequence) in the cytoplasm with accumulation in the nucleus, particularly in the nucleoli. (E & F) Exclusion of SPRYSEC effector GPLIN_001465500 from the cytoplasm but accumulation in the nucleoplasm, particularly surrounding the nucleolus. (G-H) Expression of GPLIN_000662500 (ortholog of H. glycines G20E03 sequence) in peroxisomes. (I-J) Expression of GPLIN_000457000 (ortholog of H. glycines hgsec4 sequence) in the peroxisome membrane (arrows). In all images the eGFP-effector fusion is seen in green, with monomeric-red fluorescent protein (mRFP) co-label in magenta, from either transgenic expression of an mRFP-histone fusion (A-F) or co-expression of an mRFP-tagged peroxisome marker (G-J), and autofluorescent chloroplasts displayed in blue. Scale bars represent 50 μm for A, D, 20 μm for C, 10 μm for G, H, and 5 μm for E, F, I and J. Localisation of Agrobacterium-mediated transient expression of eGFP-effector fusions was observed 48 h post inoculation by confocal microscopy.

Globodera pallida induces profound gene expression changes in its host when establishing a feeding site. It is possible that effectors which target the nuclei may manipulate gene expression directly either as a part of the process of establishing the feeding site or in order to suppress defence signalling pathways. Alternatively, such effectors may act indirectly, targeting host transcription factors or other nuclear proteins, which may interfere with nuclear functions. Effectors from a range of cyst [32] and root-knot nematodes [33] that localise to the nucleus have been identified, demonstrating the importance of targeting this structure for a range of biotrophic PPN groups (reviewed in [36]), as is the case for other plant pathogens. By contrast, this is the first report of any nematode effector that targets the peroxisome. Indeed, no direct evidence for targeting of this organelle by any plant pathogen has been reported previously, although a peroxisome targeting signal has been identified in the DspA/E effector of Erwinia amylovora[37]. Plant peroxisomes are involved in a variety of metabolic processes, including responses to abiotic stress and production of auxins [38], both of which are clearly of relevance to the G. pallida life cycle. For example, part of the biosynthetic pathway for jasmonic acid occurs in the peroxisomes and a mutant line lacking one of the enzymes involved in this process in Solanum lycopersicum is impaired in its ability to mount defence responses against some insects [39]. The peroxisome is also an important site for production of hydrogen peroxide [40], which is known to be deployed as part of the defence response against PPN [41]. Less is known about the nature and roles of proteins present in the peroxisomal membrane, although several are known to control the import/export of other proteins and metabolites into/from the peroxisome (reviewed in [42]. Information on the effectors from G. pallida that target the peroxisomes will be useful in terms of framing future functional studies.

Novel effector candidates

The analysis performed as part of the G. pallida genome project [27], identified 117 potential novel effectors (Additional file 3: Table S5). These predicted proteins have a signal peptide, lack a transmembrane domain and were upregulated at J2 (vs egg) or at early parasitic stages (vs J2). We have subsequently performed further analysis of these 117 proteins. We first examined the expression profiles of each of the candidates across the G. pallida life cycle as only the expression at two life stages was considered during the identification process. This analysis showed that, as for the orthologues of effectors from other species (above), the expression profiles formed five clusters (Figure 5): J2 only (28 sequences), J2 and parasitic (46 sequences), J2 and male (8 sequences), parasitic (4 sequences) and parasitic and male (31 sequences). The candidate effectors were also analysed for the presence of Pfam domains. Twenty seven potential Pfam domains were identified for 11 of the 117 predicted proteins (Table 2) while the other 106 sequences contained no known domains and, since these also have no BLAST matches against NR using the stringency requirements described here, they therefore represent novel proteins. A BLAST search of assembled G. rostochiensis transcripts showed that all but two of the G. pallida novel effector candidates had strong matches in G. rostochiensis (not shown). Since the transcriptome assembly (generated from J2s and parasitic G. rostochiensis 14 days after infection) is likely to be incomplete, compared to the full G. pallida genome, it is reasonable to conclude that the “novel” effector candidates do not include any G. pallida specific sequences. By contrast, only 17 of the sequences had matches in the M. incognita proteins predicted from the genome sequence, suggesting that the dataset is enriched in proteins specific to cyst nematodes.

Figure 5

Expression profiles of “novel” effector candidates across the G. pallida life cycle. The Y-axis represents reads per kilobase per million. Clusters represent: (A) J2 specific; (B) J2 and parasitic; (C) J2 and male; (D) parasitic only; (E) parasitic and male.

Table 2 PFam analysis of “novel” candidate effectors

In order to confirm that at least some of the novel effector candidates represented genuine effectors we performed in situ hybridisation to determine the site of expression of some of the candidates. In situ hybridisation with nematodes is a labour intensive process meaning that it was not feasible to analyse most of the sequences. However, two of the “novel” candidate effectors that were upregulated at J2 (GPLIN_000333000 and GPLIN_000834600) were expressed in the subventral and dorsal gland cells respectively (Figure 6A and B). GPLIN_000333000 is part of a small family (~7 sequences) of related proteins, each of which has a signal peptide. These data therefore suggest that it represents a novel effector gene family in G. pallida.

Figure 6

Localisation of “novel” G. pallida candidate effectors expression by in situ hybridisation to preparasitic second stage juveniles (J2s). Sections of nematodes were incubated with antisense probes designed based on DNA coding sequence for the following gene loci (A) GPLIN_000333000; (B) GPLIN_000834600; dg, dorsal gland cell; m, metacorpal bulb; svg, subventral gland cell.

Genomic organisation of effectors and diversification by alternative splicing

Several of the G. pallida effector gene families are considerably expanded and this may be a reflection of selection pressure to avoid detection or maintain function within a host that is under selection pressure to evade infection. For example, one G. pallida effector gene family, similar to previously identified effectors from H. glycines (4D06 and several other related sequences) consists of 39 members in G. pallida. However, the most striking examples of expansion are found in the SPRY domain proteins and a family of proteins similar to a sequence annotated as “H. avenae dorsal gland protein” in the database. These families consist of 299 and 295 members respectively. We analysed clustering of these sequences in the G. pallida genome and found that 180/295 of the sequences similar to “H. avenae dorsal gland protein” are clustered (i.e. present on the same G. pallida scaffold) with the largest cluster consisting of 10 similar genes in a group of 14 genes on scaffold 299 (Additional file 4: Figure S1). Similarly, 152/299 of the SPRY domain proteins are located in clusters present on individual scaffolds, with the largest cluster consisting of 7 genes on scaffold 141 (Additional file 4: Figure S1). These estimates of clustering are likely to be significant underestimates, given the fragmented nature of the G. pallida genome sequence. These data suggest that the SPRY domain and “dorsal gland protein” gene families are likely to have expanded, at least in part, by a process of tandem duplication followed by diversification.

It has recently been shown that for several plant pathogens, including the late blight pathogen Phytophthora infestans, effector sequences are more likely to be found in repeat-rich, gene sparse regions of the genome [43]. This property has been used to identify novel candidate virulence factors [44]. Calculating the distance between each gene and its nearest neighbours (the flanking intergenic regions) and comparing these values for all genes in a genome, can be used as a way of calculating whether a gene resides in a gene rich or gene sparse environment [45]. We therefore examined whether or not this pattern is also true in the G. pallida genome with a particular focus on the SPRYSECs as the nature of this substantial gene family and the absence of a similarly expanded family in other cyst nematodes suggested that it may be in the process of rapid evolution. Although this type of analysis is challenging for genomes as fragmented as that of G. pallida, with two data points being lost at the end of each contig, the pattern reported for other plant pathogens did not seem to be followed in G. pallida, with the distribution patterns of the effectors not significantly different from those of other genes (Additional file 5: Figure S2).

Alternative splicing generates multiple mRNA transcripts from a single precursor and can lead to changes in RNA or protein levels, or can give rise to different forms of related proteins from a single gene. These may differ in terms of function or cellular localisation (reviewed in [46]). In addition, effectors are in direct contact with the host and are therefore likely to be under strong diversifying selection pressure in order to produce variants that evade detection. Alternative splicing is one mechanism by which this may be achieved and has been described previously for one cyst nematode effector [47]. We analysed the G. pallida effectors similar to those from other nematode species (orthologues of other previously characterised nematode effectors as described above) for evidence of alternative splicing. Mapping of transcripts to the predicted genome sequence for these genes indicated that at least 38% undergo alternative splicing. Some effector gene families showed evidence for more frequent alternative splicing: Nineteen of 39 members of the G. pallida gene family similar to the H. glycines “4D06” genes and 3 of 4 G. pallida orthologues of the H. glycines G20E03 sequence showed evidence for alternative splicing. Evidence was obtained for alternative splicing occurring within a single life stage (Additional file 6: Figure S3A). Intriguingly, evidence supporting the generation of different mRNAs at different life stages was also obtained (Additional file 6: Figure S3B). Although functional studies are required in order to assess the relevance of this observation, this raises the possibility that alternative splicing generates forms of a protein with different functional properties at different life stages, or that this may be another mechanism by which effector activity is regulated. In addition, these data suggest that alternative splicing may be relatively widespread within effector sequences of cyst nematodes.

The data presented here demonstrate that G. pallida deploys a significant number of effectors in order to promote its biotrophic interactions with the host. These effectors are under precise temporal regulation and target a variety of host subcellular structures and, most likely, a variety of host processes. The availability of these sequences will greatly facilitate future functional studies on this, and related cyst nematodes and enormously expand the range of potential targets for novel control strategies based, for example, on RNA-interference [48]. A proof of concept study targeting an effector of root-knot nematode [49] demonstrated the feasibility of such approach in the model plant Arabidopsis, where it conferred broad resistance to several Meloidogyne species.


  • The genome of G. pallida has several hundred genes that could encode putative effectors. Further diversity in effectors is generated through alternative splicing. In the present G. pallida assembly there was no evidence to suggest that effectors are preferentially localised in gene sparse regions of the genome.

  • A bioinformatic pipeline combining analysis of the presence of a signal peptide, absence of a transmembrane domain and expression profiling can be used to identify novel effectors.

  • Effectors display distinct temporal expression profiles across the G. pallida life cycle, suggesting that different effectors are deployed in order to support different phases of the parasitic process.

  • Effectors of G. pallida target a range of host structures including the nucleus and peroxisomes.


G. pallidaresources & alternative splicing analysis

The biological material used for both genome and transcriptome analysis originate from a standard Pa2/3 pathotype of G. pallida population “Lindley” which is held at the James Hutton Institute, Dundee, UK [50]. The G. pallida predicted protein set version 1.0 (16th May 2012) was used for identification of effectors. This protein set is available at and was used for the detailed analysis of the G. pallida genome [27]. Expression profiles of effectors across the life cycle were determined analysing the RNAseq information available for G. pallida[27]; replicated RNAseq datasets from eggs (containing unhatched J2), invasive stage J2, parasitic nematodes at 7, 14, 21, 28 and 35 days post infection (dpi) and adult males. Distribution of effectors in gene sparse versus gene rich regions was analysed as described in [45]. Scaffold drawings indicating the direction and location of genes were produced using a combination of Biopython [51] and Genome Diagram [52].

For analysis of alternative splicing, genomic regions of the effectors present in other nematode species (referred to as “orthologue effectors” below) were first extracted from the genome sequence and compared to the orthologues using BLAT. Life stage specific de novo transcriptome assemblies for the J2 and 7dpi RNAseq data were generated using Trinity as previously described [53] and compared to the genomic regions using BLAT with 90% similarity cut off. All tracks were visualised in IGV (Integrative Genome Viewer [54]. Potential alternative splicing events, where more than one transcript mapped to the same genomic copy, were checked manually to assess the impact of putative alternative splice events.

Identification of G. pallidaorthologues of previously characterised effectors

A list of known effectors from other PPN was collated using data from H. glycines gland cell ESTs [10, 55], microarray analysis [56], effectors identified from cDNA-AFLP analysis on G. rostochiensis[9], G. rostochiensis and G. pallida ESTs [57] and effectors identified from M. incognita[58]. The list also included effectors that had previously been identified from G. pallida[59, 60]. In addition, a list of G. rostochiensis effectors was provided by Dr. G. Smant (Wageningen University, Netherlands). The collated effector list, consisting of 133 G. rostochiensis sequences, 53 H. glycines sequences, three G. pallida sequences and 35 M. incognita sequences, was subjected to a local, command line BLAST [61] against the G. pallida genome sequence. This search used an E-value threshold of 10−5 with low complexity filtering turned off.

CAZymes and other cell wall modifying proteins

The CAZymes Analysis Toolkit (CAT) [62]; was used to identify putative carbohydrate active enzymes (CAZymes) with a predefined CAZyme database on the G. pallida predicted protein set. Putative CAZymes were manually annotated using a combination of BLASTP Vs NR database, NCBI’s Conserved Domain Database service [63] and InterProScan [64] to determine to presence of the catalytic domains. Genes of interest were identified by parsing the CAT output files. Databases of expansins and carbohydrate binding module (CBM) genes from other PPN were used for BLASTP searches against the G. pallida genome.

Identification of novel candidate effectors

The predicted G. pallida protein set was first analysed using a standard secretory protein identification protocol. Proteins that had a predicted signal peptide and no transmembrane domain were identified using SignalP 3.0 [65] followed by TMHMM [66], based on the methodology used in [8] using the Galaxy tools and workflow described in [67]. Expression profiles of the genes that passed these filters were then analysed using DESeq ( [68] in order to identify genes that were significantly upregulated in J2 compared to eggs or upregulated at 7 dpi compared to J2. Genes that passed this expression profiling filter were then BLAST searched against the NR database and those that obviously had functions unrelated to parasitism (e.g. collagens, digestive proteinases) were manually removed. In some cases the results of this BLAST searching provided functional information about the novel putative effectors. The putative effector list was thus analysed for any known domains using Pfam rules defined in (July 2012), using HMMER [69]. All sequences passing these filters were then BLAST searched against the proteins predicted in the M. incognita genome sequence [21] and against an assembled transcriptome derived from an RNAseq dataset of G. rostochiensis J2 and parasitic nematodes at 14 dpi (S. Eves van den Akker, unpublished). The sequences were also BLAST searched against dbEST in order to identify matches in other plant parasitic nematodes.

Analysis of effector expression profiles

The expression profiles of putative effectors were analysed using the normalised RNAseq data generated as part of the G. pallida genome project. The MBClusterseq program ( was used to separate the effectors with similar expression profiles into clusters. Inspection of the results of this analysis revealed that some clusters showed very similar patterns and genes in such clusters were subsequently merged into the same cluster.

In situhybridisation

The spatial expression patterns of some candidate effectors were examined by in situ hybridisation as previously described [60]. The sequences of primers used to amplify fragments of these genes for probe synthesis are provided in Additional file 7: Table S1.

Cloning and characterisation of effectors

Messenger RNAs were isolated from J2 or parasitic stage G. pallida (population Lindley) using a Dynabeads mRNA Direct Micro kit (Invitrogen) and treated with RQ1 DNase (Promega). cDNA was synthesised from approximately 400 ng purified mRNA using the Superscript III system (Invitrogen) with poly(dT) primers following the manufacturer’s instructions. For cloning, the coding sequences of selected effector candidates were amplified by PCR from cDNA, excluding the predicted signal peptide sequence but with the ACCATG leader sequence and a stop codon in the forward and reverse primer respectively (Additional file 8: Table S2). PCR was performed using the proof reading KOD DNA polymerase (Novagen) and products were resolved on 1.5% (w/v) agarose gels. Amplification products of the expected size were purified from gels using the QIAquick Gel Extraction Kit (QIAGEN) and inserted into the pCR8/GW/TOPO Gateway ENTRY vector by TA cloning following the manufacturer’s instructions (Invitrogen). Clones were subsequently recombined into the binary pK7WGF2 expression vector [70] for fusion with the enhanced green fluorescent protein (eGFP) tag using LR clonase (Invitrogen) following the manufacturer’s instructions. The integrity of the effector sequence in both ENTRY clones and in the destination vectors, as well as the fusion with the eGFP were confirmed by sequencing. For Agrobacterium-mediated transient expression assays, the eGFP-fusion expression vectors (Spectinomycin selection) were transferred by electroporation to Agrobacterium tumefaciens strain GV3101 that contains a helper vector encoding virGN54D (Gentamycin selection) [71].

Transient expression and analysis of subcellular localisation in planta

For subcellular localisation of the eGFP-effector fusions in planta, the constructs were transiently expressed in leaves of 4-week-old Nicotiana benthamiana using Agrobacterium-mediated transformation. Agrobacterium clones were grown overnight at 28°C in 5 mL Luria Bertani (LB) medium containing 25 μg/L Gentamycin and 100 μg/mL Spectinomycin. Bacterial cells were pelleted by centrifugation, rinsed and resuspended in infiltration buffer containing 10 mM MgCl2, 10 mM MES (2-[N-Morpholino] ethane sulfonic acid), and 200 μM acetosyringone, and adjusted to an optical density at 600 nm (OD600nm) of 1. Bacteria were then incubated for at least 3 h in the dark at room temperature prior to further dilution in infiltration buffer to OD600nm of 0.02 per construct and infiltration on the abaxial side of the leaves using a 1-mL needleless syringe.

For co-localisation analysis, bacteria were either infiltrated in leaves of transgenic N. benthamiana line (CB157) expressing a nuclear histone marker fused to mRFP (mRFP-H2B) [72] or co-infiltrated into wild-type plants together with an A. tumefaciens clone containing a peroxisome marker fused to mRFP (PfluB4; 50 μg/mL Kanamycin selection) [73]. Localisations were imaged 48 h post inoculation using either a Zeiss LSM 710 or a Leica SP2 confocal laser-scanning microscope. eGFP was imaged with an excitation wavelength (λ) of 488 nm and emission at λ495-530 nm (λ505-530 nm for SP2). Autofluorescence from chlorophyll generated by excitation at this wavelength was collected at λ657-737 nm (SP2 λ650-700 nm). mRFP was imaged sequentially with an excitation at λ561 nm and emission at λ592-632 nm (SP2 λ580-610 nm). Western blots were performed using standard protocols to analyse the size of expressed proteins – in all cases where a signal was detected the size of the band was in agreement with the predicted size of the eGFP-effector fusion protein with no evidence for degradation (not shown).

Availability of supporting data

Sequence data analysed in this paper were part of the G. pallida genome project. All data generated in this project have been submitted to the Genbank database under the accession number PRJEB123. Data and annotation have been submitted to Wormbase and are available at The G. pallida genome assembly and functional annotation is available from and via GeneDB at


  1. 1.

    Nicol JM, Turner SJ, Coyne DL, den Nijs L, Hockland S, Maafi ZT: Current Nematode Threats to World Agriculture. Genomics and Molecular Genetics of Plant-Nematode Interactions. Edited by: Jones JT, Gheysen G, Fenoll C. 2011, Dordrecht: Springer, 21-43.

    Chapter  Google Scholar 

  2. 2.

    Jones JT, Haegeman A, Danchin EGJ, Gaur HS, Helder J, Jones MGK, Kikuchi T, Manzanilla-Lopez R, Palomares-Rius J, Wesemael WML, Perry RN: Top 10 plant parasitic nematodes in molecular plant pathology. Mol Plant Pathol. 2013, 14: 946-961. 10.1111/mpp.12057.

    PubMed  Article  Google Scholar 

  3. 3.

    Evans K, Franco J, Descurrah MM: Distribution of species of potato cyst nematodes in South America. Nematologica. 1975, 21: 365-369. 10.1163/187529275X00103.

    Article  Google Scholar 

  4. 4.

    Hockland S, Niere B, Grenier E, Blok V, Phillips M, Den Nijs L, Anthoine G, Pickup J, Viane N: An evaluation of the implications of virulence in non-European populations of Globodera pallida and G. rostochiensis for potato cultivation in Europe. Nematology. 2012, 14: 1-13. 10.1163/138855411X587112.

    Article  Google Scholar 

  5. 5.

    Haegeman A, Jones JT, Danchin E: Horizontal gene transfer in nematodes: a catalyst for plant parasitism?. Mol Plant-Microbe Interact. 2011, 24: 879-887. 10.1094/MPMI-03-11-0055.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Sobczak M, Golinowski W: Cyst Nematodes and Syncytia. Genomics and Molecular Genetics of Plant-Nematode Interactions. Edited by: Jones JT, Gheysen G, Fenoll C. 2011, Dordrecht: Springer, 61-82.

    Chapter  Google Scholar 

  7. 7.

    Hussey RS, Mimms CW: Ultrastructure of esophogeal glands and their secretory granuls in the root-knot nematode Meloidogyne incognita. Protoplasma. 1990, 162: 99-107.

    Article  Google Scholar 

  8. 8.

    Jones JT, Kumar A, Pylypenko LA, Thirugnanasambandam A, Castelli L, Chapman S, Cock PJ, Grenier E, Lilley CJ, Phillips MS, Blok VC: Identification and functional characterization of effectors in expressed sequence tags from various life cycle stages of the potato cyst nematode Globodera pallida. Mol Plant Pathol. 2009, 10: 815-828. 10.1111/j.1364-3703.2009.00585.x.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Qin L, Overmars B, Helder J, Popeijus H, van der Voort JR, Groenink W, van Koert P, Schots A, Bakker J, Smant G: An efficient cDNA-AFLP-based strategy for the identification of putative pathogenicity factors from the potato cyst nematode Globodera rostochiensis. Mol Plant-Microbe Interact. 2000, 13: 830-836. 10.1094/MPMI.2000.13.8.830.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Gao BL, Allen R, Maier T, Davis EL, Baum TJ, Hussey RS: The parasitome of the phytonematode Heterodera glycines. Mol Plant-Microbe Interact. 2003, 16: 720-726. 10.1094/MPMI.2003.16.8.720.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Maier TR, Hewezi T, Peng JQ, Baum TJ: Isolation of whole esophageal gland cells from plant-parasitic nematodes for transcriptome analyses and effector identification. Mol Plant-Microbe Interact. 2013, 26: 31-35. 10.1094/MPMI-05-12-0121-FI.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    De Boer JM, Yan Y, Smant G, Davis EL, Baum TJ: In-situ hybridization to messenger RNA in Heterodera glycines. J Nematol. 1998, 30: 309-312.

    CAS  PubMed Central  PubMed  Google Scholar 

  13. 13.

    Haegeman A, Mantelin S, Jones JT, Gheysen G: Secretions of plant parasitic nematodes. Gene. 2012, 492: 19-31. 10.1016/j.gene.2011.10.040.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Jones JDG, Dangl JL: The plant immune system. Nature. 2006, 444: 323-329. 10.1038/nature05286.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Jaouannet M, Magliano M, Arguel MJ, Gourges M, Evangelisti E, Abad P, Rosso MN: The root-knot nematode calreticulin Mi-CRT is a key effector in plant defense suppression. Mol Plant-Microbe Interact. 2013, 26: 97-105. 10.1094/MPMI-05-12-0130-R.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Postma WJ, Slootweg EJ, Rehman S, Finkers-Tomczak A, Tytgat TOG, van Gelderen K, Lozaon-Torres JL, Roosien J, Pomp R, van Schaik C, Bakker J, Goverse A, Smant G: The effector SPRYSEC-19 of Globodera rostochiensis suppresses CC-NB-LRR mediated disease resistance in plants. Plant Physiol. 2012, 160: 944-954. 10.1104/pp.112.200188.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  17. 17.

    Lee C, Chronis D, Kenning C, Peret B, Hewezi T, Davis EL, Baum TJ, Hussey RS, Bennet M, Mitchum MG: The novel cyst nematode effector protein 19C07 interacts with the Arabidoipsis auxin influx transporter LAX3 to control feeding site development. Plant Physiol. 2011, 155: 866-880. 10.1104/pp.110.167197.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  18. 18.

    Wang X, Mitchum MG, Gao B, Li C, Diab H, Baum TJ, Hussey RS, Davis EL: A parasitism gene from a plant-parasitic nematode with function similar to CLAVATA3/ESR (CLE) of Arabidopsis thaliana. Mol Plant Pathol. 2005, 6: 187-191. 10.1111/j.1364-3703.2005.00270.x.

    PubMed  Article  Google Scholar 

  19. 19.

    Hewezi T, Howe P, Maier T, Hussey RS, Goellner-Mitchum M, Davis EL, Baum TJ: Cellulose binding protein from the parasitic nematode Heterodera schachtii interacts with Arabidopsis pectin methylesterase: cooperative cell wall modification during parasitism. Plant Cell. 2008, 20: 3080-3093. 10.1105/tpc.108.063065.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  20. 20.

    Hewezi T, Howe P, Maier T, Hussey RS, Goellner-Mitchum M, Davis EL, Baum TJ: Arabidopsis spermidine synthase is targeted by an effector protein of the cyst nematode Heterodera schachtii. Plant Physiol. 2010, 152: 968-984. 10.1104/pp.109.150557.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  21. 21.

    Abad P, Gouzy J, Aury J-M, Castagnone-Sereno P, Danchin EGJ, Deleury E, Perfus-Barbeoch L, Anthouard V, Artiguenave F, Blok VC, Caillaud M-C, Coutinho PM, Dasilva C, De Luca F, Deau F, Esquibet M, Flutre T, Goldstone JV, Hamamouch N, Hewezi T, Jaillon O, Jubin C, Leonetti P, Magliano M, Maier TR, Markov GV, McVeigh P, Pesole G, Poulain J, Robinson-Rechavi M, et al: Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita. Nat Biotechnol. 2008, 26: 909-915. 10.1038/nbt.1482.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Opperman CH, Bird DM, Williamson VM, Rokhsar DS, Burke M, Cohn J, Cromer J, Diener S, Gajan J, Graham S, Houfek TD, Liu Q, Mitros T, Schaff J, Schaffer R, Scholl E, Sosinski BR, Thomas VP, Windham E: Sequence and genetic map of Meloidogyne hapla: a compact nematode genome for plant parasitism. Proc Natl Acad Sci U S A. 2008, 105: 14802-14807. 10.1073/pnas.0805946105.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  23. 23.

    Bellafiore S, Shen Z, Rosso M-N, Abad P, Shih P, Briggs SP: Direct identification of the Meloidogyne incognita secretome reveals proteins with host cell reprogramming potential. PLoS Pathog. 2008, 4 (10): e1000192-10.1371/journal.ppat.1000192.

    PubMed Central  PubMed  Article  Google Scholar 

  24. 24.

    Wang X, Moreno YA, Wu H, Ma C, Li Y, Zhang J, Yang C, Sun S, Ma W, Geary TG: Proteomic profiles of soluble proteins from the esophageal gland in female Meloidogyne incognita. Int J Parasitol. 2012, 42: 1177-1183. 10.1016/j.ijpara.2012.10.008.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Kikuchi T, Cotton JA, Dalzell JJ, Hasegawa K, Kanzaki N, McVeigh P, Takanashi T, Tsai IJ, Assefa SA, Cock PJA, Otto TD, Hunt M, Reid AJ, Sanchez-Flores A, Tsuchihara K, Yokoi T, Larsson MC, Miwa J, Maule AG, Sahashi N, Jones JT, Berriman M: Genomic insights into the origin of parasitism in the emerging plant pathogen Bursaphelenchus xylophilus. PLoS Pathog. 2011, 7: e1002219-10.1371/journal.ppat.1002219.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  26. 26.

    Shinya R, Morisaka H, Kikuchi T, Takeuchi Y, Ueda M, Futai K: Secretome analysis of the pine wood nematode Bursaphelenchus xylophilus reveals the tangled roots of parasitism and its potential for molecular mimicry. PLoS ONE. 2013, 8 (6): e67377-10.1371/journal.pone.0067377.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  27. 27.

    Cotton JA, Lilley CJ, Jones LM, Kikuchi T, Reid AJ, Thorpe P, Tsai IJ, Beasley H, Blok VC, Cock PJA, Eves-van den Akker S, Holroyd N, Hunt M, Mantelin S, Naghra H, Pain A, Palomares-Rius JE, Zarowiecki M, Berriman M, Jones JT, Urwin PE: The genome and life-stage specific transcriptomes of Globodera pallida elucidate key aspects of plant parasitism by a cyst nematode. Genome Biol. 2014, 15: R43-10.1186/gb-2014-15-3-r43.

    PubMed Central  PubMed  Article  Google Scholar 

  28. 28.

    Smant G, Stokkermans JP, Yan Y, de Boer JM, Baum TJ, Wang X, Hussey RS, Gommers FJ, Henrissat B, Davis EL, Helder J, Schots A, Bakker J: Endogenous cellulases in animals: isolation of beta-1, 4-endoglucanase genes from two species of plant-parasitic cyst nematodes. Proc Natl Acad Sci U S A. 1998, 95: 4906-4911. 10.1073/pnas.95.9.4906.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  29. 29.

    Meyer WE, Schuster LN, Bartelmes G, Dieterich C, Sommer RJ: Horizontal gene transfer of microbial cellulases into nematode genomes is associated with functional assimilation and gene turnover. BMC Evolutionary Biol. 2011, 11: 13-10.1186/1471-2148-11-13.

    Article  Google Scholar 

  30. 30.

    Goellner M, Wang X, Davis EL: Endo-ß-1,4-glucanase expression in compatible plant-nematode interactions. Plant Cell. 2001, 13: 2241-2255.

    CAS  PubMed Central  PubMed  Google Scholar 

  31. 31.

    Sobczak M, Fudali S, Wieczorek K: Cell Wall Modifications Induced by Nematodes. Genomics and Molecular Genetics of Plant-Nematode Interactions. Edited by: Jones JT, Gheysen G, Fenoll C. 2011, Dordrecht: Springer, 395-422.

    Chapter  Google Scholar 

  32. 32.

    Elling AA, Davis EL, Hussey RS, Baum TJ: Active uptake of cyst nematode parasitism proteins into the plant cell nucleus. Int J Parasitol. 2007, 37: 1269-1279. 10.1016/j.ijpara.2007.03.012.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Jaouannet M, Perfus-Barbeoch L, Deleury E, Magliano M, Engler G, Vieira P, Danchin EGJ, Da Rocha M, Coquillard P, Abad P, Rosso MN: A root-knot nematode-secreted protein is injected into giant cells and targeted to the nuclei. New Phytol. 2012, 194: 924-931. 10.1111/j.1469-8137.2012.04164.x.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Lin B, Zhuo K, Wu P, Cui R, Zhang L-H, Liao J: A novel effector protein, MJ-NULG1a, targeted to giant cell nuclei plays a role in Meloidogyne javanica parasitism. Mol Plant-Microbe Interact. 2013, 26: 55-66. 10.1094/MPMI-05-12-0114-FI.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Wright KM, Cowan GH, Lukhovitskaya NI, Tilsner J, Roberts AG, Savenkov EI, Torrance L: The N-terminal domain of PMTV TGB1 movement protein is required for nuclear localisation, microtubule association and long-distance movement. Mol Plant-Microbe Interact. 2010, 23: 1486-1497. 10.1094/MPMI-05-10-0105.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Quentin M, Abad P, Favery B: Plant parasitic nematode effectors target host defense and nuclear functions to establish feeding cells. Front Plant Sci. 2013, 4: doi:10.3389/fpls.2013.00053

    Google Scholar 

  37. 37.

    Siamer S, Gaubert S, Boureau T, Brisset M-N, Barny MA: Mutational analysis of a predicted double β-propeller domain of the DspA/E effector of Erwinia amylovora. FEMS Microbiol Lett. 2013, 342: 54-61. 10.1111/1574-6968.12108.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Hu J, Baker A, Bartel B, Linka L, Mullen RT, Reumann S, Zolman BK: Plant peroxisomes: biogenesis and function. Plant Cell. 2012, 24: 2279-2303. 10.1105/tpc.112.096586.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  39. 39.

    Li C, Schilmiller AL, Liu G, Lee GI, Jayanty S, Sageman C, Vrebalov J, Giovannoni JJ, Yagi K, Kobayashi Y, Howe GA: Role of beta-oxidation in jasmonate biosynthesis and systemic wound signaling in tomato. Plant Cell. 2005, 17: 971-986. 10.1105/tpc.104.029108.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  40. 40.

    Mhamdi A, Noctor G, Baker A: Plant catalases: peroxisomal redox guardians. Arch Biochem Biophys. 2012, 525: 181-194. 10.1016/

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Waetzig GH, Sobczak M, Grundler FMW: Localization of hydrogen peroxide during the defence response of Arabidopsis thaliana against the plant parasitic nematode Heterodera glycines. Nematology. 1999, 1: 681-686. 10.1163/156854199508702.

    CAS  Article  Google Scholar 

  42. 42.

    Theodoulou FL, Berhardt K, Linka N, Baker A: Peroxisome membrane proteins: multiple trafficking routes and multiple functions?. Biochem J. 2013, 451: 345-352. 10.1042/BJ20130078.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Haas BJ, Kamoun S, Zody MC, Jiang RHY, Handsaker RE, Cano LM, Grabherr M, Kodira CD, Raffaele S, Torto-Alalibo T, Bozkurt TO, Ah-Fong AMV, Alvarado L, Anderson VL, Armstrong MR, Avrova A, Baxter L, Beynon J, Boevink PC, Bollmann SR, Bos JIB, Bulone V, Cai G, Cakir C, Carrington JC, Chawner M, Conti L, Costanzo S, Ewan R, Fahlgren N, et al: Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature. 2009, 461: 393-398. 10.1038/nature08358.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Raffaele S, Win J, Cano LM, Kamoun S: Analyses of genome architecture and gene expression reveal novel candidate virulence factors in the secretome of Phytophthora infestans. BMC Genomics. 2010, 11: 637-10.1186/1471-2164-11-637.

    PubMed Central  PubMed  Article  Google Scholar 

  45. 45.

    Saunders DGO, Win J, Kamoun S, Raffaele S: Two dimensional data binning for the analysis of genome architecture in filamentous plant pathogens and othe eukaryotes. Meth Mol Biol. 2014, in press

    Google Scholar 

  46. 46.

    Staiger D, Brown JWS: Alternative splicing at the intersection of biological timing, development and stress responses. Plant Cell. 2013, 25: 3640-3656. 10.1105/tpc.113.113803.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  47. 47.

    Lu SW, Tian DH, Borchardt-Wier HB, Wang XH: Alternative splicing: a novel mechanism of regulation identified in the chorismate mutase gene of the potato cyst nematode Globodera rostochiensis. Mol Biochem Parasitol. 2008, 162: 1-15. 10.1016/j.molbiopara.2008.06.002.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Atkinson HJ, Lilley CJ, Urwin PE: Strategies for transgenic nematode control in developed and developing world crops. Curr Opin Biotech. 2012, 23: 251-256. 10.1016/j.copbio.2011.09.004.

    CAS  PubMed  Article  Google Scholar 

  49. 49.

    Huang G, Allen R, Davis EL, Baum TJ, Hussey RS: Engineering broad root-knot resistance in transgenic plants by RNAi silencing of a conserved and essential root-knot nematode parasitism gene. Proc Natl Acad Sci U S A. 2006, 103: 14302-14306. 10.1073/pnas.0604698103.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  50. 50.

    Phillips MS, Trudgill DL: Variation of virulence, in terms of quantitative reproduction of Globodera pallida populations, from Europe and South America, in relation to resistance from Solanum vernei and S. tuberosum ssp andigena CPC 2802. Nematologica. 1998, 44: 409-423. 10.1163/005525998X00070.

    Article  Google Scholar 

  51. 51.

    Cock PJ, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, Friedberg I, Hamelryck T, Kauff F, Wilczynski B, de Hoon MJ: Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009, 25: 1422-1423. 10.1093/bioinformatics/btp163.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  52. 52.

    Pritchard L, White JA, Birch PRJ, Toth IK: GenomeDiagram: a python package for the visualization of large-scale genomic data. Bioinformatics. 2006, 22: 616-617. 10.1093/bioinformatics/btk021.

    CAS  PubMed  Article  Google Scholar 

  53. 53.

    Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, MacManes MD, Ott M, Orvis J, Pochet N, Strozzzi F, Weeks N, Westerman R, William T, Dewey CN, Henschel R, Leduc RD, Freidman N, Regev A: De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013, 8: 1494-1512. 10.1038/nprot.2013.084.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Merisov JP: Integrative genomics viewer. Nat Biotechnol. 2011, 29: 24-26. 10.1038/nbt.1754.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  55. 55.

    Wang X, Allen R, Ding X, Goellner M, Maier T, De Boer JM, Baum TJ, Hussey RS, Davis EL: Signal peptide-selection of cDNA cloned directly from the esophageal gland cells of the soybean cyst nematode Heterodera glycines. Mol Plant-Microbe Interact. 2001, 14: 536-544. 10.1094/MPMI.2001.14.4.536.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    De Boer JM, Mcdermott JP, Wang X, Maier T, Qui F, Hussey RS, Davis EL, Baum TJ: The use of DNA microarrays for the developmental expression analysis of cDNAs from the oesophageal gland cell region of Heterodera glycines. Mol Plant Pathol. 2002, 3: 261-270. 10.1046/j.1364-3703.2002.00122.x.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Popeijus H, Blok VC, Cardle L, Bakker E, Phillips MS, Helder J, Smant G, Jones JT: Analysis of genes expressed in second stage juveniles of the potato cyst nematodes Globodera rostochiensis and G. pallida using the expressed sequence tag approach. Nematology. 2000, 2: 567-574. 10.1163/156854100509358.

    CAS  Article  Google Scholar 

  58. 58.

    Huang GZ, Gao BL, Maier T, Allen R, Davis EL, Baum TJ, Hussey RS: A profile of putative parasitism genes expressed in the esophageal gland cells of the root-knot nematode Meloidogyne incognita. Mol Plant-Microbe Interact. 2003, 16: 376-381. 10.1094/MPMI.2003.16.5.376.

    CAS  PubMed  Article  Google Scholar 

  59. 59.

    Blanchard A, Fouville D, Esquibet M, Mugniery D, Grenier E: Sequence polymorphism of 2 pioneer genes expressed in phytoparasitic nematodes showing different host ranges. J Heredity. 2007, 98: 611-619. 10.1093/jhered/esm050.

    CAS  Article  Google Scholar 

  60. 60.

    Jones JT, Furlanetto C, Bakker E, Banks B, Blok VC, Chen Q, Prior A: Characterisation of a chorismate mutase from the potato cyst nematode Globodera pallida. Mol Plant Pathol. 2003, 4: 43-50. 10.1046/j.1364-3703.2003.00140.x.

    CAS  PubMed  Article  Google Scholar 

  61. 61.

    Altschul SF, Madden TL, Schañffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nuc Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    CAS  Article  Google Scholar 

  62. 62.

    Park BH, Karpinets TV, Syed MH, Leuze MR, Uberbacher EC: CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. Glycobiology. 2010, 20: 1574-1584. 10.1093/glycob/cwq106.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke ZX, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thank N, Yamashita RA, Zhang DC, Zhang NG, Zheng C, Bryant SH: CDD: A Conserved Domain Database for the functional annotation of proteins. Nuc Acids Res. 2011, 39: D225-D229. 10.1093/nar/gkq1189.

    CAS  Article  Google Scholar 

  64. 64.

    Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nuc Acid Res. 2005, 33: W116-W120. 10.1093/nar/gki442.

    CAS  Article  Google Scholar 

  65. 65.

    Dyrl Bendtsen J, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.

    Article  Google Scholar 

  66. 66.

    Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.

    CAS  PubMed  Article  Google Scholar 

  67. 67.

    Cock PJA, Grüning BA, Paszkiewicz KH, Pritchard L: Galaxy tools and workflows for sequence analysis with applications in molecular plant pathology. PeerJ. 2013, 1: e167-

    PubMed Central  PubMed  Article  Google Scholar 

  68. 68.

    Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): 106-10.1186/gb-2010-11-10-r106.

    Article  Google Scholar 

  69. 69.

    Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nuc Acids Res. 2011, 39: W29-W37. 10.1093/nar/gkr367.

    CAS  Article  Google Scholar 

  70. 70.

    Karimi M, Inze D, Depicker A: Gateway vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 2002, 7: 193-195. 10.1016/S1360-1385(02)02251-3.

    CAS  PubMed  Article  Google Scholar 

  71. 71.

    van der Fits L, Deakin EA, Hoge JHC, Memelink J: The ternary transformation system: constitutive virG on a compatible plasmid dramatically increases Agrobacterium-mediated plant transformation. Plant Mol Biol. 2000, 43: 495-502. 10.1023/A:1006440221718.

    CAS  PubMed  Article  Google Scholar 

  72. 72.

    Martin K, Kopperud K, Chakrabarty R, Banerjee R, Brooks R, Goodin MM: Transient expression in Nicotiana benthamiana fluorescent marker lines provides enhanced definition of protein localization, movement and interactions in planta. Plant J. 2009, 59: 150-162. 10.1111/j.1365-313X.2009.03850.x.

    CAS  PubMed  Article  Google Scholar 

  73. 73.

    McLellan H, Boevink PC, Armstrong MR, Pritchard L, Gomez S, Morales J, Beynon JL, Whisson SC, Birch PRJ: An RxLR effector from Phytophthora infestans prevents re-localisation of two plant NAC transcription factors from the endoplasmic reticulum to the nucleus. PLoS Pathog. 2013, doi:10.1371/journal.ppat.1003670

    Google Scholar 

Download references


This project was funded through BBSRC grants BB/F00334X/1 and BB/H000801/1, Scottish Government Contract Research Fund project CRF0926 and through a BBSRC CASE award for which the industrial partner was Mylnefield Research Services. The James Hutton Institute receives funding from the Scottish Government. EG was funded through ERASMUS MUNDUS programme 2008–102 (EUMAINE). The authors thank Prof. B. Braeckman and Dr F. Matthijssens for assistance provided to EG, James Cotton for help with analysis of the G. pallida genome and Dr L. Pritchard for assistance with the analysis of gene density. The support of Ailsa Smith and Anne Holt is gratefully acknowledged.

Author information



Corresponding authors

Correspondence to Peter E Urwin or John T Jones.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PT identified the effector complement of G. pallida, analysed genome distributions, assisted in cloning of the effectors and assisted in drafting the MS, SM was responsible for cloning of the effector sequences and assisted in the localisation experiments, PJAC assisted with all aspects of the bioinformatic analysis, VCB assisted in the genome and transcriptome analysis, MCC undertook most of the in situ hybridisation work, SEvdA assisted with in situ hybridisation and undertook the analysis of alternative splicing, EG assisted with cloning of effectors and analysis of subcellular localisations, CJL assisted in the design and coordination of the study, assisted in the transciptome analysis and assisted in drafting the MS, GS assisted in the bioinformatic identification of effectors, AJR undertook the transcriptome analysis, KMW assisted with subcellular localisation studies, PEU assisted in the design and coordination of the study and in drafting the MS, JTJ assisted in the design and coordination of the study, undertook some of the in situ hybridisation analysis and took primary responsibility for drafting the MS. All authors have critically read and revised the MS and all agree to be accountable for all aspects of the work.

Electronic supplementary material


Additional file 1: Table S3: G. pallida genes similar to effectors characterised from other nematodes (excluding cell wall degrading and modifying enzymes). Expression profiles of each gene as inferred from RNAseq analysis are indicated. (DOCX 23 KB)


Additional file 2: Table S4: Previously identified effectors from other cyst nematodes absent from the current G. pallida genome assembly. (DOCX 15 KB)


Additional file 3: Table S5: Globodera pallida secreted proteins up-regulated in J2 or early parasitic stages that may represent novel effector candidates. (DOCX 16 KB)


Additional file 4: Figure S1: Locations of exons (coloured arrows) in scaffolds 299 and 141 (B). Exons of genes encoding proteins similar to “H. avenae dorsal gland cell protein” are indicated in purple on panel A and exons of genes encoding proteins that contain a SPRY domain are indicated in yellow on panel B. Other predicted genes are indicated in blue or green (alternating genes). Direction of arrows indicates orientation of predicted open reading frames. Grey shading indicates unsequenced regions of scaffolds. (PNG 495 KB)


Additional file 5: Figure S2: Distribution of various effector groups across gene sparse and gene rich regions of the Globodera pallida genome. Heat map generated reflecting gene density and the distribution of three classes of effector candidates: the SPRYSECs, G. pallida orthologs of effectors from other PPN and “novel” effectors. (PNG 12 KB)

Figure S3A:

Additional file 6: Figure S3: Example of potential alternative splicing events in G. pallida effectors Figure S3A: Alternative splicing of GPLIN_000359000 within one life stage. Black bar indicates predicted sequence from genome with bars showing predicted coding region from gene model; introns are shown as lines. Red bars indicate de novo assembled transcripts from RNA extracted from parasitic nematodes 7dpi, red lines indicate gaps compared to genome sequence. Figure S3B: Alternative splicing of GPLIN_000243800 between life stages. Black bar indicates predicted sequence from genome with bars showing predicted coding region from gene model; introns are shown as lines. Red and blue bars/lines indicate de novo assembled transcripts from parasitic nematodes 7dpi and J2s respectively. (TIFF 2 MB)


Additional file 7: Table S1: Primers used to generate DNA fragments used for synthesis of probes for in situ hybridisation. (DOCX 14 KB)


Additional file 8: Table S2: Sequences of primers used for cloning full length effectors (without signal peptides). (DOCX 14 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Thorpe, P., Mantelin, S., Cock, P.J. et al. Genomic characterisation of the effector complement of the potato cyst nematode Globodera pallida. BMC Genomics 15, 923 (2014).

Download citation


  • Plant parasitic nematode
  • Effector
  • Genome
  • Transcriptome