The cathepsin S cysteine proteinase of the burrowing nematode Radopholus similis is essential for the reproduction and invasion

Advances in Migratory Plant Endoparasitic Nematode Effectors

Unlike sedentary plant-parasitic nematodes, migratory plant endoparasitic nematodes (MPENs) are unable to establish permanent feeding sites, and all developmental stages (except eggs) can invade and feed on plant tissues and can be easily overlooked because of the unspecific symptoms. They cause numerous economic losses in agriculture, forestry, and horticulture. In order to understand the pathogenetic mechanism of MPENs, here we describe research on functions and host targets focused on currently identified effectors from six MPENs, namely Radopholus similis, Pratylenchus spp., Ditylenchus destructor, Bursaphelenchus xylophilus, Aphelenchoides besseyi, and Hirschmanniella oryzae. This information will provide valuable insights into understanding MPEN effectors and for future fostering advancements in plant protection.

Single Nematode Transcriptomic Analysis, Using Long-Read Technology, Reveals Two Novel Virulence Gene Candidates in the Soybean Cyst Nematode, Heterodera glycines

The soybean cyst nematode (Heterodera glycines, SCN), is the most damaging disease of soybean in North America. While management of this pest using resistant soybean is generally still effective, prolonged exposure to cultivars derived from the same source of resistance (PI 88788) has led to the emergence of virulence. Currently, the underlying mechanisms responsible for resistance breakdown remain unknown. In this study, we combined a single nematode transcriptomic profiling approach with long-read sequencing to reannotate the SCN genome. This resulted in the annotation of 1932 novel transcripts and 281 novel gene features. Using a transcript-level quantification approach, we identified eight novel effector candidates overexpressed in PI 88788 virulent nematodes in the late infection stage. Among these were the novel gene Hg-CPZ-1 and a pioneer effector transcript generated through the alternative splicing of the non-effector gene Hetgly21698. While our results demonstrate that alternative splicing in effectors does occur, we found limited evidence of direct involvement in the breakdown of resistance. However, our analysis highlighted a distinct pattern of effector upregulation in response to PI 88788 resistance indicative of a possible adaptation process by SCN to host resistance.

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.

Comparative transcriptomic analysis of candidate effectors to explore the infection and survival strategy of Bursaphelenchus xylophilus during different interaction stages with pine trees

BACKGROUND

The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages.

RESULTS

In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana.

CONCLUSIONS

The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors.

A Venom Allergen-Like Protein, RsVAP, the First Discovered Effector Protein of Radopholus similis That Inhibits Plant Defense and Facilitates Parasitism

Radopholus similis is a migratory endoparasitic nematode that is extremely harmful to host plants. Venom allergen-like proteins (VAPs) are members of the cysteine-rich secretory protein family that are widely present in plants and animals. In this study, we cloned a VAP gene from R. similis, designated as RsVAP. RsVAP contains an open reading frame of 1089 bp encoding 362 amino acids. RsVAP is specifically expressed in the esophageal gland, and the expression levels of RsVAP are significantly higher in juveniles than in other life stages of R. similis. This expression pattern of RsVAP was consistent with the biological characteristics of juveniles of R. similis, which have the ability of infection and are the main infection stages of R. similis. The pathogenicity and reproduction rate of R. similis in tomato was significantly attenuated after RsVAP was silenced. In tobacco leaves transiently expressing RsVAP, the pathogen-associated molecular pattern-triggered immunity (PTI) induced by a bacterial flagellin fragment (flg22) was inhibited, while the cell death induced by two sets of immune elicitors (BAX and Gpa2/RBP-1) was repressed. The RsVAP-interacting, ras-related protein RABA1d (LeRabA1d) was identified in tomato hosts by yeast two-hybrid and co-immunoprecipitation assays. RsVAP may interact with LeRabA1d to affect the host defense response, which in turn facilitates nematode infection. This study provides the first evidence for the inhibition of plant defense response by a VAP from migratory plant-parasitic nematodes, and, for the first time, the target protein of R. similis in its host was identified.

