CN113755497B - 芋球茎发育过程中内参基因筛选及其应用 - Google Patents

芋球茎发育过程中内参基因筛选及其应用 Download PDF

Info

Publication number
CN113755497B
CN113755497B CN202111121263.4A CN202111121263A CN113755497B CN 113755497 B CN113755497 B CN 113755497B CN 202111121263 A CN202111121263 A CN 202111121263A CN 113755497 B CN113755497 B CN 113755497B
Authority
CN
China
Prior art keywords
taro
gene
genes
reference gene
sequence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111121263.4A
Other languages
English (en)
Other versions
CN113755497A (zh
Inventor
何芳练
董伟清
蒋慧萍
黄诗宇
唐娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangxi Zhuang Nationality Autonomous Region Academy of Agricultural Sciences
Original Assignee
Guangxi Zhuang Nationality Autonomous Region Academy of Agricultural Sciences
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangxi Zhuang Nationality Autonomous Region Academy of Agricultural Sciences filed Critical Guangxi Zhuang Nationality Autonomous Region Academy of Agricultural Sciences
Priority to CN202111121263.4A priority Critical patent/CN113755497B/zh
Publication of CN113755497A publication Critical patent/CN113755497A/zh
Application granted granted Critical
Publication of CN113755497B publication Critical patent/CN113755497B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869—Methods for sequencing
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/13—Plant traits
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/166—Oligonucleotides used as internal standards, controls or normalisation probes
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y02A50/30—Against vector-borne diseases, e.g. mosquito-borne, fly-borne, tick-borne or waterborne diseases whose impact is exacerbated by climate change

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Biotechnology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Botany (AREA)
  • Mycology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明涉及植物基因技术领域，特别涉及芋球茎发育过程中内参基因筛选及其应用，本发明对芋球茎发育期的样本进行RNA‑Seq测序，筛选表达量中等、稳定的基因作为候选内参基因，利用软件评估并分析各个候选内参基因在球茎的整个发育阶段的表达稳定性，并利用目的基因验证，最终筛选出比常用的管家基因更为稳定的内参基因gene‑Taro_030192，最大程度地减小了由于内参基因的选择造成样品之间和样品内部的差异，为后期利用q RT‑PCR探索芋球茎发育相关关键基因奠定了坚实的基础。

Description

芋球茎发育过程中内参基因筛选及其应用
【技术领域】
本发明涉及植物基因技术领域，特别涉及芋球茎发育过程中内参基因筛选及其应用。
【背景技术】
芋(Colocasia esculenta)又称芋头、芋艿、毛芋，是天南星科芋属的多年生草本植物。芋最初产于东南亚，之后向东传播到远东亚洲和大洋洲，向西传播到地中海地区，目前在全球范围内广泛种植，是很多非洲和亚洲太平洋岛国居民的主要粮食，与其它块茎和块根植物相比，芋具有更高的营养价值。芋球茎的碳水化合物含量几乎是马铃薯的两倍，蛋白质含量也比山药、木薯、红薯高11％，由于存在如此丰富的营养物质含量，了解这些营养物质组成的遗传基础很重要。
但是目前对该作物的研究不多，近年来，随着高通量测序和分子生物学研究的不断发展，基因的表达分析也逐渐应用于揭示芋基因表达和调控机理的研究。随着芋基因组的测序，更加加快了芋相关性状的遗传机理解析。
实时荧光定量PCR(quantitative Real-time PCR,qRT-PCR)是目前检测基因表达水平的主要分析手段。qRT-PCR结果的可靠性受到初始样本量、RNA完整性、cDNA质量和扩增效率等因素的影响。因此，表达稳定的内参基因对比较不同样本目的基因表达水平至关重要，正确选择合适的内参基因进行校正和标准化，才能获得准确的qRT-PCR结果。目前，植物中常用来作为内参基因的基因有肌动蛋白(actin)、微管蛋白基因(β-tubulin,TUB)、转录延伸因子基因(elongation factor,EF1A和EF1B)、真核起始因子(eukaryotic initiationfactor,eIF)、泛素结合酶基因(ubiquitin-conjugating enzyme,UBC)、组蛋白(Histone,H3-1)、甘油醛-3-磷酸脱氢酶(glyceraldehyde-3-phosphate dehydrogenase,GAPDH)等。但是随着研究的深入，以往认为的很多内参基因(如actin等)在一定的实验条件下其表达量并不稳定。因此，在实际研究中，需要针对样本的不同组织或者不同生长发育时间段进行取材来筛选出稳定表达的内参基因。
本申请基于课题组前期对芋球茎不同发育阶段的样本进行转录本测序的数据，利用qRT-PCR分析候选基因在芋球茎8组不同发育时期表达稳定性，利用三种内参基因稳定性分析软件geNorm、NormFinder和BestKeeper分析候选内参基因在不同发育时期的稳定性并进行综合分析排序，为今后能更好开展芋球茎发育相关功能基因表达研究提供了参考依据。
【发明内容】
鉴于上述内容，有必要挖掘芋的内参基因，可以很好的对芋的表达分析提供重要参考。
为达到上述目的，本发明所采用的技术方案是：
芋内参基因gene-Taro_030192，所述内参基因gene-Taro_030192的核苷酸序列如序列表SEQ ID NO.1所示。
进一步的，所述引物对其上游序列如序列表SEQ ID NO.2所示，下游序列如序列表SEQ ID NO.3所示。
本发明还包括所述的芋内参基因gene-Taro_030192在芋球茎发育不同时期进行荧光定量表达分析的应用。
进一步的，所述芋球茎发育不同时期分别为：播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎。
进一步的，所述荧光定量表达分析优选的引物对其上游序列如序列表SEQ IDNO.2所示，下游序列如序列表SEQ ID NO.3所示。
进一步的，所述荧光定量表达分析的PCR反应体系共20μL为：2X SYBR GreenMaster Mix 10μl、10μM上游引物0.4μl、10μM上游引物0.4μl、RNase-free water 8.2μL、cDNA1μl；PCR反应程序为95℃预变性3min；95℃变性10s，57℃退火15min，72℃延伸20s总共45个循环；采集荧光信号。
本发明的芋内参基因gene-Taro_030192，筛选方法为：1、对芋球茎整个发育期连续取样，每隔30d进行规律性取样，即分别取播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎，共8组样本，标记为YTS1、YTS2、YTS3、YTS4、YTS5、YTS6、YTS7、YTS8为芋球茎材料，提取芋球茎RNA；2、反转录合成相应的芋球茎cDNA；3、根据不同基因设计不同引物并在荧光定量PCR仪上进行反应，绘制溶解曲线，进行稳定性分析，获得稳定性好、表达效果好的内参基因。
本发明具有如下有益效果：
本发明对芋球茎整个发育期连续取样，每隔30d进行规律性取样，即分别取播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎，共8组样本，标记为YTS1、YTS2、YTS3、YTS4、YTS5、YTS6、YTS7、YTS8为芋球茎材料，根据RNA-Seq测序数据分析，利用软件评估并分析各个候选内参基因在球茎的整个发育阶段中的表达稳定性，并利用目的基因验证，最终筛选出比常用的管家基因更为稳定的内参基因gene-Taro_030192，最大程度地减小了由于内参基因的选择造成样品之间和样品内部的差异，为后期利用q RT-PCR探索芋球茎发育相关关键基因奠定了坚实的基础。
【附图说明】
图1是球茎样品RNA琼脂糖凝胶电泳图，图中泳道1为YTS1、2为YTS2、3为YTS3、4为YTS4、5为YTS5、6为YTS6、7为YTS7、8为YTS8；
图2为候选内参基因gene-Taro_030192在球茎发育的PCR检测电泳图；
图3为候选内参基因gene-Taro_005478在球茎发育的PCR检测电泳图；
图4为候选内参基因gene-Taro_012907在球茎发育的PCR检测电泳图；
图5为候选内参基因gene-Taro_016279在球茎发育的PCR检测电泳图；
图6为候选内参基因gene-Taro_031479在球茎发育的PCR检测电泳图；
图2-6中，M为Marker DL2000；泳道1-3为芋头YTS1的3个平行样；泳道4-6为芋头YTS2的3个平行样；泳道7-9为芋头YTS3的3个平行样；泳道10-12为芋头YTS4的3个平行样；泳道13-15为芋头YTS5的3个平行样；泳道16-18为芋头YTS6的3个平行样；泳道19-21为芋头YTS7的3个平行样；泳道22-24为芋头YTS8的3个平行样。
图7为候选内参基因gene-Taro_030192基因的qPCR的扩增曲线图；
图8为候选内参基因gene-Taro_005478基因的qPCR的扩增曲线图；
图9为候选内参基因gene-Taro_012907基因的qPCR的扩增曲线图；
图10为候选内参基因gene-Taro_016279基因的qPCR的扩增曲线图；
图11为候选内参基因gene-Taro_031479基因的qPCR的扩增曲线图；
图12为候选内参基因gene-Taro_030192基因的qPCR的溶解曲线图；
图13为候选内参基因gene-Taro_005478基因的qPCR的溶解曲线图；
图14为候选内参基因gene-Taro_012907基因的qPCR的溶解曲线图；
图15为候选内参基因gene-Taro_016279基因的qPCR的溶解曲线图；
图16为候选内参基因gene-Taro_031479基因的qPCR的溶解曲线图；
图17为五个候选内参基因在芋球茎不同发育时期的Ct值箱线图；
图18为geNorm软件对候选内参基因的表达稳定排序；
图19为geNorm软件分析内参基因变异系数柱状图；
图20为gene-Taro_030192基因在芋头球茎中做为内参基因的表达情况图。
【具体实施方式】
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施的限制。
实施例1：
本实施例为内参基因gene-Taro_030192的获得，其通过如下方法得到：
一、本实施例的试验材料：
本实施例的芋材料为广西农业科学院生物技术研究所选育的芋新品种桂芋2号。2020年3月份种植于广西农业科学院里建科研基地。球茎整个发育期连续取样，每隔30d进行规律性取样。分别取播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎，共8组样本，标记为YTS1、YTS2、YTS3、YTS4、YTS5、YTS6、YTS7、YTS8。每组样本三个生物学重复。将采集的球茎削皮，切成2cm左右小块，放入50ml的无菌冷冻管中，液氮速冻后放入-80℃冰箱中备用。
二、总RNA提取、检测与c DNA合成：
利用美基生物难提植物总RNA小提试剂盒(R4165)提取球茎总RNA，操作步骤参照试剂盒说明书操作。使用1.5％琼脂糖凝胶电泳检测提取RNA的完整性，同时利用Nanodrop2100检测RNA的浓度和质量。使用诺唯赞公司的HiScript II Q RT SuperMix forqPCR(+gDNA wiper)(R223-01)的反转录试剂盒，将1μg的总RNA反转录成cDNA。总反应体系8μL，4×gDNA wiper Mix 2μL，总RNA1μg，加RNase-free water纯水补足至8μL。轻柔混匀后进行反转录反应，反应条件设置为：50℃、15min，85℃、5s，随即将cDNA 储存于-20℃冰箱备用。
三、内参基因的选择与特异性引物的设计：
根据芋球茎不同发育阶段的转录组测序数据，挑选表达量值相对稳定，使用变异系数(CV)进行评估，选择CV值小于0.2，且cDNA≧1000bp的基因作为候选内参基因。从候选内参基因中，根据芋基因的表达量值，选取了5个FPKM在25～35的基因，分别为gene-Taro_005478，gene-Taro_012907，gene-Taro_016279，gene-Taro_030192，gene-Taro_031479(表1)。利用Primer Premier 6软件设计5个候选内参基因的特异性引物，引物长度在146～185bp，引物由北京擎科生物技术有限公司(南宁)进行合成(表2)。
表1候选内参基因在球茎发育过程中的FPKM值
S1 S2 S3 S4 S5 S6 平均值
gene-Taro_016279 25.923763 27.793249 24.374704 23.942102 25.453297 28.683028 26.028357
gene-Taro_031479 27.445554 29.943567 26.516506 27.26381 24.116269 27.436972 27.120446
gene-Taro_012907 28.687375 29.151117 33.726094 26.461276 26.208292 25.275093 28.251541
gene-Taro_005478 34.746262 33.329018 32.09345 25.122231 30.080798 35.480905 31.808777
gene-Taro_030192 35.327808 38.398286 30.189635 28.885954 27.971349 30.10928 31.813719
表2芋5个候选内参基因的引物
Figure BDA0003277269230000051
四、内参基因的qPCR：
将cDNA样品稀释10倍作为模板上机检测。定量PCR实验，采用的仪器为：AceQ_qPCR_SYBR Green_Master_Mix定量PCR仪：德国AnalytikJena qTOWERE2.2。
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、cDNA1μl；
PCR反应程序为95℃预变性3min；95℃变性10s，57℃退火15min，72℃延伸20s总共45个循环；采集荧光信号，每组样品各设三次生物重复和三次技术重复。
上述PCR反应过程中各基因所用到的引物如上表2所示：
五、数据处理与分析
利用geNorm、NormFinder和BestKeeper三个算法不相同的软件单独分析候选内参基因的稳定性高低。geNorm和NormFinder是利用2-△CT值表示其内参基因的表达稳定性。而BestKeeper根据各候选内参基因的Ct值评价其稳定性。最后通过几何平均值法进行综合排序。
六、结果与分析
总RNA样品与引物特异性检测
提取球茎发育8个阶段的RNA后对其质量进行评估。如电泳结果所示(图1)，显示有两条完整条带，且28S亮度高于18S，说明RNA完整性较好，无降解。NanoDrop2100测得RNA浓度为161.5～339.0ng/μl，OD260/280均在1.91～2.11之间，表示RNA纯度较高。结果如图1所示，图1为各个样品的RNA凝胶电泳图；图中从左到右的泳道依次为1为YTS1、2为YTS2、3为YTS3、4为YTS4、5为YTS5、6为YTS6、7为YTS7、8为YTS8。从图可见，上述8个条带均表现出图谱清晰，条带亮度高，总RNA完整性好，符合试验要求。
候选内参基因的特异性分析
内参基因的PCR扩增产物电泳检测结果(如图2)所示，其余4个对比基因的PCR扩增产物电泳检测结果(如图3-6)显示；图2-6中M为Marker DL2000；泳道1-3为芋头YTS1的3个平行样；泳道4-6为芋头YTS2的3个平行样；泳道7-9为芋头YTS3的3个平行样；泳道10-12为芋头YTS4的3个平行样；泳道13-15为芋头YTS5的3个平行样；泳道16-18为芋头YTS6的3个平行样；泳道19-21为芋头YTS7的3个平行样；泳道22-24为芋头YTS8的3个平行样。从图中可见，5个候选内参基因均可见与预期产物长度大小相一致的单一条带，不存在引物二聚体。
对qRT-PCR结果中5个候选基因的扩增曲线进行分析，结果如图7-11所示；图7为gene-Taro_030192的扩增曲线、图8为gene-Taro_005478的扩增曲线、图9为gene-Taro_012907的扩增曲线、图10为gene-Taro_016279的扩增曲线、图11为gene-Taro_031479的扩增曲线；由图可知，5个候选内参基因的扩增曲线均很好。
对5个候选基因的溶解曲线进行分析，结果如图12-16所示，图12为gene-Taro_030192基因的溶解曲线、图13为gene-Taro_005478基因的溶解曲线、图14为gene-Taro_012907基因的溶解曲线、图15为gene-Taro_016279基因的溶解曲线、图16为gene-Taro_031479基因的溶解曲线；
从图12-图16的溶解曲线分析发现，溶解曲线均只有一个单峰，表明扩增条带单一，特异性强，没有非特异性扩增出现。因此，所设计的实时荧光定量PCR引物特异性强，扩增效率高，能够用于芋球茎发育的荧光定量PCR的内参引物试验。
候选基因的Ct值分析
Ct值是表达基因丰富度的重要衡量标准。Ct值与表达丰富程度呈反比关系。对5个内参基因在芋球茎不同发育时期进行表达丰度分析发现(结果如图17所示)，5个候选内参基因的Ct值全部介于26.38～32.68之间，表明这些基因的整体表达量比较接近。其中，gene-Taro_005478在球茎不同发育阶段的Ct值为27.18～28.59，最大值和最小值相差1.46个循环；gene-Taro_012907在球茎不同发育阶段的Ct值为27.86～29.62，最大值和最小值相差1.76个循环；gene-Taro_016279在球茎不同发育阶段的Ct值为26.38～27.42，最大值和最小值相差1.03个循环；gene-Taro_030192在球茎不同发育阶段的Ct值为27.73～29.35，最大值和最小值相差1.62个循环；gene-Taro_031479在球茎不同发育阶段的Ct值为30.12～32.68，最大值和最小值相差2.55个循环。
内参基因的稳定性分析
本研究用内参分析软件geNorm、Normfinder及BestKeeper软件作为辅助，对上述5个候选内参基因从多方面一致评估分析。
geNorm软件对候选内参基因的稳定性评估
geNorm程序通过计算出每个内参基因稳定性的M值来筛选出稳定性较好的内参基因，判定标准为M值越小内参基因稳定性越好，反之，则稳定性越差，若M值大于1.5，则不考虑为内参基因。按照表达稳定性从高到低，依次为gene-Taro_030192(0.386)﹥gene-Taro_012907(0.390)﹥gene-Taro_005478(0.455)﹥gene-Taro_016279(0.512)﹥gene-Taro_031479(0.621)。所有候选内参基因M值都小于1.5，均有可能是适合的内参基因。其中以gene-Taro_030192和gene-Taro_012907的M值最小，稳定性也越好(图18)。
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.112)﹤0.15(图19)，内参基因组合最佳数目是2个。因此，芋球茎不同发育时期最合适的内参基因有gene-Taro_030192和gene-Taro_012907。
NormFinder软件对候选内参基因的稳定性评估
NormFinder软件程序算法与geNorm程序类似，也是先获得内参基因表达稳定值，再根据稳定值大小来筛选最合适的内参基因，判定标准为表达稳定值最小的候选内参基因为合适内参基因。NormFinder程序不但可以比较候选内参基因的表达差异，还可以汁算样品组间的变异，但该程序只能筛选出一个最合适的内参基因。
NormFinder软件对5个候选内参基因的评估结果显示，内参基因表达稳定值gene-Taro_030192(0.045)﹥gene-Taro_012907(0.057)﹥gene-Taro_005478(0.225)﹥gene-Taro_016279(0.289)﹥gene-Taro_031479(0.399)，结果与geNorm软件评估结果一致。在芋球茎不同发育过程最合适内参基因为gene-Taro_030192。
BestKeeper对候选内参基因的分析结果
BestKeeper软件通过计算样品数据，获得基因间相关配对系数(r)，标准偏差(SD)和变异相关系数(CV)，利用这3项参数判定内参的稳定性。当相关系数r越高、变异系数(CV)和标准偏差(SD)越低，则内参基因就越稳定。另外，当标准偏差SD大于1时，认为该候选内参基因表达不稳定，具体如表3。
根据SD值排序，gene-Taro_016279﹥gene-Taro_005478﹥gene-Taro_030192﹥gene-Taro_012907﹥gene-Taro_031479，但所有候选的5个内参基因SD值均小于1，最合适的为gene-Taro_016279，gene-Taro_005478，gene-Taro_030192。
表3 BestKeeper软件分析内参基因稳定值
Figure BDA0003277269230000081
Figure BDA0003277269230000091
内参基因的稳定性分析
通过对三种内参基因软件进行比较和分析，发现BestKeeper软件的分析结果与geNorm、NormFinder的分析结果有差异。由于geNorm和NormFinder统计算法是以Ct值为原始数据，求2-ΔCt值，得出的结果一致。在Excel中对3款软件的分析结果进行综合性评价；几何平均值越小，内参基因则越稳定。从表4可以看出gene-Taro_030192在三款软件的几何平均值为1.44，在5个候选基因中几何平均值最小，为最稳定的内参基因。因此，5个候选基因中，gene-Taro_030192为最合适的内参基因。
表4 geNorm、NormFinder和BestKeeper三款软件综合分析结果
geNorm排序 NomFinder排序 BestKeeper排序 几何平均值
gene-Taro_005478 3 3 2 2.62
gene-Taro_012907 2 2 4 2.52
gene-Taro_016279 4 4 1 2.52
gene-Taro_030192 1 1 3 1.44
gene-Taro_031479 5 5 5 5.00
实施例2：
根据实施例1筛选出来的稳定性最佳的基因gene-Taro_030192进行进一步的验证，保证得出真实可靠的定量数据。以gene-Taro_030192作为内参基因，选择球茎发育过程中淀粉合成关键酶基因CeAGPS1基因作为目的基因，分析目的基因在芋球茎发育的8个时期表达，具体如下：
分别提取不同时期：即芋整个发育期，即播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎(标记为YTS1、YTS2、YTS3、YTS4、YTS5、YTS6、YTS7、YTS8)时期的cDNA为作为定量PCR的模板，以gene-Taro_030192基因为内参基因，以ADP-葡萄糖焦磷酸化酶小亚基基因CeAGPS1为目的基因设计引物，引物见表5：
表5内参基因和目的基因的引物序列
Figure BDA0003277269230000092
在实时荧光定量PCR检测仪上进行荧光定量反应，反应的
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、cDNA1μl；
PCR反应程序为95℃预变性3min；95℃变性10s，57℃退火15min，72℃延伸20s总共45个循环；采集荧光信号。
反应结束后使用2-ΔΔCT法所有的数据均表示为平均值±标准偏差。实时荧光定量PCR数据分析用SPSS 17.0和Microsoft Office Excel 2007进行统计分析，采用SPSS17.0软件方差分析，采用Excel 2010软件制图。
得到的反应结果如图20所示，表明：内参基因gene-Taro_030192在芋整个发育期，即播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎均能稳定的表达，能很好反映CeAGPS1在球茎发育过程的表达量变化。
综上所述，基因gene-Taro_030192可在芋球茎发育过程中进行稳定表达，能做为芋球茎发育不同时期荧光定量的内参基因，能为今后对芋重要功能性状基因的研究奠定基础。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明的保护范围应以所附权利要求为准。
序列表
<110> 广西壮族自治区农业科学院
<120> 芋球茎发育过程中内参基因筛选及其应用
<160> 3
<170> SIPOSequenceListing 1.0
<210> 1
<211> 1540
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 1
ctcttctcca gcgagccaga aagggggtgg gaagaggaga gagataactt tttagagaga 60
gaaagggagc gagaacaaca gtggaaacaa gagagataga gagggggcgg aatcgaagag 120
acgcagacgc gggcgacccc ccttctcttc tctctctctc tctctctctc cccagcgacg 180
cgatctcctt cgcgccgctc ccctcaccat ctctgtcccc cggcgaccgc ttgccgcaaa 240
gcttcccaga tgtacttcct ctctgctgag tgatctccgc tgcttgcttt cccagaatga 300
gggtttggga gcagtcgtga ggcaccgtcc ggttcagcgt gggttgcgga ggatgtctgc 360
caagccgcaa gtcaccatca ccctcgggcg tagcggccaa gttgtcaaga gagctggggg 420
catttcagaa ggtgctcatt ctgattatgg gtccttggca ggaggcaaga gacctatcaa 480
ggatagattg aggaataact cagatggtta tggtggccag tcccagaatc ctcgaattgg 540
tcagaatgac ctcaggttca aacttatgcg gaaaaacttg cccagaccaa tgaatatgga 600
ctctaagatg catgatgtgg accttcgcga taggctttca cgaactgccc aagtgccttt 660
gagatttgag gtacgcaaac atcaaatgga agcaaagacc tcggggtttc ctggtcggat 720
tccttctaca agaagtgcag atgatttgct gcgtatggat tctcgaagaa gatcatattc 780
tcctgggcct tccgatagat tcagacacag atccccagac aggcttttaa gccattctag 840
gggactgtcc ccaccaagaa gttatgatga cttgcgacag cattcgtcga taaggtctgt 900
cgatcctcca aggtcgtcat cctatatgac caaaactatg gttgatgctc ccaggccctc 960
cttcacagca aaggctcctt atcctgttga accttctaaa caagctatgc ggcttcctac 1020
agccaatggt acactgaaga gctcatatat ggctgaggag ccaaccagtg ttagtaattt 1080
gctgcattct ctgggtttgg gaaaatatgc cattcttttc caagcagaag aagtggacat 1140
ggctgcacta aagcaaatgg gagacaatga tctgaaggag ttggggatac caatgagcag 1200
tgccagagcc cgaggatgtt acggtgtcac gagcactggg gttgtggata gaggggaggg 1260
accagggaac atgggacggc agatgaggaa gacacaaagg tggaggcagt ggtcagatct 1320
ggcccggcat tatatgacct gggtgctgat ctccttttcg gcctctgtgc ctcttccctg 1380
gtcaaatccc catcatagat ctcgtcagat ctactacaga ctgccattgt tcaccaatat 1440
ctctctggcc ttgcagatat gtcctccaga tgctgcagta gatactgtca aggccggctg 1500
cttccatgaa ggagacagtg cacgccacct gactagctga 1540
<210> 2
<211> 22
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 2
cctggtcgga ttccttctac aa 22
<210> 3
<211> 22
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 3
gctgtcgcaa gtcatcataa ct 22

