一种泡桐丛枝病相关基因网络分析方法与流程

1.本发明属于分子生物学技术领域，具体涉及一种泡桐丛枝病相关基因网络分析方法。


背景技术：

2.泡桐，树皮灰色、灰褐色或灰黑色，幼时平滑，老时纵裂；假二杈分枝，单叶，对生，叶大，卵形，全缘或有浅裂，具长柄，柄上有绒毛；喜光，较耐阴，喜温暖气候，耐寒性不强；其生长速度快，耐瘠薄，材性优良，在防风固沙、保障粮食安全和木材供应等方面起着重要的作用。
3.在生产中，由植原体引起的泡桐丛枝病是制约泡桐产业发展的重要因素，由于泡桐丛枝植原体不能体外培养，严重限制了泡桐丛枝病发病机理的阐明。泡桐丛枝病的发生、发展是一个非常复杂的过程，目前泡桐丛枝病发病机理仍未阐明，丛枝病的防治的有效方法仍未建立，因此，还需要开辟新的途径进行研究。
4.因此，如何提供一种泡桐丛枝病的分析方法，为泡桐丛枝病的防治提供依据是本领域亟待解决的问题。


技术实现要素：

5.本发明公开了一种泡桐丛枝病相关基因网络分析方法，可用于实现相关编码rna和非编码rna的联合分析、cerna等，从而快速全面准确地获得与特泡桐丛枝病相关的rna数据信息，进而将调控机制研究延伸到“网、多层面”结合的立体模式，有助于泡桐丛枝病的防治研究。
6.为了实现上述目的，本发明采用如下技术方案：
7.一种泡桐丛枝病相关基因网络分析方法，包括下述步骤：
8.(1)建立泡桐丛枝病发生模拟系统
9.利用不同浓度的甲基磺酸甲酯mms和/或利福平rif处理泡桐丛枝病苗，模拟泡桐丛枝病恢复和发病的过程，确定测序取样时间点；
10.(2)全转录组测序及文库构建
11.提取泡桐总rna，进行总rna样品检测，构建全转录组文库；
12.所述全转录组文库包括小rna文库和去除核糖体rna的链特异性文库；
13.所述小rna文库用于获得mirna序列信息；
14.所述去除核糖体rna的链特异性文库用于获得mrna、lncrna和circrna的序列信息；
15.(3)生物信息学分析
16.分别对mirna、mrna、lncrna和circrna数据进行数据分析，以及cerna分析，筛选泡桐丛枝病相关调控rna及cerna网络。
17.mms和/或rif处理可有效模拟泡桐丛枝病恢复和发病的过程，进而为调控rna的筛
选提供有利条件，保证mirna、mrna、lncrna和circrna分析以及cerna分析的可靠性。
18.优选地，步骤(1)中，
19.60mg
·
l
‑1mms或100mg
·
l
‑1rif处理白花泡桐病苗，用于使病苗恢复健康；
20.20mg
·
l
‑1mms或30mg
·
l
‑1rif处理白花泡桐病苗，用于使病苗恢复健康，但此时，相当于无症状感染者，将恢复健康的病苗继代后，感病状态又出现。
21.优选地，步骤(1)中，
22.60mg
·
l
‑1mms或100mg
·
l
‑1rif处理白花泡桐病苗后，取样时间为处理后5
‑
20d多时间点取样；
23.20mg
·
l
‑1mms或30mg
·
l
‑1rif处理白花泡桐病苗后，取样时间为处理后10
‑
30d至少两次取样，继代后20
‑
40d至少两次取样。
24.优选地，步骤(2)中，
25.小rna文库测序策略为50se；
26.去除核糖体rna的链特异性文库测序策略为150pe。
27.优选地，步骤(3)中，
28.数据分析包括基本分析、mirna、mrna、lncrna、circrna的表达水平分析、差异表达分析、go富集分析、kegg富集分析和cerna分析。
29.优选地，步骤(3)中，
30.筛选泡桐丛枝病相关调控rna时，包括如下步骤：
31.1)筛选出泡桐丛枝病恢复过程表达量升高或下降的rna，以及泡桐丛枝病发病过程表达量下降或升高的rna，
32.2)比较获得泡桐丛枝病恢复过程与桐丛枝病发病过程表达量变化趋势相反的rna；
33.3)对于同种处理试剂的不同剂量组，筛选出共同确定的rna。
34.优选地，步骤(3)中，
35.对于不同处理试剂组，分别按照1)
‑
3)筛选出调控rna，再比较筛选出不同处理试剂组共同确定的rna。
36.优选地，步骤(3)中，
37.筛选出的泡桐丛枝病相关mirna包括：pf
‑
mir122、pf
‑
mir156f
‑
5p、pf
‑
mir159a
‑
3p、pf
‑
mir169f；
38.筛选出的泡桐丛枝病相关mrna包括：多聚半乳糖醛酸酶、cullin 1、泛素结合酶e2、squamosa启动子结合蛋白3、squamosa启动子结合蛋白6、脱落酸受体pyr1、肌醇加氧酶、葡糖醛酸激酶、延伸因子1
‑
α、mads box蛋白、9
‑
顺式
‑
环氧类胡萝卜素双加氧酶、泛素化结构域蛋白dsk2a、质体移动受损蛋白2、蓝光下的叶绿体运动蛋白1、泛素化受体rad23c、叶绿体积累响应所需的j结构域蛋白、branched1
‑
like相关mrna；
39.筛选出的泡桐丛枝病相关lncrna包括：mstrg.18394.1、mstrg.22708.1、mstrg.29927.1、mstrg.33514.1；
40.筛选出的泡桐丛枝病相关circrna包括：circrna7714、circrna315、circrna6856、circrna466、circrna2705、circrna6438。
41.优选地，步骤(3)中，
42.分析得到的泡桐丛枝病相关cerna包括circrna7714
‑
mir156f
‑
5p
‑
spl3/6。
43.综上所述，本发明通过建立泡桐丛枝病发生模拟系统、取样全转录组测序、构建文库、对mirna、mrna、lncrna、circrna进行生物信息学分析，筛选出筛选泡桐丛枝病相关调控rna及cerna网络，为泡桐丛枝病的防治提供有力依据。
附图说明
44.附图1所示为泡桐丛枝病相关rna表达量变化趋势图；其中，
45.a.60mg
·
l
‑1mms模拟泡桐丛枝病恢复过程中rna表达量呈下降趋势；
46.b.60mg
·
l
‑1mms模拟泡桐丛枝病恢复过程中rna表达量呈升高趋势；
47.c.20mg
·
l
‑1mms模拟泡桐丛枝病恢复过程中rna表达量呈下降趋势；
48.d.20mg
·
l
‑1mms模拟泡桐丛枝病恢复过程中rna表达量呈上升趋势；
49.e.20mg
·
l
‑1mms模拟泡桐丛枝病发病过程中rna表达量呈下降趋势；
50.f.20mg
·
l
‑1mms模拟泡桐丛枝病发病过程中rna表达量呈上升趋势
51.图2所示为泡桐丛枝病相关cerna网络。
52.图3所示为qrt
‑
pcr验证结果。
53.图4所示为毛果杨中mir156f
‑
5p过表达验证结果。
具体实施方式
54.下面对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
55.实施例1
56.一种泡桐丛枝病相关基因网络分析方法，包括下述步骤：
57.1.建立泡桐丛枝病发生模拟系统
58.利用不同浓度的甲基磺酸甲酯(mms)和利福平(rif)处理白花泡桐丛枝病苗，模拟泡桐丛枝病恢复和发病的过程，确定测序取样时间点。
59.60mg
·
l
‑1mms或100mg
·
l
‑1rif处理白花泡桐病苗后，可使病苗恢复健康；取样时间为处理后5、10、15、20d。
60.20mg
·
l
‑1mms或30mg
·
l
‑1rif处理白花泡桐病苗后，也可使病苗恢复健康，但此时，相当于无症状感染者，将恢复健康的病苗继代后，感病状态又出现。取样时间为处理后10、30d，继代后20、40d。
61.在上述确定的取样时间点，取长势均匀的泡桐顶芽，液氮速冻保存于
‑
80℃冰箱，为后续试验做准备。并且设置健康苗组pf和患病苗组pfi作为对照。
62.每个样品三个生物学重复，各样品分组如表1所示。
63.表1
64.[0065][0066]
2.全转录组测序及文库构建
[0067]
提取各泡桐样品总rna，进行总rna样品检测，构建全转录组文库。
[0068]
全转录组文库包括小rna文库和去除核糖体rna的链特异性文库；小rna文库用于获得mirna序列信息；去除核糖体rna的链特异性文库用于获得mrna、lncrna和circrna的序列信息。
[0069]
采用50se测序策略构建小rna文库：
[0070]
样品总rna经检测合格后，采用truseq small rna sample prep kits(illumina,san diego,usa)试剂盒制备小rna测序文库。首先连接3’和5’端rna接头序列，再使用试剂盒进行逆转录反应，将带有接头的srna片段反转录成cdna，然后进行pcr扩增。扩增产物经过page胶电泳后，回收片段构建文库。库制备工作完成后，对构建好的文库使用illumina hiseq2500进行测序，测序读长为单端1
×
50bp。
[0071]
采用150pe测序策略构建去除核糖体rna的链特异性文库：
[0072]
样品总rna经检测合格后，采用illumina公司ribo
‑
zero gold试剂盒去除核糖体rna，纯化收回剩余rna。纯化回收的rna被片段化试剂(franmentation buffer)随机打断成
短片段，以此为模板，用六碱基随机引物(random hexamers)合成第一链cdna，随后加入缓冲液、dntp和dnapolymerase i进行第二链cdna合成。ampure xpbeads纯化双链产物，然后将dna的粘性末端修复为平末端，3’末端加接头，ampure xpbeads进行片段选择，最后进行pcr扩增获得最终测序文库，文库质检合格后采用illuminahiseq4000进行测序，测序读长为2
×
150bp(pe150)。
[0073]
3.mirna分析
[0074]
mirna测序获得原始数据705,982,419个reads。过滤掉除3
′
接头序列，碱基长度小于18nt的序列，富含(80％)a或c或g或t，只有a、c没有g、t，或只有g、t没有a、c，或者是连续的核苷酸二聚体，三聚体序列。然后测得序列与mrna，rfam(包含rrna、trna、snrna和snorna等)和repbase数据库进行比对，结果如表2所示。
[0075]
表2 mirna文库比对统计结果
[0076]
[0077][0078]
过滤后得到的cleanreads有401,203,924个，clean reads用于mirna的比对鉴定和预测分析。通过与白花泡桐基因组和mirbase的比对，共鉴定到668个mirna，其中427个是已知mirna，241个是新mirna。
[0079]
响应丛枝病mirna的筛选：
[0080]
mirna表达量以归一化后的nor值进行分析。
[0081]
在mms处理组中，60mg
·
l
‑1mms组模拟泡桐丛枝病恢复的过程，找出表达量呈升高(图1b)或下降(图1a)的mirna；20mg
·
l
‑1mms组模拟泡桐丛枝病恢复及发病的过程，在恢复过程中，鉴定出表达量升高(图1d)或下降(图1c)的mirna，发病过程中，找出表达量下降(图1e)或升高(图1f)的mirna，再找出恢复和发病过程中表达量变化趋势相反的转录本。
[0082]
为了消除mms的影响，找出60mg
·
l
‑1mms组和20mg
·
l
‑1mms组共同确定到的泡桐丛枝病相关mirna。
[0083]
对于rif处理组，100mg
·
l
‑1rif处理组参考60mg
·
l
‑1mms组，30mg
·
l
‑1rif处理组参考20mg
·
l
‑1mms组，再找出共同确定到的泡桐丛枝病相关mirna。
[0084]
最终，为了消除试剂的影响，找出mms处理组和rif处理组共同确定到的响应泡桐丛枝病mirna，如表3所示。
[0085]
表3 泡桐丛枝病相关mirna
[0086][0087]
4.mrna分析
[0088]
采用去除rrna的建库方式进行测序，整个测序过程中既能检测lncrna和mrna的表达水平，也能检测circrna的表达水平。结果如表4所示，本次测序得到12,423,427,214个clean reads，共鉴定到28,514个转录本。
[0089]
表4 去除核糖体rna的链特异性文库与白花泡桐参考基因组比对统计结果
[0090][0091][0092]
注：valid reads：质控后的数据；mapped reads：能比对到基因组上的reads数；
per million fragments))为单位，具体计算公式如下:
[0102][0103]
设fpkm(a)为基因a的表达量，则c为唯一比对到基因a的fragments数量，n为唯一比对到参考基因的总fragments数，l为基因a编码区的碱基数。
[0104]
参考响应泡桐丛枝病mirna的筛选方法，得到的泡桐丛枝病相关lncrna如表6所示。
[0105]
表6 泡桐丛枝病相关lncrna
[0106][0107]
6.circrna分析
[0108]
共鉴定到5,118个circrna，circrna表达定量的结果以fpkm(fragments per kb permillion fragments)为单位，具体计算公式如下:
[0109][0110]
设fpkm(a)为基因a的表达量，则c为唯一比对到基因a的fragments数量，n为唯一比对到参考基因的总fragments数，l为基因a编码区的碱基数。参考响应泡桐丛枝病mirna的筛选方法，得到的泡桐丛枝病相关circrna如表7所示。
[0111]
表7 泡桐丛枝病相关circrna
[0112][0113]
7.cerna分析
[0114]
鉴定到的28,514个转录本、5,118个lncrna、9,927个circrna和668个mirna构建circrna
‑
mirna
‑
mrna和lncrna
‑
mirna
‑
mrna调控网络，结果显示共获得到36,531条lncrna
‑
mirna
‑
mrna和84,445条circrna
‑
mirna
‑
mrna互作关系。重点关注鉴定到的与丛枝病有关的17个转录本，4个lncrna，6个circrna，4个mirna涉及的调控关系，结果见图2，从中发现了circrna7714
‑
pf
‑
mir156f
‑
5p
‑
spl3/6是这一关键的cerna。
[0115]
8.circrna7714
‑
pf
‑
mir156f
‑
5p
‑
spl3/6调控研究
[0116]
首先，进行qrt
‑
pcr验证，结果见图3，从中发现circrna7714、spl3、spl6感病后下调，mir156f
‑
5p上调，与测序结果一致。
[0117]
进一步地，在毛果杨(populus trichocarpa)中进行mir156f
‑
5p过表达分析，结果显示mir156f
‑
5p过表达的转基因苗中呈现了分枝增多，植株矮小，叶片减小现象(图4)。
[0118]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
[0119]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对上述实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。