Next Article in Journal
Thrombin Preconditioning Boosts Biogenesis of Extracellular Vesicles from Mesenchymal Stem Cells and Enriches Their Cargo Contents via Protease-Activated Receptor-Mediated Signaling Pathways
Previous Article in Journal
Role of Leukotriene B4 Receptor-2 in Mast Cells in Allergic Airway Inflammation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 12
10.3390/ijms20122898
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus)

by
Bowen Wang
1,
Xin Hao
1,
Jiayao Xu
1,
Yan Ma
2,* and
Ling Ma
1,*
1
College of Forestry, Northeast Forestry University, Harbin 150040, China
2
College of Management, Harbin University of Commerce, Harbin 150028, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(12), 2898; https://doi.org/10.3390/ijms20122898
Submission received: 9 April 2019 / Revised: 6 June 2019 / Accepted: 11 June 2019 / Published: 14 June 2019
(This article belongs to the Section Biochemistry)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: The causal agent of pine wilt disease is the pine wood nematode (PWN) (Bursaphelenchus xylophilus), whose ability to adapt different ecological niches is a crucial determinant of their invasion to colder regions. To discover the molecular mechanism of low temperature response mechanism, we attempted to study the molecular response patterns under low temperature from B. xylophilus with a comprehensive RNA sequencing analysis and validated the differentially expressed genes (DEGs) with quantitative real-time polymerase chain reaction (qRT-PCR). Bioinformatic software was utilized to isolate and identify the low-temperature-related BxGPCR genes. Transcript abundance of six low-temperature-related BxGPCR genes and function of one of the BxGPCR genes are studied by qRT-PCR and RNA interference. Results: The results showed that we detected 432 DEGs through RNA sequencing between low-temperature-treated and ambient-temperature-treated groups nematodes. The transcript level of 6 low-temperature-related BxGPCR genes increased at low temperature. And, the survival rates of BxGPCR17454 silenced B. xylophilus revealed a significant decrease at low temperature. Conclusion: in conclusion, this transcriptome-based study revealed a crucial role of BxGPCR17454 in low temperature response process of pine wood nematode. These discoveries would assist the development of management and methods for efficient control of this devastating pine tree pest.

Graphical Abstract

1. Introduction

As one of the most dangerous plant pests in the world, pine wood nematode (Bursaphelenchus xylophilus, PWN) causes devastating pine wilt diseases to the pine trees in Asia, Europe, and North America [1]. Unfortunately, the infestation area of PWN, in all probability, will continue expanding to colder regions of Asia and Europe [2,3,4,5]. For many parasites, PWN included, their ability to adapt different ecological niches is a crucial determinant of their invasion to new regions [6]. The formation of low-temperature-induced diapause stage of PWN is one of the most pivotal factors for PWN’s survival and expansion. In the autumn, PWN will gradually stop their growth and gradually enter winter diapause in response to the decreased environmental temperature. Additionally, the second-stage propagative juveniles (J2) gradually turn into specialized third-stage dauer larva (DL3), so that they can enter winter diapause together with the adults [7,8,9]. All the low-temperature-induced nematodes revealed similar characteristics: extended lifespan and accumulated fat content [8].
Although very little is known about the low temperature-sensing mechanism in PWN, much effort has been invested in probing the functions of some genes that respond to the low temperature in the model nematode Caenorhabditis elegans. The low temperature response process of C. elegans is regulated by many genes and groups including Cyclic guanosine monophosphate (cGMP) pathway, insulin-signaling pathway [10], G protein signaling [11], and phospholipid saturation [12], among others. The analysis of this widely studied model nematode inspired our analysis for the mechanism study of low temperature sensation in PWN. Some efforts have been made in the study of molecular mechanism of PWN’s response to low temperature in our previous research: we proved the key role of patched-related protein gene Bx-DAF6 [13], cGMP pathway [14], and stearoyl-coA desaturase gene Bx-SCD [15]. However, the mechanism of low-temperature-induced lifespan extension of PWN remains elusive.
In this study, we performed RNA sequencing on mixed-stage PWN between ambient temperature treated and low temperature treated groups. The results of RNA sequencing illustrated an essential role of G protein signaling during PWN’s early response to the low temperature. G protein coupled receptors (GPCRs) are a functionally diverse superfamily which are known as widely related to the environmental signal detection process [16,17,18]. Thus, we identified six GPCRs genes in G protein signaling from PWN and studied their response pattern to low temperature. We also validated the function of one GPCRs gene in low temperature-sensing with RNAi method. The results reveal that BxGPCR17454 plays a key role in PWN’s response to low temperature. This means the BxGPCR17454 in PWN may have the potential to be treated as a target for the management of this dangerous pest.

2. Results

2.1. RNA Sequencing and General Transcription Patterns

We performed RNA sequencing and obtained 312,138,316 high quality clean reads in all six samples after removal of rRNA and low-quality reads (Table 1). The read counts were converted to FPKM. Dendrogram clustering of the biological duplicates for each RNA sequencing condition reveals conserved alignment between replicates. (Figure 1A). The proportion of total reads that mapped to the reference genome ranged from 70.47% to 75.61%, and correlation values were significantly higher between the duplicated samples than among the treatments (Figure 1B).

2.2. B. xylophilus Genes Differentially Expressed in Response to Low Temperature

Between low-temperature-treated and CK nematodes, we identified 432 DEGs of which 314 were up-regulated, while 118 DEGs were down-regulated (Figure 2A, Table S3). Among all the DEGs, the gene expression patterns were similar between the biological duplicates but significantly different between the low-temperature-treated and CK group. (Figure 2B). GO enrichment showed that the DEGs mainly distributed in the single-organism process, cellular process, metabolic process, developmental process, and response to stimulus in biological process terms and binding and catalytic activity in molecular function terms (Figure 2C). Our results also indicated that low temperature influenced the expression of several gene groups with the putative function of environmental signal detection and stress adaptation (Table S3). These genes groups include four cytochrome P450s (BXY_0803400, BXY_0817800, BXY_1185500 and BXY_1697600), six heat shock proteins (BXY_0165400, BXY_0586000, BXY_0640100, BXY_0768000, BXY_1274600 and BXY_1563600) and 16 GPCRs (Table S4), among others.

