De novo transcriptome sequencing and analysis of the cereal cyst nematode, Heterodera avenae

Chromosome-level genome assembly of the cereal cyst nematode Heterodera flipjevi

As an economically important plant parasitic nematode (PPN), Heterodera filipjevi causes great damage on wheat, and now it was widely recorded in many countries. While multiple genomes of PPNs have been published, high-quality genome assembly and annotation on H. filipjevi have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for H. filipjevi, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The genome consists of 9 pseudo-chromosomes that contain 134.19 Mb of sequence, with a scaffold N50 length of 11.88 Mb. In total, 10,036 genes were annotated, representing 75.20% of the total predicted protein-coding genes. Our study provides the first chromosome-scale genome for H. filipjevi, which is also the inaugural high-quality genome of cereal cyst nematodes (CCNs). It provides a valuable genomic resource for further biological research and pest management of cereal cyst nematodes disease.

Identification of genes governing resistance to PCN (Globodera rostochiensis) through transcriptome analysis in Solanum tuberosum

Potato cyst nematodes (PCNs) are major pests worldwide that affect potato production. The molecular changes happening in the roots upon PCN infection are still unknown. Identification of transcripts and genes governing PCN resistance will help in the development of resistant varieties. Hence, differential gene expression of compatible (Kufri Jyoti) and incompatible (JEX/A-267) potato genotypes was studied before (0 DAI) and after (10 DAI) inoculation of Globodera rostochiensis J2s through RNA sequencing (RNA-Seq). Total sequencing reads generated ranged between 33 and 37 million per sample, with a read mapping of 48-84% to the potato reference genome. In the infected roots of the resistant genotype JEX/A-267, 516 genes were downregulated, and 566 were upregulated. In comparison, in the susceptible genotype Kufri Jyoti, 316 and 554 genes were downregulated and upregulated, respectively. Genes encoding cell wall proteins, zinc finger protein, WRKY transcription factors, MYB transcription factors, disease resistance proteins, and pathogenesis-related proteins were found to be majorly involved in the incompatible reaction after PCN infection in the resistant genotype, JEX/A-267. Furthermore, RNA-Seq results were validated through quantitative real-time PCR (qRT-PCR), and it was observed that ATP, FLAVO, CYTO, and GP genes were upregulated at 5 DAI, which was subsequently downregulated at 10 DAI. The genes encoding ATP, FLAVO, LBR, and GP were present in > 1.5 fold before infection in JEX-A/267 and upregulated 7.9- to 27.6-fold after 5 DAI; subsequently, most of these genes were downregulated to 0.9- to 2.8-fold, except LBR, which was again upregulated to 44.4-fold at 10 DAI.

PAL-mediated SA biosynthesis pathway contributes to nematode resistance in wheat

The pathogen cereal cyst nematode (CCN) is deleterious to Triticeae crops and is a threat to the global crop yield. Accession no. 1 of Aegilops variabilis, a relative of Triticum aestivum (bread wheat), is highly resistant to CCN. Our previous study demonstrated that the expression of the phenylalanine ammonia lyase (PAL) gene AevPAL1 in Ae. variabilis is strongly induced by CCN. PAL, the first enzyme of phenylpropanoid metabolism, is involved in abiotic and biotic stress responses. However, its role in plant-CCN interaction remains unknown. In the present study, we proved that AevPAL1 helps to confer CCN resistance through affecting the synthesis of salicylic acid (SA) and downstream secondary metabolites. The silencing of AevPAL1 increased the incidence of CCN infection in roots and decreased the accumulation of SA and phenylalanine (Phe)-derived specialized metabolites. The exogenous pre-application of SA also improved CCN resistance. Additionally, the functions of PAL in phenylpropanoid metabolism correlated with tryptophan decarboxylase (TDC) functioning in tryptophan metabolism pathways. The silencing of either AevPAL1 or AevTDC1 exhibited a concomitant reduction in the expression of both genes and the contents of metabolites downstream of PAL and TDC. These results suggested that AevPAL1, possibly in coordination with AevTDC1, positively contributes to CCN resistance by altering the downstream secondary metabolites and SA content in Ae. variabilis. Moreover, AevPAL1 overexpression significantly enhanced CCN resistance in bread wheat and did not exhibit significant negative effects on yield-related traits, suggesting that AevPAL1 is valuable for the genetic improvement of CCN resistance in bread wheat.

A Revised Adaptation of the Smart-Seq2 Protocol for Single-Nematode RNA-Seq

The advancement of transcriptomic studies in plant parasitic nematodes will greatly benefit from the development of single-nematode RNA-seq methods. Since many plant parasitic nematodes are obligate parasites, it is often difficult to efficiently obtain sufficient amounts of nematodes for transcriptomic studies. Here we have adapted SMART-Seq2 for single-nematode RNA-seq requiring only an individual nematode for a sample replicate. This protocol provides a detailed step-by-step procedure of the RNA-seq workflow starting from lysis of the nematode to quantification of transcripts using a user-friendly online platform.

Transcriptomic changes in the pre-parasitic juveniles of Meloidogyne incognita induced by silencing of effectors Mi-msp-1 and Mi-msp-20

Plant-parasitic root-knot nematode Meloidogyne incognita uses an array of effector proteins to establish successful plant infections. Mi-msp-1 and Mi-msp-20 are two known effectors secreted from nematode subventral oesophageal glands; Mi-msp-1 being a putative secretory venom allergen AG5-like protein, whereas Mi-msp-20 is a pioneer gene with a coiled-coil motif. Expression of specific effector is known to cause disturbances in the expression of other effectors. Here, we used RNA-Seq to investigate the pleiotropic effects of silencing Mi-msp-1 and Mi-msp-20. A total of 25.1-51.9 million HQ reads generated from Mi-msp-1 and Mi-msp-20 silenced second-stage juveniles (J2s) along with freshly hatched J2s were mapped to an already annotated M. incognita proteome to understand the impact on various nematode pathways. As compared to control, silencing of Mi-msp-1 caused differential expression of 29 transcripts, while Mi-msp-20 silencing resulted in differential expression of a broader set of 409 transcripts. In the Mi-msp-1 silenced J2s, cytoplasm (GO:0005737) was the most enriched gene ontology (GO) term, whereas in the Mi-msp-20 silenced worms, embryo development (GO:0009792), reproduction (GO:0000003) and nematode larval development (GO:0002119) were the most enriched terms. Limited crosstalk was observed between these two effectors as a sheer 5.9% of the up-regulated transcripts were common between Mi-msp-1 and Mi-msp-20 silenced nematodes. Our results suggest that in addition to the direct knock-down caused by silencing of Mi-msp-1 and Mi-msp-20, the cascading effect on other genes might also be contributing to a reduction in nematode's parasitic abilities.

Signatures of adaptation to a monocot host in the plant-parasitic cyst nematode Heterodera sacchari

