高通量筛选玉米内参基因的方法
【专利摘要】本发明公开了一种高通量筛选玉米内参基因的方法，该方法利用源于具有高度多样性的不同玉米自交系、不同组织、不同发育时期以及不同胁迫处理后构建的高通量转录组测序(RNA-seq)数据库，大规模分析整个转录组的表达情况，筛选在不同发育时期、不同条件下表达最恒定的基因作为内参基因，为转录组表达研究提供准确的对照。本方法利用整个转录组数据，分析全基因组所有表达的基因，系统性地筛选玉米内参基因，相对传统筛选方法更加准确、可靠。
【专利说明】高通量筛选玉米内参基因的方法
【技术领域】
[0001]本发明涉及分子遗传学领域，尤其涉及一种高通量筛选玉米内参基因的方法。
【背景技术】
[0002]基因表达(gene express1n)是指细胞在生命过程中，把储存在脱氧核糖核酸(DNA)序列中的遗传信息经过转录和翻译，转变成具有生物活性的蛋白质分子。该过程决定了细胞的分化及形态的发生。每个基因转录产生信使核糖核酸(mRNA)的量，受到时空等多种因素的调控，最终影响个体的生长发育、形态结构及生物学功能。基因调控是现代分子生物学研究的中心课题之一。为研究基因表达调控的机制就必须了解不同基因在特定条件下的表达水平。在检测基因表达量时取样、上样等实验操作会对结果产生影响。为消除实验误差，就必须使用内源性参照基因作为参照物。在理论上，内参基因在各组织和细胞中的表达恒定，在检测基因的表达水平变化时作为参照物，用于校正实验误差，保证结果的准确性。
[0003]常用的内源性参照基因多为已知持家基因，这类基因是指所有细胞中均表达的一类基因，其产物对维持细胞基本生命活动是必需的，如微管蛋白基因、糖酵解酶系基因与核糖体蛋白基因等。早期人们认为持家基因在不同组织和不同发育时期的表达量也是一致的，但随着研究的深入，人们发现不同持家基因在不同条件下的表达量也有差异。因此，发掘表达恒定的内参基因对基因表达调控的研究至关重要。
[0004]目前，玉米内参基因的筛选主要是先选择几个或十几个候选持家基因，然后分析比较它们在不同条件下的表达变异，选择相对表达最稳定的一个作为内参基因。需要说明的是，这种方法筛选内参基因的覆盖面非常低，其是根据之前界定的部分持家基因，仅仅比较部分持家基因的表达稳定性，数据量太小，不足以筛选出真正表达稳定的内参基因。
【发明内容】
[0005]本发明要解决的技术问题是解决现有技术中筛选内参基因的方法覆盖面较低，数据量较小，不足以选择出真正表达稳定的内参基因的缺陷。
[0006]为了解决上述技术问题，本发明采用的技术方案如下:一种高通量筛选玉米内参基因的方法，选择源于具有高度多样性的不同玉米自交系、不同组织、不同发育时期以及不同胁迫处理后构建的高通量转录组测序(RNA-seq)数据库，其筛选步骤如下:
[0007](I)利用SolexaQA软件将RNA-seq的数据进行过滤，除去比较差的数据，留下质量较闻的序列；
[0008](2)利用Bowtie软件将得到的序列锚定到参考基因组上，当一个序列能够锚定到基因组上的不同位点，则除去该序列，只选择在基因组中有单一位点的序列；
[0009] (3)使用 Cuff links 软件计算各转录本的 FPKM 值(fragments per kilobase oftranscript per mill1n mapped reads),去掉FPKM小于10的转录本,余下的转录本用于后续分析；
[0010](4)利用NormFinder软件计算所有基因在不同条件下的表达稳定性系数，选取不同RNA-seq数据库中都能达到前20%稳定表达的基因，并整合分析不同品种、不同组织、不同发育时期以及不同胁迫处理后得到的所有RNA-seq数据库，筛选在不同RNA-seq数据库中都能稳定表达的基因作为候选内参基因；
[0011](5)选取不同品种进行胁迫处理，在不同发育时期取样提取RNA，根据候选内参基因的序列设计特异引物，通过qPCR检测候选内参基因及现有内参基因的表达水平，利用NormFinder计算各内参基因在不同条件下的表达稳定性系数，比较候选内参基因与现有内参基因在不同材料及不同处理条件下的表达稳定性，选择最优内参基因。
[0012]进一步的，所述步骤(1)中比较差的数据为Phred值低于20、长度小于20bp的数据。
[0013]本发明的高通量筛选玉米内参基因的方法，其有益效果包括:
[0014](I)使用高通量转录组测序(RNA-seq)的数据，大规模分析整个转录组的表达情况，筛选在不同发育时期、不同条件下表达最恒定的基因作为内参基因，系统性筛选地内参基因相对现有内参基因更加准确、可靠。
[0015](2)使用高通量转录组测序(RNA-seq)的数据，能够同时筛选整个玉米基因组所有的表达基因，大大提闻了选择的效率。
[0016](3)目前已知源于不同品种、不同组织、不同发育时期以及不同胁迫处理的RNA-seq数据库已经超过50个，属于公开资源，可直接在网站上下载使用，通过数据分析选择出候选内参基因后，再通过qPCR验证其表达稳定性，大大减少了实验工作量，降低了实验成本。
【专利附图】
【附图说明】
[0017]图1为本实施例高通量筛选玉米内参基因方法的流程图；
[0018]图2为利用本实施例高通量筛选玉米内参基因方法获得的内参基因(DPP9、DUF)与现有内参基因(ACT和GAPDH)的表达稳定性示意图。
【具体实施方式】
[0019]下述实施例中实验方法如无特殊说明，均为常规实验方法。下述实施例中所述的实验试剂及耗材如无特殊说明，均来自常规生化试剂公司。
[0020]下面结合附图与【具体实施方式】对本发明的高通量筛选玉米内参基因的方法作进一步的详细说明。
[0021]本实施例为高通量筛选玉米内参基因的方法，选择源于具有高度多样性的不同玉米自交系、不同组织、不同发育时期以及不同胁迫处理后构建的高通量转录组测序(RNA-seq)数据库，其筛选步骤如图1所示，包括:
[0022](I)利用SolexaQA软件将RNA-seq的数据进行过滤，除去Phred值低于20、长度小于20bp的质量比较差的数据，留下质量较高的序列；
[0023](2)利用Bowtie软件将测序得到的序列锚定到参考基因组上，当一个序列能够锚定到基因组上的不同位点，则除去该序列，只选择在基因组中有单一位点的序列；
[0024](3)使用 Cuff links 软件计算各转录本的 FPKM 值(fragments per kilobase oftranscript per mill1n mapped reads),去掉FPKM小于10的转录本,余下的转录本用于后续分析；
[0025](4)利用NormFinder软件计算所有基因在不同条件下的表达稳定性系数，选取不同数据库中都能达到前20%稳定表达的基因，整合分析不同品种、不同组织、不同发育时期以及不同胁迫处理后得到的所有RNA-seq数据库，筛选在不同数据库中都能达到前20%稳定表达的基因作为候选内参基因；
[0026](5)选取 9 个不同玉米品种(B73、Mol7、Z1144、Ky21、M37W、M162W、Mol4、Mo28 以及Mo51)，种子发芽24小时后移栽至盆钵中。从三叶期开始进行氮处理，每隔两天施尿素一次，对照组只浇水。发芽20天、40天后取叶片利用TRIzol方法提取RNA，用TransScriptgDNA去除及cDNA合成试剂盒反转录为cDNA。随机选择两个候选内参基因DPP9、DUF (序列号分别为:GRMZM2G174572，GRMZM2G163888)，根据候选内参基因及现有内参基因ACT和GAPDH(序列号分别为:GRMZM2G126010，GRMZM2G046804)的序列用软件Primer (3)设计引物，其引物如表1所示。
[0027]表1通过Primer3设计的内参基因的qPCR引物序列
[0028]
【权利要求】
1.高通量筛选玉米内参基因的方法，其特征在于:选择源于具有高度多样性的不同玉米自交系、不同组织、不同发育时期以及不同胁迫处理后构建的高通量转录组测序(RNA-seq)数据库，其筛选步骤如下: (1)利用SolexaQA软件将RNA-seq的数据进行过滤，除去比较差的数据，留下质量较高的序列； (2)利用Bowtie软件将得到的序列锚定到参考基因组上，当一个序列能够锚定到基因组上的不同位点，则除去该序列，只选择在基因组中有单一位点的序列； (3)使用Cuff links 软件计算各转录本的 FPKM 值(fragments per kilobase oftranscript per mill1n mapped re ads),去掉FPKM小于10的转录本,余下的转录本用于后续分析； (4)利用NormFinder软件计算所有基因在不同条件下的表达稳定性系数，选取不同数据库中都能达到前20%稳定表达的基因，并整合分析不同品种、不同组织、不同发育时期以及不同胁迫处理后得到的所有RNA-seq数据库，筛选在不同数据库中都能稳定表达的基因作为候选内参基因； (5)选取不同品种进行胁迫处理，在不同发育时期取样提取RNA，根据候选内参基因的序列设计特异引物，通过qPCR检测候选内参基因及现有内参基因的表达水平，利用NormFinder软件计算各内参基因在不同条件下的表达稳定性系数，比较候选内参基因与现有内参基因在不同材料及不同处理条件下的表达稳定性，选择最优内参基因。
2.按照权利要求1所述的高通量筛选玉米内参基因的方法，其特征在于:所述步骤(1)中比较差的数据为Phred值低于20、长度小于20bp的数据。
【文档编号】C12Q1/68GK104131107SQ201410396506
【公开日】2014年11月5日 申请日期:2014年8月12日 优先权日:2014年8月12日
【发明者】林峰, 赵涵 申请人:江苏省农业科学院