Current Insights into Migratory Endoparasitism: Deciphering the Biology, Parasitism Mechanisms, and Management Strategies of Key Migratory Endoparasitic Phytonematodes

Despite their physiological differences, sedentary and migratory plant-parasitic nematodes (PPNs) share several commonalities. Functional characterization studies of key effectors and their targets identified in sedentary phytonematodes are broadly applied to migratory PPNs, generalizing parasitism mechanisms existing in distinct lifestyles. Despite their economic significance, host-pathogen interaction studies of migratory endoparasitic nematodes are limited; they have received little attention when compared to their sedentary counterparts. Because several migratory PPNs form disease complexes with other plant-pathogens, it is important to understand multiple factors regulating their feeding behavior and lifecycle. Here, we provide current insights into the biology, parasitism mechanism, and management strategies of the four-key migratory endoparasitic PPN genera, namely Pratylenchus, Radopholus, Ditylenchus, and Bursaphelenchus. Although this review focuses on these four genera, many facets of feeding mechanisms and management are common across all migratory PPNs and hence can be applied across a broad genera of migratory phytonematodes.

Recombinant expression, characterization and phylogenetic studies of novels cystatins-like proteins of sweet orange (Citrus sinensis) and clementine (Citrus clementina)

Phytocystatins are plant cystatins that are related to several physiological processes regulating endogenous cysteine proteases involved in seed development and germination, programmed cell death and response to stress conditions. In addition, phytocystatins can act in plant defense against exogenous peptidases from herbivorous insects, pathogens and nematodes. Considering that Citrus fruits are important to human nutrition and represent a high value crop in worldwide agriculture, in the present work, we performed the identification of putative cystatins from Citrus sinensis and from Citrus clementine and submitted them to phylogenetic analysis. Six cystatins from each species were identified as orthologous and classified into three well supported phylogenetic groups. Five cystatins representative of the phylogenetic groups were recombinantly expressed and the in vitro studies revealed them to be potent inhibitors against the cysteine peptidases papain, legumain, human cathepsins (B, L, S, K) and a cathepsin B-like from Diaphorina citri (the Asian Citrus psyllid). Our findings provide the C. clementina and C. sinensis cystatins classification and an enzyme-inhibitor interactions profile, which may reflect an evolutionary process of Citrus cystatins related to gene functions as initial germination rates and seedlings development as well associated to plant defense against pathogens, as insects and nematodes.

Identification and Characterization of a Fatty Acid- and Retinoid-Binding Protein Gene (Ar-far-1) from the Chrysanthemum Foliar Nematode, Aphelenchoides ritzemabosi

The chrysanthemum foliar nematode (CFN), Aphelenchoides ritzemabosi, is a migratory, plant-parasitic nematode that is widely distributed and infects the aboveground parts of many plants. The fatty acid- and retinoid-binding proteins (FAR) are nematode-specific proteins that are involved in the development, reproduction, and infection of nematodes and are secreted into the tissues to disrupt the plant defense reaction. In this study, we obtained the full-length sequence of the FAR gene (Ar-far-1) from CFN, which is 727 bp and includes a 546 bp ORF that encodes 181 amino acids. Ar-FAR-1 from CFN has the highest sequence similarity to Ab-FAR-1 from A. besseyi, and they are located within the same branch of the phylogenetic tree. Fluorescence-based ligand-binding analysis confirmed that recombinant Ar-FAR-1 was bound to fatty acids and retinol. Ar-far-1 mRNA was expressed in the muscle layer, intestine, female genital system, and egg of CFN, and more highly expressed in females than in males among the four developmental stages of CFN. We demonstrated that the reproduction number and infection capacity of CFN decreased significantly when Ar-far-1 was effectively silenced by in vitro RNAi. Ar-far-1 plays an important role in the development, reproduction, infectivity, and pathogenesis of CFN and may be used as an effective target gene for the control of CFN. The results provide meaningful data about the parasitic and pathogenic genes of CFN to study the interaction mechanism between plant-parasitic nematodes and hosts.