Interactions between plant-parasitic nematodes and their hosts are mediated by effectors, i.e. secreted proteins that manipulate the plant to the benefit of the pathogen. To understand the role of effectors in host adaptation in nematodes, we analysed the transcriptome of Heterodera sacchari, a cyst nematode parasite of rice (Oryza sativa) and sugarcane (Saccharum officinarum). A multi-gene phylogenetic analysis showed that H. sacchari and the cereal cyst nematode Heterodera avenae share a common evolutionary origin and that they evolved to parasitise monocot plants from a common dicot-parasitic ancestor. We compared the effector repertoires of H. sacchari with those of the dicot parasites Heterodera glycines and Globodera rostochiensis to understand the consequences of this transition. While, in general, effector repertoires are similar between the species, comparing effectors and non-effectors of H. sacchari and G. rostochiensis shows that effectors have accumulated more mutations than non-effectors. Although most effectors show conserved spatiotemporal expression profiles and likely function, some H. sacchari effectors are adapted to monocots. This is exemplified by the plant-peptide hormone mimics, the CLAVATA3/EMBRYO SURROUNDING REGION-like (CLE) effectors. Peptide hormones encoded by H. sacchari CLE effectors are more similar to those from rice than those from other plants, or those from other plant-parasitic nematodes. We experimentally validated the functional significance of these observations by demonstrating that CLE peptides encoded by H. sacchari induce a short root phenotype in rice, whereas those from a related dicot parasite do not. These data provide a functional example of effector evolution that co-occurred with the transition from a dicot-parasitic to a monocot-parasitic lifestyle.

Colonization of Beauveria bassiana 08F04 in root-zone soil and its biocontrol of cereal cyst nematode (Heterodera filipjevi)

Cereal cyst nematodes cause serious yield losses of wheat in Hunaghuai winter wheat growing region in China. Beauveria bassiana 08F04 isolated from the surface of cysts is a promising biological control agent for cereal cyst nematodes. As the colonization capacity is a crucial criteria to assess biocontrol effectiveness for a microbial agent candidate, we aimed to label B. bassiana 08F04 for efficient monitoring of colonization in the soil. The binary pCAM-gfp plasmid containing sgfp and hph was integrated into B. bassiana 08F04 using the Agrobacterium tumefaciens-mediated transformation. The transformation caused a significant change in mycelial and conidial yields, and in extracellular chitinase activity in some transformants. The cultural filtrates of some transformants also decreased acetylcholinesterase activity and the survival of Heterodera filipjevi second-stage juveniles relative to the wild-type strain. One transformant (G10) had a growth rate and biocontrol efficacy similar to the wild-type strain, so it was used for a pilot study of B. bassiana colonization conducted over 13 weeks. Real-time PCR results and CFU counts revealed that the population of G10 increased quickly over the first 3 weeks, then decreased slowly over the following 4 weeks before stabilizing. In addition, the application of wild-type B. bassiana 08F04 and transformant G10 significantly reduced the number of H. filipjevi females in roots by 64.4% and 60.2%, respectively. The results of this study have practical applications for ecological, biological and functional studies of B. bassiana 08F04 and for bionematicide registration.

Transcriptome and Parasitome Analysis of Beet Cyst Nematode Heterodera schachtii

Beet cyst nematodes depend on a set of secretory proteins (effectors) for the induction and maintenance of their syncytial feeding sites in plant roots. In order to understand the relationship between the beet cyst nematode H. schachtii and its host, identification of H. schachtii effectors is crucial and to this end, we sequenced a whole animal pre-infective J2-stage transcriptome in addition to pre- and post-infective J2 gland cell transcriptome using Next Generation Sequencing (NGS) and identified a subset of sequences representing putative effectors. Comparison between the transcriptome of H. schachtii and previously reported related cyst nematodes and root-knot nematodes revealed a subset of esophageal gland related sequences and putative effectors in common across the tested species. Structural and functional annotation of H. schachtii transcriptome led to the identification of nearly 200 putative effectors. Six putative effector expressions were quantified using qPCR and three of them were functionally analyzed using RNAi. Phenotyping of the RNAi nematodes indicated that all tested genes decrease the level of nematodes pathogenicity and/or the average female size, thereby regulating cyst nematode parasitism. These discoveries contribute to further understanding of the cyst nematode parasitism.

A secreted metal-binding protein protects necrotrophic phytopathogens from reactive oxygen species

Few secreted proteins involved in plant infection common to necrotrophic bacteria, fungi and oomycetes have been identified except for plant cell wall-degrading enzymes. Here we study a family of iron-binding proteins that is present in Gram-negative and Gram-positive bacteria, fungi, oomycetes and some animals. Homolog proteins in the phytopathogenic bacterium Dickeya dadantii (IbpS) and the fungal necrotroph Botrytis cinerea (BcIbp) are involved in plant infection. IbpS is secreted, can bind iron and copper, and protects the bacteria against H₂O₂-induced death. Its 1.7 Å crystal structure reveals a classical Venus Fly trap fold that forms dimers in solution and in the crystal. We propose that secreted Ibp proteins binds exogenous metals and thus limit intracellular metal accumulation and ROS formation in the microorganisms.

The authors identify a family of iron-binding proteins that is present in phytopathogenic bacteria, fungi and oomycetes. Some of these proteins are secreted, bind metals, protect the pathogen from H₂O₂-induced death, and are involved in plant infection.

De novo transcriptome assembly of a facultative parasitic nematode Pelodera (syn. Rhabditis) strongyloides

Pelodera strongyloides is a generally free-living gonochoristic facultative nematode. The whole genomic sequence of P. strongyloides remains unknown but 4 small subunit ribosomal RNA (ssrRNA) gene sequences are available. This project launched a de novo transcriptome assembly with 100 bp paired-end RNA-seq reads from normal, starved and wet-plate cultured animals. Trinity assembly tool generated 104,634 transcript contigs with N50 contig being 2195 bp and average contig length at 1103 bp. Transcriptome BLASTX matching results of five nematodes (C. elegans, Strongyloides stercoralis, Necator americanus, Trichuris trichiura, and Pristionchus pacificus) were consistent with their evolutionary relationships. Sixteen genes were identified to be homologous to key elements of the C. elegans RNA interference system, such as Dicer, Argonaute, RNA-dependent RNA polymerase and double strand RNA transport proteins. In starved samples, we observed up-regulation of cuticle related genes and 3 dauer formation genes. Dauer morphology was captured with enlarged phasmid under light microscopy, and dauer and normal larvae counts in clumps had a Pearson's product-moment correlation of 0.805 with P-value = 0.0088. Our results demonstrate that P. strongyloides could be used for studying nematode-related human or pet parasitic diseases. The sequenced assembled transcriptome reported here may be useful to understand the evolution of parasitism in Nematoda.

Plant Proteins and Processes Targeted by Parasitic Nematode Effectors

Sedentary endoparasitic nematodes, such as root-knot nematodes (RKN; Meloidogyne spp.) and cyst nematodes (CN; Heterodera spp. and Globodera spp.) cause considerable damage to agricultural crops. RKN and CN spend most of their life cycle in plant roots, in which they induce the formation of multinucleate hypertrophied feeding cells, called "giant cells" and "syncytia," respectively. The giant cells result from nuclear divisions of vascular cells without cytokinesis. They are surrounded by small dividing cells and they form a new organ within the root known as a root knot or gall. CN infection leads to the fusion of several root cells into a unique syncytium. These dramatically modified host cells act as metabolic sinks from which the nematode withdraws nutrients throughout its life, and they are thus essential for nematode development. Both RKN and CN secrete effector proteins that are synthesized in the oesophageal glands and delivered to the appropriate cell in the host plant via a syringe-like stylet, triggering the ontogenesis of the feeding structures. Within the plant cell or in the apoplast, effectors associate with specific host proteins, enabling them to hijack important processes for cell morphogenesis and physiology or immunity. Here, we review recent findings on the identification and functional characterization of plant targets of RKN and CN effectors. A better understanding of the molecular determinants of these biotrophic relationships would enable us to improve the yields of crops infected with parasitic nematodes and to expand our comprehension of root development.

Biocontrol of Cereal Crop Diseases Using Streptomycetes