Claims (6)

1.芋内参基因gene-Taro_030192，其特征在于，所述内参基因gene-Taro_030192的核苷酸序列如序列表SEQ ID NO.1所示。
2.检测如权利要求1所述芋内参基因gene-Taro_030192的引物对，其特征在于，所述引物对其上游序列如序列表SEQ ID NO.2所示，下游序列如序列表SEQ ID NO.3所示。
3.如权利要求1所述的芋内参基因gene-Taro_030192在芋球茎发育不同时期进行荧光定量表达分析的应用。
4.根据权利要求3所述的应用，其特征在于，所述芋球茎发育不同时期分别为：播种后30d、60d、90d、120d、150d、180d、210d、240d的球茎。
5.根据权利要求3所述的应用，其特征在于，所述荧光定量表达分析采用的引物对其上游序列如序列表SEQ ID NO.2所示，下游序列如序列表SEQ ID NO.3所示。
6.根据权利要求3所述的应用，其特征在于，所述荧光定量表达分析的PCR反应体系共20μL为：2×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℃变性10s，57℃退火15min，72℃延伸20s总共45个循环；采集荧光信号。
CN202111121263.4A 2021-09-24 2021-09-24 芋球茎发育过程中内参基因筛选及其应用 Active CN113755497B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111121263.4A CN113755497B (zh) 2021-09-24 2021-09-24 芋球茎发育过程中内参基因筛选及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111121263.4A CN113755497B (zh) 2021-09-24 2021-09-24 芋球茎发育过程中内参基因筛选及其应用

