A novel sugar beet cyst nematode effector 2D01 targets the Arabidopsis HAESA receptor-like kinase

Plant-parasitic cyst nematodes use a stylet to deliver effector proteins produced in oesophageal gland cells into root cells to cause disease in plants. These effectors are deployed to modulate plant defence responses and developmental programmes for the formation of a specialized feeding site called a syncytium. The Hg2D01 effector gene, coding for a novel 185-amino-acid secreted protein, was previously shown to be up-regulated in the dorsal gland of parasitic juveniles of the soybean cyst nematode Heterodera glycines, but its function has remained unknown. Genome analyses revealed that Hg2D01 belongs to a highly diversified effector gene family in the genomes of H. glycines and the sugar beet cyst nematode Heterodera schachtii. For functional studies using the model Arabidopsis thaliana-H. schachtii pathosystem, we cloned the orthologous Hs2D01 sequence from H. schachtii. We demonstrate that Hs2D01 is a cytoplasmic effector that interacts with the intracellular kinase domain of HAESA (HAE), a cell surface-associated leucine-rich repeat (LRR) receptor-like kinase (RLK) involved in signalling the activation of cell wall-remodelling enzymes important for cell separation during abscission and lateral root emergence. Furthermore, we show that AtHAE is expressed in the syncytium and, therefore, could serve as a viable host target for Hs2D01. Infective juveniles effectively penetrated the roots of HAE and HAESA-LIKE2 (HSL2) double mutant plants; however, fewer nematodes developed on the roots, consistent with a role for this receptor family in nematode infection. Taken together, our results suggest that the Hs2D01-AtHAE interaction may play an important role in sugar beet cyst nematode parasitism.

The genome and lifestage-specific transcriptomes of a plant-parasitic nematode and its host reveal susceptibility genes involved in trans-kingdom synthesis of vitamin B5

Plant-parasitic nematodes are a major threat to crop production in all agricultural systems. The scarcity of classical resistance genes highlights a pressing need to find new ways to develop nematode-resistant germplasm. Here, we sequence and assemble a high-quality phased genome of the model cyst nematode Heterodera schachtii to provide a platform for the first system-wide dual analysis of host and parasite gene expression over time, covering all major parasitism stages. Analysis of the hologenome of the plant-nematode infection site identified metabolic pathways that were incomplete in the parasite but complemented by the host. Using a combination of bioinformatic, genetic, and biochemical approaches, we show that a highly atypical completion of vitamin B5 biosynthesis by the parasitic animal, putatively enabled by a horizontal gene transfer from a bacterium, is required for full pathogenicity. Knockout of either plant-encoded or now nematode-encoded steps in the pathway significantly reduces parasitic success. Our experiments establish a reference for cyst nematodes, further our understanding of the evolution of plant-parasitism by nematodes, and show that congruent differential expression of metabolic pathways in the infection hologenome represents a new way to find nematode susceptibility genes. The approach identifies genome-editing-amenable targets for future development of nematode-resistant crops.

miR778 mediates gene expression, histone modification, and DNA methylation during cyst nematode parasitism

Despite the known critical regulatory functions of microRNAs, histone modifications, and DNA methylation in reprograming plant epigenomes in response to pathogen infection, the molecular mechanisms underlying the tight coordination of these components remain poorly understood. Here, we show how Arabidopsis (Arabidopsis thaliana) miR778 coordinately modulates the root transcriptome, histone methylation, and DNA methylation via post-transcriptional regulation of the H3K9 methyltransferases SU(var)3-9 homolog 5 (SUVH5) and SUVH6 upon infection by the beet cyst nematode Heterodera schachtii. miR778 post-transcriptionally silences SUVH5 and SUVH6 upon nematode infection. Manipulation of the expression of miR778 and its two target genes significantly altered plant susceptibility to H. schachtii. RNA-seq analysis revealed a key role of SUVH5 and SUVH6 in reprograming the transcriptome of Arabidopsis roots upon H. schachtii infection. In addition, chromatin immunoprecipitation (ChIP)-seq analysis established SUVH5 and SUVH6 as the main enzymes mediating H3K9me2 deposition in Arabidopsis roots in response to nematode infection. ChIP-seq analysis also showed that these methyltransferases possess distinct DNA binding preferences in that they are targeting transposable elements under noninfected conditions and protein-coding genes in infected plants. Further analyses indicated that H3K9me2 deposition directed by SUVH5 and SUVH6 contributes to gene expression changes both in roots and in nematode feeding sites and preferentially associates with CG DNA methylation. Together, our results uncovered multi-layered epigenetic regulatory mechanisms coordinated by miR778 during Arabidopsis-H. schachtii interactions.

Targeted transcriptomics reveals signatures of large-scale independent origins and concerted regulation of effector genes in Radopholus similis

The burrowing nematode, Radopholus similis, is an economically important plant-parasitic nematode that inflicts damage and yield loss to a wide range of crops. This migratory endoparasite is widely distributed in warmer regions and causes extensive destruction to the root systems of important food crops (e.g., citrus, banana). Despite the economic importance of this nematode, little is known about the repertoire of effectors owned by this species. Here we combined spatially and temporally resolved next-generation sequencing datasets of R. similis to select a list of candidates for the identification of effector genes for this species. We confirmed spatial expression of transcripts of 30 new candidate effectors within the esophageal glands of R. similis by in situ hybridization, revealing a large number of pioneer genes specific to this nematode. We identify a gland promoter motif specifically associated with the subventral glands (named Rs-SUG box), a putative hallmark of spatial and concerted regulation of these effectors. Nematode transcriptome analyses confirmed the expression of these effectors during the interaction with the host, with a large number of pioneer genes being especially abundant. Our data revealed that R. similis holds a diverse and emergent repertoire of effectors, which has been shaped by various evolutionary events, including neofunctionalization, horizontal gene transfer, and possibly by de novo gene birth. In addition, we also report the first GH62 gene so far discovered for any metazoan and putatively acquired by lateral gene transfer from a bacterial donor. Considering the economic damage caused by R. similis, this information provides valuable data to elucidate the mode of parasitism of this nematode.

Esophageal Gland RNA-Seq Resource of a Virulent and Avirulent Population of the Soybean Cyst Nematode Heterodera glycines

The soybean cyst nematode Heterodera glycines is the most economically devastating pathogen of soybean in the United States and threatens to become even more damaging through the selection of virulent nematode populations in the field that can overcome natural resistance mechanisms in soybean cultivars. This pathogen, therefore, demands intense transcriptomic/genomic research inquiries into the biology of its parasitic mechanisms. H. glycines delivers effector proteins that are produced in specialized gland cells into the soybean root to enable infection. The study of effector proteins, thus, is particularly promising when exploring novel management options against this pathogen. Here, we announce the availability of a gland cell-specific RNA-seq resource. These data represent an expression snapshot of gland cell activity during early soybean infection of a virulent and an avirulent H. glycines population, providing a unique and highly valuable resource for scientists examining effector biology and nematode virulence.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Toward genetic modification of plant-parasitic nematodes: delivery of macromolecules to adults and expression of exogenous mRNA in second stage juveniles

Plant-parasitic nematodes are a continuing threat to food security, causing an estimated 100 billion USD in crop losses each year. The most problematic are the obligate sedentary endoparasites (primarily root knot nematodes and cyst nematodes). Progress in understanding their biology is held back by a lack of tools for functional genetics: forward genetics is largely restricted to studies of natural variation in populations and reverse genetics is entirely reliant on RNA interference. There is an expectation that the development of functional genetic tools would accelerate the progress of research on plant-parasitic nematodes, and hence the development of novel control solutions. Here, we develop some of the foundational biology required to deliver a functional genetic tool kit in plant-parasitic nematodes. We characterize the gonads of male Heterodera schachtii and Meloidogyne hapla in the context of spermatogenesis. We test and optimize various methods for the delivery, expression, and/or detection of exogenous nucleic acids in plant-parasitic nematodes. We demonstrate that delivery of macromolecules to cyst and root knot nematode male germlines is difficult, but possible. Similarly, we demonstrate the delivery of oligonucleotides to root knot nematode gametes. Finally, we develop a transient expression system in plant-parasitic nematodes by demonstrating the delivery and expression of exogenous mRNA encoding various reporter genes throughout the body of H. schachtii juveniles using lipofectamine-based transfection. We anticipate these developments to be independently useful, will expedite the development of genetic modification tools for plant-parasitic nematodes, and ultimately catalyze research on a group of nematodes that threaten global food security.