The genome of the migratory nematode, Radopholus similis, reveals signatures of close association to the sedentary cyst nematodes

Radopholus similis, commonly known as the burrowing nematode, is an important pest of myriad crops and ornamentals including banana (Musa spp.) and Citrus spp. In order to characterize the potential role of putative effectors encoded by R. similis genes we compared predicted proteins from a draft R. similis genome with other plant-parasitic nematodes in order to define the suite of excreted/secreted proteins that enable it to function as a parasite and to ascertain the phylogenetic position of R. similis in the Tylenchida order. Identification and analysis of candidate genes encoding for key plant cell-wall degrading enzymes including GH5 cellulases, PL3 pectate lyases and GH28 polygalactouranase revealed a pattern of occurrence similar to other PPNs, although with closest phylogenetic associations to the sedentary cyst nematodes. We also observed the absence of a suite of effectors essential for feeding site formation in the cyst nematodes. Clustering of various orthologous genes shared by R. similis with other nematodes showed higher overlap with the cyst nematodes than with the root-knot or other migratory endoparasitic nematodes. The data presented here support the hypothesis that R. similis is evolutionarily closer to the cyst nematodes, however, differences in the effector repertoire delineate ancient divergence of parasitism, probably as a consequence of niche specialization. These similarities and differences further underscore distinct evolutionary relationships during the evolution of parasitism in this group of nematodes.

Transcriptome analysis of Globodera pallida from the susceptible host Solanum tuberosum or the resistant plant Solanum sisymbriifolium

A transcriptome analysis of G. pallida juveniles collected from S. tuberosum or S. sisymbriifolium 24 h post infestation was performed to provide insights into the parasitic process of this nematode. A total of 41 G. pallida genes were found to be significantly differentially expressed when parasitizing the two plant species. Among this set, 12 were overexpressed when G. pallida was parasitizing S. tuberosum and 29 were overexpressed when parasitizing S. sisymbriifolium. Out of the 12 genes, three code for secretory proteins; one is homologous to effector gene Rbp-4, the second is an uncharacterized protein with a signal peptide sequence, and the third is an ortholog of a Globodera rostochiensis effector belonging to the 1106 effector family. Other overexpressed genes from G. pallida when parasitizing S. tuberosum were either unknown, associated with a stress or defense response, or associated with sex differentiation. Effector genes namely Eng-1, Cathepsin S-like cysteine protease, cellulase, and two unknown genes with secretory characteristics were over expressed when G. pallida was parasitizing S. sisymbriifolium relative to expression from S. tuberosum. Our findings provide insight into gene regulation of G. pallida while infecting either the trap crop S. sisymbriifolium or the susceptible host, S. tuberosum.

Screening of reference genes in real-time PCR for Radopholus similis

Six candidate reference genes were chosen from the transcriptome database of Radopholus similis using the bioinformatics method, including four conventional reference genes (actin, Eukaryotic translation initiation factor 5A (eIF5A), Tubulin alpha (a-tubulin), ubiquitin (UBI)) and two new candidate reference genes (Ribosomal protein S21 (Rps21) and Serine/threonine protein phosphatase PP1-β catalytic subunit (β-PP1)). In addition, a traditional reference gene 18S ribosomal RNA (18S rRNA) obtained from NCBI databases was also added to the analysis. Real-time PCR was used to detect the expression of seven candidate reference genes in six populations of R. similis and four developmental stages (female, male, larva and egg) of a population. The stability of the expression of candidate genes was evaluated by three software programs, BestKeeper, geNorm and NormFinder. The results showed that eIF5A is the most suitable reference gene for gene functional research of different populations, while both Rps21 and eIF5A are the most suitable reference genes for different developmental stages of a population. Therefore, eIF5A is the best reference gene for studying R. similis. However, one defect of this study is that only seven candidate reference genes were analyzed; ideally, more genes should be tested.

Molecular identification and functional characterization of the cathepsin B gene (Ab-cb-1) in the plant parasitic nematode Aphelenchoides besseyi

