Figures
Abstract
As the causal agent of pine wilt disease (PWD), the pine wood nematode (PWN), Bursaphelenchus xylophilus, causes huge economic losses by devastating pine forests worldwide. However, the pathogenesis-related genes of B. xylophilus are not well characterized. Thus, DNA microarrays were used to investigate differential gene expression in PWN where Pinus thunbergii was inoculated with nematodes, compared with those cultured on Botrytis cinerea. The microarrays comprised 31121 probes, 1310 (4.2%) of which were differentially regulated (changes of >2-fold, P < 0.01) in the two growth conditions. Of these 1310 genes, 633 genes were upregulated, whereas 677 genes were downregulated. Gene Ontology (GO) categories were assigned to the classes Cellular Component, Molecular Function, and Biological Process. The comparative gene expression analysis showed that a large number of the pathogenesis-related genes of B. xylophilus, such as pectate lyase genes, cytochrome P450s, UGTs, and ABC transporter genes, were highly expressed when B. xylophilus infected P. thunbergii. Annotation analysis indicated that these genes contributed to cell wall degradation, detoxification, and the reproduction process. The microarray results were validated using quantitative RT-PCR (qRT-PCR). The microarray data confirmed the specific expression of B. xylophilus genes during infection of P. thunbergii, which provides basic information that facilitates a better understanding of the molecular mechanism of PWD.
Citation: Qiu X, Wu X, Huang L, Tian M, Ye J (2013) Specifically Expressed Genes of the Nematode Bursaphelenchus Xylophilus Involved with Early Interactions with Pine Trees. PLoS ONE 8(10): e78063. https://doi.org/10.1371/journal.pone.0078063
Editor: Zhanjiang Liu, Auburn University, United States of America
Received: July 13, 2013; Accepted: September 16, 2013; Published: October 14, 2013
Copyright: © 2013 Ye et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the State Forest Administration of China (number 201204501), the project of National Science Research of Jiangsu University in China (number 11KJA220002), National Nature Science Foundation of China (number 31000303), the Natural Science Foundation of Jiangsu Province (number BK2010566). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Pine wilt disease (PWD) is one of the most serious diseases [1,2] that affects coniferous forests around the world and it is considered to be caused by pine wood nematodes (PWNs), Bursaphelenchus xylophilus (Steiner & Buhrer) Nickle [3,4]. In Asia, PWD was introduced into Japan during the early 20th century [5] and it spread subsequently to other countries, including China [6] and Korea [7], where pine trees have no natural resistance to B. xylophilus. Later, PWD spread to European countries, including Finland [8], Portugal [9], and Russia [10]. B. xylophilus is considered to be a native of North America [11], where the local forests are generally resistant or tolerant [12]. B. xylophilus causes severe economic, environmental, and social impacts [13,14] in non-indigenous areas and it has attracted attention because of its high pathogenicity in host trees. As the main pathogenic agent of PWD [15], PWNs lead to the death of pine trees within several months, especially between May and September. B. xylophilus is vectored by longhorn beetles in the genus Monochamus [16–19] and is transmitted via an unusual tree to tree infection route [20]. It is widely known that the introduction and expansion of B. xylophilus into non-native areas was mediated mainly by the international wood trade via the short- or long-distance transportation of pine wood, as well as packing materials and wood products infected by PWN [21–23]. Therefore, there is an urgent need to develop effective management strategies to control PWD [24].
PWN is a serious threat to forests in Asia and Europe, and the risk of this problem is likely to increase due to climate change. However, the pathogenic mechanism of PWD is not clear at present [25]. Little was known about the molecular pathogenicity of B. xylophilus until ten years ago, but much progress has been made recently since the development of biotechnology techniques [21]. High-throughput sequencing, which is a powerful method for gene research, has been used to screen resistant genes in trees infected by B. xylophilus [26]. Microsatellite markers have also been used to study the genetic diversity of PWN to understand its invasion route and the host colonization process [24].
Information about the pathogenesis-related genes of B. xylophilus is essential for understanding the pathogenic mechanism of PWD, but there has been little research in this area. Thus, we compared the differential gene expression of B. xylophilus in two growth conditions: growth on Botrytis cinerea and after inoculating Pinus thunbergii with PWNs. The goal of this study was to analyze the specifically expressed PWN genes involved with the early interactions between B. xylophilus and P. thunbergii, and to screen the pathogenesis-related genes of B. xylophilus using DNA microarrays. Quantitative RT-PCR (qRT-PCR) was used to validate the results obtained using the microarray assays.
Materials and Methods
B. Xylophilus growth conditions and sampling
P. thunbergii seedlings (2 years old) obtained from the greenhouse at Nanjing Forest University were transplanted into pots (30 cm in diameter, 25 cm in height) and maintained with a relative humidity of 70%. The seedlings (height c. 80 cm) were watered every other day and maintained with a photoperiod of 14 h day (25°C) and 10 h night (20°C). The highly virulent AMA3cl strain of PWN was used in the experiment, which was maintained by Lihua Zhu in our laboratory [27]. Two treatments were applied: (i) a suspension of 5000 nematodes (a mixture of adults and juveniles) was used to inoculate the fungus B. cinerea (the fungus had been incubated at 25°C for 6 days) on potato dextrose agar medium and grown at 25°C for a further 7 days; (ii) the same amount of nematodes was pipetted into wounds (2 cm in length) in P. thunbergii seedlings at about 50 cm above the soil level. The inoculated seedlings were cultivated in the greenhouse for 7 days at 25°C during the daytime and 20°C at night, with 70% humidity. The nematodes were separated from B. cinerea and P. thunbergii seedlings using a Baermann funnel. Next, B. xylophilus was collected by centrifugation at 3000 rpm for 1 min and frozen in liquid nitrogen, before further RNA isolation.
RNA isolation
The RNA was extracted from frozen nematodes using an RNAprep Kit (Tiangen, China) and purified further with an RNAclean Kit (Tiangen, China), according to the manufacturer’s protocol. The RNA was quantified at 260 nm using a spectrophotometer and examined by electrophoresis on a 1.5% agarose gel.
Microarray construction, hybridization, and data analysis
The microarray experiments were performed with the help of Shanghai Biotechnology Corporation. In total, 31121 oligonucleotides were synthesized based on the whole genome sequences (obtained from http://www.ncbi.nlm.nih.gov/genome) of B. xylophilus. A silane mixture was used to prepare hydrophobic glass slides, which contained exposed hydroxyl groups to facilitate nucleotide combining. The synthesized oligonucleotides were spotted robotically onto the hydrophobic surfaces of the glass slides. The first-strand cDNA was synthesized using a T7-oligo (dT) promoter primer. Subsequently, double-stranded cDNA was produced using a combination of Affinity Script RNase Block Mix, dNTP mix, DTT, and First Strand Buffer. The cRNA was generated from cDNAs using a Genechip IVT Labeling Kit (Affymetrix, USA) and purified with an RNeasy Mini Kit (Qiagen, Germany). The cRNAs were then hybridized at 65°C for 17 h with a microarray that contained 31121 probe sets. After hybridization, the slides were washed with GE wash buffer, according to the manufacturer’s instructions, and scanned using an Agilent microarray scanner (Cat#G2565CA, Agilent Technologies, Santa Clara, CA, USA). The data generated by the scanner were normalized using the Quantile algorithm, Gene Spring software 11.0 (Agilent Technologies, Santa Clara, CA, USA). The gene expression data were analyzed using a t-test (P < 0.01). Only genes with significantly different expression levels were screened for further analysis. The gene data analysis and functional annotation were performed using SBC Analysis System (http://www.sas.ebioservice.com). The results of microarray analysis were deposited in the NCBI database (Gene Expression Omnibus) and the accession number is GSE50481.
Quantitative real-time RT-PCR
Real-time quantitative RT-PCR (qRT-PCR) was used to verify the microarray results. Gene-specific primers (Table S1) were designed using Primer Premier 5 software. The actin gene of B. xylophilus was selected as the internal control. The first-strand cDNA was synthesized using a Prime Script 1st strand cDNA synthesis Kit (TaKaRa, Japan). The cDNA samples were diluted to 20 ng/µl. Real-time PCR was conducted using a 20 µl reaction volume, which contained 2 µl of template, 10 µl SYBR Premix Ex Tap, 0.4 µl ROX Reference Dye II, 0.4 µl forward primer, 0.4 µl reverse primer, and 6.8 µl ddH2O. PCR amplification was carried out using the following conditions: denaturation at 95°C for 30 s, followed by 40 cycles of amplification where each cycle comprised denaturation at 95°C for 5 s, and annealing and extension at 60°C for 34 s. The quantitative variations in the selected genes in the two growth conditions were evaluated using the relative quantification method (ΔΔ CT) [28].
Results
Screening of the specifically expressed genes of B. Xylophilus when cultured on B. cinerea and inoculated into P. thunbergii
B. xylophilus was separated from B. cinerea and P. thunbergii seedlings using Baermann funnels. To reduce the likelihood of biological errors, nematodes were collected from five inoculated P. thunbergii seedlings that exhibited similar symptoms and combined as one sample. There were three replicates for each treatment. Differences in the transcripts of B. xylophilus in the two different growth conditions were detected by genome-wide expression profiling using the Agilent B. xylophilus GeneChip microarray containing 31121 probe sets. Over 31000 probe signals were detected from the microarrays and 1310 (4.2%) B. xylophilus genes were differentially expressed (Figure 1). Of these 1310 genes, 633 genes were upregulated, whereas 677 genes were downregulated. Of the 633 upregulated genes, 569 genes were upregulated by 2- to 5-fold and 64 genes were upregulated >5-fold. Among the specifically expressed genes, we focused on those related to cell wall degradation, detoxification, and reproduction to identify pathogenesis-related genes.
The expression data for the 1310 genes that changed significantly were clustered using the SBC Analysis System (http://www.sas.ebioservice.com). The heat map shows the genes that were upregulated (red), downregulated (green), and no change (black).
Functional annotation
Functional annotation was performed to assign the genes of B. xylophilus with Gene Ontology (GO) terms. The main GO categories included Cellular Component, Molecular Function, and Biological Process. Based on the annotations (Figure 2) for Biological Process, 13.4% and 12.9% of the assignments belonged to the categories “Metabolic Process” (GO: 0008152) and “Cellular Process” (GO: 0009987), respectively, followed by “Developmental Process” (GO: 0032502, 10.5%) and “Multicellular Organismal Process” (GO: 0032501, 10.5%). In addition, the “Binding” (GO: 0005488) and “Catalytic Activity” (GO: 0003824) categories were prominent among the Molecular Function terms with 41.5% and 37.7% of the assignments, respectively, followed by the categories “Transporter Activity” (GO: 0005215, 6.2%), “Transcription Regulator Activity” (GO: 0030528, 3.5%), and “Electron Carrier Activity” (GO: 0009055, 3.5%). Furthermore, the categories of “Cell Part” (GO: 0044464) and “Cell” (GO: 0005623) had the same proportion, followed by “Organelle” (GO: 0043226) and “Macromolecular Complex” (GO: 0032991), which accounted for 9.1% and 3.4% of assignments. Interestingly, a recent study reported a similar classification of the annotated amino acid sequences in Pinus pinaster and Pinus pinea after inoculation with B. xylophilus [28]. In general, the aforementioned GO classes accounted for the majority of the specifically expressed genes.
The genes that changed significantly were divided into three functional sub-categories based on the Cellular Component (A), Molecular Function (B), and Biological Process (C) categories, according to the Gene Ontology (GO) principles.
The gene expression levels of B. xylophilus changed dramatically when the nematode was used to inoculate P. thunbergii compared with those grown on B. cinerea. We found that 520 genes were upregulated at least 3-fold, including 176 annotated genes. The GO analysis showed that these highly expressed genes (changes of >3-fold) were related mainly to metal ion binding, transferase activity, protein binding, nucleic acid binding, nucleotide binding, oxidoreductase activity, and hydrolase activity (Figure 3). The upregulated metal ion binding-related genes may contribute significantly to the enhancement of various enzyme activities and increase the adaptability of B. xylophilus. Interestingly, transferases related to metal ions were also upregulated. The protein binding factors associated with regulation probably activated the over-expression of the pathogenicity-related genes of B. xylophilus during the PWD process. The high expression levels of nucleic and nucleotide binding-related genes indicated that B. xylophilus can reproduce in P. thunbergii, which is the causal premise of PWD. In addition, a number of oxidoreductase and hydrolase genes, such as glucose oxidase, pectate lyase, and glutathione S-transferase genes, were detected by the microarrays. According to the functional annotation, we found that some of the cell wall degradation-related genes were upregulated significantly. These genes are considered to be key factors that allow B. xylophilus to invade its host.
The genes that changed significantly were designated as those where the expression level changed by more than 3-fold when Bursaphelenchus xylophilus was used to inoculate Pinus thunbergii compared with those grown on Botrytis cinerea.
qRT-PCR analysis
To validate the microarray results, 20 genes were selected at random from the highly expressed sequences and analyzed by qRT-PCR. We designed specific-primers for nine downregulated genes (Probe name: nl-p01087, nl-p02937, nl-p03084, nl-p04066, nl-p04652, nl-p15102, nl-p17373, nl-p20437, nl-p24677) and 11 upregulated genes (Probe name: nl-p03314, nl-p04906, nl-p10686, nl-p10889, nl-p12558, nl-p13660, nl-p14243, nl-p14336, nl-p22988, nl-p27181, nl-p29107). The expression patterns of the 20 candidate genes detected by qRT-PCR were similar to those detected by the microarray (Figure 4), which demonstrated the reliability of the microarray data.
Twenty genes were chosen at random from the highly expressed sequences and their expression levels were assessed by qRT-PCR. The numbers on the horizontal axis represent the probe name in the microarray, while those on the vertical axis are the genes that were upregulated or downregulated.
Discussion
Valuable information was found by identifying specifically expressed genes based on the differential gene expression of PWNs used to inoculate P. thunbergii compared with those grown on B. cinerea [26]. The plant cell wall is the primary barrier that a pathogen must penetrate [29]. The cell wall is comprised mainly of pectin, cellulose, and hemicellulose [30]. Cellulose is a polymer of (1–4)-linked β-D-glucose that is interlocked by hemicellulose to form a strong elastic network [31,32]. The network is embedded in pectin, which is the most complex component of the cell wall polysaccharides. The plant cell wall plays a critical role in resisting invasion by B. xylophilus. To invade the host successfully, B. xylophilus must produce a range of cell wall degradation enzymes that break down this barrier. Therefore, the high level expression of cell wall degradation-related genes is essential to allow B. xylophilus to infect its plant host.
Cell wall-degrading enzymes have been proposed to determine the pathogenicity of PWN. Plant parasitic nematodes have a wide range of cell wall-degrading enzymes, such as cellulase, polygalacturonase [33], pectate lyases [34], and xylanases [35], which contribute to the nematode’s ability to feed on plant tissues. Cellulase is one of the key enzymes associated with the pathogenicity of PWN [36]. Only the GHF5 and GHF45 families of cellulases have been found in nematodes and it has been proposed that they were acquired via horizontal gene transfer (HGT) events [37]. Mayer [38] studied the evolution of the cellulase genes of ten nematode species based on a phylogenetic framework and demonstrated that nematode species with a cellulase gene acquired via HGT exhibited cellulase activities. It has been suggested that the cellulase gene was functionally integrated into the nematode’s genome. In the present study, a number of cellulose degradation-related genes were detected by our microarray analysis, such as beta-1,3-endoglucanase, beta-1,4-endoglucanase, and cellulase (Table 1). The expression of beta-1,4-endoglucanase genes did not change greatly (0.9-fold). However, the expression levels of cellulase genes were downregulated significantly. During the early stage of infection, it is possible that cell wall components other than cellulose are degraded primarily by B. xylophilus.
| Genes | Function | Fold change |
|---|---|---|
| cht-1 | Chitinase activity | 0.1354 |
| BxPel2 | Pectate lyase activity | 6.6523 |
| Bx-13g-1 | Beta-1,3-endoglucanase activity | 0.1269 |
| Bx-eng-1 | Beta-1,4-endoglucanase activity | 0.8897 |
| Bx-eng-2 | Beta-1,4-endoglucanase activity | 0.8897 |
| Bx-eng-3 | Beta-1,4-endoglucanase activity | 0.8897 |
| BX-C10 | Cellulase activity | 0.2666 |
| BX-C12 | Cellulase activity | 0.2666 |
| BX-C13 | Cellulase | 0.2666 |
Table 1. Putative cell wall degradation-related genes of Bursaphelenchus xylophilus.
Pectin-degrading enzymes, such as polygalacturonases and pectate lyases (PLs), are the primary cell-degrading enzymes secreted by pathogens, which can expose more polymers that need to be degraded by other cell-degrading enzymes [39,40]. PLs are thought to be important pathogenic factors for phytopathogens [41]. A pectate lyase 2（pel-2）gene was detected in B. xylophilus based on the microarray assays (Table 1) and its expression was upregulated >6-fold. The high expression of PLs may enhance the infection capacity of B. xylophilus and facilitate its migration through plant tissues. Kang et al. (2012) showed that the migration activity of B. xylophilus plays an important role in efficient parasitism [42]. The pel-2 gene has also been detected in potato cyst nematodes (PCNs) and experiments demonstrated that the silencing of pel-2 in PCN juveniles reduced the infection efficiency greatly [43].
As one of the primary natural biopolymers, chitin is a major component of fungal cell walls. B. xylophilus can secrete chitinase to digest fungal cell walls throughout its growth cycle [44]. The expression of a chitinase gene (cht-1) was downregulated dramatically when B. xylophilus was used to inoculate P. thunbergii compared with culture on B. cinerea, probably because of the change in the available sources that B. xylophilus could fed upon. The formation and degradation of chitin are essential for nematode oviposition and molting cycle development. Nematode metamorphosis is regulated by the release of ecdysone, which is related to the regulation of chitin synthesis. A previous report showed that ecdysone plays a role in regulating the expression levels of the CHS-1 and CHS-2 genes during Drosophila metamorphosis [45]. However, abnormal chitin synthesis could lead to disorders during nematode development. Several proteins, such as yeast CHS4p and CHS5p, are related to the regulation of chitin synthesis in yeast [46]. Interestingly, the expression of chitinase is restricted to the period of molt and pupation, and different growth conditions [47]. Our microarray results showed that the expression levels of ecdysteroid-related genes were upregulated dramatically when P. thunbergii was inoculated with B. xylophilus. For example, the ptc-1 and rack-1 genes were upregulated by 4.2-fold and 3.6-fold, respectively.
The cytochrome p450 (CYP450) family includes enzymes that play important roles during the biotransformation of secondary metabolites [48]. Studies have shown that C. elegans CYP450s catalyze a series of exogenous and endogenous substrates, which cope with variation in the growth conditions [49]. This phenomenon may also occur in B. xylophilus. Pine trees generate numerous secondary metabolites, such as terpenoids and cyclic aromatics, to combat the invasion by B. xylophilus [26]. B. xylophilus must also mobilize defensive reactions to avoid damage by these complex compounds [50,51]. The presence of chemical compounds and plant secondary metabolites in the environment can induce the expression of CYP450 genes [52,53]. Based on complete genome sequencing data, 80 and 76 CYP450 genes were detected in C. elegans and B. xylophilus, respectively [54]. As shown in Table 2, the expression levels of seven CYP450 genes changed strikingly according to the microarray analysis. Of these, four genes were upregulated by at least 3-fold, whereas three genes were downregulated. It is surprising that the CYP450 genes of B. xylophilus not only generate enzymes to utilize secondary metabolites, but they also produce toxic metabolites that damage P. thunbergii [32]. Thus, the over-expression of CYP450 genes may enhance the pathogenicity of B. xylophilus and play a crucial role in the disease process [55]. However, the pathogenic role of the CYP450 genes in B. xylophilus remains unknown [56]. The gene function annotations showed that the W01A11.1 gene of B. xylophilus, which was upregulated by 6.3-fold, was a response to toxins, thereby suggesting that P. thunbergii generates a range of secondary metabolites as toxins that resist invasion by B. xylophilus. Furthermore, the accumulation of secondary metabolites produced by P. thunbergii may induce the over-expression of B. xylophilus CYP450 genes [57,58].
| Genes | Function | Fold change |
|---|---|---|
| ugt-49 | UDP-Glucuronosyl Transferase | 3.1834 |
| ugt-47 | UDP-Glucuronosyl Transferase | 3.5447 |
| ugt-59 | UDP-Glucuronosyl Transferase | 3.3254 |
| ugt-54 | UDP-Glucuronosyl Transferase | 7.509 |
| cyp-31A3 | Monooxygenase activity | 0.3014 |
| cyp-25A5 | Monooxygenase activity | 0.1778 |
| cyp-13A11 | Monooxygenase activity | 3.0041 |
| cyp-33C1 | Monooxygenase activity | 0.3102 |
| cyp-33C4 | Monooxygenase activity | 4.4049 |
| cyp-33C9 | Monooxygenase activity | 6.2141 |
| cyp-33D3 | Monooxygenase activity | 3.2066 |
| abtm-1 | ABC Transporter | 3.8127 |
| gst-33 | Glutathione S-Transferase | 3.3241 |
| W01A11.1 | response to toxin | 6.3557 |
| Y52B11A.8 | lipid catabolic process | 5.1726 |
Table 2. Putative detoxification-related genes of Bursaphelenchus xylophilus.
At present, the detoxification process of B. xylophilus is divided into three phases [54]. First, CYPs are the primary proteins in the first phase, which provide the enzyme substrates for the next stage. Second, the glutathione S-transferases (GSTs) and UDP-glucuronosyl transferases (UGTs) are essential in the second phase. A series of detoxification reactions occurs during this stage, which produce a high level of efflux. As shown in Table 2, one GST and four UGTs were found to be upregulated by at least 3-fold. In the last phase, ATP-binding cassette (ABC) transporters are the main group responsible for the efflux of detoxified molecules. One ABC transporter was found to be highly upregulated according to the microarray analysis.
In particular, the expression levels of genes related to detoxification were upregulated significantly, which suggests that B. xylophilus enhanced its gene expression in response to the secondary metabolites produced by pine trees. One gene, Y52B11A.8, which is related to the lipid catabolism process, was detected by the microarray. The expression of Y52B11A.8 was upregulated by more than 5-fold. This indicated that lipids, which are likely to have nematocidal activity, were produced by pine trees after invasion by B. xylophilus.
The reproductive ability of B. xylophilus in pine trees is a crucial causal factor in PWD [56]. Many reports have suggested that the virulence level of B. xylophilus is associated with its reproductive ability in in vitro or in vivo conditions [59,60]. When two PWN isolates (virulent and avirulent) were inoculated into pine trees, the PWN population of the virulent isolate was greater than that of the avirulent isolate. The rate of population increase by the virulent isolate was also greater than that of the avirulent isolate in PWD-killed seedlings [60–62]. The translocation of the virulent isolate into the xylem resin canals and phloem tissues was faster than that of the avirulent isolate after they were inoculated into separate pine trees [63]. Aikawa et al. [64] reported that the number of third- and fourth-stage dispersal juveniles (JIIIs and JIVs) of a virulent isolate from adult beetles was much higher than that of the avirulent isolates. This indicated that the fertility of the virulent isolate was higher than that of the avirulent isolate. Previous studies have shown that the virulence of nematodes is closely correlated with their reproductive ability [65]. The reproductive ability of B. xylophilus in host trees may be an important marker of its virulence, which could be used to evaluate the virulence level of B. xylophilus. A series of reproduction-related genes were detected by the microarray analysis (Table 3). The functions of reproduction-related genes were involved with egg hatching, germ cell development and meiosis, and the positive regulation of the growth rate, etc. The results indicated that the expression levels of the reproduction-related genes of B. xylophilus, such as rpl-24.2, R11A8.2, set-16, and F25H5.6, were upregulated dramatically by growth in P. thunbergii compared with those when grown on B. cinerea.
| Genes | Function | Fold change |
|---|---|---|
| vps-34 | nematode larval development | 4.0836 |
| ngp-1 | genitalia development | 3.5102 |
| F25H5.6 | morphogenesis of an epithelium | 4.3557 |
| rpl-24.2 | body morphogenesis | 8.386 |
| vps-32.1 | embryonic development ending in birth or egg hatching | 3.3897 |
| ptc-1 | oviposition and molting cycle | 4.2054 |
| rsr-2 | nematode larval development | 4.0309 |
| rpl-16 | structural constituent of ribosome | 3.3006 |
| set-16 | nematode larval development | 5.4682 |
| hoe-1 | germ cell development and meiosis | 3.1255 |
| rack-1 | oviposition and genitalia development | 3.6128 |
| R11A8.2 | positive regulation of growth rate | 6.194 |
| his-66 | nucleosome assembly | 3.6435 |
| mua-1 | nematode larval development | 4.0318 |
Table 3. Putative reproduction-related genes of Bursaphelenchus xylophilus.
In the early stages of PWD, ligno-suberization and wound periderm were observed around the cortex resin canal. In the advanced stage, the number of nematodes increased dramatically and the cambium of pine trees was destroyed [25,66]. The results of previous studies have shown that B. xylophilus feeds mainly on the xylem ray parenchyma cells of pine trees [67]. The rapid reproduction rate of B. xylophilus in P. thunbergii probably leads to more serious destruction of the cortex resin canal. However, more of the cortex parenchyma was fed on by B. xylophilus, which induces rapid increases in the cavitation area. The accumulation of the toxic substances produced by B. xylophilus also aggravated the wilting of needles. Thus, a high frequency of nematode reproduction is probably an important pathogenicity factor related to PWD.
PWD is a complex disease and its mechanism remains unclear. The specific genes related to cellulose degradation, hydrolase, and detoxification by B. xylophilus may play key roles in this disease. Thus, the mechanism of PWD may be involved with variation in the transcript levels of B. xylophilus genes. Our microarray analysis demonstrated that the expression levels of B. xylophilus genes changed significantly after it was used to inoculate P. thunbergii seedlings.
The present study investigated the differential gene expression of B. xylophilus cultured on B. cinerea and when used to inoculate P. thunbergii. The results showed that the gene expression of B. xylophilus changed significantly and a number of genes were upregulated dramatically. We found that some of the upregulated genes were related to cell wall degradation and reproduction, which suggests that B. xylophilus not only succeeded in invading P. thunbergii, but it also reproduced within P. thunbergii. In addition, we detected several toxin-related genes, which indicate that B. xylophilus may produce toxins and utilize the secondary metabolites of its host, and this process is possibly involved with the wilting of pine trees. These data facilitate a better understanding of the molecular mechanism of PWD and they may help to develop effective control strategies to combat B. xylophilus. However, further investigations are required to understand the functions of these specifically expressed genes.
Supporting Information
Table S1.
Primers used in the quantitative real-time PCR (qRT-PCR) analysis.
https://doi.org/10.1371/journal.pone.0078063.s001
(DOCX)
Acknowledgments
We are grateful to Professor Li-Hua Zhu for kindly providing the B. xylophilus isolate and we thank Professor De-Wei Li for critically reviewing the manuscript. We would also like to thank Shanghai Biotechnology Corporation for their helpful technical support.
Author Contributions
Conceived and designed the experiments: JRY LH. Performed the experiments: XWQ LH XQW. Analyzed the data: XWQ. Contributed reagents/materials/analysis tools: JRY XWQ MQT. Wrote the manuscript: XWQ LH.
References
- 1. Yang BJ, Wang Q (1989) Distribution of the pinewood nematode in China and susceptibility of some Chinese and exotic pines to the nematode. Can J Forest Res 19: 1527-1530. doi:https://doi.org/10.1139/x89-232.
- 2. Zhang F, Kajiwara J, Mori Y, Ohira M, Tsutsumi Y et al. (2013) Metabolites from resistant and susceptible Pinus thunbergii after inoculation with pine wood nematode. Am J Plant Sci 4: 512-518.
- 3. Steiner G, Buhrer EM (1934) Aphelenchoides xylophilus n. sp., a nematode associated with blue-stain and other fungi in timber. J Agr Res 48: 949-951.
- 4. Valadas V, Laranjo M, Barbosa P, Espada M, Mota M et al. (2012) The pine wood nematode, Bursaphelenchus xylophilus, in Portugal: possible introductions and spread routes of a serious biological invasion revealed by molecular methods. Nematology 14: 899-911. doi:https://doi.org/10.1163/156854112X632673.
- 5. Mamiya Y, Kiyohara T (1972) Description of Bursaphelenchus lignicolus n. sp. (Nematoda: Aphelenchoididae) from pine wood and histopathology of nematode-infested trees. Nematologica 18: 120-124. doi:https://doi.org/10.1163/187529272X00296.
- 6. Cheng H, Lin M, Li W, Fang Z (1983) The occurrence of a pine wilting disease caused by a nematode found in Nanjing. Forest Pest Dis 4: 1-5.
- 7. Yi CK, Byun B, Park J, Yang S, Chang K (1989) First finding of the pine wood nematode, Bursaphelenchus xylophilus (Steiner et Buhrer) Nickle and its insect vector in Korea. Res Rep Forest Res I 38: 141-149.
- 8. Tomminen J, Nuorteva M, Pulkkinen M, Väkevä J (1989) Occurrence of the nematode Bursaphelenchus mucronatus Mamiya & Enda 1979 (Nematoda: Aphelenchoididae) in Finland. Silva Fenn 23: 271-277.
- 9. Mota MM, Braasch H, Bravo MA, Penas AC, Burgermeister W et al. (1999) First report of Bursaphelenchus xylophilus in Portugal and in Europe. Nematology 1: 7-8.
- 10. Braasch H, Tomiczek C, Metge K, Hoyer U, Burgermeister W et al. (2001) Records of Bursaphelenchus spp. (Nematoda, Parasitaphelenchidae) in coniferous timber imported from the Asian part of Russia. Forest Pathol 31: 129-140. doi:https://doi.org/10.1046/j.1439-0329.2001.00233.x.
- 11. Kanzaki N, Futai K (2002) A PCR primer set for determination of phylogenetic relationships of Bursaphelenchus species within the xylophilus group. Nematology 4: 35-41. doi:https://doi.org/10.1163/156854102760082186.
- 12. Rutherford T, Mamiya Y, Webster J (1990) Nematode-induced pine wilt disease: factors influencing its occurrence and distribution. Forest Sci 36: 145-155.
- 13. Taga Y, Goto S, Matsunaga K, Togashi K (2011) Temporal changes in characteristics of populations originating from interbreeding between Bursaphelenchus xylophilus and B. mucronatus. Nematology 13: 701-712. doi:https://doi.org/10.1163/138855410X543166.
- 14. Soliman T, Mourits MC, van der Werf W, Hengeveld GM, Robinet C et al. (2012) Framework for modelling economic impacts of invasive species, applied to pine wood nematode in Europe. PLOS ONE 7: e45505. doi:https://doi.org/10.1371/journal.pone.0045505. PubMed: 23029059.
- 15. Kiyohara T, Tokushige Y (1971) Inoculation experiments of a nematode, Bursaphelenchus sp., onto pine trees. J Jpn Forestr Soc 53: 210-218.
- 16. Linit MJ (1988) Nemtaode-vector relationships in the pine wilt disease system. J Nematol 20: 227–235. PubMed: 19290206.
- 17. Sousa E, Bravo MA, Pires J, Naves P, Penas AC et al. (2001) Bursaphelenchus xylophilus (Nematoda; Aphelenchoididae) associated with Monochamus galloprovincialis (Coleoptera; Cerambycidae) in Portugal. Nematology 3: 89-91. doi:https://doi.org/10.1163/156854101300106937.
- 18. Naves P, Kenis M, Sousa E (2005) Parasitoids associated with Monochamus galloprovincialis (Oliv.) (Coleoptera: Cerambycidae) within the pine wilt nematode-affected zone in Portugal. J Pest Sci 78: 57-62. doi:https://doi.org/10.1007/s10340-004-0068-z.
- 19. Togashi K, Shigesada N (2006) Spread of the pinewood nematode vectored by the Japanese pine sawyer: modeling and analytical approaches. Popul Ecol 48: 271-283. doi:https://doi.org/10.1007/s10144-006-0011-7.
- 20. Li H, Shen P, Fu P, Lin M, Moens M (2007) Characteristics of the emergence of Monochamus alternatus, the vector of Bursaphelenchus xylophilus (Nematoda: Aphelenchoididae), from Pinus thunbergii logs in Nanjing, China, and of the transmission of the nematodes through feeding wounds. Nematology 9: 807-816.
- 21. Jones JT, Moens M, Mota M, Li H, Kikuchi T (2008) Bursaphelenchus xylophilus: opportunities in comparative genomics and molecular host–parasite interactions. Mol Plant Pathol 9: 357-368. doi:https://doi.org/10.1111/j.1364-3703.2007.00461.x. PubMed: 18705876.
- 22. Robinet C, Roques A, Pan H, Fang G, Ye J et al. (2009) Role of human-mediated dispersal in the spread of the pinewood nematode in China. PLOS ONE 4: e4646. doi:https://doi.org/10.1371/journal.pone.0004646. PubMed: 19247498.
- 23. Robinet C, Van Opstal N, Baker R, Roques A (2011) Applying a spread model to identify the entry points from which the pine wood nematode, the vector of pine wilt disease, would spread most rapidly across Europe. Biol Invasions 13: 2981-2995. doi:https://doi.org/10.1007/s10530-011-9983-0.
- 24. Mallez S, Castagnone C, Espada M, Vieira P, Eisenback JD et al. (2013) First insights into the genetic diversity of the pinewood nematode in its native area using new polymorphic microsatellite loci. PLOS ONE 8: e59165. doi:https://doi.org/10.1371/journal.pone.0059165. PubMed: 23554990.
- 25. Wang Z, Wang CY, Fang ZM, Zhang DL, Liu L et al. (2010) Advances in research of pathogenic mechanism of pine wilt disease. Afr J Microbiol Res 4: 437-442.
- 26. Santos CS, Pinheiro M, Silva AI, Egas C, Vasconcelos MW (2012) Searching for resistance genes to Bursaphelenchus xylophilus using high throughput screening. BMC Genomics 13: 599. doi:https://doi.org/10.1186/1471-2164-13-599. PubMed: 23134679.
- 27. Zhu L-h, Ye J, Negi S, Xu X-l, Wang Z-l et al. (2012) Pathogenicity of aseptic Bursaphelenchus xylophilus. PLOS ONE 7: e38095. doi:https://doi.org/10.1371/journal.pone.0038095. PubMed: 22662271.
- 28. Revel AT, Talaat AM, Norgard MV (2002) DNA microarray analysis of differential gene expression in Borrelia burgdorferi, the Lyme disease spirochete. Proc Natl Acad Sci USA 99: 1562-1567. doi:https://doi.org/10.1073/pnas.032667699. PubMed: 11830671.
- 29. Taiz L (1984) Plant cell expansion: regulation of cell wall mechanical properties. Annu Rev Plant Physiol 35: 585–657. doi:https://doi.org/10.1146/annurev.arplant.35.1.585.
- 30. Gilbert HJ (2010) The biochemistry and structural biology of plant cell wall deconstruction. Plant Physiol 153: 444-455. doi:https://doi.org/10.1104/pp.110.156646. PubMed: 20406913.
- 31. Carpita NC, Gibeaut DM (1993) Structural models of primary cell walls in flowering plants: consistency of molecular structure with the physical properties of the walls during growth. Plant J 3: 1-30. doi:https://doi.org/10.1046/j.1365-313X.1993.00999.x. PubMed: 8401598.
- 32. Cosgrove DJ (2005) Growth of the plant cell wall. Nat Rev Mol Cell Biol 6: 850-861. doi:https://doi.org/10.1038/nrm1746. PubMed: 16261190.
- 33. Jaubert S, Laffaire J-B, Abad P, Rosso M-N (2002) A polygalacturonase of animal origin isolated from the root-knot nematode Meloidogyne incognita. FEBS Lett 522: 109-112. doi:https://doi.org/10.1016/S0014-5793(02)02906-X. PubMed: 12095628.
- 34. Popeijus H, Overmars H, Jones J, Blok V, Goverse A et al. (2000) Enzymology: degradation of plant cell walls by a nematode. Nature 406: 36-37. doi:https://doi.org/10.1038/35017641. PubMed: 10894530.
- 35.
Mitreva-Dautova M, Roze E, Overmars H, de Graaff L, Schots A et al. (2006) A symbiont-independent endo-1, 4-β-xylanase from the plant-parasitic nematode Meloidogyne incognita. Mol Plant Microbe In 19: 521-529.
- 36. Tardy F, Nasser W, Robert-Baudouy J, Hugouvieux-Cotte-Pattat N (1997) Comparative analysis of the five major Erwinia chrysanthemi pectate lyases: enzyme characteristics and potential inhibitors. J Bacteriol 179: 2503-2511.
- 37. Kikuchi T, Jones JT, Aikawa T, Kosaka H, Ogura N (2004) A family of glycosyl hydrolase family 45 cellulases from the pine wood nematode Bursaphelenchus xylophilus. FEBS Lett 572: 201-205. doi:https://doi.org/10.1016/j.febslet.2004.07.039. PubMed: 15304348.
- 38. Mayer WE, Schuster LN, Bartelmes G, Dieterich C, Sommer RJ (2011) Horizontal gene transfer of microbial cellulases into nematode genomes is associated with functional assimilation and gene turnover. BMC Evol Biol 11: 13. doi:https://doi.org/10.1186/1471-2148-11-13. PubMed: 21232122.
- 39. Jones T, Anderson A, Albersheim P (1972) Host-pathogen interactions IV. Studies on the polysaccharide-degrading enzymes secreted by Fusarium oxysporum f. sp. lycopersici. Physiol Plant Pathol 2: 153-166. doi:https://doi.org/10.1016/0048-4059(72)90023-9.
- 40. Tomassini A, Sella L, Raiola A, D’Ovidio R, Favaron F (2009) Characterization and expression of Fusarium graminearum endopolygalacturonases in vitro and during wheat infection. Plant Pathol 58: 556-564. doi:https://doi.org/10.1111/j.1365-3059.2008.02019.x.
- 41. Yakoby N, Beno-Moualem D, Keen NT, Dinoor A, Pines O et al. (2001) Colletotrichum gloeosporioides pelB is an important virulence factor in avocado fruit-fungus interaction. Mol Plant Microbe In 14: 988-995.
- 42. Kang JS, Koh YH, Moon YS, Lee SH (2012) Molecular properties of a venom allergen-like protein suggest a parasitic function in the pinewood nematode Bursaphelenchus xylophilus. Int J Parasitol 42: 63-70. doi:https://doi.org/10.1016/j.ijpara.2011.10.006. PubMed: 22142561.
- 43. Vanholme B, Van Thuyne W, Vanhouteghem K, De Meutter J, Cannoot B et al. (2007) Molecular characterization and functional importance of pectate lyase secreted by the cyst nematode Heterodera schachtii. Mol Plant Pathol 8: 267-278. doi:https://doi.org/10.1111/j.1364-3703.2007.00392.x. PubMed: 20507498.
- 44. Merzendorfer H, Zimoch L (2003) Chitin metabolism in insects: structure, function and regulation of chitin synthases and chitinases. J Exp Biol 206: 4393-4412. doi:https://doi.org/10.1242/jeb.00709. PubMed: 14610026.
- 45. Gagou ME, Kapsetaki M, Turberg A, Kafetzopoulos D (2002) Stage-specific expression of the chitin synthase DmeChSA and DmeChSB genes during the onset of Drosophila metamorphosis. Insect Biochem Mol Biol 32: 141-146. doi:https://doi.org/10.1016/S0965-1748(01)00101-1. PubMed: 11755055.
- 46. Ono N, Yabe T, Sudoh M, Nakajima T, Yamada-Okabe T et al. (2000) The yeast Chs4 protein stimulates the trypsin-sensitive activity of chitin synthase 3 through an apparent protein–protein interaction. Microbiology 146: 385-391. PubMed: 10708377.
- 47. Lemos FJ, Secundino NF, Páscoa V, Pereira ST, Pimenta PF (2002) Presence of chitinase and beta-N-acetylglucosaminidase in the Aedes aegypti: a chitinolytic system involving peritrophic matrix formation and degradation. Insect Biochem Mol Biol 32: 1723-1729. doi:https://doi.org/10.1016/S0965-1748(02)00112-1. PubMed: 12429124.
- 48. Urlacher VB, Girhard M (2012) Cytochrome P450 monooxygenases: an update on perspectives for synthetic application. Trends Biotechnol 30: 26-36. doi:https://doi.org/10.1016/j.tibtech.2011.06.012. PubMed: 21782265.
- 49. Nebert DW, Russell DW (2002) Clinical importance of the cytochromes P450. Lancet 360: 1155-1162. doi:https://doi.org/10.1016/S0140-6736(02)11203-7. PubMed: 12387968.
- 50. Langenheim JH (1994) Higher plant terpenoids: a phytocentric overview of their ecological roles. J Chem Ecol 20: 1223-1280. doi:https://doi.org/10.1007/BF02059809.
- 51. Takeuchi Y, Kanzaki N, Futai K (2006) Volatile compounds in pine stands suffering from pine wilt disease: qualitative and quantitative evaluation. Nematology 8: 869-879. doi:https://doi.org/10.1163/156854106779799295.
- 52. Gonzalez FJ, Nebert DW (1990) Evolution of the P450 gene superfamily: animal-plant ‘warfare’, molecular drive and human genetic differences in drug oxidation. Trends Genet 6: 182-186. doi:https://doi.org/10.1016/0168-9525(90)90174-5. PubMed: 2196721.
- 53. Nelson DR, Kamataki T, Waxman DJ, Guengerich FP, Estabrook RW et al. (1993) The P450 superfamily: update on new sequences, gene mapping, accession numbers, early trivial names of enzymes, and nomenclature. DNA Cell Biol 12: 1-51. doi:https://doi.org/10.1089/dna.1993.12.1. PubMed: 7678494.
- 54. Kikuchi T, Cotton JA, Dalzell JJ, Hasegawa K, Kanzaki N et al. (2011) Genomic insights into the origin of parasitism in the emerging plant pathogen Bursaphelenchus xylophilus. PLOS Pathog 7: e1002219. PubMed: 21909270.
- 55. Menzel R, Bogaert T, Achazi R (2001) A systematic gene expression screen of Caenorhabditis elegans Cytochrome P450 genes reveals CYP35 as strongly xenobiotic inducible. Arch Biochem Biophys 395: 158-168. doi:https://doi.org/10.1006/abbi.2001.2568. PubMed: 11697852.
- 56. Kulas J, Schmidt C, Rothe M, Schunck W-H, Menzel R (2008) Cytochrome P450-dependent metabolism of eicosapentaenoic acid in the nematode Caenorhabditis elegans. Arch Biochem Biophys 472: 65-75. doi:https://doi.org/10.1016/j.abb.2008.02.002. PubMed: 18282462.
- 57. Nebert DW, Dieter MZ (2000) The evolution of drug metabolism. Pharmacology 61: 124-135.
- 58. Schuetz EG (2001) Induction of cytochromes P450. Curr Drug Metab 2: 139-147. doi:https://doi.org/10.2174/1389200013338595. PubMed: 11469722.
- 59. Shinya R, Takeuchi Y, Ichimura K, Takemoto S, Futai K (2012) Establishment of a set of inbred strains of the pine wood nematode, Bursaphelenchus xylophilus (Aphelenchida: Aphelenchoididae), and evidence of their varying levels of virulence. Appl Entomol Zool 47: 341-350. doi:https://doi.org/10.1007/s13355-012-0124-8.
- 60. Aikawa T, Kikuchi T (2007) Estimation of virulence of Bursaphelenchus xylophilus (Nematoda: Aphelenchoididae) based on its reproductive ability. Nematology 9: 371-377. doi:https://doi.org/10.1163/156854107781352007.
- 61. Wang Y, Yamada T, Sakaue D, Suzuki K (2005) Variations in life history parameters and their influence on rate of population increase of different pathogenic isolates of the pine wood nematode, Bursaphelenchus xylophilus. Nematology 7: 459-467. doi:https://doi.org/10.1163/156854105774355545.
- 62. Aikawa T, Kikuchi T, Kosaka H (2006) Population structure of Bursaphelenchus xylophilus within single Pinus thunbergii trees inoculated with two nematode isolates. Forest Pathol 36: 1-13. doi:https://doi.org/10.1111/j.1439-0329.2006.00428.x.
- 63. Ichihara Y, Fukuda K, Suzuki K (2000) Early symptom development and histological changes associated with migration of Bursaphelenchus xylophilus in seedling tissues of Pinus thunbergii. Plant Dis 84: 675-680. doi:https://doi.org/10.1094/PDIS.2000.84.6.675.
- 64. Aikawa T, Togashi K, Kosaka H (2003) Different developmental responses of virulent and avirulent isolates of the pinewood nematode, Bursaphelenchus xylophilus (Nematoda: Aphelenchoididae), to the insect vector, Monochamus alternatus (Coleoptera: Cerambycidae). Environ Entomol 32: 96-102. doi:https://doi.org/10.1603/0046-225X-32.1.96.
- 65. Mota MM, Takemoto S, Takeuchi Y, Hara N, Futai K (2006) Comparative studies between Portuguese and Japanese isolates of the pinewood nematode, Bursaphelenchus xylophilus. J Nematol 38: 429–433. PubMed: 19259459.
- 66. Fukuda K (1997) Physiological process of the symptom development and resistance mechanism in pine wilt disease. J Forest Res 2: 171-181. doi:https://doi.org/10.1007/BF02348216.
- 67. Fukuda K, Utsuzawa S, Sakaue D (2007) Correlation between acoustic emission, water status and xylem embolism in pine wilt disease. Tree Physiol 27: 969-976. doi:https://doi.org/10.1093/treephys/27.7.969. PubMed: 17403649.