Phytonematode peptide effectors exploit a host post-translational trafficking mechanism to the ER using a novel translocation signal

Cyst nematodes induce a multicellular feeding site within roots called a syncytium. It remains unknown how root cells are primed for incorporation into the developing syncytium. Furthermore, it is unclear how CLAVATA3/EMBRYO SURROUNDING REGION (CLE) peptide effectors secreted into the cytoplasm of the initial feeding cell could have an effect on plant cells so distant from where the nematode is feeding as the syncytium expands. Here we describe a novel translocation signal within nematode CLE effectors that is recognized by plant cell secretory machinery to redirect these peptides from the cytoplasm to the apoplast of plant cells. We show that the translocation signal is functionally conserved across CLE effectors identified in nematode species spanning three genera and multiple plant species, operative across plant cell types, and can traffic other unrelated small peptides from the cytoplasm to the apoplast of host cells via a previously unknown post-translational mechanism of endoplasmic reticulum (ER) translocation. Our results uncover a mechanism of effector trafficking that is unprecedented in any plant pathogen to date, andthey illustrate how phytonematodes can deliver effector proteins into host cells and then hijack plant cellular processes for their export back out of the cell to function as external signaling molecules to distant cells.

A new esophageal gland transcriptome reveals signatures of large scale de novo effector birth in the root lesion nematode Pratylenchus penetrans

BACKGROUND

The root lesion nematode Pratylenchus penetrans is a migratory plant-parasitic nematode responsible for economically important losses in a wide number of crops. Despite the importance of P. penetrans, the molecular mechanisms employed by this nematode to promote virulence remain largely unknown.

RESULTS

Here we generated a new and comprehensive esophageal glands-specific transcriptome library for P. penetrans. In-depth analysis of this transcriptome enabled a robust identification of a catalogue of 30 new candidate effector genes, which were experimentally validated in the esophageal glands by in situ hybridization. We further validated the expression of a multifaceted network of candidate effectors during the interaction with different plants. To advance our understanding of the "effectorome" of P. penetrans, we adopted a phylogenetic approach and compared the expanded effector repertoire of P. penetrans to the genome/transcriptome of other nematode species with similar or contrasting parasitism strategies. Our data allowed us to infer plausible evolutionary histories that shaped the effector repertoire of P. penetrans, as well as other close and distant plant-parasitic nematodes. Two remarkable trends were apparent: 1) large scale effector birth in the Pratylenchidae in general and P. penetrans in particular, and 2) large scale effector death in sedentary (endo) plant-parasitic nematodes.

CONCLUSIONS

Our study doubles the number of validated Pratylenchus penetrans effectors reported in the literature. The dramatic effector gene gain in P. penetrans could be related to the remarkable ability of this nematode to parasitize a large number of plants. Our data provide valuable insights into nematode parasitism and contribute towards basic understating of the adaptation of P. penetrans and other root lesion nematodes to specific host plants.

Large tandem duplications affect gene expression, 3D organization, and plant-pathogen response

Rapid plant genome evolution is crucial to adapt to environmental changes. Chromosomal rearrangements and gene copy number variation (CNV) are two important tools for genome evolution and sources for the creation of new genes. However, their emergence takes many generations. In this study, we show that in Arabidopsis thaliana, a significant loss of ribosomal RNA (rRNA) genes with a past history of a mutation for the chromatin assembly factor 1 (CAF1) complex causes rapid changes in the genome structure. Using long-read sequencing and microscopic approaches, we have identified up to 15 independent large tandem duplications in direct orientation (TDDOs) ranging from 60 kb to 1.44 Mb. Our data suggest that these TDDOs appeared within a few generations, leading to the duplication of hundreds of genes. By subsequently focusing on a line only containing 20% of rRNA gene copies (20rDNA line), we investigated the impact of TDDOs on 3D genome organization, gene expression, and cytosine methylation. We found that duplicated genes often accumulate more transcripts. Among them, several are involved in plant–pathogen response, which could explain why the 20rDNA line is hyper-resistant to both bacterial and nematode infections. Finally, we show that the TDDOs create gene fusions and/or truncations and discuss their potential implications for the evolution of plant genomes.

Screening soybean cyst nematode effectors for their ability to suppress plant immunity

The soybean cyst nematode (SCN), Heterodera glycines, is one of the most destructive pathogens of soybeans. SCN is an obligate and sedentary parasite that transforms host plant root cells into an elaborate permanent feeding site, a syncytium. Formation and maintenance of a viable syncytium is an absolute requirement for nematode growth and reproduction. In turn, sensing pathogen attack, plants activate defence responses and may trigger programmed cell death at the sites of infection. For successful parasitism, H. glycines must suppress these host defence responses to establish and maintain viable syncytia. Similar to other pathogens, H. glycines engages in these molecular interactions with its host via effector proteins. The goal of this study was to conduct a comprehensive screen to identify H. glycines effectors that interfere with plant immune responses. We used Nicotiana benthamiana plants infected by Pseudomonas syringae and Pseudomonas fluorescens strains. Using these pathosystems, we screened 51 H. glycines effectors to identify candidates that could inhibit effector-triggered immunity (ETI) and/or pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI). We identified three effectors as ETI suppressors and seven effectors as PTI suppressors. We also assessed expression modulation of plant immune marker genes as a function of these suppressors.

Targeted suppression of soybean BAG6-induced cell death in yeast by soybean cyst nematode effectors

While numerous effectors that suppress plant immunity have been identified from bacteria, fungi, and oomycete pathogens, relatively little is known for nematode effectors. Several dozen effectors have been reported from the soybean cyst nematode (SCN). Previous studies suggest that a hypersensitive response-like programmed cell death is triggered at nematode feeding sites in soybean during an incompatible interaction. However, virulent SCN populations overcome this incompatibility using unknown mechanisms. A soybean BAG6 (Bcl-2 associated anthanogene 6) gene previously reported by us to be highly up-regulated in degenerating feeding sites induced by SCN in a resistant soybean line was attenuated in response to a virulent SCN population. We show that GmBAG6-1 induces cell death in yeast like its Arabidopsis homolog AtBAG6 and also in soybean. This led us to hypothesize that virulent SCN may target GmBAG6-1 as part of their strategy to overcome soybean defence responses during infection. Thus, we used a yeast viability assay to screen SCN effector candidates for their ability to specifically suppress GmBAG6-1-induced cell death. We identified several effectors that strongly suppressed cell death mediated by GmBAG6-1. Two effectors identified as suppressors showed direct interaction with GmBAG6-1 in yeast, suggesting that one mechanism of cell death suppression may occur through an interaction with this host protein.

A role for Arabidopsis growth-regulating factors 1 and 3 in growth-stress antagonism

Growth-regulating factors (GRFs) belong to a small family of transcription factors that are highly conserved in plants. GRFs regulate many developmental processes and plant responses to biotic and abiotic stimuli. Despite the importance of GRFs, a detailed mechanistic understanding of their regulatory functions is still lacking. In this study, we used ChIP sequencing (ChIP-seq) to identify genome-wide binding sites of Arabidopsis GRF1 and GRF3, and correspondingly their direct downstream target genes. RNA-sequencing (RNA-seq) analysis revealed that GRF1 and GRF3 regulate the expression of a significant number of the identified direct targets. The target genes unveiled broad regulatory functions of GRF1 and GRF3 in plant growth and development, phytohormone biosynthesis and signaling, and the cell cycle. Our analyses also revealed that clock core genes and genes with stress- and defense-related functions are most predominant among the GRF1- and GRF3-bound targets, providing insights into a possible role for these transcription factors in mediating growth-defense antagonism and integrating environmental stimuli into developmental programs. Additionally, GRF1 and GRF3 target molecular nodes of growth-defense antagonism and modulate the levels of defense- and development-related hormones in opposite directions. Taken together, our results point to GRF1 and GRF3 as potential key determinants of plant fitness under stress conditions.

Re-targeting of a plant defense protease by a cyst nematode effector