A growing world population and an increasing demand for greater food production requires that crop losses caused by pests and diseases are dramatically reduced. Concurrently, sustainability targets mean that alternatives to chemical pesticides are becoming increasingly desirable. Bacteria in the plant root microbiome can protect their plant host against pests and pathogenic infection. In particular, Streptomyces species are well-known to produce a range of secondary metabolites that can inhibit the growth of phytopathogens. Streptomyces are abundant in soils and are also enriched in the root microbiomes of many different plant species, including those grown as economically and nutritionally valuable cereal crops. In this review we discuss the potential of Streptomyces to protect against some of the most damaging cereal crop diseases, particularly those caused by fungal pathogens. We also explore factors that may improve the efficacy of these strains as biocontrol agents in situ, as well as the possibility of exploiting plant mechanisms, such as root exudation, that enable the recruitment of microbial species from the soil to the root microbiome. We argue that a greater understanding of these mechanisms may enable the development of protective plant root microbiomes with a greater abundance of beneficial bacteria, such as Streptomyces species.

Life-stage specific transcriptomes of a migratory endoparasitic plant nematode, Radopholus similis elucidate a different parasitic and life strategy of plant parasitic nematodes

Radopholus similis is an important migratory endoparasitic nematode, severely harms banana, citrus and many other commercial crops. Little is known about the molecular mechanism of infection and pathogenesis of R. similis. In this study, 64761 unigenes were generated from eggs, juveniles, females and males of R. similis. 11443 unigenes showed significant expression difference among these four life stages. Genes involved in host parasitism, anti-host defense and other biological processes were predicted. There were 86 and 102 putative genes coding for cell wall degrading enzymes and antioxidase respectively. The amount and type of putative parasitic-related genes reported in sedentary endoparasitic plant nematodes are variable from those of migratory parasitic nematodes on plant aerial portion. There were no sequences annotated to effectors in R. similis, involved in feeding site formation of sedentary endoparasites nematodes. This transcriptome data provides a new insight into the parasitic and pathogenic molecular mechanisms of the migratory endoparasitic nematodes. It also provides a broad idea for further research on R. similis.

RNAi-Based Biocontrol of Wheat Nematodes Using Natural Poly-Component Biostimulants

With the growing global demands on sustainable food production, one of the biggest challenges to agriculture is associated with crop losses due to parasitic nematodes. While chemical pesticides have been quite successful in crop protection and mitigation of damage from parasites, their potential harm to humans and environment, as well as the emergence of nematode resistance, have necessitated the development of viable alternatives to chemical pesticides. One of the most promising and targeted approaches to biocontrol of parasitic nematodes in crops is that of RNA interference (RNAi). In this study we explore the possibility of using biostimulants obtained from metabolites of soil streptomycetes to protect wheat (Triticum aestivum L.) against the cereal cyst nematode Heterodera avenae by means of inducing RNAi in wheat plants. Theoretical models of uptake of organic compounds by plants, and within-plant RNAi dynamics, have provided us with useful insights regarding the choice of routes for delivery of RNAi-inducing biostimulants into plants. We then conducted in planta experiments with several streptomycete-derived biostimulants, which have demonstrated the efficiency of these biostimulants at improving plant growth and development, as well as in providing resistance against the cereal cyst nematode. Using dot blot hybridization we demonstrate that biostimulants trigger a significant increase of the production in plant cells of si/miRNA complementary with plant and nematode mRNA. Wheat germ cell-free experiments show that these si/miRNAs are indeed very effective at silencing the translation of nematode mRNA having complementary sequences, thus reducing the level of nematode infestation and improving plant resistance to nematodes. Thus, we conclude that natural biostimulants produced from metabolites of soil streptomycetes provide an effective tool for biocontrol of wheat nematode.

Two venom allergen-like proteins, HaVAP1 and HaVAP2, are involved in the parasitism of Heterodera avenae

Despite the fact that venom allergen-like proteins (VAPs) have been identified in many animal- and plant-parasitic nematodes, studies on VAPs in Heterodera avenae, which is an important phytonematode, are still in their infancy. Here, we isolated, cloned and characterized two VAPs, named HaVAP1 and HaVAP2, from H. avenae. The two encoded proteins, HaVAP1 and HaVAP2, harbour an SCP-like domain each, but share only 38% identity with each other. HaVAP1 and HaVAP2 are expressed in subventral and dorsal oesophageal glands, respectively. HaVAP1 is expressed mainly at the early stages, whereas HaVAP2 accumulates principally at the late stages. Both HaVAP1 and HaVAP2 are secreted when expressed in Nicotiana benthamiana leaves, but HaVAP1 is delivered into chloroplasts, whereas HaVAP2 is translocated to the nucleus without signal peptides. Knocking down HaVAP1 increased the virulence of H. avenae. In contrast, silencing of HaVAP2 hampered the parasitism of H. avenae. Both HaVAP1 and HaVAP2 suppressed the cell death induced by BAX in N. benthamiana leaves. Moreover, HaVAP2 physically interacted with a CYPRO4-like protein (HvCLP) of Hordeum vulgare in the nucleus of the plant. It is reasonable to speculate that the changes in the transcript of HvCLP are associated with HaVAP2 during the parasitism of H. avenae. All results obtained in this study show that both HaVAP1 and HaVAP2 are involved in the parasitism of H. avenae, but they possess different functions, broadening our understanding of the parasitic mechanism of H. avenae.

Transcriptional profiling of wheat (Triticum aestivum L.) during a compatible interaction with the cereal cyst nematode Heterodera avenae

Cereal cyst nematode (CCN, Heterodera avenae) presents severe challenges to wheat (Triticum aestivum L.) production worldwide. An investigation of the interaction between wheat and CCN can greatly improve our understanding of how nematodes alter wheat root metabolic pathways for their development and could contribute to new control strategies against CCN. In this study, we conducted transcriptome analyses of wheat cv. Wen 19 (Wen19) by using RNA-Seq during the compatible interaction with CCN at 1, 3 and 8 days past inoculation (dpi). In total, 71,569 transcripts were identified, and 10,929 of them were examined as differentially expressed genes (DEGs) in response to CCN infection. Based on the functional annotation and orthologous findings, the protein phosphorylation, oxidation-reduction process, regulation of transcription, metabolic process, transport, and response process as well as many other pathways previously reported were enriched at the transcriptional level. Plant cell wall hydrolysis and modifying proteins, auxin biosynthesis, signalling and transporter genes were up-regulated by CCN infection to facilitate penetration, migration and syncytium establishment. Genes responding to wounding and jasmonic acid stimuli were enriched at 1 dpi. We found 16 NBS-LRR genes, 12 of which were down-regulated, indicating the repression of resistance. The expression of genes encoding antioxidant enzymes, glutathione S-transferases and UDP-glucosyltransferase was significantly up-regulated during CCN infection, indicating that they may play key roles in the compatible interaction of wheat with CCN. Taken together, the results obtained from the transcriptome analyses indicate that the genes involved in oxidation-reduction processes, induction and suppression of resistance, metabolism, transport and syncytium establishment may be involved in the compatible interaction of Wen 19 with CCN. This study provides new insights into the responses of wheat to CCN infection. These insights could facilitate the elucidation of the potential mechanisms of wheat responses to CCN.

Resistance to Cereal Cyst Nematodes in Wheat and Barley: An Emphasis on Classical and Modern Approaches

