Application of RNA Interference in the Pinewood Nematode, Bursaphelenchus xylophilus

摘要

Here, we introduce a detailed soaking method of RNA interference in Bursaphelenchus xylophilus to facilitate the study of gene functions.

摘要

The pinewood nematode, Bursaphelenchus xylophilus, is one of the most destructive invasive species worldwide, causing the wilting and eventual death of pine trees. Despite the recognition of their economic and environmental significance, it has thus far been impossible to study the detailed gene functions of plant-parasitic nematodes (PPNs) using conventional forward genetics and transgenic methods. However, as a reverse genetics technology, RNA interference (RNAi) facilitates the study of the functional genes of nematodes, including B. xylophilus.

This paper outlines a new protocol for RNAi of the ppm-1 gene in B. xylophilus, which has been reported to play crucial roles in the development and reproduction of other pathogenic nematodes. For RNAi, the T7 promoter was linked to the 5′-terminal of the target fragment by polymerase chain reaction (PCR), and double-stranded RNA (dsRNA) was synthesized by in vitro transcription. Subsequently, dsRNA delivery was accomplished by soaking the nematodes in a dsRNA solution mixed with synthetic neurostimulants. Synchronized juveniles of B. xylophilus (approximately 20,000 individuals) were washed and soaked in dsRNA (0.8 µg/μL) in the soaking buffer for 24 h in the dark at 25 °C.

The same quantity of nematodes was placed in a soaking buffer without dsRNA as a control. Meanwhile, another identical quantity of nematodes was placed in a soaking buffer with green fluorescent protein (gfp) gene dsRNA as a control. After soaking, the expression level of the target transcripts was determined using real-time quantitative PCR. The effects of RNAi were then confirmed using microscopic observation of the phenotypes and a comparison of the body size of the adults among the groups. The current protocol can help advance research to better understand the functions of the genes of B. xylophilus and other parasitic nematodes toward developing control strategies through genetic engineering.

引言

Plant-parasitic nematodes (PPNs) are a continuing threat to food security and forest ecosystems. They cause an estimated 100 billion USD in economic losses each year¹, the most problematic of which are primarily root-knot nematodes, cyst nematodes, and pinewood nematodes. The pinewood nematode, Bursaphelenchus xylophilus, is a migratory, endoparasitic nematode, which is the causal pathogen of pine wilt disease². It has caused great harm to pine forests worldwide³. Using the terminology of Van Megen et al.⁴, B. xylophilus is a member of the Parasitaphelenchidae and belongs to clade 10, whereas most other major plant parasites belong to clade 12.

As an independent and recently evolved plant parasite, B. xylophilus is an attractive model for comparative studies. To date, there has been substantial research on root-knot nematodes and cyst nematodes belonging to clade 12, which are obligate, sedentary endoparasites and are some of the most intensely studied nematodes. However, conducting further research in this important area comes with a major challenge: the function of parasitism genes is a research bottleneck. Functional studies generally include ectopic expression and knockdown/out experiments but rely on effective genetic transformation protocols for the nematode. As a result, reverse genetics in PPNs almost exclusively relies on gene silencing by RNAi.

RNAi, a mechanism widely present in eukaryotic cells, silences gene expression by introducing double-stranded RNA (dsRNA)⁵. To date, the posttranscriptional gene-silencing mechanism induced by dsRNA has been found in all studied eukaryotes, and RNAi technology, as a tool of functional genomics research and other applications, has developed rapidly in many organisms. Since the discovery of the RNAi machinery in Caenorhabditis elegans in 1998⁶, RNAi techniques have become effective methods for identifying the gene function of nematodes and are proposed as a new way to effectively control pathogenic nematodes⁷.

RNAi is technically facile-soaking the juveniles in dsRNA can suffice; however, the efficacy and reproducibility of this approach vary widely with the nematode species and the target gene⁸. The silencing of 20 genes involved in the RNAi pathways of the root-knot nematode, Meloidogyne incognita, was investigated using long dsRNAs as triggers, resulting in diverse responses, including an increase and no change in the expression of some genes⁹. These results show that target genes may respond to RNAi knockdown differently, necessitating an exhaustive assessment of their suitability as targets for nematode control via RNAi. However, there is currently a paucity of research on the developmental and reproductive biology of B. xylophilus.

As a continuation of previous work^10,11,12,13, we describe here a protocol for applying RNAi to study the function of the ppm-1 gene of B. xylophilus, including the synthesis of dsRNA, synthetic neurostimulant soaking, and quantitative polymerase chain reaction (qPCR) detection. The knowledge gained from this experimental approach will likely contribute markedly to understanding basic biological systems and preventing pine wilt disease.