Plants mount defense responses during pathogen attacks, and robust host defense suppression by pathogen effector proteins is essential for infection success. 4E02 is an effector of the sugar beet cyst nematode Heterodera schachtii. Arabidopsis thaliana lines expressing the effector-coding sequence showed altered expression levels of defense response genes, as well as higher susceptibility to both the biotroph H. schachtii and the necrotroph Botrytis cinerea, indicating a potential suppression of defenses by 4E02. Yeast two-hybrid analyses showed that 4E02 targets A. thaliana vacuolar papain-like cysteine protease (PLCP) 'Responsive to Dehydration 21A' (RD21A), which has been shown to function in the plant defense response. Activity-based protein profiling analyses documented that the in planta presence of 4E02 does not impede enzymatic activity of RD21A. Instead, 4E02 mediates a re-localization of this protease from the vacuole to the nucleus and cytoplasm, which is likely to prevent the protease from performing its defense function and at the same time, brings it in contact with novel substrates. Yeast two-hybrid analyses showed that RD21A interacts with multiple host proteins including enzymes involved in defense responses as well as carbohydrate metabolism. In support of a role in carbohydrate metabolism of RD21A after its effector-mediated re-localization, we observed cell wall compositional changes in 4E02 expressing A. thaliana lines. Collectively, our study shows that 4E02 removes RD21A from its defense-inducing pathway and repurposes this enzyme by targeting the active protease to different cell compartments.

Homeostasis in the soybean miRNA396-GRF network is essential for productive soybean cyst nematode infections

Heterodera glycines, the soybean cyst nematode, penetrates soybean roots and migrates to the vascular cylinder where it forms a feeding site called the syncytium. MiRNA396 (miR396) targets growth-regulating factor (GRF) genes, and the miR396-GRF1/3 module is a master regulator of syncytium development in model cyst nematode H. schachtii infection of Arabidopsis. Here, we investigated whether this regulatory system operates similarly in soybean roots and is likewise important for H. glycines infection. We found that a network involving nine MIR396 and 23 GRF genes is important for normal development of soybean roots and that GRF function is specified in the root apical meristem by miR396. All MIR396 genes are down-regulated in the syncytium during its formation phase while, specifically, 11 different GRF genes are up-regulated. The switch to the syncytium maintenance phase coincides with up-regulation of MIR396 and down-regulation of the 11 GRF genes specifically via post-transcriptional regulation by miR396. Furthermore, interference with the miR396-GRF6/8-13/15-17/19 regulatory network, through either overexpression or knockdown experiments, does not affect the number of H. glycines juveniles that enter the vascular cylinder to initiate syncytia, but specifically inhibits efficient H. glycines development to adult females. Therefore, homeostasis in the miR396-GRF6/8-13/15-17/19 regulatory network is essential for productive H. glycines infections.

The genome of the soybean cyst nematode (Heterodera glycines) reveals complex patterns of duplications involved in the evolution of parasitism genes

BACKGROUND

Heterodera glycines, commonly referred to as the soybean cyst nematode (SCN), is an obligatory and sedentary plant parasite that causes over a billion-dollar yield loss to soybean production annually. Although there are genetic determinants that render soybean plants resistant to certain nematode genotypes, resistant soybean cultivars are increasingly ineffective because their multi-year usage has selected for virulent H. glycines populations. The parasitic success of H. glycines relies on the comprehensive re-engineering of an infection site into a syncytium, as well as the long-term suppression of host defense to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms of effector acquisition, diversification, and selection need to be understood before effective control strategies can be developed, but the lack of an annotated genome has been a major roadblock.

RESULTS

Here, we use PacBio long-read technology to assemble a H. glycines genome of 738 contigs into 123 Mb with annotations for 29,769 genes. The genome contains significant numbers of repeats (34%), tandem duplicates (18.7 Mb), and horizontal gene transfer events (151 genes). A large number of putative effectors (431 genes) were identified in the genome, many of which were found in transposons.

CONCLUSIONS

This advance provides a glimpse into the host and parasite interplay by revealing a diversity of mechanisms that give rise to virulence genes in the soybean cyst nematode, including: tandem duplications containing over a fifth of the total gene count, virulence genes hitchhiking in transposons, and 107 horizontal gene transfers not reported in other plant parasitic nematodes thus far. Through extensive characterization of the H. glycines genome, we provide new insights into H. glycines biology and shed light onto the mystery underlying complex host-parasite interactions. This genome sequence is an important prerequisite to enable work towards generating new resistance or control measures against H. glycines.

Heterodera glycines utilizes promiscuous spliced leaders and demonstrates a unique preference for a species-specific spliced leader over C. elegans SL1

Spliced leader trans-splicing (SLTS) plays a part in the maturation of pre-mRNAs in select species across multiple phyla but is particularly prevalent in Nematoda. The role of spliced leaders (SL) within the cell is unclear and an accurate assessment of SL occurrence within an organism is possible only after extensive sequencing data are available, which is not currently the case for many nematode species. SL discovery is further complicated by an absence of SL sequences from high-throughput sequencing results due to incomplete sequencing of the 5'-ends of transcripts during RNA-seq library preparation, known as 5'-bias. Existing datasets and novel methodology were used to identify both conserved SLs and unique hypervariable SLs within Heterodera glycines, the soybean cyst nematode. In H. glycines, twenty-one distinct SL sequences were found on 2,532 unique H. glycines transcripts. The SL sequences identified on the H. glycines transcripts demonstrated a high level of promiscuity, meaning that some transcripts produced as many as nine different individual SL-transcript combinations. Most uniquely, transcriptome analysis revealed that H. glycines is the first nematode to demonstrate a higher SL trans-splicing rate using a species-specific SL over well-conserved Caenorhabditis elegans SL-like sequences.

SCNBase: a genomics portal for the soybean cyst nematode (Heterodera glycines)

Soybean is an important worldwide crop, and farmers continue to experience significant yield loss due to the soybean cyst nematode (SCN), Heterodera glycines. This soil-borne roundworm parasite is rated the most important pathogen problem in soybean production. The infective nematodes enter into complex interactions with their host plant by inducing the development of specialized plant feeding cells that provide the parasites with nourishment. Addressing the SCN problem will require the development of genomic resources and a global collaboration of scientists to analyze and use these resources. SCNBase.org was designed as a collaborative hub for the SCN genome. All data and analyses are downloadable and can be analyzed with three integrated genomic tools: JBrowse, Feature Search and BLAST. At the time of this writing, a number of genomic and transcriptomic data sets are already available, with 43 JBrowse tracks and 21 category pages describing SCN genomic analyses on gene predictions, transcriptome and read alignments, effector-like genes, expansion and contraction of genomic repeats, orthology and synteny with related nematode species, Single Nucleotide Polymorphism (SNPs) from 15 SCN populations and novel splice sites. Standard functional gene annotations were supplemented with orthologous gene annotations using a comparison to nine related plant-parasitic nematodes, thereby enabling functional annotations for 85% of genes. These annotations led to a greater grasp on the SCN effectorome, which include over 3324 putative effector genes. By designing SCNBase as a hub, future research findings and genomic resources can easily be uploaded and made available for use by others with minimal needs for further curation. By providing these resources to nematode research community, scientists will be empowered to develop novel, more effective SCN management tools.

An Effector from the Cyst Nematode Heterodera schachtii Derepresses Host rRNA Genes by Altering Histone Acetylation

Cyst nematodes are plant-pathogenic animals that secrete effector proteins into plant root cells to alter host gene expression and reprogram these cells to form specialized feeding sites, known as syncytia. The molecular mechanisms of these effectors are mostly unknown. We determined that the sugar beet cyst nematode (Heterodera schachtii) 32E03 effector protein strongly inhibits the activities of Arabidopsis thaliana histone deacetylases including the HDT1 enzyme, which has a known function in the regulation of rRNA gene expression through chromatin modifications. We determined that plants expressing the 32E03 coding sequence exhibited increased acetylation of histone H3 along the rDNA chromatin. At low 32E03 expression levels, these chromatin changes triggered the derepression of a subset of rRNA genes, which were conducive to H. schachtii parasitism. By contrast, high levels of 32E03 caused profound bidirectional transcription along the rDNA, which triggered rDNA-specific small RNA production leading to RNA-directed DNA methylation and silencing of rDNA, which inhibited nematode development. Our data show that the 32E03 effector alters plant rRNA gene expression by modulating rDNA chromatin in a dose-dependent manner. Thus, the 32E03 effector epigenetically regulates plant gene expression to promote cyst nematode parasitism.

The plant-parasitic cyst nematode effector GLAND4 is a DNA-binding protein