Cereal cyst nematodes (CCNs) are among the most important nematode pests that limit production of small grain cereals like wheat and barley. These nematodes alone are estimated to reduce production of crops by 10% globally. This necessitates a huge enhancement of nematode resistance in cereal crops against CCNs. Nematode resistance in wheat and barley in combination with higher grain yields has been a preferential research area for cereal nematologists. This usually involved the targeted genetic exploitations through natural means of classical selection breeding of resistant genotypes and finding quantitative trait luci (QTLs) associated with resistance genes. These improvements were based on available genetic diversity among the crop plants. Recently, genome-wide association studies have widely been exploited to associate nematode resistance or susceptibility with particular regions of the genome. Use of biotechnological tools through the application of various transgenic strategies for enhancement of nematode resistance in various crop plants including wheat and barley had also been an important area of research. These modern approaches primarily include the use of gene silencing, exploitation of nematode effector genes, proteinase inhibitors, chemodisruptive peptides and a combination of one or more of these approaches. Furthermore, the perspective genome editing technologies including CRISPR-Cas9 could also be helpful for improving CCN resistance in wheat and barley. The information provided in this review will be helpful to enhance resistance against CCNs and will attract the attention of the scientific community towards this neglected area.

Effector gene birth in plant parasitic nematodes: Neofunctionalization of a housekeeping glutathione synthetase gene

Plant pathogens and parasites are a major threat to global food security. Plant parasitism has arisen four times independently within the phylum Nematoda, resulting in at least one parasite of every major food crop in the world. Some species within the most economically important order (Tylenchida) secrete proteins termed effectors into their host during infection to re-programme host development and immunity. The precise detail of how nematodes evolve new effectors is not clear. Here we reconstruct the evolutionary history of a novel effector gene family. We show that during the evolution of plant parasitism in the Tylenchida, the housekeeping glutathione synthetase (GS) gene was extensively replicated. New GS paralogues acquired multiple dorsal gland promoter elements, altered spatial expression to the secretory dorsal gland, altered temporal expression to primarily parasitic stages, and gained a signal peptide for secretion. The gene products are delivered into the host plant cell during infection, giving rise to "GS-like effectors". Remarkably, by solving the structure of GS-like effectors we show that during this process they have also diversified in biochemical activity, and likely represent the founding members of a novel class of GS-like enzyme. Our results demonstrate the re-purposing of an endogenous housekeeping gene to form a family of effectors with modified functions. We anticipate that our discovery will be a blueprint to understand the evolution of other plant-parasitic nematode effectors, and the foundation to uncover a novel enzymatic function.

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.

Characterization of Putative Effectors from the Cereal Cyst Nematode Heterodera avenae

Few molecular details of effectors of Heterodera avenae parasitism are known. We performed a high-throughput sequencing analysis of the H. avenae transcriptome at five developmental stages. A total of 82,549 unigenes were ultimately obtained, and 747 transcripts showed best hits to genes putatively encoding carbohydrate-active enzymes in plant-parasitic nematodes that play an important role in the invasion process. A total of 1,480 unigenes were homologous to known phytonematode effectors, and 63 putative novel effectors were identified in the H. avenae transcriptomes. Twenty-three unigenes were analyzed by qRT-PCR and confirmed to be highly expressed during at least one developmental stage. For in situ hybridization, 17 of the 22 tested putative effectors were specifically expressed and located in the subventral gland cells, and five putative novel effectors were specifically expressed in the dorsal gland. Furthermore, 115 transcripts were found to have putative lethal RNA interference (RNAi) phenotypes. Three target genes with lethal RNAi phenotypes and two of the four tested putative effectors were associated with a decrease in the number of cysts through in vitro RNAi technology. These transcriptomic data lay a foundation for further studies of interactions of H. avenae with cereal and H. avenae parasitic control.

Transcriptional responses of wheat and the cereal cyst nematode Heterodera avenae during their early contact stage

Cereal cyst nematode (Heterodera avenae) is attracted to and aggregated around wheat roots to initiate infection, but this interaction between wheat and the nematode is not fully understood. The transcriptional responses of both wheat and H. avenae were examined during their early contact stage by mRNA sequencing analysis; certain numbers of the differentially expressed genes (DEGs) were validated using quantitative real-time PCR. The immobile host wheat root only had 93 DEGs (27 up-regulated and 66 down-regulated), while the mobile plant parasitic nematode reacted much more actively with 879 DEGs (867 up-regulated and 12 down-regulated). Among them, a number of wheat DEGs (mostly down-regulated) were involved in biotic stress pathways, while several putative effector genes were up-regulated in the nematode DEGs. One putative chitinase-like effector gene of H. avenae was able to suppress BAX-triggered programmed cell death in Nicotiana benthamiana. Results of these experiments demonstrated that nematode responded more actively than wheat during the contact stage of parasitism. The parasite's responses mainly involved up-regulation of genes including at least one anti-plant-defence effector gene, whereas the host responses mainly involved down-regulation of certain defence-related genes.

Cereal Cyst Nematodes: A Complex and Destructive Group of Heterodera Species

Small grain cereals have served as the basis for staple foods, beverages, and animal feed for thousands of years. Wheat, barley, oats, rye, triticale, rice, and others are rich in calories, proteins, carbohydrates, vitamins, and minerals. These cereals supply 20% of the calories consumed by people worldwide and are therefore a primary source of energy for humans and play a vital role in global food and nutrition security. Global production of small grains increased linearly from 1960 to 2005, and then began to decline. Further decline in production is projected to continue through 2050 while global demand for these grains is projected to increase by 1% per annum. Currently, wheat, barley, and oat production exceeds consumption in developed countries, while in developing countries the consumption rate is higher than production. An increasing demand for meat and livestock products is likely to compound the demand for cereals in developing countries. Current production levels and trends will not be sufficient to fulfill the projected global demand generated by increased populations. For wheat, global production will need to be increased by 60% to fulfill the estimated demand in 2050. Until recently, global wheat production increased mostly in response to development of improved cultivars and farming practices and technologies. Production is now limited by biotic and abiotic constraints, including diseases, nematodes, insect pests, weeds, and climate. Among these constraints, plant-parasitic nematodes alone are estimated to reduce production of all world crops by 10%. Cereal cyst nematodes (CCNs) are among the most important nematode pests that limit production of small grain cereals. Heavily invaded young plants are stunted and their lower leaves are often chlorotic, forming pale green patches in the field. Mature plants are also stunted, have a reduced number of tillers, and the roots are shallow and have a "bushy-knotted" appearance. CCNs comprise a number of closely-related species and are found in most regions where cereals are produced.

Genome Evolution of Plant-Parasitic Nematodes

Plant parasitism has evolved independently on at least four separate occasions in the phylum Nematoda. The application of next-generation sequencing (NGS) to plant-parasitic nematodes has allowed a wide range of genome- or transcriptome-level comparisons, and these have identified genome adaptations that enable parasitism of plants. Current genome data suggest that horizontal gene transfer, gene family expansions, evolution of new genes that mediate interactions with the host, and parasitism-specific gene regulation are important adaptations that allow nematodes to parasitize plants. Sequencing of a larger number of nematode genomes, including plant parasites that show different modes of parasitism or that have evolved in currently unsampled clades, and using free-living taxa as comparators would allow more detailed analysis and a better understanding of the organization of key genes within the genomes. This would facilitate a more complete understanding of the way in which parasitism has shaped the genomes of plant-parasitic nematodes.

Heterodera schachtii Tyrosinase-like protein - a novel nematode effector modulating plant hormone homeostasis

The beet cyst nematode Heterodera schachtii causes major yield losses in sugar beet. Understanding the interaction between H. schachtii and its host plant is important for developing a sustainable management system. Nematode effectors play a crucial role in initializing and sustaining successful parasitism. In our study, we identified a gene (Hs-Tyr) encoding a tyrosinase functional domain (PF00264). We describe Hs-Tyr as a novel nematode effector. Hs-Tyr is localized in the nematode esophageal gland. Up-regulation of its expression coincided with the parasitic developmental stages of the nematode. Silencing Hs-Tyr by RNA interference made the treated nematodes less virulent. When RNAi-treated nematodes succeeded in infecting the plant, developing females and their associated syncytial nurse cells were significantly smaller than in control plants. Ectopically expressing the Hs-Tyr effector in Arabidopsis increased plant susceptibility to H. schachtii, but not to the root-knot nematode Meloidogyne incognita. Interestingly, Hs-Tyr in the plant promoted plant growth and changed the root architecture. Additionally, the expression of Hs-Tyr in Arabidopsis caused changes in the homeostasis of several plant hormones especially auxin and the ethylene precursor aminocyclopropane-carboxylic acid.

Comparative analysis of pre- and post-parasitic transcriptomes and mining pioneer effectors of Heterodera avenae

BACKGROUND

The cereal cyst nematode (CCN, Heterodera avenae) is a devastating pathogen of wheat and barley crops in many countries. We aimed to prioritize genetic and molecular targets for H. avenae control via the powerful and integrative bioinformatics platform.

RESULTS

Here, we sequenced mRNA isolated from Chinese H. avenae at pre-parasitic (consisting of egg, J1 and hatched-J2) stages and post-parasitic (consisting of parasitic-J2, J3, J4 and adults) stages. Total 1,066,719 reads of whole life cycle transcriptomes were assembled into 10,811 contigs with N50 length of 1754 bp and 71,401 singletons. Comparative analyses of orthologous among H. avenae and 7 other nematodes with various life-styles revealed the significance and peculiarity of neurological system for sedentary phytonematode. KEGG pathway enrichment demonstrated active crosstalk events of nervous system at pre-parasitic stages, and 6 FMRFamide-like neuropeptides were verified to display an expression peak at the hatched-J2 stage in H. avenae. Furthermore, multiple approaches were undertaken to mine putative effectors and parasitism-specific genes. Notably, H. avenae might represent the first phytonematode reported to possess the pioneer effectors with RxLR motif and potential effectors with homologies to Ant-5/Ant-34.

CONCLUSION

Our work provides valuable resources for in-depth understanding the parasitism and pathogenicity of H. avenae, as well as developing new targets-oriented strategies on effective managements.

The Distribution of Lectins across the Phylum Nematoda: A Genome-Wide Search

Nematodes are a very diverse phylum that has adapted to nearly every ecosystem. They have developed specialized lifestyles, dividing the phylum into free-living, animal, and plant parasitic species. Their sheer abundance in numbers and presence in nearly every ecosystem make them the most prevalent animals on earth. In this research nematode-specific profiles were designed to retrieve predicted lectin-like domains from the sequence data of nematode genomes and transcriptomes. Lectins are carbohydrate-binding proteins that play numerous roles inside and outside the cell depending on their sugar specificity and associated protein domains. The sugar-binding properties of the retrieved lectin-like proteins were predicted in silico. Although most research has focused on C-type lectin-like, galectin-like, and calreticulin-like proteins in nematodes, we show that the lectin-like repertoire in nematodes is far more diverse. We focused on C-type lectins, which are abundantly present in all investigated nematode species, but seem to be far more abundant in free-living species. Although C-type lectin-like proteins are omnipresent in nematodes, we have shown that only a small part possesses the residues that are thought to be essential for carbohydrate binding. Curiously, hevein, a typical plant lectin domain not reported in animals before, was found in some nematode species.

Large-Scale Identification and Characterization of Heterodera avenae Putative Effectors Suppressing or Inducing Cell Death in Nicotiana benthamiana

Heterodera avenae is one of the most important plant pathogens and causes vast losses in cereal crops. As a sedentary endoparasitic nematode, H. avenae secretes effectors that modify plant defenses and promote its biotrophic infection of its hosts. However, the number of effectors involved in the interaction between H. avenae and host defenses remains unclear. Here, we report the identification of putative effectors in H. avenae that regulate plant defenses on a large scale. Our results showed that 78 of the 95 putative effectors suppressed programmed cell death (PCD) triggered by BAX and that 7 of the putative effectors themselves caused cell death in Nicotiana benthamiana. Among the cell-death-inducing effectors, three were found to be dependent on their specific domains to trigger cell death and to be expressed in esophageal gland cells by in situ hybridization. Ten candidate effectors that suppressed BAX-triggered PCD also suppressed PCD triggered by the elicitor PsojNIP and at least one R-protein/cognate effector pair, suggesting that they are active in suppressing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). Notably, with the exception of isotig16060, these putative effectors could also suppress PCD triggered by cell-death-inducing effectors from H. avenae, indicating that those effectors may cooperate to promote nematode parasitism. Collectively, our results indicate that the majority of the tested effectors of H. avenae may play important roles in suppressing cell death induced by different elicitors in N. benthamiana.

Transcriptome Analysis of the Chrysanthemum Foliar Nematode, Aphelenchoides ritzemabosi (Aphelenchida: Aphelenchoididae)

The chrysanthemum foliar nematode (CFN), Aphelenchoides ritzemabosi, is a plant parasitic nematode that attacks many plants. In this study, a transcriptomes of mixed-stage population of CFN was sequenced on the Illumina HiSeq 2000 platform. 68.10 million Illumina high quality paired end reads were obtained which generated 26,817 transcripts with a mean length of 1,032 bp and an N50 of 1,672 bp, of which 16,467 transcripts were annotated against six databases. In total, 20,311 coding region sequences (CDS), 495 simple sequence repeats (SSRs) and 8,353 single-nucleotide polymorphisms (SNPs) were predicted, respectively. The CFN with the most shared sequences was B. xylophilus with 16,846 (62.82%) common transcripts and 10,543 (39.31%) CFN transcripts matched sequences of all of four plant parasitic nematodes compared. A total of 111 CFN transcripts were predicted as homologues of 7 types of carbohydrate-active enzymes (CAZymes) with plant/fungal cell wall-degrading activities, fewer transcripts were predicted as homologues of plant cell wall-degrading enzymes than fungal cell wall-degrading enzymes. The phylogenetic analysis of GH5, GH16, GH43 and GH45 proteins between CFN and other organisms showed CFN and other nematodes have a closer phylogenetic relationship. In the CFN transcriptome, sixteen types of genes orthologues with seven classes of protein families involved in the RNAi pathway in C. elegans were predicted. This research provides comprehensive gene expression information at the transcriptional level, which will facilitate the elucidation of the molecular mechanisms of CFN and the distribution of gene functions at the macro level, potentially revealing improved methods for controlling CFN.

Molecular Characterization of A Novel Effector Expansin-like Protein from Heterodera avenae that Induces Cell Death in Nicotiana benthamiana

Cereal cyst nematodes are sedentary biotrophic endoparasites that maintain a complex interaction with their host plants. Nematode effector proteins are synthesized in the oesophageal glands and are secreted into plant tissues through the stylet. To understand the function of nematode effectors in parasitic plants, we cloned predicted effectors genes from Heterodera avenae and transiently expressed them in Nicotiana benthamiana. Infiltration assays showed that HaEXPB2, a predicted expansin-like protein, caused cell death in N. benthamiana. In situ hybridization showed that HaEXPB2 transcripts were localised within the subventral gland cells of the pre-parasitic second-stage nematode. HaEXPB2 had the highest expression levels in parasitic second-stage juveniles. Subcellular localization assays revealed that HaEXPB2 could be localized in the plant cell wall after H. avenae infection.This The cell wall localization was likely affected by its N-terminal and C-terminal regions. In addition, we found that HaEXPB2 bound to cellulose and its carbohydrate-binding domain was required for this binding. The infectivity of H. avenae was significantly reduced when HaEXPB2 was knocked down by RNA interference in vitro. This study indicates that HaEXPB2 may play an important role in the parasitism of H. avenae through targeting the host cell wall.

