Exploring the transcriptome of the burrowing nematode Radopholus similis

Transgenic East African Highland Banana Plants Are Protected against Radopholus similis through Host-Delivered RNAi

The burrowing nematode Radopholus similis is considered a major problem of intensive banana cultivation. It can cause extensive root damage resulting in the toppling disease of banana, which means that plants fall to the ground. Soaking R. similis in double-stranded (ds) RNA of the nematode genes Rps13, chitin synthase (Chs-2), Unc-87, Pat-10 or beta-1,4-endoglucanase (Eng1a) suppressed reproduction on carrot discs, from 2.8-fold (Chs-2) to 7-fold (Rps13). The East African Highland Banana cultivar Nakitembe was then transformed with constructs for expression of dsRNA against the same genes, and for each construct, 30 independent transformants were tested with nematode infection. Four months after transfer from in vitro culture to the greenhouse, the banana plants were transferred to a screenhouse and inoculated with 2000 nematodes per plant, and thirteen weeks later, they were analyzed for several parameters including plant growth, root necrosis and final nematode population. Plants with dsRNA constructs against the nematode genes were on average showing lower nematode multiplication and root damage than the nontransformed controls or the banana plants expressing dsRNA against the nonendogenous gene. In conclusion, RNAi seems to efficiently protect banana against damage caused by R. similis, opening perspectives to control this pest.

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.

The Wolbachia Symbiont: Here, There and Everywhere

Wolbachia symbionts, first observed in the 1920s, are now known to be present in about 30-70% of tested arthropod species, in about half of tested filarial nematodes (including the majority of human filarial nematodes), and some plant-parasitic nematodes. In arthropods, they are generally viewed as parasites while in nematodes they appear to be mutualists although this demarcation is not absolute. Their presence in arthropods generally leads to reproductive anomalies, while in nematodes, they are generally required for worm development and reproduction. In mosquitos, Wolbachia inhibit RNA viral infections, leading to populational reductions in human RNA virus pathogens, whereas in filarial nematodes, their requirement for worm fertility and survival has been channeled into their use as drug targets for filariasis control. While much more research on these ubiquitous symbionts is needed, they are viewed as playing significant roles in biological processes, ranging from arthropod speciation to human health.

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.

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.

A Draft Genome Sequence of the Burrowing Nematode Radopholus similis

Radopholus similis also known as the burrowing nematode is a devastating pest of banana (Musa spp.) and many economically important crops and ornamentals. In this publication, we present the genome assembly of R. similis.

Genome Announcement: The Draft Genomes of Two Radopholus similis populations from Costa Rica

Radopholus similis is an economically important pest of both banana and citrus in tropical regions. Here we present draft genomes from two populations of R. similis from Costa Rica that were created and assembled using short read libraries from Illumina HiSeq technology.

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.

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.

Identification of candidate effector genes of Pratylenchus penetrans

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

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.

The mitochondrial genome of Globodera ellingtonae is composed of two circles with segregated gene content and differential copy numbers

Background

The evolution of animal mitochondrial (mt) genomes has resulted in a highly conserved structure: a single compact circular chromosome approximately 14 to 20 kb long. Within the last two decades exceptions to this conserved structure, such as the division of the genome into multiple chromosomes, have been reported in a diverse set of metazoans. We report on the two circle multipartite mt genome of a newly described cyst nematode, Globodera ellingtonae.

Results

The G. ellingtonae mt genome was found to be comprised of two circles, each larger than any other multipartite circular mt chromosome yet reported, and both were larger than the single mt circle of the model nematode Caenorhabditis elegans. The genetic content of the genome was disproportionately divided between the two circles, although they shared a ~6.5 kb non-coding region. The 17.8 kb circle (mtDNA-I) contained ten protein-coding genes and two tRNA genes, whereas the 14.4 kb circle (mtDNA-II) contained two protein-coding genes, 20 tRNA genes and both rRNA genes. Perhaps correlated with this division of genetic content, the copy number of mtDNA-II was more than four-fold that of mtDNA-I in individual nematodes. The difference in copy number increased between second-stage and fourth-stage juveniles.

Conclusions

The segregation of gene types to different mt circles in G. ellingtonae could provide benefit by localizing gene functional types to independent transcriptional units. This is the first report of both two-circle and several-circle mt genomes within a single genus. The differential copy number associated with this multipartite mt organization could provide a model system for deconstructing mechanisms regulating mtDNA copy number both in somatic cells and during germline development.

Electronic supplementary material

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

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 of glutathione S-transferases from the pine wood nematode, Bursaphelenchus xylophilus

We have previously identified two secreted glutathione S-transferases (GST) expressed in the pharyngeal gland cell of Bursaphelenchus xylophilus, which are upregulated post infection of the host. This study examines the functional role of GSTs in B. xylophilus biology. We analysed the expression profiles of all predicted GSTs in the genome and the results showed that they belong to kappa and cytosolic subfamilies and the majority are upregulated post infection of the host. A small percentage is potentially secreted and none is downregulated post infection of the host. One secreted protein was confirmed as a functional GST and is within a cluster that showed the highest expression fold change in infection. This enzyme has a protective activity that may involve host defences, namely in the presence of terpenoid compounds and peroxide products. These results suggest that GSTs secreted into the host participate in the detoxification of host-derived defence compounds and enable successful parasitism.