Cyst nematodes are plant pathogens that infect a wide range of economically important crops. One parasitic mechanism employed by cyst nematodes is the production and in planta delivery of effector proteins to modify plant cells and suppress defences to favour parasitism. This study focuses on GLAND4, an effector of Heterodera glycines and H. schachtii, the soybean and sugar beet cyst nematodes, respectively. We show that GLAND4 is recognized by the plant cellular machinery and is transported to the plant nucleus, an organelle for which little is known about plant nematode effector functions. We show that GLAND4 has DNA-binding ability and represses reporter gene expression in a plant transcriptional assay. One DNA fragment that binds to GLAND4 is localized in an Arabidopsis chromosomal region associated with the promoters of two lipid transfer protein genes (LTP). These LTPs have known defence functions and are down-regulated in the nematode feeding site. When expressed in Arabidopsis, the presence of GLAND4 causes the down-regulation of the two LTP genes in question, which is also associated with increased susceptibility to the plant-pathogenic bacterium Pseudomonas syringae. Furthermore, overexpression of one of the LTP genes reduces plant susceptibility to H. schachtii and P. syringae, confirming that LTP repression probably suppresses plant defences. This study makes GLAND4 one of a small subset of characterized plant nematode nuclear effectors and identifies GLAND4 as the first DNA-binding, plant-parasitic nematode effector.

The novel cyst nematode effector protein 30D08 targets host nuclear functions to alter gene expression in feeding sites

Cyst nematodes deliver effector proteins into host cells to manipulate cellular processes and establish a metabolically hyperactive feeding site. The novel 30D08 effector protein is produced in the dorsal gland of parasitic juveniles, but its function has remained unknown. We demonstrate that expression of 30D08 contributes to nematode parasitism, the protein is packaged into secretory granules and it is targeted to the plant nucleus where it interacts with SMU2 (homolog of suppressor of mec-8 and unc-52 2), an auxiliary spliceosomal protein. We show that SMU2 is expressed in feeding sites and an smu2 mutant is less susceptible to nematode infection. In Arabidopsis expressing 30D08 under the SMU2 promoter, several genes were found to be alternatively spliced and the most abundant functional classes represented among differentially expressed genes were involved in RNA processing, transcription and binding, as well as in development, and hormone and secondary metabolism, representing key cellular processes known to be important for feeding site formation. In conclusion, we demonstrated that the 30D08 effector is secreted from the nematode and targeted to the plant nucleus where its interaction with a host auxiliary spliceosomal protein may alter the pre-mRNA splicing and expression of a subset of genes important for feeding site formation.

Identification of candidate effector genes of Pratylenchus penetrans

Pratylenchus penetrans is one of the most important species of root lesion nematodes (RLNs) because of its detrimental and economic impact in a wide range of crops. Similar to other plant-parasitic nematodes (PPNs), P. penetrans harbours a significant number of secreted proteins that play key roles during parasitism. Here, we combined spatially and temporally resolved next-generation sequencing datasets of P. penetrans to select a list of candidate genes aimed at the identification of a panel of effector genes for this species. We determined the spatial expression of transcripts of 22 candidate effectors within the oesophageal glands of P. penetrans by in situ hybridization. These comprised homologues of known effectors of other PPNs with diverse putative functions, as well as novel pioneer effectors specific to RLNs. It is noteworthy that five of the pioneer effectors encode extremely proline-rich proteins. We then combined in situ localization of effectors with available genomic data to identify a non-coding motif enriched in promoter regions of a subset of P. penetrans effectors, and thus a putative hallmark of spatial expression. Expression profiling analyses of a subset of candidate effectors confirmed their expression during plant infection. Our current results provide the most comprehensive panel of effectors found for RLNs. Considering the damage caused by P. penetrans, this information provides valuable data to elucidate the mode of parasitism of this nematode and offers useful suggestions regarding the potential use of P. penetrans-specific target effector genes to control this important pathogen.

Novel global effector mining from the transcriptome of early life stages of the soybean cyst nematode Heterodera glycines

Soybean cyst nematode (SCN) Heterodera glycines is an obligate parasite that relies on the secretion of effector proteins to manipulate host cellular processes that favor the formation of a feeding site within host roots to ensure its survival. The sequence complexity and co-evolutionary forces acting upon these effectors remain unknown. Here we generated a de novo transcriptome assembly representing the early life stages of SCN in both a compatible and an incompatible host interaction to facilitate global effector mining efforts in the absence of an available annotated SCN genome. We then employed a dual effector prediction strategy coupling a newly developed nematode effector prediction tool, N-Preffector, with a traditional secreted protein prediction pipeline to uncover a suite of novel effector candidates. Our analysis distinguished between effectors that co-evolve with the host genotype and those conserved by the pathogen to maintain a core function in parasitism and demonstrated that alternative splicing is one mechanism used to diversify the effector pool. In addition, we confirmed the presence of viral and microbial inhabitants with molecular sequence information. This transcriptome represents the most comprehensive whole-nematode sequence currently available for SCN and can be used as a tool for annotation of expected genome assemblies.

Suppression or Activation of Immune Responses by Predicted Secreted Proteins of the Soybean Rust Pathogen Phakopsora pachyrhizi

Rust fungi, such as the soybean rust pathogen Phakopsora pachyrhizi, are major threats to crop production. They form specialized haustoria that are hyphal structures intimately associated with host-plant cell membranes. These haustoria have roles in acquiring nutrients and secreting effector proteins that manipulate host immune systems. Functional characterization of effector proteins of rust fungi is important for understanding mechanisms that underlie their virulence and pathogenicity. Hundreds of candidate effector proteins have been predicted for rust pathogens, but it is not clear how to prioritize these effector candidates for further characterization. There is a need for high-throughput approaches for screening effector candidates to obtain experimental evidence for effector-like functions, such as the manipulation of host immune systems. We have focused on identifying effector candidates with immune-related functions in the soybean rust fungus P. pachyrhizi. To facilitate the screening of many P. pachyrhizi effector candidates (named PpECs), we used heterologous expression systems, including the bacterial type III secretion system, Agrobacterium infiltration, a plant virus, and a yeast strain, to establish an experimental pipeline for identifying PpECs with immune-related functions and establishing their subcellular localizations. Several PpECs were identified that could suppress or activate immune responses in nonhost Nicotiana benthamiana, N. tabacum, Arabidopsis, tomato, or pepper plants.

Cooperative Regulatory Functions of miR858 and MYB83 during Cyst Nematode Parasitism

MicroRNAs (miRNAs) recently have been established as key regulators of transcriptome reprogramming that define cell function and identity. Nevertheless, the molecular functions of the greatest number of miRNA genes remain to be determined. Here, we report cooperative regulatory functions of miR858 and its MYB83 transcription factor target gene in transcriptome reprogramming during Heterodera cyst nematode parasitism of Arabidopsis (Arabidopsis thaliana). Gene expression analyses and promoter-GUS fusion assays documented a role of miR858 in posttranscriptional regulation of MYB83 in the Heterodera schachtii-induced feeding sites, the syncytia. Constitutive overexpression of miR858 interfered with H. schachtii parasitism of Arabidopsis, leading to reduced susceptibility, while reduced miR858 abundance enhanced plant susceptibility. Similarly, MYB83 expression increases were conducive to nematode infection because overexpression of a noncleavable coding sequence of MYB83 significantly increased plant susceptibility, whereas a myb83 mutation rendered the plants less susceptible. In addition, RNA-seq analysis revealed that genes involved in hormone signaling pathways, defense response, glucosinolate biosynthesis, cell wall modification, sugar transport, and transcriptional control are the key etiological factors by which MYB83 facilitates nematode parasitism of Arabidopsis. Furthermore, we discovered that miR858-mediated silencing of MYB83 is tightly regulated through a feedback loop that might contribute to fine-tuning the expression of more than a thousand of MYB83-regulated genes in the H. schachtii-induced syncytium. Together, our results suggest a role of the miR858-MYB83 regulatory system in finely balancing gene expression patterns during H. schachtii parasitism of Arabidopsis to ensure optimal cellular function.

Arabidopsis miR827 mediates post-transcriptional gene silencing of its ubiquitin E3 ligase target gene in the syncytium of the cyst nematode Heterodera schachtii to enhance susceptibility