Parasitic Nematode Immunomodulatory Strategies: Recent Advances and Perspectives

More than half of the described species of the phylum Nematoda are considered parasitic, making them one of the most successful groups of parasites. Nematodes are capable of inhabiting a wide variety of niches. A vast array of vertebrate animals, insects, and plants are all identified as potential hosts for nematode parasitization. To invade these hosts successfully, parasitic nematodes must be able to protect themselves from the efficiency and potency of the host immune system. Innate immunity comprises the first wave of the host immune response, and in vertebrate animals it leads to the induction of the adaptive immune response. Nematodes have evolved elegant strategies that allow them to evade, suppress, or modulate host immune responses in order to persist and spread in the host. Nematode immunomodulation involves the secretion of molecules that are capable of suppressing various aspects of the host immune response in order to promote nematode invasion. Immunomodulatory mechanisms can be identified in parasitic nematodes infecting insects, plants, and mammals and vary greatly in the specific tactics by which the parasites modify the host immune response. Nematode-derived immunomodulatory effects have also been shown to affect, negatively or positively, the outcome of some concurrent diseases suffered by the host. Understanding nematode immunomodulatory actions will potentially reveal novel targets that will in turn lead to the development of effective means for the control of destructive nematode parasites.

De Novo Analysis of the Transcriptome of Meloidogyne enterolobii to Uncover Potential Target Genes for Biological Control

Meloidogyne enterolobii is one of the obligate biotrophic root-knot nematodes that has the ability to reproduce on many economically-important crops. We carried out de novo sequencing of the transcriptome of M. enterolobii using Roche GS FLX and obtained 408,663 good quality reads that were assembled into 8193 contigs and 31,860 singletons. We compared the transcripts in different nematodes that were potential targets for biological control. These included the transcripts that putatively coded for CAZymes, kinases, neuropeptide genes and secretory proteins and those that were involved in the RNAi pathway and immune signaling. Typically, 75 non-membrane secretory proteins with signal peptides secreted from esophageal gland cells were identified as putative effectors, three of which were preliminarily examined using a PVX (pGR107)-based high-throughput transient plant expression system in Nicotiana benthamiana (N. benthamiana). Results showed that these candidate proteins suppressed the programmed cell death (PCD) triggered by the pro-apoptosis protein BAX, and one protein also caused necrosis, suggesting that they might suppress plant immune responses to promote pathogenicity. In conclusion, the current study provides comprehensive insight into the transcriptome of M. enterolobii for the first time and lays a foundation for further investigation and biological control strategies.

Advances in Understanding the Molecular Mechanisms of Root Lesion Nematode Host Interactions

Root lesion nematodes (RLNs) are one of the most economically important groups of plant nematodes. As migratory endoparasites, their presence in roots is less obvious than infestations of sedentary endoparasites; nevertheless, in many instances, they are the major crop pests. With increasing molecular information on nematode parasitism, available data now reflect the differences and, in particular, similarities in lifestyle between migratory and sedentary endoparasites. Far from being unsophisticated compared with sedentary endoparasites, migratory endoparasites are exquisitely suited to their parasitic lifestyle. What they lack in effectors required for induction of permanent feeding sites, they make up for with their versatile host range and their ability to move and feed from new host roots and survive adverse conditions. In this review, we summarize the current molecular data available for RLNs and highlight differences and similarities in effectors and molecular mechanisms between migratory and sedentary endoparasitic nematodes.

Characterization and analysis of a de novo transcriptome from the pygmy grasshopper Tetrix japonica

The pygmy grasshopper Tetrix japonica is a common insect distributed throughout the world, and it has the potential for use in studies of body colour polymorphism, genomics and the biology of Tetrigoidea (Insecta: Orthoptera). However, limited biological information is available for this insect. Here, we conducted a de novo transcriptome study of adult and larval T. japonica to provide a better understanding of its gene expression and develop genomic resources for future work. We sequenced and explored the characteristics of the de novo transcriptome of T. japonica using Illumina HiSeq 2000 platform. A total of 107 608 206 paired-end clean reads were assembled into 61 141 unigenes using the trinity software; the mean unigene size was 771 bp, and the N50 length was 1238 bp. A total of 29 225 unigenes were functionally annotated to the NCBI nonredundant protein sequences (Nr), NCBI nonredundant nucleotide sequences (Nt), a manually annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of putative genes that are potentially involved in pigment pathways, juvenile hormone (JH) metabolism and signalling pathways were identified in the T. japonica transcriptome. Additionally, 165 769 and 156 796 putative single nucleotide polymorphisms occurred in the adult and larvae transcriptomes, respectively, and a total of 3162 simple sequence repeats were detected in this assembly. This comprehensive transcriptomic data for T. japonica will provide a usable resource for gene predictions, signalling pathway investigations and molecular marker development for this species and other pygmy grasshoppers.

De novo analysis of the transcriptome of Pratylenchus zeae to identify transcripts for proteins required for structural integrity, sensation, locomotion and parasitism

The root lesion nematode Pratylenchus zeae, a migratory endoparasite, is an economically important pest of major crop plants (e.g. cereals, sugarcane). It enters host roots, migrates through root tissues and feeds from cortical cells, and defends itself against biotic and abiotic stresses in the soil and in host tissues. We report de novo sequencing of the P. zeae transcriptome using 454 FLX, and the identification of putative transcripts encoding proteins required for movement, response to stimuli, feeding and parasitism. Sequencing generated 347,443 good quality reads which were assembled into 10,163 contigs and 139,104 singletons: 65% of contigs and 28% of singletons matched sequences of free-living and parasitic nematodes. Three-quarters of the annotated transcripts were common to reference nematodes, mainly representing genes encoding proteins for structural integrity and fundamental biochemical processes. Over 15,000 transcripts were similar to Caenorhabditis elegans genes encoding proteins with roles in mechanical and neural control of movement, responses to chemicals, mechanical and thermal stresses. Notably, 766 transcripts matched parasitism genes employed by both migratory and sedentary endoparasites in host interactions, three of which hybridized to the gland cell region, suggesting that they might be secreted. Conversely, transcripts for effectors reported to be involved in feeding site formation by sedentary endoparasites were conspicuously absent. Transcripts similar to those encoding some secretory-excretory products at the host interface of Brugia malayi, the secretome of Meloidogyne incognita and products of gland cells of Heterodera glycines were also identified. This P. zeae transcriptome provides new information for genome annotation and functional analysis of possible targets for control of pratylenchid nematodes.

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.

De novo transcriptomic analysis of the female and male adults of the blood fluke Schistosoma turkestanicum

Background

Schistosoma turkestanicum is a parasite of considerable veterinary importance as an agent of animal schistosomiasis in many countries, including China. The S. turkestanicum cercariae can also infect humans, causing cercarial dermatitis in many countries and regions of the world. In spite of its significance as a pathogen of animals and humans, there is little transcriptomic and genomic data in the public databases.

Methods

Herein, we performed the transcriptome Illumina RNA sequencing (RNA-seq) of adult males and females of S. turkestanicum and de novo transcriptome assembly.

Results