2.3. Validation of DEGs by Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

12 DEGs were selected randomly for qRT-PCR validation. In general, the patterns of up-regulation and down-regulation were consistent with those obtained from the RNA sequencing analysis. Thus, qRT-PCR analysis confirmed that the changes detected in the RNA sequencing analysis were reliable (Figure 3, Table S2).

2.4. Identification and Transcript Abundance Analysis of 6 Low-temperature-related BxGPCRs

GPCRs are a group of genes which has been reported as closely related to the environmental signal detection process [16,17,18]. For this reason, we studied the relationship of low temperature stress and BxGPCRs expression. Structurally, GPCRs are also known as 7TM receptors because they have seven transmembrane (TM) domains [18,19]. Thus, we then used conserved domain analysis and transmembrane region prediction to filter the 16 differently expressed genes with GPCRs putative function description. BxGPCRs with incomplete conserved domain (Table S5) or with numbers of TM not equal to seven (Figure S1) are not further considered. There are six low-temperature-related BxGPCRs left after this filtering. According to these six genes’ ID BXY_1375800, BXY_1672700, BXY_0593200, BXY_0521000, BXY_1391500, and BXY_1745400, we named them BxGPCR13758, BxGPCR16727, BxGPCR05932, BxGPCR05210, BxGPCR13915, and BxGPCR17454, respectively. Blastp results showed that the deduced amino acid sequence of these six low-temperature-related BxGPCRs have a relatively high level of identity with other GPCR proteins of several nematodes. On this basis, we downloaded selected homologous amino acid sequences from NCBI (Table S6), performed alignment of the deduced amino acid sequences (Figure S2) and constructed phylogenetic tree (Figure 4A).
To validate the transcript pattern of BxGPCRs under low temperature, we measured transcript levels of six BxGPCRs at 5 °C for 1, 3, 5, and 7 days, and 25 °C for 1, 3, 5, and 7 days, respectively, with qRT-PCR (Figure 4B). The result showed that all of 6 BxGPCRs genes revealed higher transcript levels at low temperature than regular temperature over 7 days. The gene expression culminates at 3 days.

2.5. RNAi Validation of Low-temperature-related BxGPCR

The FITC treated nematodes reflected green fluorescence after exposed under ultraviolet light (Figure 5A), which indicated that the dsRNA can be fully absorbed by nematodes with soaking method. According to the result of expression changes between BxGPCR17454 dsRNA-treated and dsRNA-free nematodes, BxGPCR17454 gene was significantly silenced after dsRNA treatment. The mean expression level of BxGPCR17454 was decreased to 11.85% compared to the dsRNA-free nematodes. However, the dsRNA had no obvious effect on the transcript level of another internal control gene β-actin (Figure 5B). The alignment of six BxGPCRs nucleic sequences also illustrated that RNAi primer sequences is unique for BxGPCR17454 (Figure S3). These results indicated that the BxGPCR17454 RNAi was potent and specific in this study.
The results of the survival rates calculation showed that B. xylophilus revealed a significantly decreased survival rate after RNAi of BxGPCR17454 under low temperature compared to dsRNA-free group over 30 days (Figure 5C). This indicated that gene silencing of BxGPCR17454 can significantly decrease the survival rate of B. xylophilus at low temperature. BxGPCRs plays a key role in the low temperature response process of B. xylophilus. Although the transcript level of BxGPCR17454 was evaluated in the first 7 days according to the qRT-PCR results, there are no differences between the survival rate of dsRNA-free group and dsRNA-treated groups in the first 15 days. This may be because there are other downstream genes are influenced by the RNAi of BxGPCR17454. All these genes together with BxGPCR17454 influenced the survival rate of B. xylophilus under low temperature. Further studies need to be carried out to detect the downstream genes of BxGPCR17454 in the future.

3. Discussion

It is commonly known that low temperature response processes of both homeotherms and poikilotherms are not passive thermodynamic process but active genetic-promoted processes [20]. In this study, our results coincided with this theory. Transcriptome analysis revealed that low temperature response process of this plant parasite nematode is also a genetic process regulated by many genes. The past 15 years has witnessed a rapid progress in the understanding of how low temperature influenced the genetic response patterns of many plant pests with RNA sequencing method [21,22,23,24,25,26]. However, little is known about how B. xylophilus, which is also one of the most destructive pine tree pests, genetically responded to low temperature.
The identification and characterization of genes involved in low temperature responses is essential to elucidate low temperature defense mechanisms and develop effective control strategies. Changes in the transcriptome of B. xylophilus during low temperature directly reflect the impact of low temperature on the genetic activities. In this article, we conducted a comprehensive transcriptome analysis and characterized the gene expression profiles of B. xylophilus under low temperatures. Through the analysis of DEGs sets, we identified transcriptome changes in B. xylophilus in response to low temperature and validated the transcriptome result with qRT-PCR. The DEG analysis revealed that many gene groups including heat shock proteins, cytochrome P450s and GPCRs, among others, may potentially have played a crucial role in low temperature response process of B. xylophilus. On this basis, six low-temperature-related BxGPCRs genes were identified by phylogenetic analysis. In addition, we found these six low-temperature-related BxGPCRs genes revealed higher transcript abundance under low temperature than ambient temperature. Furthermore, the RNAi method was utilized to study the functions of one of the low-temperature-related BxGPCRs which we named BxGPCR17454. The results indicated that interference of BxGPCR17454 can significantly reduce the survival rate of B. xylophilus under low temperature. Taking all these factors into consideration, we hypothesized that G protein signaling may play a crucial role for the low temperature response process in B. xylophilus. This is widely consistent with previous research results in other organisms such as nematodes [10,11,27], bacterium [28], plants [29,30], fish [31], and mammals [32,33,34]. Interference of BxGPCR17451 may represent a novel approach for the prevention of this dangerous plant pest. Apart from G protein signaling, others DEGs such as heat shock proteins, are widely investigated as a big gene group regarding their role in thermal stress response, development and lifespan regulation [35,36,37,38]. Cytochrome P450s are not only reported as an import detoxification gene group [39], but also thermal-stress-related genes [40,41]. Further validations are expected to be made in these gene groups in the future. Other gene groups and pathways that have been reported to have a low temperature response function in model organism C. elegans include TRP channels [20,42,43], the cAMP pathway [44], and fatty acid metabolism pathway [45,46]; their orthologous genes in B. xylophilus may also share similar low temperature response functions. Further efforts are also expected to be made in the relevant fields in the future.
This study focused on the molecular response pattern of B. xylophilus under low temperature, identification and response pattern validation of six low-temperature-related BxGPCRs genes as well as functional validation of BxGPCR17454. The results shaped an important role of a G protein signaling gene BxGPCR17454 and potential value of other low-temperature-related genes in the low temperature response process of B. xylophilus. These discoveries would assist the development of management and methods for efficient control of this pine tree pest.