MicroRNAs (miRNAs) are a major class of small non-coding RNAs with emerging functions in biotic and abiotic interactions. Here, we report on a new functional role of Arabidopsis miR827 and its NITROGEN LIMITATION ADAPTATION (NLA) target gene in mediating plant susceptibility to the beet cyst nematode Heterodera schachtii. Cyst nematodes are sedentary endoparasites that induce the formation of multinucleated feeding structures termed syncytia in the roots of host plants. Using promoter:GUS fusion assays we established that miR827 was activated in the initial feeding cells and this activation was maintained in the syncytium during all sedentary stages of nematode development. Meanwhile, the NLA target gene, which encodes an ubiquitin E3 ligase enzyme, was post-transcriptionally silenced in the syncytium to permanently suppress its activity during all nematode parasitic stages. Overexpression of miR827 in Arabidopsis resulted in hyper-susceptibility to H. schachtii. In contrast, inactivation of miR827 activity through target mimicry or by overexpression a miR827-resistant cDNA of NLA produced the opposite phenotype of reduced plant susceptibility to H. schachtii. Gene expression analysis of several pathogenesis-related genes together with Agrobacterium-mediated transient expression in Nicotiana benthamiana provided strong evidence that miR827-mediated downregulation of NLA to suppress basal defense pathways. In addition, using yeast two-hybrid screens we identified several candidates of NLA-interacting proteins that are involved in a wide range of biological processes and molecular functions, including three pathogenesis-related proteins. Taken together, we conclude that nematode-activated miR827 in the syncytium is necessary to suppress immune responses in order to establish infection and cause disease.

A Plasmodium-like virulence effector of the soybean cyst nematode suppresses plant innate immunity

Heterodera glycines, the soybean cyst nematode, delivers effector proteins into soybean roots to initiate and maintain an obligate parasitic relationship. HgGLAND18 encodes a candidate H. glycines effector and is expressed throughout the infection process. We used a combination of molecular, genetic, bioinformatic and phylogenetic analyses to determine the role of HgGLAND18 during H. glycines infection. HgGLAND18 is necessary for pathogenicity in compatible interactions with soybean. The encoded effector strongly suppresses both basal and hypersensitive cell death innate immune responses, and immunosuppression requires the presence and coordination between multiple protein domains. The N-terminal domain in HgGLAND18 contains unique sequence similarity to domains of an immunosuppressive effector of Plasmodium spp., the malaria parasites. The Plasmodium effector domains functionally complement the loss of the N-terminal domain from HgGLAND18. In-depth sequence searches and phylogenetic analyses demonstrate convergent evolution between effectors from divergent parasites of plants and animals as the cause of sequence and functional similarity.

A Small Cysteine-Rich Protein from the Asian Soybean Rust Fungus, Phakopsora pachyrhizi, Suppresses Plant Immunity

The Asian soybean rust fungus, Phakopsora pachyrhizi, is an obligate biotrophic pathogen causing severe soybean disease epidemics. Molecular mechanisms by which P. pachyrhizi and other rust fungi interact with their host plants are poorly understood. The genomes of all rust fungi encode many small, secreted cysteine-rich proteins (SSCRP). While these proteins are thought to function within the host, their roles are completely unknown. Here, we present the characterization of P. pachyrhizi effector candidate 23 (PpEC23), a SSCRP that we show to suppress plant immunity. Furthermore, we show that PpEC23 interacts with soybean transcription factor GmSPL12l and that soybean plants in which GmSPL12l is silenced have constitutively active immunity, thereby identifying GmSPL12l as a negative regulator of soybean defenses. Collectively, our data present evidence for a virulence function of a rust SSCRP and suggest that PpEC23 is able to suppress soybean immune responses and physically interact with soybean transcription factor GmSPL12l, a negative immune regulator.

A cyst nematode effector binds to diverse plant proteins, increases nematode susceptibility and affects root morphology

Cyst nematodes are plant-parasitic roundworms that are of significance in many cropping systems around the world. Cyst nematode infection is facilitated by effector proteins secreted from the nematode into the plant host. The cDNAs of the 25A01-like effector family are novel sequences that were isolated from the oesophageal gland cells of the soybean cyst nematode (Heterodera glycines). To aid functional characterization, we identified an orthologous member of this protein family (Hs25A01) from the closely related sugar beet cyst nematode H. schachtii, which infects Arabidopsis. Constitutive expression of the Hs25A01 CDS in Arabidopsis plants caused a small increase in root length, accompanied by up to a 22% increase in susceptibility to H. schachtii. A plant-expressed RNA interference (RNAi) construct targeting Hs25A01 transcripts in invading nematodes significantly reduced host susceptibility to H. schachtii. These data document that Hs25A01 has physiological functions in planta and a role in cyst nematode parasitism. In vivo and in vitro binding assays confirmed the specific interactions of Hs25A01 with an Arabidopsis F-box-containing protein, a chalcone synthase and the translation initiation factor eIF-2 β subunit (eIF-2bs), making these proteins probable candidates for involvement in the observed changes in plant growth and parasitism. A role of eIF-2bs in the mediation of Hs25A01 virulence function is further supported by the observation that two independent eIF-2bs Arabidopsis knock-out lines were significantly more susceptible to H. schachtii.

Horizontal gene transfer of acetyltransferases, invertases and chorismate mutases from different bacteria to diverse recipients

BACKGROUND

Hoplolaimina plant-parasitic nematodes (PPN) are a lineage of animals with many documented cases of horizontal gene transfer (HGT). In a recent study, we reported on three likely HGT candidate genes in the soybean cyst nematode Heterodera glycines, all of which encode secreted candidate effectors with putative functions in the host plant. Hg-GLAND1 is a putative GCN5-related N-acetyltransferase (GNAT), Hg-GLAND13 is a putative invertase (INV), and Hg-GLAND16 is a putative chorismate mutase (CM), and blastp searches of the non-redundant database resulted in highest similarity to bacterial sequences. Here, we searched nematode and non-nematode sequence databases to identify all the nematodes possible that contain these three genes, and to formulate hypotheses about when they most likely appeared in the phylum Nematoda. We then performed phylogenetic analyses combined with model selection tests of alternative models of sequence evolution to determine whether these genes were horizontally acquired from bacteria.

RESULTS

Mining of nematode sequence databases determined that GNATs appeared in Hoplolaimina PPN late in evolution, while both INVs and CMs appeared before the radiation of the Hoplolaimina suborder. Also, Hoplolaimina GNATs, INVs and CMs formed well-supported clusters with different rhizosphere bacteria in the phylogenetic trees, and the model selection tests greatly supported models of HGT over descent via common ancestry. Surprisingly, the phylogenetic trees also revealed additional, well-supported clusters of bacterial GNATs, INVs and CMs with diverse eukaryotes and archaea. There were at least eleven and eight well-supported clusters of GNATs and INVs, respectively, from different bacteria with diverse eukaryotes and archaea. Though less frequent, CMs from different bacteria formed supported clusters with multiple different eukaryotes. Moreover, almost all individual clusters containing bacteria and eukaryotes or archaea contained species that inhabit very similar niches.

CONCLUSIONS

GNATs were horizontally acquired late in Hoplolaimina PPN evolution from bacteria most similar to the saprophytic and plant-pathogenic actinomycetes. INVs and CMs were horizontally acquired from bacteria most similar to rhizobacteria and Burkholderia soil bacteria, respectively, before the radiation of Hoplolaimina. Also, these three gene groups appear to have been frequent subjects of HGT from different bacteria to numerous, diverse lineages of eukaryotes and archaea, which suggests that these genes may confer important evolutionary advantages to many taxa. In the case of Hoplolaimina PPN, this advantage likely was an improved ability to parasitize plants.

Eighteen New Candidate Effectors of the Phytonematode Heterodera glycines Produced Specifically in the Secretory Esophageal Gland Cells During Parasitism

Heterodera glycines, the soybean cyst nematode, is the number one pathogen of soybean (Glycine max). This nematode infects soybean roots and forms an elaborate feeding site in the vascular cylinder. H. glycines produces an arsenal of effector proteins in the secretory esophageal gland cells. More than 60 H. glycines candidate effectors were identified in previous gland-cell-mining projects. However, it is likely that additional candidate effectors remained unidentified. With the goal of identifying remaining H. glycines candidate effectors, we constructed and sequenced a large gland cell cDNA library resulting in 11,814 expressed sequence tags. After bioinformatic filtering for candidate effectors using a number of criteria, in situ hybridizations were performed in H. glycines whole-mount specimens to identify candidate effectors whose mRNA exclusively accumulated in the esophageal gland cells, which is a hallmark of many nematode effectors. This approach resulted in the identification of 18 new H. glycines esophageal gland-cell-specific candidate effectors. Of these candidate effectors, 11 sequences were pioneers without similarities to known proteins while 7 sequences had similarities to functionally annotated proteins in databases. These putative homologies provided the bases for the development of hypotheses about potential functions in the parasitism process.

Sequence and Spatiotemporal Expression Analysis of CLE-Motif Containing Genes from the Reniform Nematode (Rotylenchulus reniformis Linford & Oliveira)