Publications (2)

Publication Number Publication Date
CN113755497A CN113755497A (zh) 2021-12-07
CN113755497B true CN113755497B (zh) 2023-03-14

Family

ID=78797297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111121263.4A Active CN113755497B (zh) 2021-09-24 2021-09-24 芋球茎发育过程中内参基因筛选及其应用

Country Status (1)

Country Link
CN (1) CN113755497B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114438243B (zh) * 2022-01-27 2023-06-27 广西壮族自治区农业科学院 芋球茎发育内参基因CeRPL28的筛选及其应用

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110982923B (zh) * 2019-12-25 2020-08-21 广西壮族自治区农业科学院 西番莲内参基因及其筛选方法和应用
CN110982925B (zh) * 2019-12-27 2022-09-13 广东省农业科学院环境园艺研究所 白掌佛焰苞不同发育时期荧光定量内参基因及应用
CN112458099B (zh) * 2020-12-09 2021-06-25 广西壮族自治区农业科学院 香葱内参基因及其筛选方法和应用
CN113046367B (zh) * 2021-04-21 2022-04-19 福建农林大学 慈姑球茎发育过程中内参基因筛选及其应用
CN113174445B (zh) * 2021-04-21 2023-03-21 广西壮族自治区农业科学院 分析慈姑不同组织的内参基因及其筛选方法和应用

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GenBank.Colocasia esculenta cultivar Longxiangyu ecotype Jiangsu isolate LXY-2018 chromosome 3, whole genome shotgun sequence.2020,Accession: CM024755.1. *

Also Published As

Publication number Publication date
CN113755497A (zh) 2021-12-07

Similar Documents

Publication Publication Date Title
Chen et al. A survey of quantitative real-time polymerase chain reaction internal reference genes for expression studies in Brassica napus
CN112575010B (zh) 山药不同组织荧光定量的内参基因及其引物与应用
CN105132417A (zh) 茶树低温胁迫下miRNA荧光定量PCR内参基因及其筛选方法和应用
CN112877465B (zh) 柳杉不同组织的荧光定量内参基因及其专用引物和应用
CN111676230A (zh) 一种丝瓜内参基因EF-1α及其引物和应用
CN113755497B (zh) 芋球茎发育过程中内参基因筛选及其应用
Salcedo et al. Comparative expression analysis of reference genes in field-grown cassava
CN114196772B (zh) 南京椴不同组织及胁迫处理条件下实时定量pcr内参基因及其筛选方法和应用
CN109207634B (zh) 海州常山干旱胁迫下荧光定量内参基因及其专用引物和应用
CN113774065B (zh) 美国白蛾不同性别成虫荧光定量内参基因及其引物和应用
CN113846108B (zh) 芋高表达内参基因Ce047468的筛选及其应用
CN113174445A (zh) 分析慈姑不同组织的内参基因及其筛选方法和应用
CN108085409B (zh) 不同组织中杉木内参基因的筛选方法和筛选基因作为内参基因的应用
CN111575401A (zh) 一种丝瓜内参基因ubq的引物及应用
CN114457187B (zh) 梅花实时荧光定量pcr分析中内参基因的筛选方法和应用
CN106754965B (zh) 一种与杨树抗逆基因表达调控相关的内参基因及其应用
CN113151505B (zh) 一种美国白蛾不同种群密度的荧光定量内参基因及其引物和应用
CN114438243B (zh) 芋球茎发育内参基因CeRPL28的筛选及其应用
CN112481404B (zh) 盐碱栽培条件下猴樟内参基因、引物及筛选方法
CN111690762B (zh) 用于多花黄精块茎基因表达分析的内参基因组合及其应用
CN115927296A (zh) 菜心高温胁迫内参基因及其筛选方法和应用
CN112980989A (zh) 柳杉不同组织的miRNA荧光定量内参基因及其引物和应用
CN112813181A (zh) 适用于雪茄烟不同组织基因表达分析的内参基因及其应用
CN105274219A (zh) CL5547.Contig2基因在南瓜基因表达实时荧光定量PCR分析中作为内参基因的用途
CN116287366A (zh) 芋内参基因Ce049358的筛选及其应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant