芋球茎发育内参基因CeRPL28的筛选及其应用

芋球茎发育内参基因cerpl28的筛选及其应用
【技术领域】
1.本发明涉及植物基因技术领域，特别涉及芋球茎发育内参基因cerpl28的筛选及其应用。


背景技术：

2.芋(colocasia esculenta)又称芋头、芋艿、毛芋，是天南星科芋属的多年生草本植物。芋最初产于东南亚，之后向东传播到远东亚洲和大洋洲，向西传播到地中海地区，目前在全球范围内广泛种植，种植最多的地域为中国、日本和太平洋诸岛。芋是很多非洲和亚洲太平洋岛国居民的主要粮食，与其它块茎和块根植物相比，芋具有更高的营养价值。芋不仅可食用，也可药用，它味甘辛、平、性滑，其叶、花、球茎均可入药，具有宽肠胃、补脾胃、腹中癖块、消痨散结等作用，主治肿毒、牛皮癣、烫火伤等症。芋可食部分是球茎，是蛋白质、碳水化合物、脂肪、粗纤维、维生素c、硫胺素、核黄素、烟酸等的来源，淀粉是芋球茎最重要的成分，芋球茎的碳水化合物含量几乎是马铃薯的两倍，蛋白质含量也比山药、木薯、红薯高11％，由于存在如此丰富的营养物质含量，了解这些营养物质组成的遗传基础很重要。
3.近年来，随着高通量测序和分子生物学研究的不断发展，基因的表达分析也逐渐应用于揭示芋基因表达和调控机理的研究。随着芋基因组的测序，更加加快了芋相关性状的遗传机理解析。
4.实时荧光定量pcr(quantitative real-time pcr,qrt-pcr)是目前检测基因表达水平的主要分析手段。qrt-pcr结果的可靠性受到初始样本量、rna完整性、cdna质量和扩增效率等因素的影响。因此，表达稳定的内参基因对比较不同样本目的基因表达水平至关重要，正确选择合适的内参基因进行校正和标准化，才能获得准确的qrt-pcr结果。目前，植物中常用来作为内参基因的基因有肌动蛋白(actin)、微管蛋白基因(β-tubulin,tub)、转录延伸因子基因(elongation factor,ef1a和ef1b)、真核起始因子(eukaryotic initiation factor,eif)、泛素结合酶基因(ubiquitin-conjugating enzyme,ubc)、组蛋白(histone,h3-1)、甘油醛-3-磷酸脱氢酶(glyceraldehyde-3-phosphate dehydrogenase,gapdh)等。但是随着研究的深入，以往认为的很多内参基因(如actin等)在一定的实验条件下其表达量并不稳定。因此，在实际研究中，需要针对样本的不同组织或者不同生长发育时间段进行取材来筛选出稳定表达的内参基因。
5.本研究基于课题组前期对芋球茎不同发育阶段的样本进行转录本测序的数据，以60s核糖体蛋白基因家族基因表达量fpkm值的稳定性，使用变异系数cv进行评估，选择cv值≦0.2，的基因作为候选内参基因。利用qrt-pcr分析5个候选基因在芋球茎8组不同发育时期表达稳定性，利用三种内参基因稳定性分析软件genorm、normfinder和bestkeeper分析这5个候选内参基因在不同发育时期的稳定性并进行综合分析排序，为后续开展芋球茎发育相关功能基因表达研究提供了参考依据。


技术实现要素：

6.鉴于上述内容，有必要提供另一芋的内参基因，该基因可以很好的对芋的不同生长发育时期进行表达分析。
7.为达到上述目的，本发明所采用的技术方案是：
8.芋球茎发育内参基因cerpl28，所述内参基因cerpl28的核苷酸序列如序列表seq id no.1所示。
9.本发明还包括检测所述芋球茎发育内参基因cerpl28的引物对，所述引物对其上游序列如序列表seq id no.2所示，下游序列如序列表seq id no.3所示。
10.本发明还包括所述的芋球茎发育内参基因cerpl28在芋球茎发育不同时期进行荧光定量表达分析中的应用。
11.进一步的，所述芋球茎发育不同时期分别为：播种后30d、60d、90d、120d、150d、180d、210d和/或240d的球茎。
12.进一步的，所述荧光定量表达分析优选的引物对其上游序列如序列表seq id no.2所示，下游序列如序列表seq id no.3所示。
13.进一步的，所述荧光定量表达分析的pcr反应体系共20μl为：2x sybr green master mix 10μl、10μm上游引物0.4μl、10μm上游引物0.4μl、rnase-free water 8.2μl、cdna 1μl；pcr反应程序为95℃预变性3min；95℃变性15s，57℃退火15s，72℃延伸20s总共45个循环；采集荧光信号。
14.本发明的芋内参基因cerpl28，筛选方法为：
15.(1)对芋球茎整个发育期连续取样，每隔30d进行规律性取样，即分别取播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎，共8组样本，标记为s1、s2、s3、s4、s5、s6、s7、s8为芋球茎材料，提取芋球茎rna；
16.(2)反转录合成相应的芋球茎cdna；
17.(3)使用芋球茎6个发育阶段的转录组数据，以60s核糖体蛋白基因家族基因表达量fpkm值的稳定性，使用变异系数cv进行评估，选择cv值≦0.2的基因作为候选内参基因，根据不同候选基因设计不同引物并在荧光定量pcr仪上进行反应，通过bestkeeper、genorm和normfinder软件进行稳定性分析，获得稳定性好、表达效果好的候选基因为内参基因。
18.本发明具有如下有益效果：
19.1、本发明内参基因cerpl28是参考了2021年9月申请人在先申请《芋球茎发育过程中内参基因筛选及其应用》的分析方法进行的分析，该分析方法针对芋球茎的整个发育期进行连续取样，每隔30d进行规律性取样，即分别取播种后30d、60d、90d、120d、150d、180d、210d和240d的球茎，以该芋球茎材料，根据rna-seq测序数据分析，利用软件评估并分析各个候选内参基因的表达稳定性，获取一个新的内参基因cerpl28，为后期利用q rt-pcr探索芋球茎发育相关关键基因提供了一个新内参基因，为后期的芋发育研究提供更多的内参基因和更强有力的理论依据。
【附图说明】
20.图1-5分别为候选内参基因：cerpp3a、cerpl22、cerpl28、cerpl18和cerpl35的pcr扩增产物电泳检测结果图；其中，m为marker dl2000；泳道1-8分别为芋头s1、s2、s3、s4、s5、
s6、s7和s8。
21.图6-10分别为候选内参基因：cerpp3a、cerpl22、cerpl28、cerpl18和cerpl35的扩增曲线图；
22.图11-15分别为候选内参基因：cerpp3a、cerpl22、cerpl28、cerpl18和cerpl35的的溶解曲线图；
23.图16为五个候选内参基因在芋球茎不同发育时期的ct值箱线图；
24.图17为genorm软件对候选内参基因的表达稳定排序；
25.图18为genorm软件分析内参基因变异系数柱状图；
26.图19为cerpl28基因在芋头球茎中做为内参基因的表达情况图。
【具体实施方式】
27.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施的限制。
28.实施例1：
29.本实施例为内参基因cerpl28的获得，其通过如下方法得到：
30.一、本实施例的试验材料：
31.芋材料为广西农业科学院生物技术研究所选育的芋新品种桂芋2号。2020年3月份种植于广西农业科学院里建科研基地。球茎整个发育期连续取样，每隔30d进行规律性取样。分别取播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎，共8组样本，标记为s1、s2、s3、s4、s5、s6、s7、s8。每组样本三个生物学重复。将采集的球茎削皮，切成2cm左右的小块，放入50ml的无菌冷冻管中，液氮速冻后放入-80℃冰箱中备用。
32.二、总rna提取、检测与c dna合成：
33.使用试剂盒提取芋头球茎的总rna，使用1.5％琼脂糖凝胶电泳检测提取rna的完整性，利用nanodrop2100检测rna的浓度和质量。使用反转录试剂盒，将总rna反转录成cdna。总反应体系为8μl，4
×
gdna wiper mix 2μl，芋球茎总rna 1μg，加rnase-free water纯水补足至8μl。轻柔混匀后进行反转录反应，反应条件设置为：50℃、15min，85℃、5s，随即将cdna储存于-20℃冰箱备用。
34.三、内参基因的选择与特异性引物的设计：
35.使用芋球茎6个发育阶段的转录组数据(genbank登录号：prjna756190)，以60s核糖体蛋白基因家族基因表达量fpkm值的稳定性，使用变异系数cv进行评估，选择cv值≦0.2，的基因作为候选内参基因。从候选内参基因中，选取了5个fpkm在22.63～224.34的基因(表1)。利用primer premier 6软件设计5个候选内参基因的特异性引物，引物长度在101～178bp，具体引物如表2所示。
36.表1候选内参基因在球茎发育过程中的fpkm值
[0037][0038]
表2芋5个候选内参基因的引物
[0039][0040]
四、内参基因的qpcr：
[0041]
对内参基因进行荧光定量pcr检测，其中，荧光定量pcr的反应体系共20μl，上游引物(f)0.4μl(10μmol/l)、下游引物(r)0.4μl(10μmol/l)，cdna模板1μl，sybr green master mix(2x)10μl，rnase-free water 8.2μl补足。反应程序：95℃预变性、3min，95℃、15s；57℃退火15s；72℃延伸20s；共45个循环。上述引物选用上表2所述的引物。
[0042]
五、数据处理与分析
[0043]
利用genorm、normfinder和bestkeeper三个算法不相同的软件单独分析候选内参基因的稳定性高低。genorm和normfinder是利用2
‑△
ct
值表示其内参基因的表达稳定性。而bestkeeper根据各候选内参基因的ct值评价其稳定性。最后通过几何平均值法进行综合排序。
[0044]
六、结果与分析
[0045]
(一)茎总rna质量检测
[0046]
提取球茎发育8个阶段的rna后对其质量进行评估。电泳结果显示有两条完整条带，且28s亮度高于18s，说明rna完整性较好，无降解。nanodrop2100测得rna浓度为161.5～339.0ng/μl，od260/280均在1.91～2.11之间，表示rna纯度较高。
[0047]
候选内参基因的特异性分析
[0048]
对表1的5个候选基因采用表2的引物进行pcr凝胶电泳扩增检测，结果如图1-5所示，图中，m为marker dl2000，泳道1-8分别为候选内参基因：cerpp3a、cerpl22、cerpl28、cerpl18和cerpl35的pcr扩增产物电泳检测结果图；其中，m为marker dl2000；泳道1-8分别为芋头s1、s2、s3、s4、s5、s6、s7和s8。从图中可见，5个候选内参基因均可见与预期产物长度大小相一致的单一条带，不存在引物二聚体。
[0049]
同时对上述内参基因进行qrt-pcr的扩增曲线如图6-10所示，图中表明，扩增曲线均保持一定的曲度，扩增效果良好；上述内参基因qrt-pcr的溶解曲线参见图11-15所示，溶
解曲线中，只有一个单峰，表明扩增条带单一，特异性强，没有非特异性扩增出现。
[0050]
由此说明，上述实时荧光定量pcr引物特异性强，扩增效率高，能够用于芋球茎发育的荧光定量pcr的内参引物试验。
[0051]
候选基因的ct值分析
[0052]
ct值是表达基因丰富度的重要衡量标准，ct值与表达丰富程度呈反比关系。对5个内参基因在芋球茎不同发育时期进行表达丰度分析结果如图16所示，5个候选内参基因的平均ct值全部介于21.90～25.52之间，除了cerpl18表达量较低外，其他4个基因的整体表达量比较接近。其中，cerpp3a在球茎不同发育阶段的ct值为22.25～23.80，最大值和最小值相差1.55个循环；cerpl22在球茎不同发育阶段的ct值为21.73～23.61，最大值和最小值相差1.88个循环；cerpl28在球茎不同发育阶段的ct值为21.00～22.51，最大值和最小值相差1.51个循环；cerpl18在球茎不同发育阶段的ct值为24.92～26.41，最大值和最小值相差1.49个循环；cerpl35在球茎不同发育阶段的ct值为20.73～22.65，最大值和最小值相差1.92个循环
[0053]
内参基因的稳定性分析
[0054]
采用分析软件genorm、normfinder及bestkeeper软件作为辅助，对上述5个候选内参基因从多方面一致评估分析。结果如下：
[0055]
genorm软件对候选内参基因的稳定性评估
[0056]
genorm程序通过计算出每个内参基因稳定性的m值来筛选出稳定性较好的内参基因，判定标准为m值越小内参基因稳定性越好，反之，则稳定性越差，若m值大于1.5，则不考虑为内参基因。按照表达稳定性从高到低，依次为cerpp3a＝cerpl28(0.2999)﹥cerpl22(0.3156)﹥cerpl35(0.3284)﹥cerpl18(0.4128)。所有候选内参基因m值都小于1.5，均有可能是适合的内参基因。其中以cerpp3a和cerpl28的m值最小，稳定性也越好(图17)
[0057]
genorm软件还可计算引入1个新的内参基因后标准化因子的配对变异v值，被称作：“配对变量(vn/n+1)”。依据公式可判断内参基因的具体个数。为确保单一内参基因不引起偏差和波动，genorm软件基于计算得到配对差值vn/n+1来保证该条件下内参基因的最合适的数量。该公式中临界值定为0.15，若vn/n+1大于0.15，则选择的内参基因个数要满足n+1个；若vn/n+1小于0.15，则选择内参基因为n个时，就可满足该软件实验要求。genorm柱状图分析表明当选择v2/3值(0.096)﹤0.15(图18)，内参基因组合最佳数目是2个。因此，芋球茎不同发育时期最合适的内参基因有cerpp3a和cerpl28。
[0058]
normfinder软件对候选内参基因的稳定性评估
[0059]
normfinder软件程序算法与genorm程序类似，也是先获得内参基因表达稳定值，再根据稳定值大小来筛选最合适的内参基因，判定标准为表达稳定值最小的候选内参基因为合适内参基因。normfinder程序不但可以比较候选内参基因的表达差异，还可以计算样品组间的变异，但该程序只能筛选出一个最合适的内参基因。
[0060]
normfinder软件对5个候选内参基因的评估结果显示，内参基因表达稳定值cerpl28(0.070)﹥cerpp3a(0.097)﹥cerpl35(0.211)﹥cerpl22(0.214)﹥cerpl18(0.344)，在芋球茎不同发育过程最合适内参基因为cerpl28。结果与genorm软件评估结果基本一致。
[0061]
bestkeeper对候选内参基因的分析结果
[0062]
bestkeeper软件通过计算样品数据，获得基因间相关配对系数(r)，标准偏差(sd)
和变异相关系数(cv)，利用这3项参数判定内参的稳定性。当相关系数(r)越高、变异系数(cv)和标准偏差(sd)越低，则内参基因就越稳定。另外，当标准偏差sd大于1时，认为该候选内参基因表达不稳定。
[0063]
根据相关系数(r)排序，候选内参基因的稳定性排序为cerpl22(0.9799)﹥cerpl28(0.977)﹥cerpp3a(0.970)﹥cerpl35(0.948)﹥cerpl18(0.786)。根据sd值进行稳定性排序，具体如下：cerpl18(0.320)﹥cerpl28(0.500)﹥cerpp3a(0.550)＝cerpl35(0.550)﹥cerpl22(0.710)，但所有候选的5个内参基因sd值均小于1。根据变异系数(cv)值进行稳定性排序，具体如下：cerpl18(1.250)﹥cerpl28(2.270)﹥cerpp3a(2.380)﹥cerpl35(2.510)﹥cerpl22(3.130)。综合相关系数(r)越高、变异系数(cv)和标准偏差(sd)三项数值，cerpl28相关系数0.977，标准偏差0.500，变异系数2.270(具体参见表3)，对bestkeeper软件的分析结果进行综合排序，结果如表4所示。
[0064]
表3 bestkeeper软件分析内参基因稳定值
[0065]
参数cerpp3acerpl22cerpl28cerpl18cerpl35相关系数(r)0.9700.9790.9770.7860.948标准偏差(sd)0.5500.7100.5000.3200.550变异系数(cv)2.3803.1302.2701.2502.510
[0066]
表4 r、sd和cv值的综合分析
[0067] 相关系数(r)排序标准偏差(sd)排序变异系数(cv)排序几何平均值cerpp3a3433.30cerpl282222.00cerpl221352.47cerpl354544.31cerpl185111.71
[0068]
由表4可知，5个候选基因的几何平均值稳定性排序：cerpl18﹥cerpl28﹥cerpl22﹥cerpp3a﹥cerpl35。
[0069]
内参基因的稳定性分析
[0070]
通过对三种内参基因软件进行比较和分析，发现bestkeeper、genorm、normfinder的分析结果略有差异。因此，对3款软件的分析结果进行综合性评价；几何平均值越小，内参基因则越稳定。从表5可以看出cerpl28在三款软件的几何平均值最小，为最稳定的内参基因。因此，5个候选基因中，cerpl28为最合适的内参基因。
[0071]
表5 genorm、normfinder和bestkeeper三款软件综合分析结果
[0072] genorm排序nomfinder排序bestkeeper排序几何平均值cerpp3a1242.00cerpl281121.26cerpl222432.88cerpl353353.56cerpl184512.71
[0073]
综合genorm、normfinder和bestkeeper三款软件的分析结果，5个内参基因中，
cerpl28均为最合适的内参基因。
[0074]
实施例2：
[0075]
为验证本研究所筛选的cerpl28基因的稳定性，挑选了1个与球茎淀粉合成相关的基因淀粉分支酶(starch branching enzyme，cesbe)，对筛选的内参基因稳定性进行验证。以cerpl28作为内参基因，以cesbe为目的基因设计引物(表6)，对芋球茎8个发育阶段：播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎(标记为s1、s2、s3、s4、s5、s6、s7、s8)cdna为作为定量pcr的模板，进行qrt-pcr表达分析：
[0076]
表6内参基因和目的基因的引物序列
[0077][0078]
上述，pcr反应体系共20μl：2x sybr green master mix 10μl、10μm上游引物0.4μl、10μm上游引物0.4μl、rnase-free water 8.2μl、cdna 1μl；
[0079]
pcr反应程序为95℃预变性3min；95℃变性15s，57℃退火15s，72℃延伸20s总共45个循环；采集荧光信号。
[0080]
反应结束后对荧光定量pcr数据分析用spss 17.0和microsoft office excel 2007进行统计分析，采用spss17.0软件方差分析，采用excel 2010软件制图，结果如图19所示：cerpl28基因可在芋球茎发育过程中进行稳定表达；说明cerpl28基因可做为芋球茎发育不同时期荧光定量的内参基因。
[0081]
综上所述，基因cerpl28可在芋球茎发育过程中进行稳定表达，能做为芋球茎发育不同时期荧光定量的内参基因，能为今后对芋重要功能性状基因的研究奠定基础。
[0082]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明的保护范围应以所附权利要求为准。