The reniform nematode, Rotylenchulus reniformis, is a sedentary semi-endoparasitic species with a host range that encompasses more than 77 plant families. Nematode effector proteins containing plant-ligand motifs similar to CLAVATA3/ESR (CLE) peptides have been identified in the Heterodera, Globodera, and Meloidogyne genera of sedentary endoparasites. Here, we describe the isolation, sequence analysis, and spatiotemporal expression of three R. reniformis genes encoding putative CLE motifs named Rr-cle-1, Rr-cle-2, and Rr-cle-3. The Rr-cle cDNAs showed >98% identity with each other and the predicted peptides were identical with the exception of a short stretch of residues at the carboxy(C)-terminus of the variable domain (VD). Each RrCLE peptide possessed an amino-terminal signal peptide for secretion and a single C-terminal CLE motif that was most similar to Heterodera CLE motifs. Aligning the Rr-cle cDNAs with their corresponding genomic sequences showed three exons with an intron separating the signal peptide from the VD and a second intron separating the VD from the CLE motif. An alignment of the RrCLE1 peptide with Heterodera glycines and Heterodera schachtii CLE proteins revealed a high level of homology within the VD region associated with regulating in planta trafficking of the processed CLE peptide. Quantitative RT-PCR (qRT-PCR) showed similar expression profiles for each Rr-cle transcript across the R. reniformis life-cycle with the greatest transcript abundance being in sedentary parasitic female nematodes. In situ hybridization showed specific Rr-cle expression within the dorsal esophageal gland cell of sedentary parasitic females.

The cyst nematode effector protein 10A07 targets and recruits host posttranslational machinery to mediate its nuclear trafficking and to promote parasitism in Arabidopsis

Plant-parasitic cyst nematodes synthesize and secrete effector proteins that are essential for parasitism. One such protein is the 10A07 effector from the sugar beet cyst nematode, Heterodera schachtii, which is exclusively expressed in the nematode dorsal gland cell during all nematode parasitic stages. Overexpression of H. schachtii 10A07 in Arabidopsis thaliana produced a hypersusceptible phenotype in response to H. schachtii infection along with developmental changes reminiscent of auxin effects. The 10A07 protein physically associates with a plant kinase and the IAA16 transcription factor in the cytoplasm and nucleus, respectively. The interacting plant kinase (IPK) phosphorylates 10A07 at Ser-144 and Ser-231 and mediates its trafficking from the cytoplasm to the nucleus. Translocation to the nucleus is phosphorylation dependent since substitution of Ser-144 and Ser-231 by alanine resulted in exclusive cytoplasmic accumulation of 10A07. IPK and IAA16 are highly upregulated in the nematode-induced syncytium (feeding cells), and deliberate manipulations of their expression significantly alter plant susceptibility to H. schachtii in an additive fashion. An inactive variant of IPK functioned antagonistically to the wild-type IPK and caused a dominant-negative phenotype of reduced plant susceptibility. Thus, exploitation of host processes to the advantage of the parasites is one mechanism by which cyst nematodes promote parasitism of host plants.

Mining novel effector proteins from the esophageal gland cells of Meloidogyne incognita

Meloidogyne incognita is one of the most economically damaging plant pathogens in agriculture and horticulture. Identifying and characterizing the effector proteins which M. incognita secretes into its host plants during infection is an important step toward finding new ways to manage this pest. In this study, we have identified the cDNAs for 18 putative effectors (i.e., proteins that have the potential to facilitate M. incognita parasitism of host plants). These putative effectors are secretory proteins that do not contain transmembrane domains and whose genes are specifically expressed in the secretory gland cells of the nematode, indicating that they are likely secreted from the nematode through its stylet. We have determined that, in the plant cells, these putative effectors are likely to localize to the cytoplasm. Furthermore, the transcripts of many of these novel effectors are specifically upregulated during different stages of the nematode's life cycle, indicating that they function at specific stages during M. incognita parasitism. The predicted proteins showed little to no homology to known proteins from free-living nematode species, suggesting that they evolved recently to support the parasitic lifestyle. On the other hand, several of the effectors are part of gene families within the M. incognita genome as well as that of M. hapla, which points to an important role that these putative effectors are playing in both parasites. With the discovery of these putative effectors, we have increased our knowledge of the effector repertoire utilized by root-knot nematodes to infect, feed on, and reproduce on their host plants. Future studies investigating the roles that these proteins play in planta will help mitigate the effects of this damaging pest.

Members of the Meloidogyne avirulence protein family contain multiple plant ligand-like motifs

Sedentary plant-parasitic nematodes engage in complex interactions with their host plants by secreting effector proteins. Some effectors of both root-knot nematodes (Meloidogyne spp.) and cyst nematodes (Heterodera and Globodera spp.) mimic plant ligand proteins. Most prominently, cyst nematodes secrete effectors that mimic plant CLAVATA3/ESR-related (CLE) ligand proteins. However, only cyst nematodes have been shown to secrete such effectors and to utilize CLE ligand mimicry in their interactions with host plants. Here, we document the presence of ligand-like motifs in bona fide root-knot nematode effectors that are most similar to CLE peptides from plants and cyst nematodes. We have identified multiple tandem CLE-like motifs conserved within the previously identified Meloidogyne avirulence protein (MAP) family that are secreted from root-knot nematodes and have been shown to function in planta. By searching all 12 MAP family members from multiple Meloidogyne spp., we identified 43 repetitive CLE-like motifs composing 14 unique variants. At least one CLE-like motif was conserved in each MAP family member. Furthermore, we documented the presence of other conserved sequences that resemble the variable domains described in Heterodera and Globodera CLE effectors. These findings document that root-knot nematodes appear to use CLE ligand mimicry and point toward a common host node targeted by two evolutionarily diverse groups of nematodes. As a consequence, it is likely that CLE signaling pathways are important in other phytonematode pathosystems as well.

Synchronization of developmental processes and defense signaling by growth regulating transcription factors

Growth regulating factors (GRFs) are a conserved class of transcription factor in seed plants. GRFs are involved in various aspects of tissue differentiation and organ development. The implication of GRFs in biotic stress response has also been recently reported, suggesting a role of these transcription factors in coordinating the interaction between developmental processes and defense dynamics. However, the molecular mechanisms by which GRFs mediate the overlaps between defense signaling and developmental pathways are elusive. Here, we report large scale identification of putative target candidates of Arabidopsis GRF1 and GRF3 by comparing mRNA profiles of the grf1/grf2/grf3 triple mutant and those of the transgenic plants overexpressing miR396-resistant version of GRF1 or GRF3. We identified 1,098 and 600 genes as putative targets of GRF1 and GRF3, respectively. Functional classification of the potential target candidates revealed that GRF1 and GRF3 contribute to the regulation of various biological processes associated with defense response and disease resistance. GRF1 and GRF3 participate specifically in the regulation of defense-related transcription factors, cell-wall modifications, cytokinin biosynthesis and signaling, and secondary metabolites accumulation. GRF1 and GRF3 seem to fine-tune the crosstalk between miRNA signaling networks by regulating the expression of several miRNA target genes. In addition, our data suggest that GRF1 and GRF3 may function as negative regulators of gene expression through their association with other transcription factors. Collectively, our data provide new insights into how GRF1 and GRF3 might coordinate the interactions between defense signaling and plant growth and developmental pathways.

The haustorial transcriptomes of Uromyces appendiculatus and Phakopsora pachyrhizi and their candidate effector families

Haustoria of biotrophic rust fungi are responsible for the uptake of nutrients from their hosts and for the production of secreted proteins, known as effectors, which modulate the host immune system. The identification of the transcriptome of haustoria and an understanding of the functions of expressed genes therefore hold essential keys for the elucidation of fungus-plant interactions and the development of novel fungal control strategies. Here, we purified haustoria from infected leaves and used 454 sequencing to examine the haustorial transcriptomes of Phakopsora pachyrhizi and Uromyces appendiculatus, the causal agents of soybean rust and common bean rust, respectively. These pathogens cause extensive yield losses in their respective legume crop hosts. A series of analyses were used to annotate expressed sequences, including transposable elements and viruses, to predict secreted proteins from the assembled sequences and to identify families of candidate effectors. This work provides a foundation for the comparative analysis of haustorial gene expression with further insights into physiology and effector evolution.

Transcriptome analyses and virus induced gene silencing identify genes in the Rpp4-mediated Asian soybean rust resistance pathway