Only a small subset of the SPRY domain gene family in Globodera pallida is likely to encode effectors, two of which suppress host defences induced by the potato resistance gene Gpa2

Analysis of the genome sequence of the potato cyst nematode, Globodera pallida, has shown that a substantial gene family (approximately 300 sequences) of proteins containing a SPRY domain is present in this species. This is a huge expansion of the gene family as compared to other organisms, including other plant-parasitic nematodes. Some SPRY domain proteins from G. pallida and G. rostochiensis have signal peptides for secretion and are deployed as effectors. One of these SPRYSEC proteins has been shown to suppress host defence responses. We describe further analysis of this gene family in G. pallida. We show that only a minority (10%) of the SPRY domain proteins in this species have a predicted signal peptide for secretion and that the presence of a signal peptide is strongly correlated with the corresponding gene being expressed at the early stages of parasitism. The data suggest that while the gene family is greatly expanded, only a minority of SPRY domain proteins in G. pallida are SPRYSEC candidate effectors. We show that several new SPRYSECs from G. pallida are expressed in the dorsal gland cell and demonstrate that some, but not all, of the SPRYSECs can suppress the hypersensitive response induced by co-expression of the resistance gene Gpa2 and its cognate avirulence factor RBP-1 in Nicotiana benthamiana.

Absence of Wolbachia endobacteria in Sri Lankan isolates of the nematode parasite of animals Setaria digitata

Setaria digitata is an animal filarial parasite with natural hosts of cattle and buffaloes that causes mild disease conditions. Infection of non-permissive hosts such as goats, sheep and horses, by this nematode can cause cerebrospinal nematodiasis that leads to lumbar paralysis and the eventual death of the animals and inflicts considerable economic losses on livestock farmers. Wolbachia are obligate mutualistic endosymbionts for some filarial nematodes and are currently being targeted for the control of diseases caused by these parasites. However, little is known about the occurrence of this endosymbiont in the Setariidae family. In this work, worms collected from infected cattle in Sri Lanka were morphologically identified as S. digitata and tested for the presence of Wolbachia by PCR screening using the WSP- and Wolbachia-specific 16S rRNA and multilocus sequence typing primers that were designed to amplify the gatB, coxA, hcpA, ftsZ and fbpA sequences of Wolbachia. The presence of endobacteria in S. digitata was also examined by whole-mount immunofluorescence staining of the parasites and transmission electron microscopic studies. These analyses did not produce evidence of presence of Wolbachia or any other endosymbiotic bacteria in S. digitata, whereas such evidence was found in Brugia malayi, which was used as a positive control in this study.

Identification and characterisation of a hyper-variable apoplastic effector gene family of the potato cyst nematodes

Sedentary endoparasitic nematodes are obligate biotrophs that modify host root tissues, using a suite of effector proteins to create and maintain a feeding site that is their sole source of nutrition. Using assumptions about the characteristics of genes involved in plant-nematode biotrophic interactions to inform the identification strategy, we provide a description and characterisation of a novel group of hyper-variable extracellular effectors termed HYP, from the potato cyst nematode Globodera pallida. HYP effectors comprise a large gene family, with a modular structure, and have unparalleled diversity between individuals of the same population: no two nematodes tested had the same genetic complement of HYP effectors. Individuals vary in the number, size, and type of effector subfamilies. HYP effectors are expressed throughout the biotrophic stages in large secretory cells associated with the amphids of parasitic stage nematodes as confirmed by in situ hybridisation. The encoded proteins are secreted into the host roots where they are detectable by immunochemistry in the apoplasm, between the anterior end of the nematode and the feeding site. We have identified HYP effectors in three genera of plant parasitic nematodes capable of infecting a broad range of mono- and dicotyledon crop species. In planta RNAi targeted to all members of the effector family causes a reduction in successful parasitism.

Sedentary plant parasitic nematodes are pathogens that invade plant roots and establish a feeding site. The feeding site is a specialist structure used by the nematode to support its development within the plant. The nematode secretes a suite of proteins, termed ‘effector proteins’ that are responsible for initiating and maintaining the feeding site. The nematode must also evade recognition by the plant defence systems throughout its lifecycle that can last for many weeks. We describe a diverse and variable effector gene family (HYP), the products of which are secreted into the plant by the nematode and are required for successful infection. The variability and modular structure of this gene family can lead to the production of a large array of effector proteins. This diversity may allow the nematodes to combat any resistance mechanisms developed by the plant. Each nematode tested within a population is genetically unique in terms of these effector genes. We found huge variation in the number, size and type of HYP effectors at the level of the individual. This may explain some of the difficulties in breeding nematode resistant plants and has profound implications for those working with other plant pathogens.

The evolution of parasitism in Nematoda