The rice white tip nematode, Aphelenchoides besseyi, is widely distributed in rice planting areas worldwide and causes serious economic losses. Cathepsin genes have been demonstrated to have importance in studying the reproduction, development, pathogenicity, and control methods of plant nematodes. In this paper, a novel cathepsin B gene, Ab-cb-1, was found and cloned. The Ab-cb-1 gene was 1347 bp in length and encodes 369 amino acids. The Ab-CB-1 protein contains characteristic occluding loops but no signal peptide. A homology analysis showed that Ab-CB-1 had the highest identity value (64%) to the known amino acid sequence of cathepsin B-like cysteine protease 6 from Toxocara canis. When Ab-cb-1 was expressed in a prokaryotic system, the protein massed approximately 45 kDa and could decompose carrot callus. Ab-cb-1 mRNA was localized in the nematode intestine. The relative expression level of Ab-cb-1 in the A. besseyi Ab-S24 population, which had high reproductivity, was approximately 6.9 times that in the Ab-N10 population, which had low reproductivity, and the difference was significant (p<0.05). The Ab-cb-1 expression level was highest in females; the expression levels in males, juveniles and eggs were 30%, 12.2% and 5% of that in females, respectively, and the differences were significant among all four stages (p<0.05). Nematodes of the Ab-S24 population were treated with Ab-cb-1 dsRNA for 12 h, 24 h, 36 h and 48 h, and their reproduction decreased with increasing time. These results demonstrated that Ab-CB-1 was a digestive enzyme with hydrolytic protease properties and that Ab-cb-1 played an important role in the reproduction of A. besseyi.

Transgenic Nicotiana benthamiana plants expressing a hairpin RNAi construct of a nematode Rs-cps gene exhibit enhanced resistance to Radopholus similis

Burrowing nematodes (Radopholus similis) cause severe harm in many agronomic and horticultural crops and are very difficult to manage. Cathepsin S is one of the most important cysteine proteinases and plays key roles in nematodes and many other parasites. To evaluate the effect of in planta RNAi on the control of this nematode, a specific fragment from the protease gene, cathepsin S (Rs-cps), was cloned into the binary vector pFGC5941 in the forward and reverse orientations to construct recombinant plant RNAi vectors. Transgenic Nicotiana benthamiana plants expressing Rs-cps dsRNA were obtained and studied. The transcript abundance of Rs-cps dsRNA appeared to be diverse in the different transgenic lines. Moreover, the bioassay results revealed that Rs-cps transgenic N. benthamiana plants were resistant to R. similis and the transcription level of Rs-cps in R. similis was drastically decreased. In addition, the reproduction and hatching rate of R. similis isolated from the Rs-cps transgenic plants were also significantly reduced. Our results suggest that Rs-cps is essential for the reproduction and pathogenicity of R. similis. This is the first study to employ in planta RNAi approach to target the Rs-cps gene for the control of plant parasitic nematodes.

Molecular Cloning, Recombinant Expression and Antifungal Activity of BnCPI, a Cystatin in Ramie (Boehmeria nivea L.)

Phytocystatins play multiple roles in plant growth, development and resistance to pests and other environmental stresses. A ramie (Boehmeria nivea L.) phytocystatin gene, designated as BnCPI, was isolated from a ramie cDNA library and its full-length cDNA was obtained by rapid amplification of cDNA ends (RACE). The full-length cDNA sequence (691 bp) consisted of a 303 bp open reading frame (ORF) encoding a protein of 100 amino acids with deduced molecular mass of 11.06 kDa and a theoretical isoelectric point (pI) of 6.0. The alignment of genome DNA (accession No. MF153097) and cDNA sequences of BnCPI showed that an intron (~104 bp) exists in the coding region. The BnCPI protein contains most of the highly conserved blocks including Gly⁵-Gly⁶ at the N-terminal, the reactive site motif QxVxG (Q⁴⁹V⁵⁰V⁵¹S⁵²G⁵³), the L⁷⁹-W⁸⁰ block and the [LVI]-[AGT]-[RKE]-[FY]-[AS]-[VI]-x-[EDQV]-[HYFQ]-N (L²²G²³R²⁴ F²⁵A²⁶V²⁷ D²⁸D²⁹H³⁰ N³¹) block that is common among plant cystatins. BLAST analysis indicated that BnCPI is similar to cystatins from Glycine max (77%), Glycine soja (76%), Hevea brasiliensis (75%) and Ricinus communis (75%). The BnCPI was subcloned into expression vector pSmart-I and then overexpressed in Escherichia coli BL21 (DE3) as a His-tagged recombinant protein. The purified reBnCPI has a molecular mass of 11.4 kDa determined by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS–PAGE). Purified reBnCPI can efficiently inhibit the protease activity of papain and ficin toward BANA (Nα-benzoyl-L-arginine-2-naphthyamide), as well as the mycelium growth of some important plant pathogenic fungi. The data further contribute to our understanding of the molecular functions of BnCPI.

Cloning and characterization of the first serine carboxypeptidase from a plant parasitic nematode, Radopholus similis

Radopholus similis is an important parasitic nematode of plants. Serine carboxypeptidases (SCPs) are peptidases that hydrolyse peptides and proteins and play critical roles in the development, invasion, and pathogenesis of certain parasitic nematodes and other animal pathogens. In this study, we obtained the full-length sequence of the SCP gene from R. similis (Rs-scp-1), which is 1665 bp long and includes a 1461-bp open reading frames encoding 486 amino acids with an 18-aa signal peptide. This gene is a double-copy gene in R. similis. Rs-scp-1 was expressed in the procorpus, esophageal glands and intestines of females and in the esophageal glands and intestines of juveniles. Rs-scp-1 expression levels were highest in females, followed by juveniles and males, and lowest in eggs. Rs-scp-1 expression levels were significantly suppressed after R. similis was soaked in Rs-scp-1 dsRNA for 12 h. Nematodes were then inoculated into Anthurium andraeanum after RNAi treatment. Compared with water treatment, R. similis treated with RNAi were reduced in number and pathogenicity. In summary, we obtained the first SCP gene from a plant parasitic nematode and confirmed its role in the parasitic process.

Comparative proteomics of hydatid fluids from two Echinococcus multilocularis isolates

The hydatid fluid (HF) that fills Echinococcus multilocularis metacestode vesicles is a complex mixture of proteins from both parasite and host origin. Here, a LC-MS/MS approach was used to compare the HF composition of E. multilocularis H95 and G8065 isolates (EmH95 and EmG8065, respectively), which present differences in terms of growth and fertility. Overall, 446 unique proteins were identified, 392 of which (88%) were from parasite origin and 54 from culture medium. At least 256 of parasite proteins were sample exclusive, and 82 of the 136 shared proteins presented differential abundance between E. multilocularis isolates. The parasite's protein repertoires in EmH95 and EmG8065 HF samples presented qualitative and quantitative differences involving antigens, signaling proteins, proteolytic enzymes, protease inhibitors and chaperones, highlighting intraspecific singularities that could be correlated to biological features of each isolate. The repertoire of medium proteins found in the HF was also differential between isolates, and the relevance of the HF exogenous protein content for the parasite's biology is discussed. The repertoires of identified proteins also provided potential molecular markers for important biological features, such as parasite growth rate and fertility, as well potential protein targets for the development of novel diagnostic and treatment strategies for alveolar echinococcosis.

BIOLOGICAL SIGNIFICANCE

E. multilocularis metacestode infection of mammal hosts involve complex interactions mediated by excretory/secretory (ES) products. The hydatid fluid (HF) that fills the E. multilocularis metacestode vesicles contains complex repertoires of parasite ES products and host proteins that mediate important molecular interactions determinant for parasite survival and development, and, consequently, to the infection outcome. HF has been also extensively reported as the main source of proteins for the immunodiagnosis of echinococcosis. The performed proteomic analysis provided a comprehensive profiling of the HF protein composition of two E. multilocularis isolates. This allowed us to identify proteins of both parasite and exogenous (medium) origin, many of which present significant differential abundances between parasite isolates and may correlate to their differential biological features, including fertility and growth rate.