Approximately 81.1 (female) and 80.5 (male) million high-quality clean reads were obtained and then 29,526 (female) and 41,346 (male) unigenes were assembled. A total of 34,624 unigenes were produced from S. turkestanicum females and males, with an average length of 878 nucleotides (nt) and N50 of 1480 nt. Of these unigenes, 25,158 (72.7 %) were annotated by blast searches against the NCBI non-redundant protein database. Among these, 21,995 (63.5 %), 22,189 (64.1 %) and 13,754 (39.7 %) of the unigenes had significant similarity in the NCBI non-redundant protein (NR), non-redundant nucleotide (NT) and Swiss-Prot databases, respectively. In addition, 3150 unigenes were identified to be expressed specifically in females and 1014 unigenes were identified to be expressed specifically in males. Interestingly, several pathways associated with gonadal development and sex maintenance were found, including the Wnt signaling pathway (103; 2 %) and progesterone-mediated oocyte maturation (77; 1.5 %).

Conclusions

The present study characterized and compared the transcriptomes of adult female and male blood fluke, S. turkestanicum. These results will not only serve as valuable resources for future functional genomics studies to understand the molecular aspects of S. turkestanicum, but also will provide essential information for ongoing whole genome sequencing efforts on this pathogenic blood fluke.

A transcriptomic insight into the infective juvenile stage of the insect parasitic nematode, Heterorhabditis indica

Background

Nematodes are the most numerous animals in the soil. Insect parasitic nematodes of the genus Heterorhabditis are capable of selectively seeking, infecting and killing their insect-hosts in the soil. The infective juvenile (IJ) stage of the Heterorhabditis nematodes is analogous to Caenorhabditis elegans dauer juvenile stage, which remains in ‘arrested development’ till it finds and infects a new insect-host in the soil. H. indica is the most prevalent species of Heterorhabditis in India. To understand the genes and molecular processes that govern the biology of the IJ stage, and to create a resource to facilitate functional genomics and genetic exploration, we sequenced the transcriptome of H. indica IJs.

Results

The de-novo sequence assembly using Velvet-Oases pipeline resulted in 13,593 unique transcripts at N50 of 1,371 bp, of which 53 % were annotated by blastx. H. indica transcripts showed higher orthology with parasitic nematodes as compared to free living nematodes. In-silico expression analysis showed 30 % of transcripts expressing with ≥100 FPKM value. All the four canonical dauer formation pathways like cGMP-PKG, insulin, dafachronic acid and TGF-β were active in the IJ stage. Several other signaling pathways were highly represented in the transcriptome. Twenty-four orthologs of C. elegans RNAi pathway effector genes were discovered in H. indica, including nrde-3 that is reported for the first time in any of the parasitic nematodes. An ortholog of C. elegans tol-1 was also identified. Further, 272 kinases belonging to 137 groups, and several previously unidentified members of important gene classes were identified.

Conclusions

We generated high-quality transcriptome sequence data from H. indica IJs for the first time. The transcripts showed high similarity with the parasitic nematodes, M. hapla, and A. suum as opposed to C. elegans, a species to which H. indica is more closely related. The high representation of transcripts from several signaling pathways in the IJs indicates that despite being a developmentally arrested stage; IJs are a hotbed of signaling and are actively interacting with their environment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2510-z) contains supplementary material, which is available to authorized users.

Analysis of the Transcriptome of the Infective Stage of the Beet Cyst Nematode, H. schachtii

The beet cyst nematode, Heterodera schachtii, is a major root pest that significantly impacts the yield of sugar beet, brassicas and related species. There has been limited molecular characterisation of this important plant pathogen: to identify target genes for its control the transcriptome of the pre-parasitic J2 stage of H. schachtii was sequenced using Roche GS FLX. Ninety seven percent of reads (i.e., 387,668) with an average PHRED score > 22 were assembled with CAP3 and CLC Genomics Workbench into 37,345 and 47,263 contigs, respectively. The transcripts were annotated by comparing with gene and genomic sequences of other nematodes and annotated proteins on public databases. The annotated transcripts were much more similar to sequences of Heterodera glycines than to those of Globodera pallida and root knot nematodes (Meloidogyne spp.). Analysis of these transcripts showed that a subset of 2,918 transcripts was common to free-living and plant parasitic nematodes suggesting that this subset is involved in general nematode metabolism and development. A set of 148 contigs and 183 singletons encoding putative homologues of effectors previously characterised for plant parasitic nematodes were also identified: these are known to be important for parasitism of host plants during migration through tissues or feeding from cells or are thought to be involved in evasion or modulation of host defences. In addition, the presence of sequences from a nematode virus is suggested. The sequencing and annotation of this transcriptome significantly adds to the genetic data available for H. schachtii, and identifies genes primed to undertake required roles in the critical pre-parasitic and early post-parasitic J2 stages. These data provide new information for identifying potential gene targets for future protection of susceptible crops against H. schachtii.

Assembly, Assessment, and Availability of De novo Generated Eukaryotic Transcriptomes

De novo assembly of a complete transcriptome without the need for a guiding reference genome is attractive, particularly where the cost and complexity of generating a eukaryote genome is prohibitive. The transcriptome should not however be seen as just a quick and cheap alternative to building a complete genome. Transcriptomics allows the understanding and comparison of spatial and temporal samples within an organism, and allows surveying of multiple individuals or closely related species. De novo assembly in theory allows the building of a complete transcriptome without any prior knowledge of the genome. It also allows the discovery of alternate splice forms of coding RNAs and also non-coding RNAs, which are often missed by proteomic approaches, or are incompletely annotated in genome studies. The limitations of the method are that the generation of a truly complete assembly is unlikely, and so we require some methods for the assessment of the quality and appropriateness of a generated transcriptome. Whilst no single consensus pipeline or tool is agreed as optimal, various algorithms, and easy to use software do exist making transcriptome generation a more common approach. With this expansion of data, questions still exist relating to how do we make these datasets fully discoverable, comparable and most useful to understand complex biological systems?

SPRYSEC Effectors: A Versatile Protein-Binding Platform to Disrupt Plant Innate Immunity

Persistent infections by sedentary plant-parasitic nematodes are a major threat to important food crops all over the world. These roundworms manipulate host plant cell morphology and physiology to establish sophisticated feeding structures. Key modifications to plant cells during their transition into feeding structures are largely attributed to the activity of effectors secreted by the nematodes. The SPRYSEC effectors were initially identified in the potato cyst nematodes Globodera rostochiensis and G. pallida, and are characterized by a single SPRY domain, a non-catalytic domain present in modular proteins with different functions. The SPRY domain is wide-spread among eukaryotes and thought to be involved in mediating protein-protein interactions. Thus far, the SPRY domain is only reported as a functional domain in effectors of plant-parasitic nematodes, but not of other plant pathogens. SPRYSEC effectors have been implicated in both suppression and activation of plant immunity, but other possible roles in nematode virulence remain undefined. Here, we review the latest reports on the structure, function, and sequence diversity of SPRYSEC effectors, which provide support for a model featuring these effectors as a versatile protein-binding platform for the nematodes to target a wide range of host proteins during parasitism.

SPRYSEC Effectors: A Versatile Protein-Binding Platform to Disrupt Plant Innate Immunity

Persistent infections by sedentary plant-parasitic nematodes are a major threat to important food crops all over the world. These roundworms manipulate host plant cell morphology and physiology to establish sophisticated feeding structures. Key modifications to plant cells during their transition into feeding structures are largely attributed to the activity of effectors secreted by the nematodes. The SPRYSEC effectors were initially identified in the potato cyst nematodes Globodera rostochiensis and G. pallida, and are characterized by a single SPRY domain, a non-catalytic domain present in modular proteins with different functions. The SPRY domain is wide-spread among eukaryotes and thought to be involved in mediating protein-protein interactions. Thus far, the SPRY domain is only reported as a functional domain in effectors of plant-parasitic nematodes, but not of other plant pathogens. SPRYSEC effectors have been implicated in both suppression and activation of plant immunity, but other possible roles in nematode virulence remain undefined. Here, we review the latest reports on the structure, function, and sequence diversity of SPRYSEC effectors, which provide support for a model featuring these effectors as a versatile protein-binding platform for the nematodes to target a wide range of host proteins during parasitism.