Nematodes are abundant and diverse, and include many parasitic species. Molecular phylogenetic analyses have shown that parasitism of plants and animals has arisen at least 15 times independently. Extant nematode species also display lifestyles that are proposed to be on the evolutionary trajectory to parasitism. Recent advances have permitted the determination of the genomes and transcriptomes of many nematode species. These new data can be used to further resolve the phylogeny of Nematoda, and identify possible genetic patterns associated with parasitism. Plant-parasitic nematode genomes show evidence of horizontal gene transfer from other members of the rhizosphere, and these genes play important roles in the parasite-host interface. Similar horizontal transfer is not evident in animal parasitic groups. Many nematodes have bacterial symbionts that can be essential for survival. Horizontal transfer from symbionts to the nematode is also common, but its biological importance is unclear. Over 100 nematode species are currently targeted for sequencing, and these data will yield important insights into the biology and evolutionary history of parasitism. It is important that these new technologies are also applied to free-living taxa, so that the pre-parasitic ground state can be inferred, and the novelties associated with parasitism isolated.

Analysis of the transcriptome of Hirschmanniella oryzae to explore potential survival strategies and host-nematode interactions

The rice root nematode Hirschmanniella oryzae is the most abundant plant-parasitic nematode in flooded rice fields and is distributed world-wide. Although it is economically less important than sedentary nematodes, it can cause severe yield reductions and economic losses in specific environmental conditions. No transcriptome data for this genus were available until now. We have performed 454 sequencing on a mixed life stages population to gain an insight into nematode-plant interactions and nematode survival strategies. The results of two assembly strategies were combined to reduce the redundancy of the data, generating a final dataset of 21 360 contigs. The data were screened for putative plant cell wall-modifying proteins, which facilitate nematode migration through host roots. A β-mannanase, previously not reported in nematodes, was detected in the dataset. The data were screened for putative effector proteins that may alter the host defence mechanism. Two enzymes, chorismate mutase and isochorismatase, thought to be involved in the salicyclic acid pathway, were identified. Experimental treatments of H. oryzae with artificial seawater showed that late embryogenesis abundant (LEA) proteins and SXP/RAL-2 are induced, suggesting that these proteins are involved in the process of anhydrobiosis. The newly generated data can highlight potential differences between sedentary and migratory nematodes, and will be useful in the further study of host-nematode interactions and the developmental biology of this nematode.

Phylogenetic analysis of entomoparasitic nematodes, potential control agents of flea populations in natural foci of plague

Entomoparasitic nematodes are natural control agents for many insect pests, including fleas that transmit Yersinia pestis, a causative agent of plague, in the natural foci of this extremely dangerous zoonosis. We examined the flea samples from the Volga-Ural natural focus of plague for their infestation with nematodes. Among the six flea species feeding on different rodent hosts (Citellus pygmaeus, Microtus socialis, and Allactaga major), the rate of infestation varied from 0 to 21%. The propagation rate of parasitic nematodes in the haemocoel of infected fleas was very high; in some cases, we observed up to 1,000 juveniles per flea specimen. Our study of morphology, life cycle, and rDNA sequences of these parasites revealed that they belong to three distinct species differing in the host specificity. On SSU and LSU rRNA phylogenies, these species representing three genera (Rubzovinema, Psyllotylenchus, and Spilotylenchus), constitute a monophyletic group close to Allantonema and Parasitylenchus, the type genera of the families Allantonematidae and Parasitylenchidae (Nematoda: Tylenchida). We discuss the SSU-ITS1-5.8S-LSU rDNA phylogeny of the Tylenchida with a special emphasis on the suborder Hexatylina.

Absence of Wolbachia endobacteria in the human parasitic nematode Dracunculus medinensis and two related Dracunculus species infecting wildlife

BACKGROUND

Wolbachia endosymbionts are a proven target for control of human disease caused by filarial nematodes. However, little is known about the occurrence of Wolbachia in taxa closely related to the superfamily Filarioidea. Our study addressed the status of Wolbachia presence in members of the superfamily Dracunculoidea by screening the human parasite Dracunculus medinensis and related species from wildlife for Wolbachia.

FINDINGS

D. medinensis, D. lutrae and D. insignis specimens were all negative for Wolbachia colonization by PCR screening for the Wolbachia ftsZ, 16S rRNA and Wolbachia surface protein (wsp) sequences. The quality and purity of the DNA preparations was confirmed by amplification of nematode 18S rRNA and cytochrome c oxidase subunit I sequences. Furthermore, Wolbachia endobacteria were not detected by whole mount fluorescence staining, or by immunohistochemistry using a Wolbachia-specific antiserum. In contrast, positive control Brugia malayi worms were shown to harbour Wolbachia by PCR, fluorescence staining and immunohistochemistry.

CONCLUSIONS

Three examined species of Dracunculus showed no evidence of Wolbachia endobacteria. This supports that members of the superfamily Dracunculoidea are free of Wolbachia. Within the order Spirurida, these endosymbionts appear restricted to the Filarioidea.

Transcriptomic analysis of the rice white tip nematode, Aphelenchoides besseyi (Nematoda: Aphelenchoididae)

Background

The rice white tip nematode Aphelenchoides besseyi, a devastating nematode whose genome has not been sequenced, is distributed widely throughout almost all the rice-growing regions of the world. The aims of the present study were to define the transcriptome of A. besseyi and to identify parasite-related, mortality-related or host resistance-overcoming genes in this nematode.

Methodology and Principal Findings

Using Solexa/Illumina sequencing, we profiled the transcriptome of mixed-stage populations of A. besseyi. A total of 51,270 transcripts without gaps were produced based on high-quality clean reads. Of all the A. besseyi transcripts, 9,132 KEGG Orthology assignments were annotated. Carbohydrate-active enzymes of glycoside hydrolases (GHs), glycosyltransferases (GTs), carbohydrate esterases (CEs) and carbohydrate-binding modules (CBMs) were identified. The presence of the A. besseyi GH45 cellulase gene was verified by in situ hybridization. Given that 13 unique A. besseyi potential effector genes were identified from 41 candidate effector homologs, further studies of these homologs are merited. Finally, comparative analyses were conducted between A. besseyi contigs and Caenorhabditis elegans genes to look for orthologs of RNAi phenotypes, neuropeptides and peptidases.

Conclusions and Significance

The present results provide comprehensive insight into the genetic makeup of A. besseyi. Many of this species' genes are parasite related, nematode mortality-related or necessary to overcome host resistance. The generated transcriptome dataset of A. besseyi reported here lays the foundation for further studies of the molecular mechanisms related to parasitism and facilitates the development of new control strategies for this species.

Top 10 plant-parasitic nematodes in molecular plant pathology

The aim of this review was to undertake a survey of researchers working with plant-parasitic nematodes in order to determine a 'top 10' list of these pathogens based on scientific and economic importance. Any such list will not be definitive as economic importance will vary depending on the region of the world in which a researcher is based. However, care was taken to include researchers from as many parts of the world as possible when carrying out the survey. The top 10 list emerging from the survey is composed of: (1) root-knot nematodes (Meloidogyne spp.); (2) cyst nematodes (Heterodera and Globodera spp.); (3) root lesion nematodes (Pratylenchus spp.); (4) the burrowing nematode Radopholus similis; (5) Ditylenchus dipsaci; (6) the pine wilt nematode Bursaphelenchus xylophilus; (7) the reniform nematode Rotylenchulus reniformis; (8) Xiphinema index (the only virus vector nematode to make the list); (9) Nacobbus aberrans; and (10) Aphelenchoides besseyi. The biology of each nematode (or nematode group) is reviewed briefly.

Exploring the host parasitism of the migratory plant-parasitic nematode Ditylenchus destuctor by expressed sequence tags analysis

The potato rot nematode, Ditylenchus destructor, is a very destructive nematode pest on many agriculturally important crops worldwide, but the molecular characterization of its parasitism of plant has been limited. The effectors involved in nematode parasitism of plant for several sedentary endo-parasitic nematodes such as Heterodera glycines, Globodera rostochiensis and Meloidogyne incognita have been identified and extensively studied over the past two decades. Ditylenchus destructor, as a migratory plant parasitic nematode, has different feeding behavior, life cycle and host response. Comparing the transcriptome and parasitome among different types of plant-parasitic nematodes is the way to understand more fully the parasitic mechanism of plant nematodes. We undertook the approach of sequencing expressed sequence tags (ESTs) derived from a mixed stage cDNA library of D. destructor. This is the first study of D. destructor ESTs. A total of 9800 ESTs were grouped into 5008 clusters including 3606 singletons and 1402 multi-member contigs, representing a catalog of D. destructor genes. Implementing a bioinformatics' workflow, we found 1391 clusters have no match in the available gene database; 31 clusters only have similarities to genes identified from D. africanus, the most closely related species to D. destructor; 1991 clusters were annotated using Gene Ontology (GO); 1550 clusters were assigned enzyme commission (EC) numbers; and 1211 clusters were mapped to 181 KEGG biochemical pathways. 22 ESTs had similarities to reported nematode effectors. Interestedly, most of the effectors identified in this study are involved in host cell wall degradation or modification, such as 1,4-beta-glucanse, 1,3-beta-glucanse, pectate lyase, chitinases and expansin, or host defense suppression such as calreticulin, annexin and venom allergen-like protein. This result implies that the migratory plant-parasitic nematode D. destructor secrets similar effectors to those of sedentary plant nematodes. Finally we further characterized the two D. destructor expansin proteins.

Transcriptome profiling of the fifth-stage larvae of Angiostrongylus cantonensis by next-generation sequencing

Angiostrongylus cantonensis is an important zoonotic nematode. It is the causative agent of eosinophilic meningitis and eosinophilic meningoencephalitis in humans. However, information of this parasite at the genomic level is very limited. In the present study, the transcriptomic profiles of the fifth-stage larvae (L5) of A. cantonensis were investigated by next-generation sequencing (NGS). In the NGS database established from the larvae isolated from the brain of Sprague–Dawley rats, 31,487 unique genes with a mean length of 617 nucleotides were assembled. These genes were found to have a 46.08 % significant similarity to Caenorhabditis elegans by BLASTx. They were then compared with the expressed sequence tags of 18 other nematodes, and significant matches of 36.09–59.12 % were found. Among these genes, 3,338 were found to participate in 124 Kyoto Encyclopedia of Genes and Genomes pathways. These pathways included 1,514 metabolisms, 846 genetic information processing, 358 environmental information processing, 264 cellular processes, and 91 organismal systems. Analysis of 30,816 sequences with the gene ontology database indicated that their annotations included 5,656 biological processes (3,364 cellular processes, 3,061 developmental processes, and 3,191 multicellular organismal processes), 7,218 molecular functions (4,597 binding and 3,084 catalytic activities), and 4,719 cellular components (4,459 cell parts and 4,466 cells). Moreover, stress-related genes (112 heat stress and 33 oxidation stress) and genes for proteases (159) were not uncommon. This study is the first NGS-based study to set up a transcriptomic database of A. cantonensis L5. The results provide new insights into the survival, development, and host–parasite interactions of this blood-feeding nematode.