4. Material and Methods

4.1. Sample Preparation

B. xylophilus, maintained in the Forestry Protection Laboratory of Northeast Forestry University, Harbin, China, were kindly provided by the Chinese Academy of Forestry, Beijing, China. The nematodes were cultured on Botrytis cinerea for 5–7 days at 25 °C. The Baermann funnel technique was used to isolate B.xylophilus from potato dextrose agar medium plates. About 60,000 nematodes were mixed evenly in 6 mL distilled water. Then, suspensions were separated into six 1.5 mL centrifuge tubes equally. Three tubes were cultured at 5 °C for 24 hours as treated group. The other three tubes were cultured at 25 °C for 24 hours as CK group. After the incubation was finished, all samples were collected in the bottom of the tubes with centrifuge. Then, the tubes were all transferred in liquid nitrogen immediately for the next analysis.

4.2. RNA Sequencing

Total RNA of above six samples was extracted with Trizol reagent (Thermo Fisher Scientific, Shanghai, China). RNA-seq library was then constructed with NEB #7530 Kit (New England Biolabs, Ipswich, MA, USA) as following method: mRNA was enriched by Oligo(dT) beads. First, the enriched mRNA was fragmented into short fragments using NEBNext First Strand Synthesis Reaction Buffer. Random primers together with ProtoScript Ⅱ Reverse Transcriptase and Marine RNase Inhibitor were then used to reverse transcript mRNA into cDNA with following procedure: 25 °C 10 min, 42 °C 15 min, 70 °C 15 min. The products were then added with Second Strand Synthesis Reaction Buffer and Second Strand Synthesis Enzyme Mix to synthesis second-strand cDNA under 16 °C 60 min. Next, the cDNA fragments were purified with 1.8X Agencourt AMPure XP Beads. Sequencing adapters were then ligated with purified second-strand cDNA with follow steps: first, purified double stranded cDNA was pre-treated with NEBNext End Repair Reaction Buffer and NEBNext End Prep Enzyme Mix under 20 °C 30 min, 65 °C30 min. Then, adapters were added in the pre-treated second-strand cDNA with Blunt/TA Ligase Master Mix and Diluted NEBNext Adaptor under 20 °C 15 min. The ligation products were size selected by agarose gel electrophoresis, PCR amplified with 12-15 cycles, and sequenced using Illumina HiSeqTM 2500. RNA extraction, RNA-seq library construction and sequencing are performed by Gene Denovo Biotechnology Co. (Guangzhou, China).

4.3. Sequencing Data Analysis