方案

The study was approved by the council for animal experimentation of Zhejiang Agricultural & Forestry University. The B. xylophilus isolate NXY61 was originally extracted from a diseased Pinus massoniana in the Ningbo area of Zhejiang province, China¹¹.

1. Gene cloning

NOTE: See the Table of Materials for details about the primers used in this protocol.

1. Collect nematodes.
   1. Culture the B. xylophilus strain on the mycelia of Botrytis cinerea on Potato Dextrose Agar (PDA) plates at 25 °C for 3-5 days.
   2. Collect the nematodes using the Bellman funnel method¹⁴.
      1. Place a clamped rubber tube below a funnel and place two layers of filter paper in the mouth of the funnel. Transfer the fungal cultures to the funnel and add water to immerse the fungal mat. Wait for 2 h, then collect the nematodes.
2. Extract the total RNA from the nematodes using a total RNA extraction reagent (see the Table of Materials)¹¹ according to the following steps.
   1. Add 500 µL of extraction reagent and 100 µL of magnetic beads to a 2 mL centrifuge tube. Aspirate 20 µL of the nematodes and transfer the sample to a grinder for grinding at 9,000 × g for 30 s. Incubate for 5 min and then centrifuge for 10 min at 12,000 × g and 4 °C.
   2. Transfer the supernatant to a new centrifuge tube. Add 100 µL of chloroform, cap the tube, and mix by inverting the tube several times. Incubate for 3 min and then centrifuge for 10 min at 12,000 × g at 4 °C.
   3. Transfer the supernatant to a new centrifuge tube. Add 250 µL of isopropyl alcohol and vortex vigorously. Centrifuge at 12,000 × g for 10 min. 
   4. Discard the supernatant. Add 500 µL of 75% ethanol to wash the RNA and then vortex the sample. Centrifuge it for 5 min at 12,000 × g and 4 °C.
   5. Air-dry the RNA pellet for 5 min.
   6. Resuspend the pellet in 30 µL of RNase-free water.
   7. Calculate the RNA concentration using the formula: A260 × dilution × 40 = µg RNA/mL. Calculate the A260/A280 ratio.
      ​NOTE: A ratio of ~2 is considered pure.
3. Perform reverse transcription of good quality RNA to obtain the cDNA template.
   1. Design and use a pair of specific primers, ppm-1-F/R (see the Table of Materials), to amplify the partial coding sequence of the Bx-ppm-1 gene in B. xylophilus (GenBank accession number QTZ96795).
   2. Clone the ppm-1 gene sequences into the pGEM-Teasy vector containing the T7 promoter following a standard cloning protocol¹¹.
      1. Set up the PCR reactions as follows: 2 µL of cDNA, 25 µL of 2x Ex Taq Polymerase Premix, 2 µL of each primer (10 pmol/l), and sterile distilled water to a final volume of 50 µL.
      2. Perform the amplification procedure as follows: 5 min at 94 °C; followed by 35 cycles of 30 s at 94 °C, 30 s at 55 °C, and 1 min at 72 °C; and a final extension step at 72 °C for 5 min.
      3. Clone the amplified products into a pGEM-T Easy vector for sequencing.

2. Synthesis of dsRNA

1. Prepare the DNA template for dsRNA synthesis using PCR with primers designed to add T7 promoter sites at both ends. Add the T7 promoter sequence to the 5' end of the primers.
2. Use the plasmid containing the ppm-1 gene fragment (894 bp) as the template for PCR and recover the fragment containing the T7 promoter¹¹. Use the PCR procedure and system described above.
3. Use an in vitro transcription kit to synthesize dsRNA¹¹.
   1. Thaw the frozen reagents on ice.
   2. Add 2 µL of 10x reaction buffer, 2 µL of enzyme mix, and 1 µg of DNA to a centrifuge tube. Add nuclease-free water to produce a standard 4 µL reaction. Then, mix equal volumes of the four ribonucleotide solutions (ATP, CTP, GTP, and UTP) together and add 8 µL of the mixture to the tube. Mix thoroughly and incubate at 37 °C for 4 h.
   3. Add 1 µL of DNase, mix well, and incubate for 15 min at 37 °C.
   4. Stop the reaction and add 30 µL of nuclease-free water and 30 µL of LiCl precipitation solution to precipitate the RNA. Mix thoroughly. Incubate at -20 °C overnight.
   5. Centrifuge for 15 min at 12,000 × g and 4 °C. Discard the supernatant.
   6. Add 1 mL of 75% ethanol to wash the RNA. Vortex the sample and centrifuge it for 10 min at 12,000 × g and 4 °C.
   7. Air-dry the RNA pellet for 3 min.
   8. Resuspend the pellet in 30 µL of RNase-free water.
   9. Analyze the quality of the dsRNA using a spectrophotometer. Pipette 1 µL of the dsRNA sample onto the measurement pedestal and set the wavelength to 340 nm. Visualize the products on a 1.0% agarose gel.

3. RNAi by soaking

1. Mix 4 µL of 5x soaking buffer (0.05% gelatin, 5.5 mM KH₂PO₄, 2.1 mM NaCl, 4.7 mM NH₄Cl, 3 mM spermidine) with the dsRNA and ddH₂O to produce a total volume of 20 µL and a final RNA concentration of 0.8 µg/μL.
2. Acquire J2 larvae.
   1. Collect the nematodes from the fungal cultures and transfer them to a glass Petri dish 6 cm in diameter. Add 10 mL of water to the dish to ensure that the nematodes can swim freely. Keep the nematodes in the dish for 30 min and wait for the eggs to adhere to the bottom.
   2. Remove the water and nematodes carefully, making sure not to disturb the eggs. Repeat the steps until all the larvae and adults are removed, leaving only the eggs in the dish.
   3. Hatch the collected eggs for 24 h in the dark at 25 °C to obtain J2 larvae. Collect the J2 larvae, place them in a tube, and wash them three times with ddH₂O for the RNAi experiment¹⁵.
   4. Transfer the J2 larvae to the 2 mL tube containing the dsRNA solution and add resorcinol solution (wrapped in tinfoil and dissolved in water) to produce a final concentration of 1.0%. Incubate the larvae with centrifugation at 15 × g on a shaking table for 24 h at 25 °C to ensure that the larvae effectively absorb the dsRNA.
   5. Soak the same quantity of nematodes in the soaking buffer without the dsRNA probe or with a GFP dsRNA probe as a control. Use the GFP gene (gfp, M62653.1) as a nonendogenous control and synthesize the dsRNA of gfp using gene-specific primers T7-GFP-F/R.

4. qPCR detection

1. Clean the J2 larvae with ddH₂O, including those with target gene interference, GFP gene interference, and the undisturbed control group. Extract the total RNAs from each group using the method described above.
2. Start the qPCR.
   1. Design the q-ppm-1-F/R primers using the desired software (see the Table of Materials).
   2. Set up the qPCR reaction in 12 µL containing 1 µL of cDNA, 6 µL of fluorescent premix, 0.4 µL of each primer (10 pmol/l), and sterile distilled water.
   3. Perform qPCR as follows: 2 min at 95 °C; followed by 40 cycles of 10 s at 95 °C, 30 s at 55 °C, and 1 min at 72 °C.
   4. Use the actb gene (GenBank accession number EU100952) and tbb-2 gene (GenBank accession number MT769316), or other genes, as internal reference genes to evaluate the changes in gene expression level after RNAi¹⁵.
3. According to the cycle threshold (Ct) value and the dissolution curve, use the 2^-ΔΔCt method to estimate the relative expression level of the target gene and verify the interference efficiency.
   1. Subtract the Ct value of the internal reference gene of each sample from the Ct value of the target gene to obtain the ΔCt value. Then, subtract the ΔCt value of the interference group from the ΔCt value of the control group to obtain the ΔΔCt value.
      ​NOTE: A ΔΔCt value greater than 0 indicates that the interference is effective.

5. Evaluate the body length of nematode adults following RNAi

1. After RNAi, culture the J2 larvae until adulthood on B. cinerea lawns in PDA plates for 60 h at 25 °C¹⁵.
2. Collect the adults using the Bellman funnel method (see step 1.1.2.1)¹⁴.
3. Acquire images of the adult nematodes under a microscope and use ImageJ software (or other measurement software) to measure the body lengths.
   1. Measure the length using ImageJ by selecting Analyze | Set Scale. If the distance is known, enter the length value of the drawn straight line. Enter the unit of length.
   2. Check the Global check box (use this standard for all pictures), and click OK to select the length to be measured with a straight line on the image to be measured.
   3. Use the command Ctrl+M (measurement) to display the results and record them in the results window.
4. Take measurements of 50 male and 50 female nematodes for statistical analysis¹².
5. Analyze the data by calculating the mean and standard deviation for each sample. Compare the means of the samples from the different groups using the Student's t-test.

结果

Analysis of ppm-1 expression of B. xylophilus after RNAi
The relative expression level of the ppm-1 gene of B. xylophilus soaked with GFP dsRNA and that soaked with target gene dsRNA was 0.92 and 0.52, respectively (the ppm-1 gene expression level of the ddH₂O-treated control group was set to 1) (Figure 1). Thus, exogenous dsRNA has no effect on the...

讨论

Although the life history and parasitic environment of B. xylophilus are different from those of other nematodes, there has been limited research on the molecular pathogenesis of this plant pathogen. Despite great progress made in the application of CRISPR/Cas9 genome editing technology in C. elegans and other nematodes, only RNAi technology applied to B. xylophilus has been published to date¹⁷. RNAi is one of the most powerful tools available to study the gene function ...

材料

本文使用的材料清单

姓名	公司	目录编号	评论
贝尔曼烟囱	无	无	分离线虫
信标设计师7.9	上海康玉生信息技术有限公司	无	设计qPCR引物
铁锈木	无	无	作为线虫的食物
Bursaphelenchus xylophilus	无	无	其编号为NXY61，最初是从中国浙江宁波的病 松松中提取的。
恒温孵化器	上海晶宏实验仪器有限公司	H1703544	与文化线虫
电泳装置	生物辐射实验室	1704466	实现电泳分析
乙醇，75%	中立药化学试剂公司	80176961	提取RNA
Ex Taq聚合酶预混剂	宝生物株式会社	RR030A	用于PCR
Ex Taq聚合酶预混剂	宝生物株式会社	RR390A	用于PCR
凝集成像仪	朗吉恩科学仪器公司	LG2020	使核酸带可见
GraphPad 棱镜8	GraphPad 棱镜	无	分析数据并制作图形
高速离心机	杭州全盛仪器有限公司	AS0813000	离心机
高通量组织研磨机	贝尔坦		提取RNA
ImageJ 软件	美国国立卫生研究院	无	测量体长
异丙醇	上海阿拉丁生化技术有限公司	L1909022	提取RNA
徕卡 DM4B 显微镜	徕卡微系统公司		观察线虫
磁珠	奥兰科技公司	150010C	提取RNA
MEGAscript T7 高产转录套件	赛莫飞世尔科学公司	AM1333	体外合成双重RNA
NanoDrop ND-2000 分光光度计	赛莫飞世尔科学公司	纳米滴落 2000/2000C	分析dsRNA的质量
PCR放大器	Bio-Rad Life Medical Products Co.	1851148	扩增核酸序列
培养皿	无	无	与文化线虫
pGEM-T 简单向量	普罗梅加公司	A1360	用于克隆
土豆葡萄糖琼脂（中等）	无	无	培养Botrytis cinerea
Prime Script RT 试剂套件带 gDNA 橡皮擦	宝生物株式会社	RR047B	转为合成cDNA
Primer Premier 5.0	PREMIER Biosoft	无	设计PCR引物
引物：PPM-1-F/R	青柯生物技术公司	无	女：5'-GATGCGAAGTTGCCAATCATCTT -3';射门：5'- CCAGATCCAGTCCACCATACACCACC -3
q-ppm-1-F/R	青柯生物技术公司	无	前锋：5'-CATCCGAATGGCAATACAG-3';R：5'-ACTATCCTCAGCGTTAGC-3'
实时热循环器及nbsp;qTOWER 2.2	分析耶拿仪器（北京）公司		对于qPCR
摇桌	上海志成分析仪器制造有限公司		浸泡线虫
立体显微镜	重庆光学仪器公司	1814120	观察线虫
T7-GFP-F/R	青柯生物技术公司	无	F： 5'-TAATACGACTCACTATAGGGAAA GGAGAAGAACTTTTCAC-3';R： 5'-TAATACGACTCACTATAGGGCTG TTACAAACCAAGAAGG-3'
T7 启动子	青柯生物技术公司	无	TAATACGACTCACTATAGGG
Takara MiniBEST 琼脂糖凝胶 DNA 提取套件	宝生物株式会社	9762	以回收DNA
TaKaRa TB 绿色预混 Ex Taq（Tli RNaseH Plus）	宝生物株式会社	RR820A	对于qPCR
三氯乙烷	上海灵风化剂公司		提取RNA
三硫试剂	赛莫飞世尔科学公司	15596026	总RNA提取试剂，用于提取RNA