Multiple ancient horizontal gene transfers and duplications in lepidopteran species

Eukaryotic horizontal gene transfer (HGT) events are increasingly being discovered yet few reports have summarized multiple occurrences in a wide range of species. We systematically investigated HGT events in the order Lepidoptera by employing a series of filters. Bombyx mori, Danaus plexippus and Heliconius melpomene had 13, 12 and 12 HGTs, respectively, from bacteria and fungi. These HGTs contributed a total of 64 predicted genes: 22 to B. mori, 22 to D. plexippus and 20 to H. melpomene. Several new genes were generated by post-transfer duplications. Post-transfer duplication of a suite of functional HGTs has rarely been reported in higher organisms. The distributional patterns of paralogues for certain genes differed in the three species, indicating potential independent duplication or loss events. All of these HGTs had homologues expressed in some other lepidopterans, indicating ancient transfer events. Most HGTs were involved in the metabolism of sugar and amino acids. These HGTs appeared to have experienced amelioration, purifying selection and accelerated evolution to adapt to the background genome of the recipient. The discovery of ancient, massive HGTs and duplications in lepidopterans and their adaptive evolution provides further insights into the evolutionary significance of the events from donors to multicellular host recipients.

Comparative transcriptomics of two pathogenic pinewood nematodes yields insights into parasitic adaptation to life on pine hosts

Bursaphelenchus xylophilus and Bursaphelenchus mucronatus are migratory endoparasitic nematodes that live in pine trees. To gain insight into their molecular similarities and differences, transcriptomes of the two nematodes were analysed. A total of 23,765 and 21,782 contigs (>300 bp) were obtained from B. xylophilus and B. mucronatus, respectively. More than 80% of the contigs could map to each other's transcriptome reciprocally. A total of 23,467 and 21,370 Open Reading Frames were predicted, respectively. Besides those known parasitism-related proteins, six new venom allergen-like proteins (VAPs) were found, which were not homologous to known VAPs. Enzymes involved in xenobiotic biodegradation were abundant in the two transcriptomes based on KEGG functional annotation. Metabolism of xenobiotics by cytochrome P450 comprised the main detoxification pathways. The mRNA expression levels of detoxification genes in nematodes living in the host were higher than those in nematodes feeding on fungus. However, there were fewer enzymes involved in the α-pinene degradation. Our results indicate that the two pinewood nematodes have evolved similar molecular mechanisms to adapt to life on pine hosts.

Genomics and transcriptomics across the diversity of the Nematoda

The diversity of biology in nematodes is reflected in the diversity of their genomes. Parasitic species in particular have evolved mechanisms to invade and outwit their hosts, and these offer opportunities for the development of control measures. Genomic analyses can reveal the molecular underpinnings of phenotypes such as parasitism and thus, initiate and support research programmes that explore the manipulation of host and parasite physiologies to achieve favourable outcomes. Wide sampling across nematode diversity allows phylogenetically informed formulation of research hypotheses, identification of core features shared by all species or important evolutionary novelties present in isolated clades. Many nematode species have been investigated through the use of the expressed sequence tag approach, which samples from the transcribed genome. Gene catalogues generated in this way can be explored to reveal the patterns of expression associated with parasitism and candidates for testing as drug targets or vaccine components. Analysis environments, such as NEMBASE facilitate exploitation of these data. The development of new high-throughput DNA-sequencing technologies has facilitated transcriptomic and genomic approaches to parasite biology. Whole genome sequencing offers more complete catalogues of genes and assists a systems approach to phenotype dissection. These efforts are being coordinated through the 959 Nematode Genomes initiative.

Functional roles of effectors of plant-parasitic nematodes

Plant pathogens have evolved a variety of different strategies that allow them to successfully infect their hosts. Plant-parasitic nematodes secrete numerous proteins into their hosts. These proteins, called effectors, have various functions in the plant cell. The most studied effectors to date are the plant cell wall degrading enzymes, which have an interesting evolutionary history since they are believed to have been acquired from bacteria or fungi by horizontal gene transfer. Extensive genome, transcriptome and proteome studies have shown that plant-parasitic nematodes secrete many additional effectors. The function of many of these is less clear although during the last decade, several research groups have determined the function of some of these effectors. Even though many effectors remain to be investigated, it has already become clear that they can have very diverse functions. Some are involved in suppression of plant defences, while others can specifically interact with plant signalling or hormone pathways to promote the formation of nematode feeding sites. In this review, the most recent progress in the understanding of the function of plant-parasitic nematode effectors is discussed.

Analysis of the transcriptome of the root lesion nematode Pratylenchus coffeae generated by 454 sequencing technology