Raw reads containing adapters, containing more than 10% of unknown nucleotides or containing more than 50% of low quality (Q-value ≤ 20) bases were removed by fastp [47] (Plant Pathology, SCRI, Invergowrie, Dundee DD2 5DA, UK. version 0.18.0) in order to get high quality clean reads. Short reads alignment tool Bowtie2 [48] was used for mapping reads to ribosome RNA (rRNA) database with default parameters. The rRNA mapped reads were removed. The remaining rRNA removed reads of each sample were then mapped to reference genome by TopHat2 [49] (Center for Bioinformatics and Computational Biology, University of Maryland, College Park, MD, USA. version 2.0.3.12). The mapping options used with TopHat are as follows: maximum read mismatch is 2, the distance between mate-pair reads is 50bp and the error of distance between mate-pair reads is ±80bp. The reference genome of B. xylophilus which was submitted directly to WormBase database by Taisei Kikuchi in 2011 [7] was downloaded from WormBase Parasite (BioProject PRJEA64437). The reference genome has 17,704 transcripts in total. We used RSEM software [50] (Department of Computer Sciences, University of Wisconsin-Madison, Madison, WI, USA) to calculate the gene expression level which was normalized by using FPKM (Fragments Per Kilobase of transcript per Million mapped reads) method. The edgeR package (Bioconductor, Roswell Park Cancer Institute, Buffalo NY, USA. http://www.r-project.org/) was used to identify differentially expressed genes with raw counts from RSEM across two groups. We identified genes with a fold change ≥2 and a false discovery rate (FDR) <0.05 in a comparison as significant differentially expressed genes (DEGs). BLASTX against the NCBInr database and Gene Ontology database was used to analysis gene description and annotations. Clustering, Pearson correlation analysis, volcano plot, heatmap, and GO categories are performed by R (Bioconductor, Roswell Park Cancer Institute, Buffalo NY, USA. version 3.2.1).

4.4. Validation of DEGs by Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

Quantitative real-time Polymerase Chain Reaction (qRT-PCR) was performed to validate the expression of DEGs. Gene-specific primers (Table S1) for 12 randomly selected DEGs were designed by Sangon Online Primer Designer (available online: https://www.sangon.com/newPrimerDesign) (Sangon Biotech Co. Shanghai, China). The Bx28S gene was used as the internal control gene (Table S1). Double-stranded cDNA was obtained using GoTaq 2-Step RT-qPCR System (Promega, Madison, WI, USA) according to instructions of the manufacturer. We used the cycle threshold data to calculate the fold of relative transcript abundance (treated group/CK group). The PCR program used in this analysis was first step: 95 for 2 min, second step: 95 for 15 s and 60 for 1 min, in 40 cycles.

4.5. Identification of Low-temperature-related BxGPCRs

According to the analyzed sequencing data, DEGs with a NR description of GPCRs were selected as candidate low-temperature-related BxGPCRs. Candidate genes with incomplete conserve domain were removed. The candidate genes whose predicted transmembrane domains unequal to seven were also removed, considering seven transmembrane domains is a typical feature of GPCRs [51,52]. Conserved domains were analyzed by NCBI Batch CD (available online: https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi) and transmembrane domains were predicted by the online tool TMHMM from the Center for Biological Sequence website (available online: http://www.cbs.dtu.dk/services/TMHMM/). BxGPCRs homologous sequences of the other organisms were obtained from BLASTP (available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi). Alignments of deduced amino acid sequences were performed by Muscle (EBI, UK) [53] and visualized by Jalview (version 2. 10. 4b 1). Phylogenetic analyses were performed with Mega 7.0 software using a maximum likelihood tree.

4.6. Transcript Abundance Analysis of 6 BxGPCRs under Low Temperature

The transcript abundance was analyzed as the method described before [14]. About 80,000 nematodes were mixed evenly in 8 mL distilled water. Then we separated the suspensions into eight 1.5 mL centrifuge tubes equally. Four tubes of nematodes were cultured at 5 °C for 1, 3, 5, and 7 days. The other four tubes of nematodes were cultured at 25 °C for 1, 3, 5, and 7 days. Total RNA of the nematodes cultured in these eight tubes were extracted and synthesized into cDNA with the method as described above after the incubation was finished. The transcript abundance and qRT-PCR Primers (Table S2) were also calculated and designed as described above. The whole transcript abundance analysis was analyzed with three replicates as three independent trials.

4.7. RNAi Validation of Low-Temperature-Related BxGPCR

We utilized RNAi method to study the function of one of the low-temperature-related BxGPCRs BxGPCR17454. The nematode RNAi was performed with soaking method as outlined in previous papers [14,39,54,55,56]. Mixed-stage B. xylophilus (male: female: juvenile ratio is approximately 1:1:2) was used in this research. MAXIscript T7/T3 RNA Synthesis Kit (Ambion, Tokyo, Japan) was used to obtain dsRNA with the following primers: i-BxGPCR17454-F (GCT AAT ACG ACT CAC TAT AGG GGG AAT CAT TGT ATA CCA GCC AAT TGC) and i-BxGPCR17454-R (AGT AAT ACG ACT CAC TAT AGG GTG GCC CGA TTG AAA TCC TGC). The size of BxGPCR17454 dsRNA for RNAi is 350 bp. The dsRNA-treated nematodes were soaked in ddH2O with final concentration of 2 mg/mL dsRNA corresponding to BxGPCR17454 sequence for 24 hours at 25 °C to uptake the dsRNA. The uptake of the dsRNA was monitored by a final nematode treatment of ddH2O containing 1 mg/mL FITC (fluorescein isothiocyanate) for 24 hours. The CK nematodes were soaked in ddH2O only. After intermittent stirring for 24 hours at 25 °C, the nematodes were washed with ddH2O to remove the external FITC and dsRNA. A fluorescence microscope was then used to take pictures of the nematode (ZEISS, Germany) to detect the uptake of FITC. Nematodes were divided into two groups after dsRNA uptake: the first group was used to examine the efficiency of RNAi with qPCR, and the second one was used to calculate survival rate under 5 °C environments. Total RNA was extracted from the CK nematodes and the dsRNA-treated nematodes after dsRNA. qRT-PCR was then performed as mentioned above with qRT-PCR primers of BxGPCR17454 and internal control primers Bx28S and β-actin [57] (Table S2). Survival rates of nematodes were calculated with 100 nematodes in 1 mL of ddH2O within each 1.5 mL centrifuge tube performed in triplicate. Then the survival rate of nematodes in each centrifuge tube was monitored every 3 days over 30 days. The survival rates of each tube were calculated as follows: number of living nematodes in each tube / number of total nematodes in each tube * 100%.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/12/2898/s1.

Author Contributions

B.W. and L.M. designed the work; B.W., X.H., and J.X. performed the experiments; B.W., J.X., Y.M., and L.M. analyzed the data; B.W. and X.H, contributed reagents/materials/analysis tools; B.W. wrote the paper.

Funding

This work was supported by grants from the Fundamental Research Funds for the Central Universities (2572018AA10); National Key Research and Development Program of China (2017YFD0600101).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mota, M.M.; Futai, K.; Vieira, P. Pine Wilt Disease and The Pinewood Nematode, Bursaphelenchus xylophilus. In Integrated Management of Fruit Crops Nematodes; Ciancio, A., Mukerji, K., Eds.; Springer: Dordrecht, The Netherlands, 2009; Volume 4, pp. 253–274. [Google Scholar]
  2. Wu, H.Y.; Tan, Q.Q.; Jiang, S.X. First report of pine wilt disease caused by Bursaphelenchus xylophilus on pinta thunbergii in the inland city of Zibo, Shandong, China. Plant Dis. 2013, 97, 1126. [Google Scholar] [CrossRef] [PubMed]
  3. Gruffudd, H.R.; Jenkins, T.A.R.; Evans, H.F. Using an evapo-transpiration model (ETpN) to predict the risk and expression of symptoms of pine wilt disease (PWD) across Europe. Boil. Invasions 2016, 18, 2823–2840. [Google Scholar] [CrossRef]
  4. Yu, Z.; Li, S.; Zhou, Y.; Zhou, R.; Wang, Y. Spatial estimation and prediction of suitable distribution of Bursaphelenchus xylophilus with different warming modes in China. J. Northeast For. Univ. 2018, 46, 85–91. [Google Scholar]
  5. Feng, Y.M.; Zhang, H.J.; Lü, Q.; Liang, J.; Zhang, X.Y. Quantification of suitability distribution region of Bursaphelenchus xylophilus in China. Scientia Silvae Sinicae 2009, 45, 65–71. [Google Scholar] [CrossRef]
  6. Devaney, E. Thermoregulation in the life cycle of nematodes. Int. J. Parasitol. 2006, 36, 641–649. [Google Scholar] [CrossRef] [PubMed]
  7. Kikuchi, T.; Cotton, J.A.; Dalzell, J.J.; Hasegawa, K.; Kanzaki, N.; McVeigh, P.; Takanashi, T.; Tsai, I.J.; Assefa, S.A.; Cock, P.J.A.; et al. Genomic Insights into the Origin of Parasitism in the Emerging Plant Pathogen Bursaphelenchus xylophilus. PLoS Pathog. 2011, 7, e1002219. [Google Scholar] [CrossRef]
  8. Zhao, L.L.; Wei, W.; Zhang, X.Y.; Kulhavy, D.; Sun, J.H. Low temperature induces two growth-arrested stages and change of secondary metabolites in Bursaphelenchus xylophilus. Nematology 2007, 9, 663–670. [Google Scholar] [CrossRef]
  9. Mamiya, Y. The life history of the pine wood nematode, Bursaphelenchus lignicolus. Jpn. J. Nematol. 1975, 5, 16–25. [Google Scholar]
  10. Ohta, A.; Ujisawa, T.; Sonoda, S.; Kuhara, A. Light and pheromone-sensing neurons regulates cold habituation through insulin signalling in Caenorhabditis elegans. Nat. Commun. 2014, 5, 4412. [Google Scholar] [CrossRef]
  11. Mertens, I.; Clinckspoor, I.; Janssen, T.; Nachman, R.; Schoofs, L. FMRFamide related peptide ligands activate the Caenorhabditis elegans orphan GPCR Y59H11AL.1. Peptides 2006, 27, 1291–1296. [Google Scholar] [CrossRef]
  12. Murray, P.; Hayward, S.A.L.; Govan, G.G.; Gracey, A.Y.; Cossins, A.R. An explicit test of the phospholipid saturation hypothesis of acquired cold tolerance in Caenorhabditis elegans. Proc. Natl. Acad. Sci. USA 2007, 104, 5489–5494. [Google Scholar] [CrossRef]
  13. Wang, F.; Li, D.L.; Ma, L.; Wang, B.W.; Chen, Q.L.; Zhang, R.Z.; Su, D.; Kang, X.Y.; Zhai, W. Identification and rnai of a Bursaphelenchus xylophilus dauer formation gene: Bx-daf6. J. Beijing For. Univ. 2016, 38, 21–27. [Google Scholar]
  14. Wang, B.; Ma, L.; Wang, F.; Wang, B.; Hao, X.; Xu, J.; Ma, Y. Low Temperature Extends the Lifespan of Bursaphelenchus xylophilus through the cGMP Pathway. Int. J. Mol. Sci. 2017, 18, 2320. [Google Scholar] [CrossRef]
  15. Wang, B.; Liu, W.; Wang, F.; Ma, L.; Li, D.; Wang, B.; Hao, X.; Northeast Forestry University. Fat accumulation in Bursaphelenchus xylophilus by positively regulating Bx-SCD under low temperature. J. Northeast For. Univ. 2017, 7. [Google Scholar]
  16. Teng, M.S.; Dekkers, M.P.J.; Ng, B.L.; Rademakers, S.; Jansen, G.; Fraser, A.G.; McCafferty, J. Expression of mammalian GPCRs in C. elegans generates novel behavioural responses to human ligands. BMC Boil. 2006, 4, 22. [Google Scholar]
  17. Cao, P.; Sun, W.; Kramp, K.; Zheng, M.; Salom, D.; Jastrzebska, B.; Jin, H.; Palczewski, K.; Feng, Z. Light-sensitive coupling of rhodopsin and melanopsin to G(i/o) and G(q) signal transduction in Caenorhabditis elegans. FASEB J. 2012, 26, 480–491. [Google Scholar] [CrossRef] [PubMed]
  18. Pierce, K.L.; Premont, R.T.; Lefkowitz, R.J. Seven-transmembrane receptors. Nat. Rev. Mol. Cell Biol. 2002, 3, 639–650. [Google Scholar] [CrossRef]
  19. Lagerström, M.C.; Schiöth, H.B. Structural diversity of G protein-coupled receptors and significance for drug discovery. Nat. Rev. Drug Discov. 2008, 7, 339–357. [Google Scholar] [CrossRef]
  20. Xiao, R.; Zhang, B.; Dong, Y.; Gong, J.; Xu, T.; Liu, J.; Xu, X.Z.S. A genetic program promotes C. elegans longevity at cold temperatures via a thermosensitive TRP channel. Cell 2013, 152, 806–817. [Google Scholar] [CrossRef] [PubMed]
  21. Huang, H.-J.; Xue, J.; Zhuo, J.-C.; Cheng, R.-L.; Xu, H.-J.; Zhang, C.-X. Comparative analysis of the transcriptional responses to low and high temperatures in three rice planthopper species. Mol. Ecol. 2017, 26, 2726–2737. [Google Scholar] [CrossRef] [PubMed]
  22. Cui, M.; Hu, P.; Wang, T.; Tao, J.; Zong, S. Differential transcriptome analysis reveals genes related to cold tolerance in seabuckthorn carpenter moth, Eogystia hippophaecolus. PLoS ONE 2017, 12, e0187105. [Google Scholar] [CrossRef] [PubMed]
  23. Schwarz, D.; Robertson, H.M.; Feder, J.L.; Varala, K.; Hudson, M.E.; Ragland, G.J.; Hahn, D.A.; Berlocher, S.H. Sympatric ecological speciation meets pyrosequencing: sampling the transcriptome of the apple maggot Rhagoletis pomonella. BMC Genom. 2009, 10, 633. [Google Scholar] [CrossRef] [PubMed]
  24. Qin, W.; Neal, S.J.; Robertson, R.M.; Westwood, J.T.; Walker, V.K. Cold hardening and transcriptional change in Drosophila melanogaster. Insect Mol. Boil. 2005, 14, 607–613. [Google Scholar] [CrossRef] [PubMed]
  25. Laayouni, H.; García-Franco, F.; Chávez-Sandoval, B.E.; Trotta, V.; Beltran, S.; Corominas, M.; Santos, M. Thermal evolution of gene expression profiles in Drosophila subobscura. BMC Evol. Boil. 2007, 7, 42. [Google Scholar] [CrossRef] [PubMed]
  26. Sørensen, J.G.; Nielsen, M.M.; Loeschcke, V. Gene expression profile analysis of Drosophila melanogaster selected for resistance to environmental stressors. J. Evol. Boil. 2007, 20, 1624–1636. [Google Scholar] [CrossRef]
  27. Ohta, A.; Kuhara, A. Molecular mechanism for trimetric G protein-coupled thermosensation and synaptic regulation in the temperature response circuit of Caenorhabditis elegans. Neurosci. Res. 2013, 76, 119–124. [Google Scholar] [CrossRef] [PubMed]
  28. Berchet, V.; Thomas, T.; Cavicchioli, R.; Russell, N.J.; Gounot, A. Structural analysis of the elongation factor G protein from the low-temperature-adapted bacterium Arthrobacter globiformis SI55. Extremophiles 2000, 4, 123–130. [Google Scholar] [CrossRef]
  29. Ma, Y.; Dai, X.; Xu, Y.; Luo, W.; Zheng, X.; Zeng, D.; Pan, Y.; Lin, X.; Liu, H.; Zhang, D.; et al. COLD1 Confers Chilling Tolerance in Rice. Cell 2015, 160, 1209–1221. [Google Scholar] [CrossRef] [Green Version]
  30. Xiong, L.; Schumaker, K.S.; Zhu, J.-K. Cell Signaling during Cold, Drought, and Salt Stress. Plant Cell 2002, 14, s165–s183. [Google Scholar] [CrossRef]
  31. Siebenaller, J.F.; Garrett, D.J. The effects of the deep-sea environment on transmembrane signaling. Comp. Biochem. Physiol. Part B Biochem. Mol. Boil. 2002, 131, 675–694. [Google Scholar] [CrossRef]
  32. Li, L.; Zhang, X. Differential inhibition of the TRPM8 ion channel by Gαq and Gα11. Channels 2013, 7, 115–118. [Google Scholar] [CrossRef] [PubMed]
  33. Clapham, D.E. TRP channels as cellular sensors. Nature 2003, 426, 517–524. [Google Scholar] [CrossRef] [PubMed]
  34. Bandell, M.; Story, G.M.; Hwang, S.W.; Viswanath, V.; Eid, S.R.; Petrus, M.J.; Earley, T.J.; Patapoutian, A. Noxious cold ion channel TRPA1 is activated by pungent compounds and bradykinin. Neuron 2004, 41, 849–857. [Google Scholar] [CrossRef]
  35. De Luca, F.; Di Vito, M.; Fanelli, E.; Reyes, A.; Greco, N.; De Giorgi, C. Characterization of the heat shock protein 90 gene in the plant parasitic nematode Meloidogyne artiellia and its expression as related to different developmental stages and temperature. Gene 2009, 440, 16–22. [Google Scholar] [CrossRef] [PubMed]
  36. Mosser, D.D.; Caron, A.W.; Bourget, L.; Denis-Larose, C.; Massie, B. Role of the human heat shock protein hsp70 in protection against stress-induced apoptosis. Mol. Cell. Boil. 1997, 17, 5317–5327. [Google Scholar] [CrossRef] [Green Version]
  37. Wentz, J.M.; Mendenhall, A.R.; Bortz, D.M. Pattern Formation in the Longevity-Related Expression of Heat Shock Protein-16.2 in Caenorhabditis elegans. Bull. Math. Boil. 2018, 80, 2669–2697. [Google Scholar] [CrossRef]
  38. Brunquell, J.; Morris, S.; Lu, Y.; Cheng, F.; Westerheide, S.D. The genome-wide role of HSF-1 in the regulation of gene expression in Caenorhabditis elegans. BMC Genom. 2016, 17, 559. [Google Scholar] [CrossRef]
  39. Xu, X.-L.; Wu, X.-Q.; Ye, J.-R.; Huang, L. Molecular Characterization and Functional Analysis of Three Pathogenesis-Related Cytochrome P450 Genes from Bursaphelenchus xylophilus (Tylenchida: Aphelenchoidoidea). Int. J. Mol. Sci. 2015, 16, 5216–5234. [Google Scholar] [CrossRef]
  40. Jia, K.; Albert, P.S.; Riddle, D.L. DAF-9, a cytochrome P450 regulating C. elegans larval development and adult longevity. Development 2002, 129, 221–231. [Google Scholar]
  41. Mak, H.Y.; Ruvkun, G. Intercellular signaling of reproductive development by the C. elegans DAF-9 cytochrome P450. Development 2004, 131, 1777–1786. [Google Scholar] [CrossRef]
  42. Chatzigeorgiou, M.; Yoo, S.; Watson, J.D.; Lee, W.-H.; Spencer, W.C.; Kindt, K.S.; Hwang, S.W.; Miller, D.M., 3rd; Treinin, M.; Driscoll, M.; et al. Specific roles for DEG/ENaC and TRP channels in touch and thermosensation in C. elegans nociceptors. Nat. Neurosci. 2010, 13, 861–868. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, B.; Xiao, R.; Ronan, E.A.; He, Y.; Hsu, A.-L.; Liu, J.; Xu, X.S. Environmental temperature differentially modulates C. elegans longevity through a thermosensitive TRP channel. Cell Rep. 2015, 11, 1414–1424. [Google Scholar] [CrossRef] [PubMed]
  44. Liu, F.; Xiao, Y.; Ji, X.-L.; Zhang, K.-Q.; Zou, C.-G. The cAMP-PKA pathway-mediated fat mobilization is required for cold tolerance in C. elegans. Sci. Rep. 2017, 7, 638. [Google Scholar] [CrossRef] [PubMed]
  45. Svensk, E.; Stahlman, M.; Andersson, C.H.; Johansson, M.; Boren, J.; Pilon, M. PAQR-2 regulates fatty acid desaturation during cold adaptation in C. elegans. PLoS Genet. 2013, 9, e1003801. [Google Scholar] [CrossRef] [PubMed]
  46. Brock, T.J.; Browse, J.; Watts, J.L. Fatty Acid Desaturation and the Regulation of Adiposity in Caenorhabditis elegans. Genetics 2007, 176, 865–875. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  48. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [Green Version]
  49. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Boil. 2013, 14, R36. [Google Scholar] [CrossRef]
  50. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef]
  51. Trzaskowski, B.; Latek, D.; Yuan, S.; Ghoshdastider, U.; Debinski, A.; Filipek, S. Action of Molecular Switches in GPCRs—Theoretical and Experimental Studies. Curr. Med. Chem. 2012, 19, 1090–1109. [Google Scholar] [CrossRef]
  52. Venkatakrishnan, A.J.; Deupi, X.; Lebon, G.; Tate, C.G.; Schertler, G.F.; Babu, M.M. Molecular signatures of G-protein-coupled receptors. Nature 2013, 494, 185–194. [Google Scholar] [CrossRef] [PubMed]
  53. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, F.; Wang, Z.; Li, D.; Chen, Q. Identification and Characterization of a Bursaphelenchus xylophilus (Aphelenchida: Aphelenchoididae) Thermotolerance-Related Gene: Bx-HSP90. Int. J. Mol. Sci. 2012, 13, 8819–8833. [Google Scholar] [CrossRef] [PubMed]
  55. Chen, Q.; Li, D.; Wang, F.; Zhang, R.; Ling, Y. Trehalose metabolism genes of Aphelenchoides besseyi (Nematoda: Aphelenchoididae) in hypertonic osmotic pressure survival. Boil. Open 2017, 6, 664–672. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, Q.; Wang, F.; Li, D.; Zhang, R.; Ling, Y. Trehalose metabolism genes render rice white tip nematode Aphelenchoides besseyi (Nematoda: Aphelenchoididae) resistant to an anaerobic environment. J. Exp. Biol. 2018, 221, jeb171413. [Google Scholar] [CrossRef] [PubMed]
  57. Li, Z.; Zhang, Q.; Zhou, X. A 2-Cys peroxiredoxin in response to oxidative stress in the pine wood nematode, Bursaphelenchus xylophilus. Sci. Rep. 2016, 6, 27438. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 02898 g001 550
Figure 1. General transcription patterns of all samples (A). Clustering of the biological duplicates for each RNA sequencing condition reveals conserved alignment between replicates. (B). Pearson correlation between low-temperature-treated and CK nematodes.
Figure 1. General transcription patterns of all samples (A). Clustering of the biological duplicates for each RNA sequencing condition reveals conserved alignment between replicates. (B). Pearson correlation between low-temperature-treated and CK nematodes.
Ijms 20 02898 g001
Ijms 20 02898 g002 550
Figure 2. Differentially expressed genes (DEGs) analysis. (A). DEGs between low-temperature-treated and CK group nematodes. Red dots represent up-regulated genes and green dots represent down-regulated genes. (B). Heat map of all the DEGs between low-temperature-treated and CK group nematodes. (C). Gene Ontology (GO) categories of DEGs between low-temperature-treated and CK group nematodes.
Figure 2. Differentially expressed genes (DEGs) analysis. (A). DEGs between low-temperature-treated and CK group nematodes. Red dots represent up-regulated genes and green dots represent down-regulated genes. (B). Heat map of all the DEGs between low-temperature-treated and CK group nematodes. (C). Gene Ontology (GO) categories of DEGs between low-temperature-treated and CK group nematodes.
Ijms 20 02898 g002
Ijms 20 02898 g003 550
Figure 3. Quantitative real-time polymerase chain reaction (qRT-PCR) identification of differentially DEGs between low-temperature-treated and CK group nematodes. The relative gene expressions were normalized with the internal control gene Bx28S. The putative function of each gene is as follows: BXY_0069900, NH127; BXY_0262500, SERA; BXY_0413400, ST2B1; BXY_0508700, LIP1; BXY_0596900, ASK12; BXY_0688500, NHR54; BXY_0794300, NPY2R; BXY_0799000, DUS9; BXY_0959300, ADH3; BXY_1220600, SOCS4; BXY_1423800, AATC; BXY_1722500, CYSK2.
Figure 3. Quantitative real-time polymerase chain reaction (qRT-PCR) identification of differentially DEGs between low-temperature-treated and CK group nematodes. The relative gene expressions were normalized with the internal control gene Bx28S. The putative function of each gene is as follows: BXY_0069900, NH127; BXY_0262500, SERA; BXY_0413400, ST2B1; BXY_0508700, LIP1; BXY_0596900, ASK12; BXY_0688500, NHR54; BXY_0794300, NPY2R; BXY_0799000, DUS9; BXY_0959300, ADH3; BXY_1220600, SOCS4; BXY_1423800, AATC; BXY_1722500, CYSK2.
Ijms 20 02898 g003
Ijms 20 02898 g004 550
Figure 4. Identification and transcript pattern analysis of six low-temperature-related BxGPCRs. (A). Analysis on phylogenetic trees of six deduced protein sequences with other organisms. Labels at the end of each branch are gene id/NCBI accessions of each gene. Different shapes stand for the Latin name of each organism. (B). Transcript abundance analysis of six BxGPCRs under low temperature.
Figure 4. Identification and transcript pattern analysis of six low-temperature-related BxGPCRs. (A). Analysis on phylogenetic trees of six deduced protein sequences with other organisms. Labels at the end of each branch are gene id/NCBI accessions of each gene. Different shapes stand for the Latin name of each organism. (B). Transcript abundance analysis of six BxGPCRs under low temperature.
Ijms 20 02898 g004
Ijms 20 02898 g005 550
Figure 5. RNAi validaton of BxGPCR17454. (A). B. xylophilus soaked in fluorescein isothiocyanate (FITC) revealed green fluorescent signal under ultraviolet light. (B). BxGPCR17454 RNAi efficiency was validated by expression changes between BxGPCR17454 dsRNA-treated and dsRNA-free nematodes. BxGPCR17454 gene was significant silenced after dsRNA treatment, while dsRNA had no obvious effect on the transcript level of internal control gene β-actin. (C). BxGPCR17454 dsRNA-treated nematodes showed significate decreased survival rate under low temperature. Data represent mean values ± SD from different repetitions. Asterisks indicates statistically significant differences (* p < 0.05, ** p < 0.01, Student’s t-test) was found between the dsRNA-treated and dsRNA-free groups.
Figure 5. RNAi validaton of BxGPCR17454. (A). B. xylophilus soaked in fluorescein isothiocyanate (FITC) revealed green fluorescent signal under ultraviolet light. (B). BxGPCR17454 RNAi efficiency was validated by expression changes between BxGPCR17454 dsRNA-treated and dsRNA-free nematodes. BxGPCR17454 gene was significant silenced after dsRNA treatment, while dsRNA had no obvious effect on the transcript level of internal control gene β-actin. (C). BxGPCR17454 dsRNA-treated nematodes showed significate decreased survival rate under low temperature. Data represent mean values ± SD from different repetitions. Asterisks indicates statistically significant differences (* p < 0.05, ** p < 0.01, Student’s t-test) was found between the dsRNA-treated and dsRNA-free groups.
Ijms 20 02898 g005
Table
Table 1. Statistics analysis of RNA sequencing data from low-temperature-treated and CK nematodes.
Table 1. Statistics analysis of RNA sequencing data from low-temperature-treated and CK nematodes.
SamplesClean ReadsHigh Quality Clean Reads (%)High Quality Clean BasesQ20 (%)Q30 (%)GC (%)Mapping Rate
Cold-151,107,56047,519,530 (92.98%)6,850,534,84895.08%86.34%48.70%75.24%
Cold-247,251,70843,777,894 (92.65%)6,230,763,13494.72%85.36%48.35%73.29%
Cold-343,275,87640,225,612 (92.95%)5,760,600,64594.92%85.93%48.62%75.36%
Warm-146,243,12242,555,170 (92.02%)6,101,120,21394.74%85.43%48.87%75.61%
Warm-245,717,82642,541,694 (93.05%)6,065,096,30094.87%85.83%48.74%75.56%
Warm-398,615,88488,413,574 (89.65%)12,105,837,90093.42%82.30%48.27%70.47%
Total332,211,976305,033,47443,113,953,040
Notes: Cold-1, Cold-2 and Cold-3 stands for low-temperature-treated nematodes sample 1, 2, and 3. Warm-1, Warm-2, and Warm-3 stands for CK group nematodes sample1, 2, and 3. Q20 stands for the percentage of bases with a Phred value >20. Q30 stands for the percentage of bases with a Phred value >30.

Share and Cite

MDPI and ACS Style

Wang, B.; Hao, X.; Xu, J.; Ma, Y.; Ma, L. Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus). Int. J. Mol. Sci. 2019, 20, 2898. https://doi.org/10.3390/ijms20122898

AMA Style

Wang B, Hao X, Xu J, Ma Y, Ma L. Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus). International Journal of Molecular Sciences. 2019; 20(12):2898. https://doi.org/10.3390/ijms20122898

Chicago/Turabian Style

Wang, Bowen, Xin Hao, Jiayao Xu, Yan Ma, and Ling Ma. 2019. "Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus)" International Journal of Molecular Sciences 20, no. 12: 2898. https://doi.org/10.3390/ijms20122898

APA Style

Wang, B., Hao, X., Xu, J., Ma, Y., & Ma, L. (2019). Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus). International Journal of Molecular Sciences, 20(12), 2898. https://doi.org/10.3390/ijms20122898

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop