- Research article
- Open Access
Comparative transcriptomic analysis of candidate effectors to explore the infection and survival strategy of Bursaphelenchus xylophilus during different interaction stages with pine trees
BMC Plant Biology volume 21, Article number: 224 (2021)
The pine wood nematode (PWN), Bursaphelenchus xylophilus, is a devastating pathogen of many Pinus species in China. The aim of this study was to understand the interactive molecular mechanism of PWN and its host by comparing differentially expressed genes and candidate effectors from three transcriptomes of B. xylophilus at different infection stages.
In total, 62, 69 and 46 candidate effectors were identified in three transcriptomes (2.5 h postinfection, 6, 12 and 24 h postinoculation and 6 and 15 d postinfection, respectively). In addition to uncharacterized pioneers, other candidate effectors were involved in the degradation of host tissues, suppression of host defenses, targeting plant signaling pathways, feeding and detoxification, which helped B. xylophilus survive successfully in the host. Seven candidate effectors were identified in both our study and the B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was identified in all three transcriptomes. These common candidate effectors were upregulated at infection stages, and one of them suppressed pathogen-associated molecular pattern (PAMP) PsXEG1-triggered cell death in Nicotiana benthamiana.
The results indicated that B. xylophilus secreted various candidate effectors, and some of them continued to function throughout all infection stages. These various candidate effectors were important to B. xylophilus infection and survival, and they functioned in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, promoting detoxification and playing virulence functions). The present results provide valuable resources for in-depth research on the pathogenesis of B. xylophilus from the perspective of effectors.
As one of the most serious coniferous forest pathogens, the pine wood nematode (PWN), Bursaphelenchus xylophilus, is an essential threat to forest ecosystems worldwide, and pine wilt disease (PWD) caused by B. xylophilus has resulted in massive economic losses in Asian and European countries, especially China and Japan [1, 2]. In China, the occurrence of PWD has been distributed in 18 provinces (No. 4 bulletin in 2020, National Forestry and Grassland Administration, China).
B. xylophilus is a migratory phytoparasitic nematode that has a unique feeding strategy, which includes phytophagous and mycetophagous stages, enabling it to reproduce and survive in host pines. Due to the complexity of the interaction between B. xylophilus and its host, the mechanisms of pathogenesis remain unclear.
Effectors are key elements in the virulence of various pathogens and parasites (including fungi, oomycetes, bacteria and plant-parasitic nematodes [PPNs]) against plants [3,4,5,6]. The effectors of PPNs are secreted into host plant tissues and facilitate invasion and migration to modulate the host immune system . To determine the functions of effectors of PPNs in the interaction process between PPNs and their hosts, it is essential to screen and identify effectors.
Compared to microarrays and expressed sequence tags (ESTs), RNA-seq allows simultaneous transcript discovery and abundance estimation as well as identification of associated molecular cellular pathways and effectors of pathogens secreted during the infection process . Many pathogen effectors have been identified by this reliable method. For example, putative effector proteins of the Gymnosporangium yamadae and G. asiaticumrice rust species, the Heterodera avenae cereal cyst nematode, the Hirschmanniella oryzae root nematode and the Pratylenchus penetrans root lesion nematode that may alter host defense mechanisms have been screened by RNA-seq [9,10,11,12]. At present, transcriptome sequencing has also been used to identify differentially expressed genes (DEGs) when B. xylophilus enters the initial phytophagous phase (2.5 h postinoculation) . Espada et al. predicted candidate B. xylophilus effectors by comparing transcriptomes between the mycetophagous and middle and later phytophagous parasitic stages (6 and 15 d postinoculation) . Moreover, in our recent study, we successfully identified and characterized three effectors (BxSapB1, BxSapB3 and Bx-FAR-1) and a novel molecular pattern, BxCDP1, from the transcriptomes of B. xylophilus during the mycetophagous and earlier infection stages (6, 12 and 24 h postinoculation) [15,16,17,18]. However, our transcriptome data were only used to screen candidate effector lists, and a large amount of data has not been thoroughly analyzed, including the types of upregulated genes, the functional annotation, the molecular cellular pathways of the involved candidate effectors and the internal relationships between candidate effectors. Moreover, previous studies have shown that the effector pathogens secreted are different in the different infection stages [10, 19].
In this study, our objective was to determine the roles of candidate effectors in B. xylophilus infection and survival by comparing the types of candidate effectors secreted by B. xylophilus during interactions with host trees at different infection stages. For this purpose, we performed a relative comparative transcriptomic analysis of B. xylophilus inoculated onto pines between our data and the above two previously reported B. xylophilus transcriptomic datasets from different stages of infection (2.5 h, 6 and 15 d postinfection). Particular emphasis was placed on the identification of DEGs, comparison of DEGs, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of candidate effectors of B. xylophilus from the three transcriptomic groups. The relative expression of five common candidate effectors from three B. xylophilus transcriptomes was detected at infection stages by quantitative real-time polymerase chain reaction (qRT-PCR). Moreover, the transient expression of the five common candidate effectors was characterized via a potato virus X (PVX) expression vector in Nicotiana benthamiana.
De novo assembly of the transcriptome and comparison of summary data from three transcriptomes for Bursaphelenchus xylophilus
In total, we generated 592,991,982 raw reads and 570,138,674 clean reads from 12 samples by Illumina sequencing, respectively, which have been shown in our previous study . The Pearson correlation between these samples was calculated (Fig. S1). In this study, we denoted our transcriptome sequencing data and two previously reported B. xylophilus transcriptomic datasets from different stages of infection at 2.5 h, 6 and 15 d postinfection as the B group, A group and C group, respectively. The comparison of summary data from the three transcriptomes (A, B and C group) for B. xylophilus is shown in Table S1. The comparative results showed that the numbers of raw reads and clean reads of samples from the A group and C group were much lower than those from our experimental data (B group).
Screening of upregulated genes during the infection stage of B. xylophilus
In our study, the genes with a P-value < 0.05 and log2 (fold change) > 1 were assigned as DEGs between the early stages of infection (6, 12 and 24 h postinfection) and the mycetophagous stage (0 h). Finally, a total of 867 DEGs were obtained. Of these DEGs, 247 were upregulated during at least one of the phytophagous time points compared with the mycetophagous stage. A Venn diagram and heatmap of DEGs at the mycetophagous stage and three early phytophagous parasitic stages are shown in Figs. S2-S3. Of the 2272 DEGs in the A group, 1143 genes were upregulated in at least one postinoculation sample. Compared to the control, 60 genes were upregulated in both postinoculation samples . In the previous literature, there was no specific indication of the number of different genes in the C group, but the top 200 sequences upregulated in the parasitic life stage of the nematode were identified .
Functional annotation of upregulated genes from B. xylophilus transcriptomes
Transcriptome annotation provides insight into the structural, functional and biological processes in which genes are involved . To determine the functional annotation of upregulated B. xylophilus genes, we used the GO database. The GO terms of the upregulated genes from the three B. xylophilus transcriptomes were compared, and the transcriptomes ranked in the top 16 terms in the molecular function category are shown in Fig. 1. The results showed that the major functions of the upregulated genes involved were biological processes of B. xylophilus, such as embryo development, positive growth rates and reproduction. This result indicated that the genes related to the growth and development of B. xylophilus played a role in the initial stage of the plant-host interaction (2.5 h postinfection) (Fig. 1a). The most highly represented GO terms in this set of 200 upregulated genes of B. xylophilus in the molecular function category were hydrolase, oxidoreductase and lyase activity in the middle and later stages of infection (6 and 15 d postinfection) (Fig. 1b). In our experimental data, the most highly represented GO terms of the 247 upregulated genes in the molecular function category were molecular function, catalytic activity, hydrolase activity and binding at the early stage of infection regardless of the time points (6, 12 and 24 h postinfection) (Fig. 1c-e). At the same time, the same GO terms were highly represented in the latter two transcriptomes in the molecular function category, such as hydrolase, oxidoreductase activity, lyase activity and ion binding. These findings indicated that B. xylophilus uses some similar responses to survive in pines.
Identification of candidate effectors in three B. xylophilus transcriptomes
In our previous experimental data, 69 out of the 247 genes were denoted as candidate effectors . According to the accepted screening criterion of candidate effectors, we screened 62 candidate effectors in the B. xylophilus transcriptome at 2.5 h postinfection. Moreover, 46 candidate effectors were identified by Espada et al. . The list of candidate effectors from the three transcriptomes is shown in Table S3. To compare the candidate effectors from different B. xylophilus transcriptomes, these candidate effectors were matched with each other by local BLAST search. The results showed that seven candidate effectors were found in both our study and B. xylophilus transcriptome at 2.5 h postinfection, and one candidate effector was found in all three B. xylophilus transcriptomes (Table 1). These results showed that some effectors of B. xylophilus played a role in the early infection stage and that some effectors played a role in the middle and later infection stages. That is, similar to other pathogens, B. xylophilus secretes different effectors at different infection stages. Nevertheless, there are also some effectors that function throughout the whole infection process.
A previous study has characterized a B. xylophilus secretome during P. thunbergii infection . To determine the number of our predicted candidate effectors also identified as secreted proteins in this study, all the candidate effectors from three transcriptomes were matched to the secretome by local BLAST search. The results showed that 59, 29 and 24 candidate effectors from the B. xylophilus transcriptome at 2.5 h postinfection, 6, 12 and 24 h postinfection and 6 d and 15 d postinfection were identified in the secretome of B. xylophilus, respectively (Table S4). The results indicated some cross-validation between these studies.
Functional annotation of candidate effectors
To determine the functional annotation of B. xylophilus candidate effectors screened from different infection stages, we also used the GO database. Because the GO annotation of the candidate effectors for 6 and 15 d of B. xylophilus infection in pine has not been published, we matched 46 candidate effectors to our transcriptome data by local BLAST search and found the corresponding gene ID numbers. Then, we searched for the GO annotation of these genes in our transcriptome GO annotation information. The GO terms ranked in the top 15 in the molecular function are shown in Fig. 2. The results showed that the most highly represented GO terms of candidate effectors from three transcriptomes in the molecular function category were all molecular function, catalytic activity and hydrolase activity and that 8 out of 15 GO terms were the same between the transcriptomes of pines infected with PWN for 2.5 h and 6, 12 and 24 h (Fig. 2a, c). Ten out of 15 GO terms were the same between the transcriptomes of pines infected with PWN for 2.5 h and 6 d and 15 d (Fig. 2a-b), and 7 out of 15 GO terms were the same between the transcriptomes of pines infected with PWN for 6, 12 and 24 h and 6 and 15 d (Fig. 2b-c). This result indicated that candidate effectors from the three transcriptome groups played similar roles in the infection stage, and their hydrolase activity and catalytic activity played an important role in degrading the plant cell wall and removing ROS by catalysis, which helped B. xylophilus to successfully infect pine.
Categories of candidate effectors from three B. xylophilus transcriptomes
The candidate effectors were searched in the SWISS-Prot database by using BLASTP. The categories of candidate effectors from the three transcriptomes are shown in Fig. 3. Forty-four out of 62 candidate effectors from the B. xylophilus transcriptome at 2.5 h postinfection were uncharacterized proteins; thus, they were denoted as pioneers. Moreover, the other candidate effectors included several classes of proteases (pepsin A and cysteine proteinase), transthyretin-like protein, ribosomal protein and lipase as well as several enzymes involved in the detoxification of xenobiotic compounds, such as epoxide hydrolase and glutathione peroxidase (Fig. 3a). Twenty-one out of 69 candidate effectors from the B. xylophilus transcriptome at 6, 12 and 24 h postinfection were pioneers. The other candidate effectors included proteases (such as cysteine proteinase B, cathepsin L and zinc metalloproteinase nas-14), carbohydrate active enzymes (CAZymes; such as pectate lyase A, pectate lyase H, pectate lyase E, and putative endoglucanase type K), lipase, probable GH family 25 lysozyme and detoxification-related proteins, including lysosomal acid phosphatase, iron/zinc purple acid phosphatase-like protein and epoxide hydrolase (Fig. 3b). Sixteen out of 46 candidate effectors from the B. xylophilus transcriptome at 6 and 15 d postinfection were pioneers. The other candidate effectors also included several classes of proteases (aspartic protease A1, cysteine proteases C1A and serine-type protease), fatty acid transport proteins, putative V5/TPx1 allergen, a lysozyme and several detoxification-related proteins, such as UDP-glucuronosyl transferase, multicopper putative acid oxidase, glutathione S-transferase, cytochrome P450, acid phosphatase and epoxide hydrolase (Fig. 3c).
Comparison of the KEGG pathway distribution for candidate effectors from the three B. xylophilus transcriptomes
The candidate effectors from three B. xylophilus transcriptomes were analyzed by BLAST using the KEGG pathway database with KOBAS software. The results showed that most candidate effectors were not assigned to known KEGG pathways (Table 2). Only one candidate effector (BXY_1074200) from the B. xylophilus transcriptomes at 2.5 h postinfection was assigned to the metabolism of xenobiotics by the cytochrome P450 pathway. At the same time, only one candidate effector (BXY_0298700) from the B. xylophilus transcriptomes at 6 and 15 d postinfection was assigned to the drug metabolism-cytochrome P450 pathway. Seven candidate effectors from the B. xylophilus transcriptomes at 6, 12 and 24 h postinfection were assigned to four pathways as follows: drug metabolism-other enzymes, metabolic pathways, lysosomes and metabolism of xenobiotics by cytochrome P450 pathways. This result indicated that the cytochrome P450 pathway is important to B. xylophilus during the entire infection stage.
Eight common candidate effectors are highly upregulated at early infection stages
According to the above transcriptome data for B. xylophilus, 8 common candidate effectors (including three known effectors, namely, BxSapB1, BxSapB3 and Bx-FAR-1) were upregulated at the infection stages. In addition to BxSapB1, BxSapB3 and Bx-FAR-1, the other five common candidate effectors were named Bx-C1, Bx-C2, Bx-C3, Bx-C4 and BxSCD5. To further confirm this finding, qRT-PCR was employed to obtain their expression profiles at the early stages of infection. The results showed that the five common candidate effectors were upregulated at early infection stages compared to the mycophagous stage. The relative expression levels of Bx-C1, Bx-C3, Bx-C4 and BxSCD5 were highest at 2.5 h compared to other infection times, and the relative expression of Bx-C2 was highest at 24 h (Fig. 4). These results suggested that these common candidate effects play important roles in the stages of B. xylophilus infection.
BxSCD5 effectively inhibits PAMP PsXEG1-triggered cell death
A previous study has shown that some Phytophthora effectors cause cell death in host and/or nonhost plants, such as the P. sojae RXLR effector, Avh241 . To test whether the 5 common candidate effectors of B. xylophilus induce cell death, they were expressed in N. benthamiana leaves using agroinfiltration. The results showed that the proportion of infiltrated sites that developed the cell death phenotype after injection of the 5 candidate effectors was almost zero, indicating that they do not trigger cell death in N. benthamiana (Fig. 5a).
PAMPs trigger strong defense responses in various plants, and many P. sojae effectors suppress immune responses, including PAMP-triggered cell death . To determine whether the 5 common candidate effectors (except BxSapB1, BxSapB3 and Bx-FAR-1) of B. xylophilus suppress P. sojae PAMP PsXEG1-triggered cell death, we expressed PsXEG1 in N. benthamiana leaves 16 h after expressing the 5 common candidate effectors or GFP (negative control) using agroinfiltration. The proportion of infiltrated sites that developed the cell death phenotype after injection of BxSCD5NSP (without SP) and the other 4 candidate effectors followed by injection of PsXEG1 was almost 100%, except for BxSCD5 (with SP). These results showed that only BxSCD5 suppressed PsXEG1-triggered cell death (Fig. 5b), indicating that BxSCD5 suppresses PsXEG1-triggered cell death in N. benthamiana when secreted into the apoplast.
It is generally accepted that the DEGs of pathogens, especially upregulated genes, are involved in the interaction between pathogens and hosts during the infection stage. In this study, the raw reads, clean reads and upregulated genes obtained by the three transcriptome groups were different, and the annotation information for the upregulated genes showed that the gene functions in B. xylophilus infection were different in different periods. In the A group (2.5 h postinfection), the upregulated genes were mainly involved in the growth and development of B. xylophilus, which may be due to adjustments made by B. xylophilus to survive and better adapt to the new environment during the sudden change from the mycetophagous stage to the phytophagous stage. In the B group (6, 12 and 24 h postinfection), the major upregulated genes were involved in molecular function, catalytic activity, hydrolase activity and binding, and the number of upregulated genes involved in these GO terms was vast. In the C group (6 and 15 d postinfection), the upregulated genes were mostly involved in hydrolase, oxidoreductase and lyase activities. Interestingly, the number of upregulated genes involved in hydrolase, oxidoreductase and lyase activities in the B group was greater than that in the C group. This result indicated that it was more important for B. xylophilus to degrade the cell wall of the host by secreting hydrolase and lyase to promote parasitism at the early infection stage. Previous studies have demonstrated that some plant cell wall-degrading enzymes (CWDEs) of B. xylophilus have roles as pathogenicity determinants, such as glycoside hydrolase family 45 cellulases, pectate lyases and b-1,3-endoglucanases [23,24,25]. In this study, several pectate lyases and endoglucanases were identified as candidate effectors in the B group, supporting the viewpoint that these CWDEs are generally regarded as “effectors”. Moreover, an oxidative burst is one of the earliest defense responses to plant pathogen attack. The transient accumulation of ROS helps to defend the host from invading pathogens and also acts as a signaling molecule to trigger various other plant defense responses . Some previous studies have shown that the capability of ROS scavenging in the host is pivotal to the high resistance of PWN and that enhancement of oxidative stress contributes to increased pathogenicity of B. xylophilus [27, 28]. The amount of ROS in the tree is critical to the survival of B. xylophilus. In addition, catalases (Bxy-ctl-1 and Bxy-ctl-2) of high virulence B. xylophilus are crucial for nematode survival under prolonged exposure to in vitro oxidative stress . Surprisingly, Bxy-ctl-1 (BUX.s00579.159) was identified in the B group (i.e., BXY_1386500) with oxidoreductase activity, participating in the oxidation-reduction process. Moreover, the significant GO terms of oxidoreductase activity and oxidation-reduction process of virulent B. xylophilus strain were enriched statistically compared to the avirulent strain . In our study, the significant GO term, oxidoreductase activity (GO:0016491), was found in both the B and C groups, and the DEG list of enriched statistically terms was mostly in the B group. Thus, we speculated that the upregulated genes functioning in catalytic and oxidoreductase activities in the B and C groups may play a key role in removing ROS by catalysis and oxidoreduction to help B. xylophilus survive in the host.
The genomic data of B. xylophilus indicate that this nematode has many specific genes . A previous study has shown that effector repertoires of sedentary PPNs (e.g., cyst and root-knot nematodes) contain hundreds of proteins implicated in the establishment of a successful plant-pathogen interaction . In this study, although the number of upregulated genes in the three transcriptomes was different, the number of candidate effectors was similar (62, 69 and 46), indicating that the number of effectors of B. xylophilus may be less than that of sedentary PPNs. Although the life cycle of B. xylophilus is different from that of other PPNs, homologous genes in different PPNs often play similar functions. Several candidate effectors screened in the present study have also been studied in other PPNs. For example, glutathione S transferase (GST) has been demonstrated to play an important role in root-knot nematodes , and GST was found in the C group in the present study. Fatty acid retinoid-binding (FAR) proteins have also been studied in several PPNs (such as Radopholus similis, Globodera pallida and P. penetran), indicating that FAR proteins are necessary for the entire nematode life cycle [33,34,35]. In the present study, B. xylophilus FAR proteins were identified in the three groups. Moreover, Bx-FAR-1 has been characterized as a key effector by Li . The VAPs of the G. rostochiensis potato cyst nematode are required for the suppression of host immunity . In the present study, one candidate effector similar to VAPs was identified in the C group.
CAZymes (such as cellulase, pectate lyase, calreticulin and expansin) are important for breaking down the polysaccharides of plant cell walls to establish infection . In the present study, some CAZymes (putative pectate lyase A, H and E as well as endoglucanase type K) were included in the candidate effectors listed in the B group, which may play roles in invasion, extension and degradation of host tissues. Fatty acid-binding proteins, retinol-binding proteins and peroxiredoxin have been reported to protect nematodes from antipathogen compounds in host defenses . In the present study, fatty acid- and retinol-binding proteins were identified in all groups. Some effectors, such as annexin, venom allergen-like protein, transthyretin-like protein and ubiquitin extension proteins, target plant signaling pathways and suppress host defenses . Here, the putative allergen, V5/TPx1, and transthyretin-like protein were identified in the C and A groups, respectively. A previous study has shown that some effectors are required for feeding efficiency, such as cathepsin. Proteolytic enzymes can be divided into four main groups as follows: cysteine, serine, aspartyl and metalloproteinases. Among them, cysteine proteinases are the most extensively studied . In our study, cysteine proteinases, including cathepsin L, were found in all three transcriptome groups. These results indicated that proteolytic enzymes, especially cysteine proteinases, of B. xylophilus play a key role in nutritional intake to promote parasitism. Some studies have also shown that cysteine proteinases play important roles in embryogenesis, development, infection, parasitism, pathogenesis and immune evasion in nematodes and many other animal parasites [39,40,41]. Moreover, a cathepsin L-like cysteine proteinase of B. xylophilus has been demonstrated to influence its development and pathogenicity . Several lipases of pathogenic bacteria and fungi have been identified as virulence factors . Here, several other digestive enzymes, such as serine-type protease, pepsin A and aspartic protease A1, were also identified in the three groups, and lipase was present in both the A and B groups. Many putative effectors have unknown functions and are called pioneers. In the present study, many pioneers were identified in three transcriptome groups, and most of them successfully matched the secretome of B. xylophilus . The functional study of these pioneers will be an important step to investigate the pathogenesis of B. xylophilus.
In addition to the above functions, we found that some effectors were related to the detoxification of xenobiotic compounds. For example, epoxide hydrolase and glutathione peroxidase were identified in the A group. Epoxide hydrolases have been characterized as virulence factors in previous studies [44, 45]. Unlike epoxide hydrolase, lysosomal acid phosphatase and iron/zinc purple acid phosphatase-like protein were included in the B group. UDP-glucuronosyl transferase, multicopper putative acid oxidase, GST, cytochrome P450 and acid phosphatase were categorized in the C group. This result suggested that both early and late infection stages of B. xylophilus require detoxification to overcome host resistance, but some detoxification substances secreted by B. xylophilus are different. Moreover, several candidate effectors in the A group, such as ribosomal proteins, were suggested to participate in RNA transcription and translation. Additionally, 17 genes encoding ribosomal proteins (RPs) have been found to be significantly altered when Meloidogyne incognita infects Solanum lycopersicum roots . These results indicated that ribosomal proteins play an important role in the interaction of phytoparasitic nematodes and their hosts.
In the process of interaction between pathogens and hosts, the secretion of PAMPs and effectors by pathogens can trigger plant defense responses, including PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). However, under the strong pressure of natural selection, pathogens will continue to produce new effectors to inhibit the PTI and ETI of the host and help pathogens escape host recognition . At the same time, synergistic cooperation between effectors can promote the successful infection of pathogens . A previous study has shown that the expression of some PR genes of pine is induced by infection with virulent isolates of B. xylophilus , indicating that B. xylophilus might secrete some PAMPs or molecular patterns, such as BxCDP1 , to induce the defense response of pines. In our study, we transiently expressed the five common candidate effectors (except the three known effectors, BxSapB1, BxSapB3 and Bx-FAR-1) in the N. benthamiana model plant for functional analysis. Among the five candidate effectors, the BxSCD5 candidate effector effectively suppressed PAMP PsXEG1-triggered cell death in N. benthamiana. Thus, we speculated that although the defense response of the host is induced after infection with B. xylophilus, there are also some effectors that inhibit host defense to some extent to help nematodes successfully infect, such as BxSCD5. Effectors of many pathogens have been identified using this strategy [49,50,51,52]. Thus, it is also feasible to identify effectors of B. xylophilus in the N. benthamiana nonhost plant. Certainly, the role effectors play in the pathogen-host interaction process still needs to be studied in their own host. Thus, the function of BxSCD5 in host pines will be studied in the future.
In the present study, common effectors (7 candidate effectors were in both the A and B groups, and 1 candidate effector was in all three groups) were identified in the three B. xylophilus transcriptomes. There are several possibilities why only a few factors were identified. First, similar to other pathogens, B. xylophilus secretes different effectors at different infection stages to cope with various survival threats, resulting in a few common effectors in three different infection stages. Second, only the top 200 upregulated genes in the parasitic life stage of the nematode were selected to predict candidate effectors in the C group (6 and 15 d postinfection) ; thus, the candidate effectors were screened from the top 200 upregulated genes, not from all of the upregulated genes, which might result in the loss of some effectors. Third, different B. xylophilus strains, pine species, tree age, inoculation environment and sequencing depth of the three transcriptomes used might lead to different DEGs and candidate effectors. Nevertheless, the BxSapB1, BxSapB3, Bx-FAR-1 and BxSCD5 key effectors were still screened from the three different B. xylophilus transcriptomes and demonstrated to play an important role in promoting nematode infection. Thus, despite its shortcomings, the relative comparative transcriptomic analysis used in the present study is feasible to some degree. To fully understand the mechanism of B. xylophilus infection, we will also study the function of noncommon candidate effectors from the three transcriptomes in the future.
The aim of the present study was to explore infection and survival strategies by comparative transcriptomic analysis of upregulated genes and candidate effectors of B. xylophilus in different infection stages. B. xylophilus secretes various candidate effectors in different infection stages with some of them continuing to function throughout all infection stages. Candidate effectors help B. xylophilus infect and survive successfully in hosts in different ways (such as breaking down host cell walls, suppressing host defenses, promoting feeding efficiency, detoxification and playing virulence functions). We provided continued evidence for the presence of ‘common’ PPN effectors and identified some novel effectors involved in the infection process of B. xylophilus. These results provide many resources for studying the pathogenesis of B. xylophilus from the perspective of effectors. Moreover, the identification and characterization of the plant targets of these effectors should also be undertaken to understand how these effectors hijack many aspects of plant cell morphology and physiology, including the immune system.
In this study, the highly virulent B. xylophilus strain AMA3 was used, which was from Anhui Province, China . The culture methods of B. xylophilus were similar to those previously reported . The culture and isolation of nematodes were performed according to our previous study .
PWN inoculation trials
Pinus thunbergii seedlings (3 years old) obtained from the experimental field of Nanjing Forestry University (Jurong Yaolingkou Forest Farm, Jiangsu, China) were used for inoculation of the PWN strain AMA3. The number of nematodes and time of inoculation have been described in our previous study . Briefly, a suspension of approximately 10,000 mixed life stage nematodes was collected from PDA as a mycetophagous control (B-0 h). The same number of nematodes (AMA3) was inoculated into the pine stems for 6 h (B-6 h), 12 h (B-12 h) and 24 h (B-24 h). The nematodes were collected and then frozen for further RNA isolation according to our previous method . Each treatment had three biological replicates.
Sample preparation and Illumina sequencing
Nematode (AMA3) RNA was extracted from the B-0 h, B-6 h, B-12 h and B-24 h samples. Then, mRNA degradation and contamination were monitored on 1% agarose gels. Then, the RNA purity, RNA concentration, RNA integrity, cDNA library preparations and sequencing were separately determined according to our previous study . The RNA-seq data for this study are submitted on the NCBI under accession number PRJNA397001.
Differential gene expression analysis
HTSeq v0.6.1 and DESeq were used to count the read numbers mapped to each gene and standardize the read count of each gene, respectively. Differential expression analysis of four treatments was performed using the DESeq R package (1.18.0), which has been previously demonstrated to be better than FPKM or RPKM . Finally, genes (with a P-value < 0.05 and log2 (fold change) > 1) were assigned as differentially expressed among four treatments. The definition of DEGs were described in our previous study .
Bursaphelenchus xylophilus transcriptomic datasets
The two previously reported B. xylophilus transcriptomic datasets from different stages of infection (2.5 h postinfection and 6 d and 15 d postinfection) [13, 14] were collected for subsequent relative comparative transcriptomic analysis with our transcriptomic experimental data. Among them, one experimental design used two postinoculation events (2.5 h postinfection in August and September) and one control (mycetophagous stage) with two biological replicates per condition, which were denoted as A-2.5 h-8, A-2.5 h-9 and A-0 h. The other dataset also used two postinoculation events (6 and 15 d postinfection) and one control (mycetophagous stage), which were denoted as C-6 d, C-15 d and C-0 d. The two postinoculation events of the C group had three biological replicates, and the control had two biological replicates.
GO enrichment analysis of differentially expressed genes
The GO enrichment analysis of B. xylophilus DEGs from the above three transcriptomic datasets was implemented by the GOseq R package, in which gene length bias was corrected . The GO terms were divided into biological processes, cellular components and molecular functions for functional categorization. GO terms with corrected P values (i.e., false discovery rate, FDR) less than 0.05 were considered significantly enriched.
Candidate effectors analysis
We screened potential effector proteins according to our previous study (the presence of an N-terminal signal peptide and the absence of a transmembrane domain) . The signal peptide and the transmembrane domain were predicted by SignalP 4.1 (http://www.cbs.dtu.dk/services/SignalP/)  and TMHMM 2.0 (http://www.cbs.dtu.dk/services/TMHMM/) , respectively. All candidate effectors were predicted functions according to our previous method .
KEGG enrichment analysis of candidate effectors of B. xylophilus
KEGG is a database resource for understanding the high-level functions and utilities of biological systems (http://www.genome.jp/kegg/) . KOBAS software was used to determine the statistical enrichment of B. xylophilus candidate effectors from the above three transcriptomic datasets in KEGG pathways. Furthermore, to improve the annotation, a BLAST search was performed against the SWISS-Prot database.
Local BLAST analysis
The amino acid sequences of candidate effectors from the above three transcriptomic datasets were matched to the secretome of B. xylophilus by a local BLAST search . The BLAST search was performed with an e-value = 1e-10, and the proteins with percentage similarity > 80% were selected.
Referring to a previously used method , approximately 10,000 AMA3 nematodes were inoculated into 3-year-old P. thunbergii seedlings, and the nematodes were then collected at 2.5, 6, 12 and 24 h after inoculation. RNA of AMA3 nematodes was extracted and reverse transcribed into cDNA. qRT-PCR assays were performed according to our previous method . The expression levels of the five common candidate effectors (except the three known effectors, BxSapB1, BxSapB3 and Bx-FAR-1) from the above three B. xylophilus transcriptomes at the early stage of host infection (2.5, 6, 12 and 24 h) were measured. Actin of B. xylophilus (GenBank EU100952) was used as a constitutively expressed endogenous control gene . All assays were performed three times. Primer sequences are provided in Table S2.
Referring to the previously used method , five common candidate effectors from the above three B. xylophilus transcriptomes were cloned with and without its native signal peptide (SP) from B. xylophilus cDNA using the specific primers listed in Table S2. Subsequently, purified PCR products were ligated into the pGR107 vector (pGR107–3*HA) using the CloneExpress II One Step Cloning Kit (Vazyme, Nanjing, China) after confirmation by sequencing.
Transient expression assays in N. benthamiana
The transformation method of the constructed plasmids, culture conditions of the Agrobacterium tumefaciens GV3101 strain and agroinfiltration assays were the same as those in our previous study . One experimental design directly infiltrated A. tumefaciens suspensions harboring five candidate effectors (with and without its native SP) into the leaves of N. benthamiana using a needleless syringe. Another experimental design infiltrated A. tumefaciens suspensions harboring five candidate effectors (with and without its native SP) 16 h prior to those infiltrated with Phytophthora sojae pathogen-associated molecular pattern (PAMP) PsXEG1. Green fluorescent protein (GFP) was used as the negative control. The phenotypes of N. benthamiana leaves were recorded at 5–7 days after infiltration. Each experiment was repeated three times, and each experiment consisted of three plants with three inoculated leaves. Finally, we counted the proportion of infiltrated sites that developed the cell death phenotype to screen the effectors that triggered cell death or inhibited PsXEG1-triggered cell death.
Availability of data and materials
The RNA-seq data in this study was available through the NCBI under accession number PRJNA397001. We have a permission to collect the plant samples we used in this study. All the data and materials that are required to reproduce these findings can be shared by contacting the corresponding author, Prof. Xiao-Qin Wu (firstname.lastname@example.org).
Pine wood nematode
Pine wilt disease
- B. xylophilus :
Expressed sequence tags
- G. asiaticumrice :
- P. penetran :
- G. rostochiensis :
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
False discovery rate
Potato virus X
Quantitative real-time polymerase chain reaction
- N. benthamiana :
Carbohydrate active enzymes
Cell wall-degrading enzymes
- P. sojae :
Potato dextrose agar
- A. tumefaciens :
Reactive oxygen species
Vicente C, Espada M, Vieira P, Mota M. Pine wilt disease: a threat to European forestry. Eur J Plant Pathol. 2011;133(1):89–99.
Jones JT, Haegeman A, Danchin EG, Gaur HS, Helder J, Jones MG, et al. Top 10 plant-parasitic nematodes in molecular plant pathology. Mol Plant Pathol. 2013;14(9):946–61. https://doi.org/10.1111/mpp.12057.
Ma Z, Song T, Zhu L, Ye W, Wang Y, Shao Y, et al. A Phytophthora sojae glycoside hydrolase 12 protein is a major virulence factor during soybean infection and is recognized as a PAMP. Plant Cell. 2015;27(7):2057–72. https://doi.org/10.1105/tpc.15.00390.
Abramovitch RB, Janjusevic R, Stebbins CE. Type III effector AvrPtoB requires intrinsic E3 ubiquitin ligase activity to suppress plant cell death and immunity. Proc Natl Acad Sci U S A. 2006;103(8):2851–6.
Dodds PN, Lawrence GJ, Catanzariti A. Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes. Proc Natl Acad Sci U S A. 2006;103(23):8888–93. https://doi.org/10.1073/pnas.0602577103.
Chen J, Lin B, Huang Q, Hu L, Zhuo K, Liao J. A novel Meloidogyne graminicola effector, MgGPP, is secreted into host cells and undergoes glycosylation in concert with proteolysis to suppress plant defenses and promote parasitism. PLoS Pathog. 2017;13(4):e1006301. https://doi.org/10.1371/journal.ppat.1006301.
Mitchum MG, Hussey RS, Baum TJ, Wang X, Elling AA, Wubben M, et al. Nematode effector proteins: an emerging paradigm of parasitism. New Phytol. 2013;199(4):879–94. https://doi.org/10.1111/nph.12323.
Sonah H, Zhang X, Deshmukh RK, Borhan MH, Fernando WG, Belanger RR. Comparative transcriptomic analysis of virulence factors in Leptosphaeria maculans during compatible and incompatible interactions with canola. Front Plant Sci. 2016;7:1784.
Bauters L, Haegeman A, Kyndt T, Gheysen G. Analysis of the transcriptome of Hirschmanniella oryzae to explore potential survival strategies and host-nematode interactions. Mol Plant Pathol. 2014;15(4):352–63. https://doi.org/10.1111/mpp.12098.
Vieira P, Maier TR, Eves-van den Akker S, Howe DK, Zasada I, Baum TJ, et al. Identification of candidate effector genes of Pratylenchus penetrans. Mol Plant Pathol. 2018;19(8):1887–907. https://doi.org/10.1111/mpp.12666.
Tao SQ, Cao B, Tian CM, Liang YM. Comparative transcriptome analysis and identification of candidate effectors in two related rust species (Gymnosporangium yamadae and Gymnosporangium asiaticum). BMC Genomics. 2017;18(1):651. https://doi.org/10.1186/s12864-017-4059-x.
Yang D, Chen C, Liu Q, Jian H. Comparative analysis of pre- and post-parasitic transcriptomes and mining pioneer effectors of Heterodera avenae. Cell Biosci. 2017;7(1):11. https://doi.org/10.1186/s13578-017-0138-6.
Tsai IJ, Tanaka R, Kanzaki N, Akiba M, Yokoi T, Espada M, et al. Transcriptional and morphological changes in the transition from mycetophagous to phytophagous phase in the plant-parasitic nematode Bursaphelenchus xylophilus. Mol Plant Pathol. 2016;17(1):77–83. https://doi.org/10.1111/mpp.12261.
Espada M, Silva AC. Eves van den Akker S, Cock PJ, Mota M, Jones JT: identification and characterization of parasitism genes from the pinewood nematode Bursaphelenchus xylophilus reveals a multilayered detoxification strategy. Mol Plant Pathol. 2016;17(2):286–95. https://doi.org/10.1111/mpp.12280.
Hu LJ, Wu XQ, Li HY, Zhao Q, Wang YC, Ye JR. An effector, BxSapB1, induces cell death and contributes to virulence in the pine wood nematode Bursaphelenchus xylophilus. Mol Plant-Microbe Interact. 2019;32(4):452–63. https://doi.org/10.1094/MPMI-10-18-0275-R.
Huang X, Hu LJ, Wu XQ. Identification of a novel effector BxSapB3 that enhances the virulence of pine wood nematode Bursaphelenchus xylophilus. Acta Biochim Biophys Sin. 2019;51(10):1071–8. https://doi.org/10.1093/abbs/gmz100.
Hu LJ, Wu XQ, Li HY, Wang YC, Huang X, Wang Y, et al. BxCDP1 from the pine wood nematode Bursaphelenchus xylophilus is recognized as a novel molecular pattern. Mol Plant Pathol. 2020;21(7):923–35. https://doi.org/10.1111/mpp.12939.
Li Y, Hu LJ, Wu XQ, Ye JR. A Bursaphelenchus xylophilus effector, Bx-FAR-1, suppresses plant defense and affects nematode infection of pine trees. Eur J Plant Pathol. 2020;153(3):637–50.
Wang QQ, Han CZ, Ferreira AO, Yu XL, Ye WW, Tripathy S, et al. Transcriptional programming and functional interactions within the Phytophthora sojae RXLR effector repertoire. Plant Cell. 2011;23(6):2064–86. https://doi.org/10.1105/tpc.111.086082.
Mutz KO, Heilkenbrinker A, Lönne M, Walter JG, Stahl F. Transcriptome analysis using next-generation sequencing. Curr Opin Biotechnol. 2013;24(1):22–30. https://doi.org/10.1016/j.copbio.2012.09.004.
Shinya R, Morisaka H, Kikuchi T, Takeuchi Y, Ueda M, Futai K. Secretome analysis of the pine wood nematode Bursaphelenchus xylophilus reveals the tangled roots of parasitism and its potential for molecular mimicry. PLoS One. 2013;8:e67377.
Yu X, Tang J, Wang Q, Ye W, Tao K, Duan S, et al. The RxLR effector Avh241 from Phytophthora sojae requires plasma membrane localization to induce plant cell death. New Phytol. 2012;196(1):247–60. https://doi.org/10.1111/j.1469-8137.2012.04241.x.
Kikuchi T, Jones JT, Aikawa T, Kosaka H, Ogura N. A family of glycosyl hydrolase family 45 cellulases from the pine wood nematode Bursaphelenchus xylophilus. FEBS Lett. 2004;572(1-3):201–5. https://doi.org/10.1016/j.febslet.2004.07.039.
Kikuchi T, Shibuya H, Jones JT. Molecular and biochemical characterization of an endo-β-1,3-glucanase from the pine wood nematode Bursaphelenchus xylophilus acquired by horizontal gene transfer from bacteria. Biochem J. 2005;389(1):117–25. https://doi.org/10.1042/BJ20042042.
Kikuchi T, Shibuya H, Aikawa T, Jones JT. Cloning and characterization of pectate lyases expressed in the esophageal gland of the pine wood nematode Bursaphelenchus xylophilus. Mol Plant-Microbe Interact. 2006;19(3):280–7. https://doi.org/10.1094/MPMI-19-0280.
Goverse A, Smant G. The activation and suppression of plant innate immunity by parasitic nematodes. Annu Rev Phytopathol. 2014;52(1):243–65. https://doi.org/10.1146/annurev-phyto-102313-050118.
Liu Q, Wei Y, Xu L, Hao Y, Chen X, Zhou Z. Transcriptomic profiling reveals differentially expressed genes associated with pine wood nematode resistance in masson pine (Pinus massoniana lamb.). Sci Rep. 2017;7(1):4693–707. https://doi.org/10.1038/s41598-017-04944-7.
Zhang W, Zhao L, Zhou J, Yu H, Zhang C, Lv Y, et al. Enhancement of oxidative stress contributes to increased pathogenicity of the invasive pine wood nematode. Philos Trans R Soc Lond Ser B Biol Sci. 2019;374(1767):1–10.
Vicente CS, Ikuyo Y, Shinya R, Mota M, Hasegawa K. Catalases induction in high virulence pine wood nematode Bursaphelenchus xylophilus under hydrogen peroxide-induced stress. PLoS One. 2015;10(4):e0123839. https://doi.org/10.1371/journal.pone.0123839.
Filipiak A, Malewski T, Matczynska E, Tomalak M. Molecular variation among virulent and avirulent strains of the quarantine nematode Bursaphelenchus xylophilus. Mol Genet Genomics. 2021;296(2):259-69.
Kikuchi T, Cotton JA, Dalzell JJ, Hasegawa K, Kanzaki N, McVeigh P, et al. Genomic insights into the origin of parasitism in the emerging plant pathogen Bursaphelenchus xylophilus. PLoS Pathog. 2011;7(9):e1002219. https://doi.org/10.1371/journal.ppat.1002219.
Dubreuil M, Magliano E, Deleury P, Abad P, Rosso MN. Transcriptome analysis of root-knot nematode functions induced in the early stages of parasitism. New Phytol. 2007;176(2):426–36. https://doi.org/10.1111/j.1469-8137.2007.02181.x.
Prior A, Jones JT, Blok VC, Beauchamp J, McDermott L, Cooper A, et al. A surface-associated retinol- and fatty acid-binding protein (Gp-FAR-1) from the potato cyst nematode Globodera pallida: lipid binding activities, structural analysis and expression pattern. Biochem J. 2001;356(2):387–94. https://doi.org/10.1042/bj3560387.
Zhang C, Xie H, Cheng X, Wang WD, Li Y, Xu CL, et al. Molecular identification and functional characterization of the fatty acid- and retinoid-binding protein gene Rs-far-1 in the burrowing nematode Radopholus similis (Tylenchida: Pratylenchidae). PLoS One. 2015;10(3):e0118414. https://doi.org/10.1371/journal.pone.0118414.
Vieira P, Kamo K, Eisenback JD. Characterization and silencing of the fatty acid- and retinol-binding Pp-far-1 gene in Pratylenchus penetrans. Plant Pathol. 2017;66(7):1214–24. https://doi.org/10.1111/ppa.12664.
Lozano-Torres JL, Wilbers RHP, Warmerdam S, Finkers-Tomczak A, Diaz-Granados A, van Schaik CC, et al. Apoplastic venom allergen-like proteins of cyst nematodes modulate the activation of basal plant innate immunity by cell surface receptors. PLoS Pathog. 2014;10(12):e1004569. https://doi.org/10.1371/journal.ppat.1004569.
Haegeman A, Mantelin S, Jones JT, Gheysen G. Functional roles of effectors of plant-parasitic nematodes. Gene. 2012;492(1):19–31. https://doi.org/10.1016/j.gene.2011.10.040.
Wang K, Li Y, Huang X, Wang DW, Xu CL, Xie H. The cathepsin S cysteine proteinase of the burrowing nematode Radopholus similis is essential for the reproduction and invasion. Cell Biosci. 2016;6(1):39. https://doi.org/10.1186/s13578-016-0107-5.
Carmona C, Dowd AJ, Smith AM, Dalton JP. Cathepsin L proteinase secreted by Fasciola hepatica in vitro prevents antibody-mediated eosinophil attachment to newly excysted juveniles. Mol Biochem Parasitol. 1993;62(1):9–17. https://doi.org/10.1016/0166-6851(93)90172-T.
Hashmi S, Britton C, Liu J, Guiliano DB, Oksov Y, Lustigman S. Cathepsin L is essential for embryogenesis and development of Caenorhabditis elegans. J Biol Chem. 2002;277(5):3477–86. https://doi.org/10.1074/jbc.M106117200.
Malagón D, Díaz-López M, Benítez R, Adroher FJ. Cathepsin B- and L-like cysteine protease activities during the in vitro development of Hysterothylacium aduncum (Nematoda: Anisakidae), a worldwide fish parasite. Parasitol Int. 2010;59(1):89–92. https://doi.org/10.1016/j.parint.2009.11.001.
Xue Q, Wu XQ, Zhang WJ, Deng LN, Wu MM. Cathepsin L-like cysteine proteinase genes are associated with the development and pathogenicity of pine wood nematode, Bursaphelenchus xylophilus. Int J Mol Sci. 2019;20(1):215. https://doi.org/10.3390/ijms20010215.
Stehr F, Kretschmar M, Kröger C, Hube B, Schäfer W. Microbial lipases as virulence factors. J Mol Catal B Enzym. 2003;22(5–6):347–55. https://doi.org/10.1016/S1381-1177(03)00049-3.
Oesch F. Mammalian epoxide hydrases: inducible enzymes catalysing the inactivation of carcinogenic and cytotoxic metabolites derived from aromatic and olefinic compounds. Xenobiotica. 1973;3(5):305–40. https://doi.org/10.3109/00498257309151525.
Hvorecny KL, Bahl CD, Kitamura S, Lee KSS, Hammock BD, Morisseau C, et al. Active-site flexibility and substrate specificity in a bacterial virulence factor: crystallographic snapshots of an epoxide hydrolase. Structure. 2017;25(5):697–707 e694. https://doi.org/10.1016/j.str.2017.03.002.
Shukla N, Yadav R, Kaur P, Rasmussen S, Goel S, Agarwal M, et al. Transcriptome analysis of root-knot nematode (Meloidogyne incognita)-infected tomato (Solanum lycopersicum) roots reveals complex gene expression profiles and metabolic networks of both host and nematode during susceptible and resistance responses. Mol Plant Pathol. 2018;19(3):615–33. https://doi.org/10.1111/mpp.12547.
Jones JD, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9. https://doi.org/10.1038/nature05286.
Hirao T, Fukatsu E, Watanabe A. Characterization of resistance to pine wood nematode infection in Pinus thunbergii using suppression subtractive hybridization. BMC Plant Biology. 2012;12(13):207-12. https://doi.org/10.1186/1471-2229-12-13.
Chen C, Liu S, Liu Q, Niu J, Liu P, Zhao J, et al. An ANNEXIN-like protein from the cereal cyst nematode Heterodera avenae suppresses plant defense. PLoS One. 2015;10(4):e0122256. https://doi.org/10.1371/journal.pone.0122256.
Zhang L, Ni H, Du X, Wang S, Ma XW, Nurnberger T, et al. The Verticillium-specific protein VdSCP7 localizes to the plant nucleus and modulates immunity to fungal infections. New Phytol. 2017;215(1):368–81. https://doi.org/10.1111/nph.14537.
Chen C, Chen Y, Jian H, Yang D, Dai Y, Pan L, et al. Large-scale identification and characterization of Heterodera avenae putative effectors suppressing or inducing cell death in Nicotiana benthamiana. Front Plant Sci. 2017;8:2062.
Li Z, Yin Z, Fan Y, Xu M, Kang Z, Huang L. Candidate effector proteins of the necrotrophic apple canker pathogen Valsa mali can suppress BAX-induced PCD. Front Plant Sci. 2015;6:579.
Ding X, Ye J, Lin S, Wu X, Li D, Nian B. Deciphering the molecular variations of pine wood nematode Bursaphelenchus xylophilus with different virulence. PLoS One. 2016;11(5):e0156040. https://doi.org/10.1371/journal.pone.0156040.
Qiu X, Wu X, Huang L, Tian M, Ye J. Specifically expressed genes of the nematode Bursaphelenchus xylophilus involved with early interactions with pine trees. PLoS One. 2013;8(10):e78063. https://doi.org/10.1371/journal.pone.0078063.
Dillies MA, Rau A, Aubert J, Hennequet-Antier C, Jeanmougin M, Servant N, et al. A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. Brief Bioinform. 2013;14(6):671–83. https://doi.org/10.1093/bib/bbs046.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biology. 2010;11(2):R14. https://doi.org/10.1186/gb-2010-11-2-r14.
Thomas NP, Soren B, Gunnar von H, Henrik N. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Sonnhammer EL, von Heijne G, Krogh A. A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998;6:175–82.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44(D1):D457–62. https://doi.org/10.1093/nar/gkv1070.
We thank Prof. Yuanchao Wang (Departments of Plant Pathology, Nanjing Agricultural University, China) for providing the vector pGR107, the seeds of Nicotiana benthamiana and technical assistant. We are also grateful to Dr. Yan Wang (Departments of Plant Pathology, Nanjing Agricultural University, China) for helpful suggestions.
XQW and LJH were the leading investigators of this research programme. LJH and XQW planned and designed the research; LJH performed the majority of experiments with the help of XLD; XQW and JRY contributed reagents, materials, and analysis tools; LJH analyzed the data; LJH wrote the paper with suggestions from XQW, XLD and JRY. All authors commented on the article before submission. The author(s) read and approved the final manuscript.
This work was supported by the National Key Research and Development Program of China (2018YFD0600203), Natural Science Foundation of Jiangsu Province (BK20200774), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), and Jiangsu Post-doctoral Research Funding Program (2020Z366). The funding organizations had no role in designing, data collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison of summary data from three transcriptomes of Bursaphelenchus xylophilus at the phytophagous and mycetophagous stages. A: Transcriptome of B. xylophilus at 2.5 h postinfection. B: Transcriptome of B. xylophilus at 6, 12 and 24 h postinfection. C: Transcriptome of B. xylophilus at 6 and 15 d postinfection.
List of primers used in this study.
List of candidate effectors in three transcriptomes of pine wood nematodes Bursaphelenchus xylophilus.
List of candidate effectors in the Bursaphelenchus xylophilus secretome during Pinus thunbergii infection.
Pearson correlation between samples. bx-A: Mycophagous stage (0 h). bx-B: Bursaphelenchus xylophilus was inoculated into each of the pine stems for 6 h. bx-C: B. xylophilus was inoculated into each of the pine stems for 12 h. bx-D: B. xylophilus was inoculated into each of the pine stems for 24 h. Each treatment had three biological replicates.
Venn diagram of differentially expressed genes (DEGs) at the mycetophagous stage and three early phytophagous parasitic stages. A total of 867 DEGs were obtained.
Heatmap of the differentially expressed genes (DEGs) at the mycetophagous stage and three early phytophagous parasitic stages.
About this article
Cite this article
Hu, LJ., Wu, XQ., Ding, XL. et al. Comparative transcriptomic analysis of candidate effectors to explore the infection and survival strategy of Bursaphelenchus xylophilus during different interaction stages with pine trees. BMC Plant Biol 21, 224 (2021). https://doi.org/10.1186/s12870-021-02993-9
- Bursaphelenchus xylophilus
- Comparative transcriptome
- Candidate effectors
- Cell death suppression