To study interactions between plants and plant-parasitic nematodes, several omics studies have nowadays become extremely useful. Since most data available so far is derived from sedentary nematodes, we decided to improve the knowledge on migratory nematodes by studying the transcriptome of the nematode Pratylenchus coffeae through generating expressed sequence tags (ESTs) on a 454 sequencing platform. In this manuscript we present the generation, assembly and annotation of over 325,000 reads from P. coffeae. After assembling these reads, 56,325 contigs and singletons with an average length of 353bp were selected for further analyses. Homology searches revealed that 25% of these sequences had significant matches to the Swiss-prot/trEMBL database and 29% had significant matches in nematode ESTs. Over 10,000 sequences were successfully annotated, corresponding to over 6000 unique Gene Ontology identifiers and 5000 KEGG orthologues. Different approaches led to the identification of different sequences putatively involved in the parasitism process. Several plant cell wall modifying enzymes were identified, including an arabinogalactan galactosidase, so far identified in cyst nematodes only. Additionally, some new putative cell wall modifying enzymes are present belonging to GHF5 and GHF16, although further functional studies are needed to determine the true role of these proteins. Furthermore, a homologue to a chorismate mutase was found, suggesting that this parasitism gene has a wider occurrence in plant-parasitic nematodes than previously assumed. Finally, the dataset was searched for orthologues against the Meloidogyne genomes and genes involved in the RNAi pathway. In conclusion, the generated transcriptome data of P. coffeae will be very useful in the future for several projects: (1) evolutionary studies of specific gene families, such as the plant cell wall modifying enzymes, (2) the identification and functional analysis of candidate effector genes, (3) the development of new control strategies, e.g. by finding new targets for RNAi and (4) the annotation of the upcoming genome sequence.

First transcriptomic analysis of the economically important parasitic nematode, Trichostrongylus colubriformis, using a next-generation sequencing approach

Trichostrongylus colubriformis (Strongylida), a small intestinal nematode of small ruminants, is a major cause of production and economic losses in many countries. The aims of the present study were to define the transcriptome of the adult stage of T. colubriformis, using 454 sequencing technology and bioinformatic analyses, and to predict the main pathways that key groups of molecules are linked to in this nematode. A total of 21,259 contigs were assembled from the sequence data produced from a normalized cDNA library; 7876 of these contigs had known orthologues in the free-living nematode Caenorhabditis elegans, and encoded, amongst others, proteins with 'transthyretin-like' (8.8%), 'RNA recognition' (8.4%) and 'metridin-like ShK toxin' (7.6%) motifs. Bioinformatic analyses inferred that relatively high proportions of the C. elegans homologues are involved in biological pathways linked to 'peptidases' (4%), 'ribosome' (3.6%) and 'oxidative phosphorylation' (3%). Highly represented were peptides predicted to be associated with the nervous system, digestion of host proteins or inhibition of host proteases. Probabilistic functional gene networking of the complement of C. elegans orthologues (n=2126) assigned significance to particular subsets of molecules, such as protein kinases and serine/threonine phosphatases. The present study represents the first, comprehensive insight into the transcriptome of adult T. colubriformis, which provides a foundation for fundamental studies of the molecular biology and biochemistry of this parasitic nematode as well as prospects for identifying targets for novel nematocides. Future investigations should focus on comparing the transcriptomes of different developmental stages, both genders and various tissues of this parasitic nematode for the prediction of essential genes/gene products that are specific to nematodes.

The Wolbachia endosymbiont as an anti-filarial nematode target

Human disease caused by parasitic filarial nematodes is a major cause of global morbidity. The parasites are transmitted by arthropod intermediate hosts and are responsible for lymphatic filariasis (elephantiasis) or onchocerciasis (river blindness). Within these filarial parasites are intracellular alpha-proteobacteria, Wolbachia, that were first observed almost 30 years ago. The obligate endosymbiont has been recognized as a target for anti-filarial nematode chemotherapy as evidenced by the loss of worm fertility and viability upon antibiotic treatment in an extensive series of human trials. While current treatments with doxycycline and rifampicin are not practical for widespread use due to the length of required treatments and contraindications, anti-Wolbachia targeting nevertheless appears a promising alternative for filariasis control in situations where current programmatic strategies fail or are unable to be delivered and it provides a superior efficacy for individual therapy. The mechanisms that underlie the symbiotic relationship between Wolbachia and its nematode hosts remain elusive. Comparative genomics, bioinfomatic and experimental analyses have identified a number of potential interactions, which may be drug targets. One candidate is de novo heme biosynthesis, due to its absence in the genome sequence of the host nematode, Brugia malayi, but presence in Wolbachia and its potential roles in worm biology. We describe this and several additional candidate targets, as well as our approaches for understanding the nature of the host-symbiont relationship.

Emerging molecular knowledge on Radopholus similis, an important nematode pest of banana

TAXONOMY

Superkingdom Eukaryota; Kingdom Metazoa; Phylum Nematoda; Class Chromadorea; Order Rhabditida; Suborder Tylenchina; Infraorder Tylenchomorpha; Superfamily Tylenchoidea; Family Pratylenchidae; Subfamily Radopholinae; Genus Radopholus.

PHYSICAL PROPERTIES

Microscopic unsegmented worm; migratory endoparasite of plants. Strong sexual dimorphism; reproduction both by amphimixis and self-fertilization.

HOSTS

Over 250 different plant species, including citrus, black pepper and banana (main host plant).

SYMPTOMS

Purple to black lesions and extensive cavities in plant roots, leading to reduced uptake of water and nutrients. In banana, this may result in poor vegetative growth, reduced bunch weight and toppling of plants.

DISEASE CONTROL

Nematicides, alternative cropping systems, nematode-free planting material, some resistant cultivars.

AGRONOMIC IMPORTANCE

Major problem in banana plantations in tropical regions worldwide.

Transcript analysis of parasitic females of the sedentary semi-endoparasitic nematode Rotylenchulus reniformis

Rotylenchulus reniformis, the reniform nematode, is a sedentary semi-endoparasitic nematode capable of infecting >300 plant species, including a large number of crops such as cotton, soybean, and pineapple. In contrast to other economically important plant-parasitic nematodes, molecular genetic data regarding the R. reniformis transcriptome is virtually nonexistant. Herein, we present a survey of R. reniformis ESTs that were sequenced from a sedentary parasitic female cDNA library. Cluster analysis of 2004 high quality ESTs produced 123 contigs and 508 singletons for a total of 631 R. reniformis unigenes. BLASTX analyses revealed that 39% of all unigenes showed similarity to known proteins (E<or=1.0e-04). R. reniformis genes homologous to known parasitism genes were identified and included beta-1,4-endoglucanase, fatty acid- and retinol-binding proteins, and an esophageal gland cell-specific gene from Heterodera glycines. Furthermore, a putative ortholog of an enzyme involved in thiamin biosynthesis, thought to exist solely in prokaryotes, fungi, and plants, was identified. Lastly, 114 R. reniformis unigenes orthologous to RNAi-lethal Caenorhabditis elegans genes were discovered. The work described here offers a glimpse into the transcriptome of a sedentary semi-endoparasitic nematode which (i) provides the transcript sequence data necessary for investigating engineered resistance against R. reniformis and (ii) hints at the existance of a thiamin biosynthesis pathway in an animal.

Identification and functional characterization of effectors in expressed sequence tags from various life cycle stages of the potato cyst nematode Globodera pallida

In this article, we describe the analysis of over 9000 expressed sequence tags (ESTs) from cDNA libraries obtained from various life cycle stages of Globodera pallida. We have identified over 50 G. pallida effectors from this dataset using bioinformatics analysis, by screening clones in order to identify secreted proteins up-regulated after the onset of parasitism and using in situ hybridization to confirm the expression in pharyngeal gland cells. A substantial gene family encoding G. pallida SPRYSEC proteins has been identified. The expression of these genes is restricted to the dorsal pharyngeal gland cell. Different members of the SPRYSEC family of proteins from G. pallida show different subcellular localization patterns in plants, with some localized to the cytoplasm and others to the nucleus and nucleolus. Differences in subcellular localization may reflect diverse functional roles for each individual protein or, more likely, variety in the compartmentalization of plant proteins targeted by the nematode. Our data are therefore consistent with the suggestion that the SPRYSEC proteins suppress host defences, as suggested previously, and that they achieve this through interaction with a range of host targets.

A unique genetic code change in the mitochondrial genome of the parasitic nematode Radopholus similis

BACKGROUND

Mitochondria (mt) contain their own autonomously replicating DNA, constituted as a small circular genome encoding essential subunits of the respiratory chain. Mt DNA is characterized by a genetic code which differs from the standard one. Interestingly, the mt genome of nematodes share some peculiar features, such as small transfer RNAs, truncated ribosomal RNAs and - in the class of Chromadorean nematodes - unidirectional transcription.

FINDINGS

We present the complete mt genomic sequence (16,791 bp) of the plant-parasitic nematode Radopholus similis (class Chromadorea). Although it has a gene content similar to most other nematodes, many idiosyncrasies characterize the extremely AT-rich mt genome of R. similis (85.4% AT). The secondary structure of the large (16S) rRNA is further reduced, the gene order is unique, the large non-coding region contains two large repeats, and most interestingly, the UAA codon is reassigned from translation termination to tyrosine. In addition, 7 out of 12 protein-coding genes lack a canonical stop codon and analysis of transcriptional data showed the absence of polyadenylation. Northern blot analysis confirmed that only one strand is transcribed and processed. Furthermore, using nucleotide content bias methods, regions for the origin of replication are suggested.

CONCLUSION

The extraordinary mt genome of R. similis with its unique genetic code appears to contain exceptional features correlated to DNA decoding. Therefore the genome may provide an incentive to further elucidate these barely understood processes in nematodes. This comprehension may eventually lead to parasitic nematode-specific control targets as healthy mitochondria are imperative for organism survival. In addition, the presented genome is an interesting exceptional event in genetic code evolution.

Characterization of a putative endoxylanase in the migratory plant-parasitic nematode Radopholus similis