Rpp4 (Resistance to Phakopsora pachyrhizi 4) confers resistance to Phakopsora pachyrhizi Sydow, the causal agent of Asian soybean rust (ASR). By combining expression profiling and virus induced gene silencing (VIGS), we are developing a genetic framework for Rpp4-mediated resistance. We measured gene expression in mock-inoculated and P. pachyrhizi-infected leaves of resistant soybean accession PI459025B (Rpp4) and the susceptible cultivar (Williams 82) across a 12-day time course. Unexpectedly, two biphasic responses were identified. In the incompatible reaction, genes induced at 12h after infection (hai) were not differentially expressed at 24 hai, but were induced at 72 hai. In contrast, genes repressed at 12 hai were not differentially expressed from 24 to 144 hai, but were repressed 216 hai and later. To differentiate between basal and resistance-gene (R-gene) mediated defence responses, we compared gene expression in Rpp4-silenced and empty vector-treated PI459025B plants 14 days after infection (dai) with P. pachyrhizi. This identified genes, including transcription factors, whose differential expression is dependent upon Rpp4. To identify differentially expressed genes conserved across multiple P. pachyrhizi resistance pathways, Rpp4 expression datasets were compared with microarray data previously generated for Rpp2 and Rpp3-mediated defence responses. Fourteen transcription factors common to all resistant and susceptible responses were identified, as well as fourteen transcription factors unique to R-gene-mediated resistance responses. These genes are targets for future P. pachyrhizi resistance research.

Nematode effector proteins: an emerging paradigm of parasitism

Phytonematodes use a stylet and secreted effectors to modify host cells and ingest nutrients to support their growth and development. The molecular function of nematode effectors is currently the subject of intense investigation. In this review, we summarize our current understanding of nematode effectors, with a particular focus on proteinaceous stylet-secreted effectors of sedentary endoparasitic phytonematodes, for which a wealth of information has surfaced in the past 10 yr. We provide an update on the effector repertoires of several of the most economically important genera of phytonematodes and discuss current approaches to dissecting their function. Lastly, we highlight the latest breakthroughs in effector discovery that promise to shed new light on effector diversity and function across the phylum Nematoda.

A virus-induced gene silencing method to study soybean cyst nematode parasitism in Glycine max

BACKGROUND

Bean pod mottle virus (BPMV) based virus-induced gene silencing (VIGS) vectors have been developed and used in soybean for the functional analysis of genes involved in disease resistance to foliar pathogens. However, BPMV-VIGS protocols for studying genes involved in disease resistance or symbiotic associations with root microbes have not been developed.

FINDINGS

Here we describe a BPMV-VIGS protocol suitable for reverse genetic studies in soybean roots. We use this method for analyzing soybean genes involved in resistance to soybean cyst nematode (SCN). A detailed SCN screening pipeline is described.

CONCLUSIONS

The VIGS method described here provides a new tool to identify genes involved in soybean-nematode interactions. This method could be adapted to study genes associated with any root pathogenic or symbiotic associations.

A ubiquitin carboxyl extension protein secreted from a plant-parasitic nematode Globodera rostochiensis is cleaved in planta to promote plant parasitism

Nematode effector proteins originating from esophageal gland cells play central roles in suppressing plant defenses and in formation of the plant feeding cells that are required for growth and development of cyst nematodes. A gene (GrUBCEP12) encoding a unique ubiquitin carboxyl extension protein (UBCEP) that consists of a signal peptide for secretion, a mono-ubiquitin domain, and a 12 amino acid carboxyl extension protein (CEP12) domain was cloned from the potato cyst nematode Globodera rostochiensis. This GrUBCEP12 gene was expressed exclusively within the nematode's dorsal esophageal gland cell, and was up-regulated in the parasitic second-stage juvenile, correlating with the time when feeding cell formation is initiated. We showed that specific GrUBCEP12 knockdown via RNA interference reduced nematode parasitic success, and that over-expression of the secreted Gr(Δ) (SP) UBCEP12 protein in potato resulted in increased nematode susceptibility, providing direct evidence that this secreted effector is involved in plant parasitism. Using transient expression assays in Nicotiana benthamiana, we found that Gr(Δ) (SP) UBCEP12 is processed into free ubiquitin and a CEP12 peptide (GrCEP12) in planta, and that GrCEP12 suppresses resistance gene-mediated cell death. A target search showed that expression of RPN2a, a gene encoding a subunit of the 26S proteasome, was dramatically suppressed in Gr(Δ) (SP) UBCEP12 but not GrCEP12 over-expression plants when compared with control plants. Together, these results suggest that, when delivered into host plant cells, Gr(Δ) (SP) UBCEP12 becomes two functional units, one acting to suppress plant immunity and the other potentially affecting the host 26S proteasome, to promote feeding cell formation.

The 8D05 parasitism gene of Meloidogyne incognita is required for successful infection of host roots

Parasitism genes encode effector proteins that are secreted through the stylet of root-knot nematodes to dramatically modify selected plant cells into giant-cells for feeding. The Mi8D05 parasitism gene previously identified was confirmed to encode a novel protein of 382 amino acids that had only one database homolog identified on contig 2374 within the Meloidogyne hapla genome. Mi8D05 expression peaked in M. incognita parasitic second-stage juveniles within host roots and its encoded protein was limited to the subventral esophageal gland cells that produce proteins secreted from the stylet. Constitutive expression of Mi8D05 in transformed Arabidopsis thaliana plants induced accelerated shoot growth and early flowering but had no visible effects on root growth. Independent lines of transgenic Arabidopsis that expressed a double-stranded RNA complementary to Mi8D05 in host-derived RNA interference (RNAi) tests had up to 90% reduction in infection by M. incognita compared with wild-type control plants, suggesting that Mi8D05 plays a critical role in parasitism by the root-knot nematode. Yeast two-hybrid experiments confirmed the specific interaction of the Mi8D05 protein with plant aquaporin tonoplast intrinsic protein 2 (TIP2) and provided evidence that the Mi8D05 effector may help regulate solute and water transport within giant-cells to promote the parasitic interaction.

Isolation of whole esophageal gland cells from plant-parasitic nematodes for transcriptome analyses and effector identification

Esophageal glands of plant-parasitic nematodes are highly specialized cells whose gene expression products include secreted effector proteins, which govern nematode parasitism of host plants. Therefore, elucidating the transcriptomes of esophageal glands with the goal of identifying nematode effectors is a promising avenue to understanding nematode parasitism and its evolutionary origins as well as to devising nematode control strategies. We have developed a method to separate and isolate individual esophageal gland cells from multiple species of plant-parasitic nematodes while preserving RNA quality. We have used such isolated gland cells for transcriptome analysis via high-throughput DNA sequencing. This method relies on the differential histochemical staining of the gland cells after homogenization of phytonematode tissues. Total RNA was extracted from whole gland cells isolated from eight different plant-parasitic nematode species. To validate this approach, the isolated RNA from three plant-parasitic nematode species-Globodera rostochiensis, Pratylenchus penetrans, and Radopholus similis-was amplified, gel purified, and used for 454 sequencing. We obtained 456,801 total reads with an average read length of 409 bp. Sequence analyses revealed the presence of homologs of previously known nematode effectors in these libraries, thus validating our approach. These data provide compelling evidence that this technical advance can be used to relatively easily and expediently discover effector repertoires of plant-parasitic nematodes.

Manipulation of plant cells by cyst and root-knot nematode effectors

A key feature of sedentary plant-parasitic nematodes is the release of effector proteins from their esophageal gland cells through their stylets into host roots. These proteinaceous stylet secretions have been shown to be crucial for successful parasitism by mediating the transition of normal root cells into specialized feeding sites and by negating plant defenses. Recent technical advances of purifying mRNA from esophageal gland cells of plant-parasitic nematodes coupled with emerging sequencing technologies is steadily expanding our knowledge of nematode effector repertoires. Host targets and biological activities of a number of nematode effectors are continuously being reported and, by now, a first picture of the complexity of sedentary nematode parasitism at the molecular level is starting to take shape. In this review, we highlight effector mechanisms that recently have been uncovered by studying the host-pathogen interaction. These mechanisms range from mediating susceptibility of host plants to the actual triggering of defense responses. In particular, we portray and discuss the mechanisms by which nematode effectors modify plant cell walls, negate host defense responses, alter auxin and polyamine signaling, mimic plant molecules, regulate stress signaling, and activate hypersensitive responses. Continuous molecular characterization of newly discovered nematode effectors will be needed to determine how these effectors orchestrate host signaling pathways and biological processes leading to successful parasitism.