RNA-Seq Based Identification of Candidate Parasitism Genes of Cereal Cyst Nematode (Heterodera avenae) during Incompatible Infection to Aegilops variabilis

One of the reasons for the progressive yield decline observed in cereals production is the rapid build-up of populations of the cereal cyst nematode (CCN, Heterodera avenae). These nematodes secrete so-call effectors into their host plant to suppress the plant defense responses, alter plant signaling pathways and then induce the formation of syncytium after infection. However, little is known about its molecular mechanism and parasitism during incompatible infection. To gain insight into its repertoire of parasitism genes, we investigated the transcriptome of the early parasitic second-stage (30 hours, 3 days and 9 days post infection) juveniles of the CCN as well as the CCN infected tissue of the host Aegilops variabilis by Illumina sequencing. Among all assembled unigenes, 681 putative genes of parasitic nematode were found, in which 56 putative effectors were identified, including novel pioneer genes and genes corresponding to previously reported effectors. All the 681 CCN unigenes were mapped to 229 GO terms and 200 KEGG pathways, including growth, development and several stimulus-related signaling pathways. Sixteen clusters were involved in the CCN unigene expression atlas at the early stages during infection process, and three of which were significantly gene-enriched. Besides, the protein-protein interaction network analysis revealed 35 node unigenes which may play an important role in the plant-CCN interaction. Moreover, in a comparison of differentially expressed genes between the pre-parasitic juveniles and the early parasitic juveniles, we found that hydrolase activity was up-regulated in pre J2s whereas binding activity was upregulated in infective J2s. RT-qPCR analysis on some selected genes showed detectable expression, indicating possible secretion of the proteins and putative role in infection. This study provided better insights into the incompatible interaction between H. avenae and the host plant Ae. varabilis. Moreover, RNAi targets with potential lethality were screened out and primarily validated, which provide candidates for engineering-based control of cereal cyst nematode in crops breeding.

Large-scale identification of wheat genes resistant to cereal cyst nematode Heterodera avenae using comparative transcriptomic analysis

Background

Cereal cyst nematode Heterodera avenae, an important soil-borne pathogen in wheat, causes numerous annual yield losses worldwide, and use of resistant cultivars is the best strategy for control. However, target genes are not readily available for breeding resistant cultivars. Therefore, comparative transcriptomic analyses were performed to identify more applicable resistance genes for cultivar breeding.

Methods

The developing nematodes within roots were stained with acid fuchsin solution. Transcriptome assemblies and redundancy filteration were obtained by Trinity, TGI Clustering Tool and BLASTN, respectively. Gene Ontology annotation was yielded by Blast2GO program, and metabolic pathways of transcripts were analyzed by Path_finder. The ROS levels were determined by luminol-chemiluminescence assay. The transcriptional gene expression profiles were obtained by quantitative RT-PCR.

Results

The RNA-sequencing was performed using an incompatible wheat cultivar VP1620 and a compatible control cultivar WEN19 infected with H. avenae at 24 h, 3 d and 8 d. Infection assays showed that VP1620 failed to block penetration of H. avenae but disturbed the transition of developmental stages, leading to a significant reduction in cyst formation. Two types of expression profiles were established to predict candidate resistance genes after developing a novel strategy to generate clean RNA-seq data by removing the transcripts of H. avenae within the raw data before assembly. Using the uncoordinated expression profiles with transcript abundance as a standard, 424 candidate resistance genes were identified, including 302 overlapping genes and 122 VP1620-specific genes. Genes with similar expression patterns were further classified according to the scales of changed transcript abundances, and 182 genes were rescued as supplementary candidate resistance genes. Functional characterizations revealed that diverse defense-related pathways were responsible for wheat resistance against H. avenae. Moreover, phospholipase was involved in many defense-related pathways and localized in the connection position. Furthermore, strong bursts of reactive oxygen species (ROS) within VP1620 roots infected with H. avenae were induced at 24 h and 3 d, and eight ROS-producing genes were significantly upregulated, including three class III peroxidase and five lipoxygenase genes.

Conclusions

Large-scale identification of wheat resistance genes were processed by comparative transcriptomic analysis. Functional characterization showed that phospholipases associated with ROS production played vital roles in early defense responses to H. avenae via involvement in diverse defense-related pathways as a hub switch. This study is the first to investigate the early defense responses of wheat against H. avenae, not only provides applicable candidate resistance genes for breeding novel wheat cultivars, but also enables a better understanding of the defense mechanisms of wheat against H. avenae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2037-8) contains supplementary material, which is available to authorized users.

Ameliorated de novo transcriptome assembly using Illumina paired end sequence data with Trinity Assembler

Advent of Next Generation Sequencing has led to possibilities of de novo transcriptome assembly of organisms without availability of complete genome sequence. Among various sequencing platforms available, Illumina is the most widely used platform based on data quality, quantity and cost. Various de novo transcriptome assemblers are also available today for construction of de novo transcriptome. In this study, we aimed at obtaining an ameliorated de novo transcriptome assembly with sequence reads obtained from Illumina platform and assembled using Trinity Assembler. We found that, primary transcriptome assembly obtained as a result of Trinity can be ameliorated on the basis of transcript length, coverage, and depth and protein homology. Our approach to ameliorate is reproducible and could enhance the sensitivity and specificity of the assembled transcriptome which could be critical for validation of the assembled transcripts and for planning various downstream biological assays.

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements

BACKGROUND

The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution.

RESULTS

The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins.

CONCLUSIONS

Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

An ANNEXIN-like protein from the cereal cyst nematode Heterodera avenae suppresses plant defense

Parasitism genes encoding secreted effector proteins of plant-parasitic nematodes play important roles in facilitating parasitism. An annexin-like gene was isolated from the cereal cyst nematode Heterodera avenae (termed Ha-annexin) and had high similarity to annexin 2, which encodes a secreted protein of Globodera pallida. Ha-annexin encodes a predicted 326 amino acid protein containing four conserved annexin domains. Southern blotting revealed that there are at least two homologies in the H. avenae genome. Ha-annexin transcripts were expressed within the subventral gland cells of the pre-parasitic second-stage juveniles by in situ hybridization. Additionally, expression of these transcripts were relatively higher in the parasitic second-stage juveniles by quantitative real-time RT-PCR analysis, coinciding with the time when feeding cell formation is initiated. Knockdown of Ha-annexin by method of barley stripe mosaic virus-based host-induced gene silencing (BSMV-HIGS) caused impaired nematode infections at 7 dpi and reduced females at 40 dpi, indicating important roles of the gene in parasitism at least in early stage in vivo. Transiently expression of Ha-ANNEXIN in onion epidermal cells and Nicotiana benthamiana leaf cells showed the whole cell-localization. Using transient expression assays in N. benthamiana, we found that Ha-ANNEXIN could suppress programmed cell death triggered by the pro-apoptotic mouse protein BAX and the induction of marker genes of PAMP-triggered immunity (PTI) in N. benthamiana. In addition, Ha-ANNEXIN targeted a point in the mitogen-activated protein kinase (MAPK) signaling pathway downstream of two kinases MKK1 and NPK1 in N. benthamiana.