Characterization and Functional Importance of Two Glycoside Hydrolase Family 16 Genes from the Rice White Tip Nematode Aphelenchoides besseyi

Simple Summary

The rice white tip nematode Aphelenchoides besseyi is a plant parasite but can also feed on fungi if this alternative nutrient source is available. Glucans are a major nutrient source found in fungi, and β-linked glucans from fungi can be hydrolyzed by β-glucanases from the glycoside hydrolase family 16 (GH16). The GH16 family is abundant in A. besseyi, but their functions have not been well studied, prompting the analysis of two GH16 members (AbGH16-1 and AbGH16-2). AbGH16-1 and AbGH16-2 are most similar to GH16s from fungi and probably originated from fungi via a horizontal gene transfer event. These two genes are important for feeding on fungi: transcript levels increased when cultured with the fungus Botrytis cinerea, and the purified AbGH16-1 and AbGH16-2 proteins inhibited the growth of B. cinerea. When AbGH16-1 and AbGH16-2 expression was silenced, the reproduction ability of A. besseyi was reduced. These findings have proved for the first time that GH16s contribute to the feeding and reproduction of A. besseyi, which thus provides novel insights into how plant-parasitic nematodes can obtain nutrition from sources other than their plant hosts.

Abstract

The glycoside hydrolase family 16 (GH16) is widely found in prokaryotes and eukaryotes, and hydrolyzes the β-1,3(4)-linkages in polysaccharides. Notably, the rice white tip nematode Aphelenchoides besseyi harbors a higher number of GH16s compared with other plant-parasitic nematodes. In this work, two GH16 genes, namely AbGH16-1 and AbGH16-2, were isolated and characterized from A. besseyi. The deduced amino acid sequences of AbGH16-1 and AbGH16-2 contained an N-terminal signal peptide and a fungal Lam16A glucanase domain. Phylogenetic analysis revealed that AbGH16-1 and AbGH16-2 clustered with ascomycete GH16s, suggesting AbGH16-1 and AbGH16-2 were acquired by horizontal gene transfer from fungi. In situ hybridization showed that both AbGH16-1 and AbGH16-2 were specifically expressed in the nematode gonads, correlating with qPCR analysis that showed the high transcript levels of the two genes in the female nematodes. AbGH16-1 and AbGH16-2 were also significantly induced in nematodes feeding on Botrytis cinerea. Characterization of the recombinant protein showed AbGH16-1 and AbGH16-2 displayed pronounced inhibition of both conidial germination and germ tube elongation of B. cinerea. In addition, silencing of AbGH16-1 and AbGH16-2 by RNA interference significantly decreased the reproduction ability of A. besseyi and had a profound impact on the development process of offspring in this nematode. These findings have firstly proved that GH16s may play important roles in A.besseyi feeding and reproduction on fungi, which thus provides novel insights into the function of GH16s in plant-parasitic nematodes.

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.

The Transcriptomes of Xiphinema index and Longidorus elongatus Suggest Independent Acquisition of Some Plant Parasitism Genes by Horizontal Gene Transfer in Early-Branching Nematodes

Nematodes have evolved the ability to parasitize plants on at least four independent occasions, with plant parasites present in Clades 1, 2, 10 and 12 of the phylum. In the case of Clades 10 and 12, horizontal gene transfer of plant cell wall degrading enzymes from bacteria and fungi has been implicated in the evolution of plant parasitism. We have used ribonucleic acid sequencing (RNAseq) to generate reference transcriptomes for two economically important nematode species, Xiphinema index and Longidorus elongatus, representative of two genera within the early-branching Clade 2 of the phylum Nematoda. We used a transcriptome-wide analysis to identify putative horizontal gene transfer events. This represents the first in-depth transcriptome analysis from any plant-parasitic nematode of this clade. For each species, we assembled ~30 million Illumina reads into a reference transcriptome. We identified 62 and 104 transcripts, from X. index and L. elongatus, respectively, that were putatively acquired via horizontal gene transfer. By cross-referencing horizontal gene transfer prediction with a phylum-wide analysis of Pfam domains, we identified Clade 2-specific events. Of these, a GH12 cellulase from X. index was analysed phylogenetically and biochemically, revealing a likely bacterial origin and canonical enzymatic function. Horizontal gene transfer was previously shown to be a phenomenon that has contributed to the evolution of plant parasitism among nematodes. Our findings underline the importance and the extensiveness of this phenomenon in the evolution of plant-parasitic life styles in this speciose and widespread animal phylum.

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.

Glycoside Hydrolase (GH) 45 and 5 Candidate Cellulases in Aphelenchoides besseyi Isolated from Bird's-Nest Fern

Five Aphelenchoides besseyi isolates collected from bird's-nest ferns or rice possess different parasitic capacities in bird's-nest fern. Two different glycoside hydrolase (GH) 45 genes were identified in the fern isolates, and only one was found in the rice isolates. A Abe GH5-1 gene containing an SCP-like family domain was found only in the fern isolates. Abe GH5-1 gene has five introns suggesting a eukaryotic origin. A maximum likelihood phylogeny revealed that Abe GH5-1 is part of the nematode monophyletic group that can be clearly distinguished from those of other eukaryotic and bacterial GH5 sequences with high bootstrap support values. The fern A. besseyi isolates were the first parasitic plant nematode found to possess both GH5 and GH45 genes. Surveying the genome of the five A. besseyi isolates by Southern blotting using an 834 bp probe targeting the GH5 domain suggests the presence of at least two copies in the fern-origin isolates but none in the rice-origin isolates. The in situ hybridization shows that the Abe GH5-1 gene is expressed in the nematode ovary and testis. Our study provides insights into the diversity of GH in isolates of plant parasitic nematodes of different host origins.

Calreticulin is required for responding to stress, foraging, and fertility in the white-tip nematode, Aphelenchoides besseyi

Calreticulin (CRT) regulates a wide array of cellular responses in physiological and pathological processes. A full-length cDNA-encoding CRT protein, namely AbCRT-1, was isolated from Aphelenchoides besseyi, an ectoparasitic plant nematode and the agent of white tip disease of rice. The deduced amino acid sequence of AbCRT-1 was highly homologous with other nematode CRTs, and showed the closest evolutionary relationship with BxCRT-1. In-situ hybridization showed that AbCRT-1 is specifically located in the oesophageal gland and gonads of A. besseyi, suggesting its potential role in parasitism and reproduction. Quantity real-time PCR analysis showed that AbCRT-1 is highly expressed in female nematodes but poorly expressed in eggs, juveniles, and male nematodes. Exposing the nematode to relatively low osmotic stress promotes the transcription of AbCRT-1 whereas extreme desiccation suppresses the transcription significantly. Nematodes in which AbCRT-1 mRNA level had been knocked down by soaking them in AbCRT-1 dsRNA solution distributed randomly and did not aggregate temporally, with a decreased capacity of food discernment. Thus the affected nematodes were markedly less fecund. These results demonstrate that AbCRT-1 is required in A. besseyi for responding to stress, foraging, and fertility.

Evolution of plant parasitism in the phylum Nematoda

Within the species-rich and trophically diverse phylum Nematoda, at least four independent major lineages of plant parasites have evolved, and in at least one of these major lineages plant parasitism arose independently multiple times. Ribosomal DNA data, sequence information from nematode-produced, plant cell wall-modifying enzymes, and the morphology and origin of the style(t), a protrusible piercing device used to penetrate the plant cell wall, all suggest that facultative and obligate plant parasites originate from fungivorous ancestors. Data on the nature and diversification of plant cell wall-modifying enzymes point at multiple horizontal gene transfer events from soil bacteria to bacterivorous nematodes resulting in several distinct lineages of fungal or oomycete-feeding nematodes. Ribosomal DNA frameworks with sequence data from more than 2,700 nematode taxa combined with detailed morphological information allow for explicit hypotheses on the origin of agronomically important plant parasites, such as root-knot, cyst, and lesion nematodes.

Distribution and evolution of glycoside hydrolase family 45 cellulases in nematodes and fungi

BACKGROUND

Horizontal gene transfer (HGT) has been suggested as the mechanism by which various plant parasitic nematode species have obtained genes important in parasitism. In particular, cellulase genes have been acquired by plant parasitic nematodes that allow them to digest plant cell walls. Unlike the typical glycoside hydrolase (GH) family 5 cellulase genes which are found in several nematode species from the order Tylenchida, members of the GH45 cellulase have only been identified in a cluster including the families Parasitaphelenchidae (with the pinewood nematode Bursaphelenchus xylophilus) and Aphelenchoididae, and their origins remain unknown.

RESULTS

In order to investigate the distribution and evolution of GH45 cellulase genes in nematodes and fungi we performed a wide ranging screen for novel putative GH45 sequences. This revealed that the sequences are widespread mainly in Ascomycetous fungi and have so far been found in a single major nematode lineage. Close relationships between the sequences from nematodes and fungi were found through our phylogenetic analyses. An intron position is shared by sequences from Bursaphelenchus nematodes and several Ascomycetous fungal species.

CONCLUSIONS

The close phylogenetic relationships and conserved gene structure between the sequences from nematodes and fungi strongly supports the hypothesis that nematode GH45 cellulase genes were acquired via HGT from fungi. The rapid duplication and turnover of these genes within Bursaphelenchus genomes demonstrate that useful sequences acquired via HGT can become established in the genomes of recipient organisms and may open novel niches for these organisms to exploit.

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.

Plant-parasitic nematode infections in rice: molecular and cellular insights

Being one of the major staple foods in the world, and an interesting model monocot plant, rice (Oryza sativa L.) has recently received attention from molecular nematologists studying the cellular and molecular aspects of the interaction between this crop and plant-parasitic nematodes. In this review, we highlight recent advances in this field, with a focus on the best-studied root-knot nematodes. Histological studies have revealed the cellular changes inside root-knot nematode-induced feeding sites, both in the compatible interaction with Oryza sativa and the incompatible interaction with the related species Oryza glaberrima. After comparing the published data from transcriptome analyses, mutant studies, and exogenous hormone applications, we provide a comprehensive model showing the role and interaction of plant hormone pathways in defense of this monocot crop against root nematodes, where jasmonate seems to play a key role. Finally, recent evidence indicates that effectors secreted from rice-infecting nematodes can suppress plant defense.