Plant-parasitic nematodes have developed an arsenal of enzymes to degrade the rigid plant cell wall. In this article, we report the presence of a putative endoxylanase in the migratory endoparasitic nematode Radopholus similis. This enzyme is thought to facilitate the migration of the nematode, as it breaks down xylan, the major component of hemicellulose. The corresponding gene (Rs-xyl1) was cloned and the sequence revealed three small introns. Interestingly, the position of all three introns was conserved in a putative endoxylanase from Meloidogyne hapla, and the position of one intron was conserved in two endoxylanases from Meloidogyne incognita, which suggests a common ancestral gene. The spatial and temporal expression of the Rs-xyl1 gene was examined by in situ hybridization and semi-quantitative reverse transcriptase-polymerase chain reaction. The putative protein consists of a signal peptide, a catalytic domain and a carbohydrate-binding module (CBM). The catalytic domain showed similarity to both glycosyl hydrolase family 5 (GHF5) and GHF30 enzymes. Using Hidden Markov Model profiles and phylogenetic analysis, we were able to show that Rs-XYL1 and its closest homologues are not members of GHF5, as suggested previously, but rather form a subclass within GHF30. Silencing the putative endoxylanase by double-stranded RNA targeting of the CBM region resulted in an average decrease in infection of 60%, indicating that the gene is important for the nematode to complete its life cycle.

Expressed sequence tags of the peanut pod nematode Ditylenchus africanus: the first transcriptome analysis of an Anguinid nematode

In this study, 4847 expressed sequenced tags (ESTs) from mixed stages of the migratory plant-parasitic nematode Ditylenchus africanus (peanut pod nematode) were investigated. It is the first molecular survey of a nematode which belongs to the family of the Anguinidae (order Rhabditida, superfamily Sphaerularioidea). The sequences were clustered into 2596 unigenes, of which 43% did not show any homology to known protein, nucleotide, nematode EST or plant-parasitic nematode genome sequences. Gene ontology mapping revealed that most putative proteins are involved in developmental and reproductive processes. In addition unigenes involved in oxidative stress as well as in anhydrobiosis, such as LEA (late embryogenesis abundant protein) and trehalose-6-phosphate synthase were identified. Other tags showed homology to genes previously described as being involved in parasitism (expansin, SEC-2, calreticulin, 14-3-3b and various allergen proteins). In situ hybridization revealed that the expression of a putative expansin and a venom allergen protein was restricted to the gland cell area of the nematode, being in agreement with their presumed role in parasitism. Furthermore, seven putative novel candidate parasitism genes were identified based on the prediction of a signal peptide in the corresponding protein sequence and homologous ESTs exclusively in parasitic nematodes. These genes are interesting for further research and functional characterization. Finally, 34 unigenes were retained as good target candidates for future RNAi experiments, because of their nematode specific nature and observed lethal phenotypes of Caenorhabditis elegans homologs.

RNAi and functional genomics in plant parasitic nematodes

Plant nematology is currently undergoing a revolution with the availability of the first genome sequences as well as comprehensive expressed sequence tag (EST) libraries from a range of nematode species. Several strategies are being used to exploit this wealth of information. Comparative genomics is being used to explore the acquisition of novel genes associated with parasitic lifestyles. Functional analyses of nematode genes are moving toward larger scale studies including global transcriptome profiling. RNA interference (RNAi) has been shown to reduce expression of a range of plant parasitic nematode genes and is a powerful tool for functional analysis of nematode genes. RNAi-mediated suppression of genes essential for nematode development, survival, or parasitism is revealing new targets for nematode control. Plant nematology in the genomics era is now facing the challenge to develop RNAi screens adequate for high-throughput functional analyses.

Evolution of GHF5 endoglucanase gene structure in plant-parasitic nematodes: no evidence for an early domain shuffling event

Background

Endo-1,4-beta-glucanases or cellulases from the glycosyl hydrolase family 5 (GHF5) have been found in numerous bacteria and fungi, and recently also in higher eukaryotes, particularly in plant-parasitic nematodes (PPN). The origin of these genes has been attributed to horizontal gene transfer from bacteria, although there still is a lot of uncertainty about the origin and structure of the ancestral GHF5 PPN endoglucanase. It is not clear whether this ancestral endoglucanase consisted of the whole gene cassette, containing a catalytic domain and a carbohydrate-binding module (CBM, type 2 in PPN and bacteria) or only of the catalytic domain while the CBM2 was retrieved by domain shuffling later in evolution. Previous studies on the evolution of these genes have focused primarily on data of sedentary nematodes, while in this study, extra data from migratory nematodes were included.

Results

Two new endoglucanases from the migratory nematodes Pratylenchus coffeae and Ditylenchus africanus were included in this study. The latter one is the first gene isolated from a PPN of a different superfamily (Sphaerularioidea); all previously known nematode endoglucanases belong to the superfamily Tylenchoidea (order Rhabditida). Phylogenetic analyses were conducted with the PPN GHF5 endoglucanases and homologous endoglucanases from bacterial and other eukaryotic lineages such as beetles, fungi and plants. No statistical incongruence between the phylogenetic trees deduced from the catalytic domain and the CBM2 was found, which could suggest that both domains have evolved together. Furthermore, based on gene structure data, we inferred a model for the evolution of the GHF5 endoglucanase gene structure in plant-parasitic nematodes. Our data confirm a close relationship between Pratylenchus spp. and the root knot nematodes, while some Radopholus similis endoglucanases are more similar to cyst nematode genes.

Conclusion

We conclude that the ancestral PPN GHF5 endoglucanase gene most probably consisted of the whole gene cassette, i.e. the GHF5 catalytic domain and the CBM2, rather than that it evolved by domain shuffling. Our evolutionary model for the gene structure in PPN GHF5 endoglucanases implies the occurrence of an early duplication event, and more recent gene duplications at genus or species level.