Comparative transcriptomic analysis of soybean cyst nematode inbred populations non-adapted or adapted on soybean rhg1-a/Rhg4-mediated resistance.

Soybean cyst nematode (SCN, Heterodera glycines) is most effectively managed through planting resistant soybean cultivars, but the repeated use of the same resistance sources has led to a widespread emergence of virulent SCN populations that can overcome soybean resistance. Resistance to SCN HG type 0 (Race 3) in soybean cultivar Forrest is mediated by an epistatic interaction between the soybean resistance genes rhg1-a and Rhg4. We previously developed two SCN inbred populations by mass-selecting SCN HG type 0 (Race 3) on susceptible and resistant recombinant inbred lines, derived from a cross between Forrest and the SCN-susceptible cultivar Essex, which differ for Rhg4. To identify SCN genes potentially involved in overcoming rhg1-a/Rhg4-mediated resistance, we conducted RNA-sequencing on early parasitic juveniles of these two SCN inbred populations infecting their respective hosts, only to discover a handful of differentially expressed genes (DEGs). However, in a comparison to early parasitic juveniles of an avirulent SCN inbred population infecting a resistant host, we discovered 59 and 171 DEGs uniquely up- or down-regulated in virulent parasitic juveniles adapted on the resistant host. Interestingly, the proteins coded by these 59 DEGs included vitamin B-associated proteins (reduced folate carrier, biotin synthase, and thiamine transporter) and nematode effectors known to play roles in plant defense suppression, suggesting that virulent SCN may exert a heightened transcriptional response to cope with enhanced plant defenses and an altered nutritional status of a resistant soybean host.

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.

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.

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.

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.

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.

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.

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.

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.

Divergent evolution of arrested development in the dauer stage of Caenorhabditis elegans and the infective stage of Heterodera glycines.

BACKGROUND: The soybean cyst nematode Heterodera glycines is the most important parasite in soybean production worldwide. A comprehensive analysis of large-scale gene expression changes throughout the development of plant-parasitic nematodes has been lacking to date. RESULTS: We report an extensive genomic analysis of H. glycines, beginning with the generation of 20,100 expressed sequence tags (ESTs). In-depth analysis of these ESTs plus approximately 1,900 previously published sequences predicted 6,860 unique H. glycines genes and allowed a classification by function using InterProScan. Expression profiling of all 6,860 genes throughout the H. glycines life cycle was undertaken using the Affymetrix Soybean Genome Array GeneChip. Our data sets and results represent a comprehensive resource for molecular studies of H. glycines. Demonstrating the power of this resource, we were able to address whether arrested development in the Caenorhabditis elegans dauer larva and the H. glycines infective second-stage juvenile (J2) exhibits shared gene expression profiles. We determined that the gene expression profiles associated with the C. elegans dauer pathway are not uniformly conserved in H. glycines and that the expression profiles of genes for metabolic enzymes of C. elegans dauer larvae and H. glycines infective J2 are dissimilar. CONCLUSION: Our results indicate that hallmark gene expression patterns and metabolism features are not shared in the developmentally arrested life stages of C. elegans and H. glycines, suggesting that developmental arrest in these two nematode species has undergone more divergent evolution than previously thought and pointing to the need for detailed genomic analyses of individual parasite species.