A soybean cyst nematode resistance gene points to a new mechanism of plant resistance to pathogens

Soybean (Glycine max (L.) Merr.) is an important crop that provides a sustainable source of protein and oil worldwide. Soybean cyst nematode (Heterodera glycines Ichinohe) is a microscopic roundworm that feeds on the roots of soybean and is a major constraint to soybean production. This nematode causes more than US$1 billion in yield losses annually in the United States alone, making it the most economically important pathogen on soybean. Although planting of resistant cultivars forms the core management strategy for this pathogen, nothing is known about the nature of resistance. Moreover, the increase in virulent populations of this parasite on most known resistance sources necessitates the development of novel approaches for control. Here we report the map-based cloning of a gene at the Rhg4 (for resistance to Heterodera glycines 4) locus, a major quantitative trait locus contributing to resistance to this pathogen. Mutation analysis, gene silencing and transgenic complementation confirm that the gene confers resistance. The gene encodes a serine hydroxymethyltransferase, an enzyme that is ubiquitous in nature and structurally conserved across kingdoms. The enzyme is responsible for interconversion of serine and glycine and is essential for cellular one-carbon metabolism. Alleles of Rhg4 conferring resistance or susceptibility differ by two genetic polymorphisms that alter a key regulatory property of the enzyme. Our discovery reveals an unprecedented plant resistance mechanism against a pathogen. The mechanistic knowledge of the resistance gene can be readily exploited to improve nematode resistance of soybean, an increasingly important global crop.

Temporal and spatial Bean pod mottle virus-induced gene silencing in soybean

Virus-induced gene silencing (VIGS) is a powerful reverse genetics tool in plant science. In this study, we investigated the temporal and spatial silencing patterns achieved by Bean pod mottle virus (BPMV)-based VIGS in soybean using virus constructs targeting green fluorescence protein (GFP). Silencing GFP enabled an in-depth analysis of silencing in soybean tissues over time in a transgenic line constitutively expressing GFP. We discovered evidence for variable GFP silencing based on insert orientation and targeted region in the coding sequence. A 3' sequence in reverse orientation produced the strongest silencing phenotypes. Furthermore, we documented that BPMV VIGS can achieve widespread silencing in a broad range of tissues, including leaves, stems, flowers and roots. Near-complete silencing was attained in leaves and flowers. Although weaker than in shoots, the observed gene silencing in soybean roots will also allow reverse genetics studies in this tissue. When GFP fluorescence was assayed in cross-sections of stems and leaf petioles, near-complete and uniform silencing was observed in all cell types. Silencing was observed from as early as 2 weeks post-virus inoculation in leaves to 7 weeks post-virus inoculation in flowers, suggesting that this system can induce and maintain silencing for significant durations.

Complex feedback regulations govern the expression of miRNA396 and its GRF target genes

The beet cyst nematode, Heterodera schachtii, is a sedentary root parasite that induces the formation of a specialized root feeding structure, the syncytium. We previously have shown that coordinated regulation of miR396 and its target genes GRF1 and GRF3 in the syncytium is required for proper formation. To gain a better understanding of this coordinated regulation, we used quantitative real-time PCR to assess the abundance of primary (pri)-miRNA396a, pri-miRNA396b and mature miRNA396 in transgenic Arabidopsis plants overexpressing either wild-type variants of the GRF1 or GRF3 coding sequences or miR396-resistant variants. We also included a grf1/grf2/grf3 triple mutant in these analyses. We observed significant decreases in the abundance of pri-miRNA396a, pri-miRNA396b and mature miR396 in the transgenic plants overexpressing GRF1 or GRF3, particularly with the miRNA396-resistant variants. In contrast, the primary transcripts and mature miRNA396 abundance were significantly increased in the grf1/grf2/grf3 triple knockout mutant. These results demonstrate that homeostasis between miR396 and the target genes GRF1 and GRF3 is established through reciprocal feedback regulation, in which GRF1/GRF3 and miR396 negatively regulate each other's expression. In addition, we found that constitutive expression of GRF1 or GRF3 decreases the mRNA abundance of other GRFs, even those that are not targeted by miR396, as well as their own endogenous transcripts, which documents further regulatory facets of this equilibrium.

The interaction of the novel 30C02 cyst nematode effector protein with a plant β-1,3-endoglucanase may suppress host defence to promote parasitism

Phytoparasitic nematodes secrete an array of effector proteins to modify selected recipient plant cells into elaborate and essential feeding sites. The biological function of the novel 30C02 effector protein of the soybean cyst nematode, Heterodera glycines, was studied using Arabidopsis thaliana as host and the beet cyst nematode, Heterodera schachtii, which contains a homologue of the 30C02 gene. Expression of Hg30C02 in Arabidopsis did not affect plant growth and development but increased plant susceptibility to infection by H. schachtii. The 30C02 protein interacted with a specific (AT4G16260) host plant β-1,3-endoglucanase in both yeast and plant cells, possibly to interfere with its role as a plant pathogenesis-related protein. Interestingly, the peak expression of 30C02 in the nematode and peak expression of At4g16260 in plant roots coincided at around 3-5 d after root infection by the nematode, after which the relative expression of At4g16260 declined significantly. An Arabidopsis At4g16260 T-DNA mutant showed increased susceptibility to cyst nematode infection, and plants that overexpressed At4g16260 were reduced in nematode susceptibility, suggesting a potential role of host β-1,3-endoglucanase in the defence response against H. schachtii infection. Arabidopsis plants that expressed dsRNA and its processed small interfering RNA complementary to the Hg30C02 sequence were not phenotypically different from non-transformed plants, but they exhibited a strong RNA interference-mediated resistance to infection by H. schachtii. The collective results suggest that, as with other pathogens, active suppression of host defence is a critical component for successful parasitism by nematodes and a vulnerable target to disrupt the parasitic cycle.

The Arabidopsis microRNA396-GRF1/GRF3 regulatory module acts as a developmental regulator in the reprogramming of root cells during cyst nematode infection

The syncytium is a unique plant root organ whose differentiation is induced by plant-parasitic cyst nematodes to create a source of nourishment. Syncytium formation involves the redifferentiation and fusion of hundreds of root cells. The underlying regulatory networks that control this unique change of plant cell fate are not understood. Here, we report that a strong down-regulation of Arabidopsis (Arabidopsis thaliana) microRNA396 (miR396) in cells giving rise to the syncytium coincides with the initiation of the syncytial induction/formation phase and that specific miR396 up-regulation in the developed syncytium marks the beginning of the maintenance phase, when no new cells are incorporated into the syncytium. In addition, our results show that miR396 in fact has a role in the transition from one phase to the other. Expression modulations of miR396 and its Growth-Regulating Factor (GRF) target genes resulted in reduced syncytium size and arrested nematode development. Furthermore, genome-wide expression profiling revealed that the miR396-GRF regulatory system can alter the expression of 44% of the more than 7,000 genes reported to change expression in the Arabidopsis syncytium. Thus, miR396 represents a key regulator for the reprogramming of root cells. As such, this regulatory unit represents a powerful molecular target for the parasitic animal to modulate plant cells and force them into novel developmental pathways.

Arabidopsis peroxidase AtPRX53 influences cell elongation and susceptibility to Heterodera schachtii

Cyst nematodes establish and maintain feeding sites (syncytia) in the roots of host plants by altering expression of host genes. Among these genes are members of the large gene family of class III peroxidases, which have reported functions in a variety of biological processes. In this study, we used Arabidopsis-Heterodera schachtii as a model system to functionally characterize peroxidase 53 (AtPRX53). Promoter assays showed that under non-infected conditions AtPRX53 is expressed mainly in the root, the hypocotyl and the base of the pistil. Under infected conditions, the AtPRX53 promoter showed upregulation at the nematode penetration sites and in their migration paths. Interestingly, strong GUS activity was observed in H. schachtii-induced syncytia during the early stage of infection and remained strong in the syncytia of third-stage juveniles. Also, AtPRX53 showed upregulation in response to wounding and jasmonic acid treatments. Manipulation of AtPRX53 expression through overexpression and knockout mutation affected both plant morphology and nematode susceptibility. While AtPRX53 overexpression lines exhibited short hypocotyls, aberrant flower development and reduced nematode susceptibility to H. schachtii, the atprx53 mutant showed long hypocotyls and a 3-carpel silique phenotype as well as a non significant increase of nematode susceptibility. Taken together these data, therefore, indicate diverse roles of AtPRX53 in the wound response, flower development and syncytium formation.