• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

巨龙竹秆形发育过程实时荧光定量PCR内参基因的筛选

郭晓娟 陈凌娜 杨汉奇

引用本文:
Citation:

巨龙竹秆形发育过程实时荧光定量PCR内参基因的筛选

    作者简介: 郭晓娟(1991-), 女, 硕士研究生, 研究方向:林木遗传育种.Email:lkygxj_417@163.com.
    通讯作者: 杨汉奇
  • 基金项目:

    云南省科学技术厅 2014HB041

    国家自然科学基金 31270662

    中央级公益性科研院所基本科研业务费专项 CAFYBB2017ZX001-8

  • 中图分类号: S718

Reference Gene Selection for Quantitative Real-Time PCR in Studying Culm Shape Development of Dendrocalamus sinicus

    Corresponding author: YANG Han-qi
  • CLC number: S718

  • 摘要: 目的 利用荧光定量PCR(qRT-PCR)和表达分析软件筛选适用于巨龙竹秆形发育研究的qRT-PCR的稳定内参基因。 方法 应用常规PCR和qRT-PCR, 分析EF-1αGAPDHActinTubulinTIP-41PP2A共6个内参基因在不同秆形巨龙竹("通直型"和"弯曲型")竹笋3个不同发育时期的表达情况。利用GeNorm、NormFinder软件对各内参基因的表达稳定性进行评估; 以筛选的内参基因对巨龙竹木质素合成中的PTAL基因的表达量进行定量分析, 验证所筛选内参基因的有效性。 结果 PCR和qRT-PCR的结果表明, 6个候选内参基因PCR目的片段电泳条带单一, 熔解曲线具有明显的单一峰。内参基因ActinEF-1αGAPDH的稳定性和相关性表现最佳; 以它们为内参, 对茎秆发育过程中PTAL基因进行定量分析, 结果显示它们均具有高效性。 结论 ActinEF-1αGAPDH可以作为巨龙竹秆形发育研究中目的基因定量分析的内参基因。为进一步分析巨龙竹秆形发育过程中关键基因表达量的变化提供研究基础。
  • 图 1  巨龙竹不同秆形各发育时期总RNA电泳图

    Figure 1.  Electrophoresis results of total RNA extracted from different types of D.sinicus at three developmental stages

    图 2  6个候选内参基因常规PCR的电泳结果

    Figure 2.  Electrophoresis results of six candidateerence genes

    图 3  6个候选内参基因的表达稳定性GeNorm分析

    Figure 3.  Expression stabilities of six candidate reference genes analyzed by GeNorm

    图 4  6个候选内参基因标准化因子的配对差异分析

    Figure 4.  Pairwise variance analyses of six candidate reference genes

    图 5  6个候选内参基因的表达稳定性NormFinder分析

    Figure 5.  Expression stabilities of six candidate reference genes analyzed by NormFinder

    图 6  不同内参基因对PTAL基因的定量分析

    Figure 6.  Quantitative analysis of different reference genes to PTAL

    表 1  采样具体信息

    Table 1.  Sampling information

    地点Site 茎秆类型Culm type 经度Longitude 纬度Latitude 海拔Altitude/m 采样时间Time
    西盟县中课镇窝笼村 通直型 99°32′06″ 22°43′48″ 1 254 2 016-07-15
    孟连县景信乡景信村 弯曲型 99°36′31″ 22°26′36″ 1 007 2 016-07-16
    下载: 导出CSV

    表 2  qRT-PCR检测中所用的6个候选内参基因的引物序列

    Table 2.  Primer sequences of six candidate reference genes for qRT-PCR analysis

    基因
    Gene name
    引物序列(5’→3’)
    Primer sequence (5’→3’)
    解链温度
    Tm/℃
    扩增产物长度Length of amplified product/bp
    PP2A F:CAGTCATCCCATGCTCTTGCT 60.6 158
    R:GAAGCTGGCTTACGTGAGATGAA 59.8
    EF-1α F:GCAACAAAACCACGCTTGAGA 59.6 186
    R:GGGTCGTGTTGAGACTGGTG 59.5
    GAPDH F:AGTCCACTGGTGTCTTCACTG 59.6 178
    R:CAATTTGTGGTGCAGCTAGCATT 60.0
    Actin F:TGCCTACATTGCCCTGGACTA 59.6 162
    R:TCCATTCCAATGAAGGATGGCT 59.8
    Tubulin F:GTCGATGCGCGAGAATACCT 59.4 186
    R:GCAGTTGCAACCATCAAGACC 59.6
    TIP-41 F:CAGCGAAATGGAAGTTCAGAAG 59.8 169
    R:TGGTCATCAGTATCCTCCCAAC 59.7
    注:正向引物(F)forwards primer,反向引物(R)reverse primer,解链温度(Tm)melting temperature。
    下载: 导出CSV
  • [1] 耿伯介, 王正平.中国植物志[M].北京:科学出版社, 1996.

    [2] 谷志佳, 杨汉奇, 孙茂盛, 等.巨龙竹资源分布特点及其开花结实现象[J].林业科学研究, 2012, 25(1):1-5. doi: 10.3969/j.issn.1001-1498.2012.01.001

    [3] 辉朝茂, 张国学, 李在留, 等.珍稀竹种巨龙竹种群特性及其保护研究[J].竹子研究会刊, 2004, 23(4):4-9.

    [4] 杨汉奇, 孙茂盛, 阮桢媛, 等.云南4种典型热带丛生竹种源分化[J].林业科学研究, 2014, 27(2):172-176.

    [5]

    Bustin S A, Beaulieu J F, Huggett J, et al. MIQE précis:Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments[J].BMC Molecular Biology, 2010, 11(1):1-5. doi: 10.1186/1471-2199-11-1
    [6]

    Chi X, Hu R, Yang Q, et al. Validation of reference genes for gene expression studies in peanut by quantitative real-time RT-PCR[J].Molecular Genetics and Genomics, 2012, 287(2):167-176. doi: 10.1007/s00438-011-0665-5
    [7] 孙美莲, 王云生, 杨冬青, 等.茶树实时荧光定量PCR分析中内参基因的选择[J].植物学报, 2010, 45(5):579-587. doi: 10.3969/j.issn.1674-3466.2010.05.007

    [8] 袁伟, 万红建, 杨悦俭.植物实时荧光定量PCR内参基因的特点及选择[J].植物学报, 2012, 47(4):427-436.

    [9] 齐飞燕, 胡陶, 彭镇华, 等.毛竹实时荧光定量PCR内参基因的筛选及成花基因PheTFL1表达分析[J].西北植物学报, 2013, 33(1):48-52. doi: 10.3969/j.issn.1000-4025.2013.01.010

    [10]

    Fan C, Ma J, Guo Q, et al. Selection of Reference Genes for Quantitative Real-Time PCR in Bamboo (Phyllostachys edulis)[J].Plos One, 2013, 8(2):e56573. doi: 10.1371/journal.pone.0056573
    [11]

    Cui K, Wang H, Liao S, et al. Transcriptome Sequencing and Analysis for Culm Elongation of the World's Largest Bamboo (Dendrocalamussinicus)[J].Plos One, 2016, 11(6):e0157362. doi: 10.1371/journal.pone.0157362
    [12] 魏洁书, 杨锦芬.应用荧光定量比较Ct值法测定基因相对表达量[J].中国科技论文在线精品论文, 2013, 6(5):390-395.

    [13]

    Fraser C M, Chapple C.The Phenylpropanoid Pathway in Arabidopsis[J].American Society of Plant Biologists, 2011, 9(1):e0152.
    [14] 吴建阳, 何冰, 杜玉洁, 等.利用geNorm、NormFinder和BestKeeper软件进行内参基因稳定性分析的方法[J].现代农业科技, 2017, (5):278-281. doi: 10.3969/j.issn.1007-5739.2017.05.174

    [15]

    Andersen C L, Jensen J L, rntoft T F.Normalization of real-time quantitative reverse transcription-PCR data:a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets[J].Cancer Research, 2004, 64(15):5245-5250. doi: 10.1158/0008-5472.CAN-04-0496
    [16]

    Zeng S, Liu Y, Wu M, et al. Identification and validation of reference genes for quantitative real-time PCR normalization and its applications in Lycium[J].Plos One, 2014, 9(5):e97039. doi: 10.1371/journal.pone.0097039
    [17] 刘洪峰, 高乐旋, 胡永红.牡丹不同发育阶段种子和花瓣组织实时荧光定量PCR中内参基因的挖掘与筛选[J].农业生物技术学报, 2015, 23(12):1639-1648.

    [18]

    Jaroová J, Kundu J K.Validation of reference genes as internal control for studying viral infections in cereals by quantitative real-time RT-PCR[J].BMC Plant Biology, 2010, 10(1):1-9. doi: 10.1186/1471-2229-10-1
    [19] 覃迎姿, 黄先益, 叶兴枝, 等.植物延伸因子eEF1A研究进展[J].南方农业学报, 2009, 40(5):472-477. doi: 10.3969/j.issn.2095-1191.2009.05.005

    [20] 李洁, 谢文光, 沈倍奋, 等.GAPDH功能多样性研究进展[J].军事医学, 2006, 30(5):483-486. doi: 10.3969/j.issn.1674-9960.2006.05.022

    [21] 郭景康, 陈青云, 戢茜, 等.拟南芥、水稻和杨树ACTIN家族全基因组分析[J].上海大学学报:自然科学版, 2009, 18(4):426-431.

    [22] 董辉, 李越中, 胡玮.γ-微管蛋白研究进展[J].生物化学与生物物理进展, 2002, 29(5):686-690. doi: 10.3321/j.issn:1000-3282.2002.05.006

    [23]

    Chao Z, Yan W, Fu M, et al. Molecular cloning, expression and functional analysis of three subunits of protein phosphatase 2A (PP2A) from black tiger shrimps (Penaeus monodon)[J].Comparative Biochemistry & Physiology Part B Biochemistry & Molecular Biology, 2017, 204(1):77-89.
    [24] 邢朝斌, 刘岩, 周秘, 等.刺五加液泡膜内在蛋白基因的克隆与表达分析[J].中草药, 2014, 45(2):250-254.

  • [1] 杨英英赵林姣杨桂娟张玉付鹏跃胡继文刘莹王楠 . ‘麦缘锦楸’叶色表型qRT-PCR内参基因筛选及验证. 林业科学研究, 2022, 35(1): 123-131. doi: 10.13275/j.cnki.lykxyj.2022.01.014
    [2] 华雅洁岳远征杨秀莲何卿 . 海州常山叶片实时荧光定量PCR的内参基因选择. 林业科学研究, 2022, 35(2): 194-202. doi: 10.13275/j.cnki.lykxyj.2022.02.023
    [3] 赵学彩郑唐春臧丽娜曲冠证 . 杨树类锌指基因ZFL的功能分析. 林业科学研究, 2013, 26(5): 562-570.
    [4] 张妍刘瀛孙丰宾戴超刘雪梅 . 白桦APETALA2(AP2)转录因子基因的分离及其表达. 林业科学研究, 2012, 25(2): 254-260.
    [5] 李红英陈萌笛王政博郝自远刘龙昌赵西平倪建伟 . 镉胁迫下构树qRT-PCR内参基因筛选及验证. 林业科学研究, 2023, 36(4): 129-138. doi: 10.12403/j.1001-1498.20220500
    [6] 李麟坤梁重钧宁楚龙张薇姜国维张鑫王利兵 . 文冠果种子发育基因表达分析内参基因的选择. 林业科学研究, 2024, 37(2): 1-7. doi: 10.12403/j.1001-1498.20230291
    [7] 周成城荣俊冬谢德金杨德明何天友郑郁善 . 福建柏实时荧光定量PCR内参基因的选择. 林业科学研究, 2021, 34(1): 137-145. doi: 10.13275/j.cnki.lykxyj.2021.01.017
    [8] 张颖陈婉婷陈冉红帅鹏李明 . 杉木实时荧光定量PCR分析中内参基因的选择. 林业科学研究, 2019, 32(2): 65-72. doi: 10.13275/j.cnki.lykxyj.2019.02.010
    [9] 吴涛韩素英张俊红李万峰杨文华齐力旺 . 日本落叶松种子萌发过程中18个miRNAs的表达变化. 林业科学研究, 2013, 26(S1): 9-17.
    [10] 谷志佳杨汉奇孙茂盛杨宇明 . 巨龙竹资源分布特点及其开花结实现象. 林业科学研究, 2012, 25(1): 1-5.
    [11] 杨汉奇孙茂盛阮桢媛董禹然梁宁 . 云南4种典型热带丛生竹的种源分化. 林业科学研究, 2014, 27(2): 168-173.
    [12] 张恺恺吕星杨立莹陈段芬邱德有杨艳芳 . 红豆杉TbAP2基因荧光定量PCR体系的建立及优化. 林业科学研究, 2019, 32(1): 39-46. doi: 10.13275/j.cnki.lykxyj.2019.01.006
    [13] 苏文会顾小平岳晋军吴晓丽朱如云林开搜 . 大木竹秆形结构的研究. 林业科学研究, 2006, 19(1): 98-101.
    [14] 吴涛李万峰张俊红韩素英杨文华齐力旺 . 落叶松实时定量PCR内参基因的筛选. 林业科学研究, 2013, 26(S1): 1-8.
    [15] 汪阳东 . 人工经营对毛竹秆形结构变异的影响. 林业科学研究, 2001, 14(3): 245-250.
    [16] 张文燕周道三马乃训叶长青曹德友张华明 . 几种造纸用丛生竹秆节育苗试验. 林业科学研究, 1997, 10(4): 425-428.
    [17] 黄焕华邵识烦陈拓童国建赖友明刘清浪 . 茶秆竹竹笋夜蛾综合防治技术研究. 林业科学研究, 2003, 16(4): 444-448.
    [18] 崔凯何彩云张建国廖声熙 . 毛竹茎秆组织速生的时空发育特征. 林业科学研究, 2012, 25(4): 425-431.
    [19] 格日乐图吴志民杨校生傅懋毅冯志龙 . 广宁茶秆竹地上生物量分布特征研究. 林业科学研究, 2011, 24(1): 127-131.
    [20] 周益权顾小平吴晓丽郑仁红朱维双 . 丛生竹秆基各笋目的出笋成竹生物学特性研究. 林业科学研究, 2011, 24(3): 314-320.
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  4608
  • HTML全文浏览量:  1720
  • PDF下载量:  451
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-07-03
  • 刊出日期:  2018-04-01

巨龙竹秆形发育过程实时荧光定量PCR内参基因的筛选

    通讯作者: 杨汉奇
    作者简介: 郭晓娟(1991-), 女, 硕士研究生, 研究方向:林木遗传育种.Email:lkygxj_417@163.com
  • 中国林业科学研究院资源昆虫研究所, 云南昆明 650224
基金项目:  云南省科学技术厅 2014HB041国家自然科学基金 31270662中央级公益性科研院所基本科研业务费专项 CAFYBB2017ZX001-8

摘要:  目的 利用荧光定量PCR(qRT-PCR)和表达分析软件筛选适用于巨龙竹秆形发育研究的qRT-PCR的稳定内参基因。 方法 应用常规PCR和qRT-PCR, 分析EF-1αGAPDHActinTubulinTIP-41PP2A共6个内参基因在不同秆形巨龙竹("通直型"和"弯曲型")竹笋3个不同发育时期的表达情况。利用GeNorm、NormFinder软件对各内参基因的表达稳定性进行评估; 以筛选的内参基因对巨龙竹木质素合成中的PTAL基因的表达量进行定量分析, 验证所筛选内参基因的有效性。 结果 PCR和qRT-PCR的结果表明, 6个候选内参基因PCR目的片段电泳条带单一, 熔解曲线具有明显的单一峰。内参基因ActinEF-1αGAPDH的稳定性和相关性表现最佳; 以它们为内参, 对茎秆发育过程中PTAL基因进行定量分析, 结果显示它们均具有高效性。 结论 ActinEF-1αGAPDH可以作为巨龙竹秆形发育研究中目的基因定量分析的内参基因。为进一步分析巨龙竹秆形发育过程中关键基因表达量的变化提供研究基础。

English Abstract

  • 竹亚科牡竹属的巨龙竹(Dendrocalamus sinicus Chia et J. L. Sun)为我国云南特有大型丛生竹种,其胸径可达30 cm,秆高可达30 m以上,是目前已知世界上最高大粗壮的竹种[1]。巨龙竹秆形在自然分布区中存在两种遗传稳定的变型,即“通直型”和“弯曲型”[2]。弯曲型竹秆基部常畸变,呈短缩、鼓胀、歪扭等形状,不利于竹板材加工;通直型巨龙竹秆形通直,平均单株质量100150 kg,单位面积年均产材量比我国主要经济竹种毛竹(Phyllostachys edulis(Carr.) H. de Lehaie.)高58倍[3]。通直型巨龙竹是培育特大型工业用材的优异竹种之一,也是我国南部热区竹产业发展和生态环境建设中最具发展潜力的优良经济竹种之一[3-4]。秆形作为巨龙竹资源培育与利用的关键性状之一,其发育机制尚不清楚,这给巨龙竹优良种源的选择和培育带来了障碍。同时,“通直型”和“弯曲型”巨龙竹作为一对自然的变异类型,也是木本竹种秆形发育研究中理想的对照系。目前,作者正在开展巨龙竹秆形发育的转录组分析,试图解析不同秆形发育过程中关键基因的时空表达模式。

    实时荧光定量PCR(quantitative real-time PCR,qRT-PCR)技术因其特异性强、灵敏度高、快速准确、重复性好等[5-6]优点被广泛用于基因表达分析。在植物学研究中,常选用稳定表达的内参基因转录延伸因子-(EF-1α)、甘油醛-3-磷酸-脱氢酶(GAPDH)、肌动蛋白基因(Actin)、微管蛋白基因(Tubilin)为内参基因[7-8]。随着基因芯片和基因组学的发展,又挖掘了新的内参基因蛋白磷酸酶2A(PP2A)和液泡膜内嵌蛋白基因(TIP-41)等[8]。各内参基因在不同组织、不同发育时期、不同处理中的稳定性也存在差异[8]。近期有关竹亚科植物内参基因的筛选已有一些相关报道,主要集中在毛竹等温带散生竹中,筛选出PP2ATIP41等作为毛竹不同组织基因表达定量分析的内参基因[9-11],但这些内参基因在不同发育时期的表达稳定性是未知的。针对巨龙竹不同秆形发育特点,本研究用qRT-PCR评价EF-1αGAPDHActinTubulinTIP-41PP2A 6个内参基因在巨龙竹茎秆发育过程中的表达稳定性,筛选在巨龙竹不同生长阶段、不同秆形变异种群中最佳内参基因,为下一步解析参与巨龙竹秆形发育关键基因的表达模式提供技术基础。

    • 以云南省西盟县中课种群(通直型)和孟连县景信种群(弯曲型)2个种群内01个月的笋(或幼竹)作为研究对象。巨龙竹发笋期(5月下旬至8月下旬)的定点观察结果显示,秆形变化主要出现在竹秆基部13 m内的竹节,1个月内秆形基本定型。因此,从每个类型种群中选择3个竹丛,每个竹丛选取出笋后5天、15天、30天的竹笋,分别代表秆形发育的初期、中期和后期。样品取自竹笋基部竹节,用干净的刀片切成小块立即放入液氮中,液氮处理后于-80℃贮存备用。

      表 1  采样具体信息

      Table 1.  Sampling information

      地点Site 茎秆类型Culm type 经度Longitude 纬度Latitude 海拔Altitude/m 采样时间Time
      西盟县中课镇窝笼村 通直型 99°32′06″ 22°43′48″ 1 254 2 016-07-15
      孟连县景信乡景信村 弯曲型 99°36′31″ 22°26′36″ 1 007 2 016-07-16
    • RNA提取试剂盒(TransZol Up Plus RNA Kit),反转录试剂盒(TransScriptⅡ One-Step gDNA Removal and cDNA Synthesis SuperMix),qRT-PCR荧光染料(TransStart Top Green qPCR SuperMix)均购自北京全式金生物技术有限公司。

    • 每个组织取3个不同的生物学重复。总RNA的提取参照TransZol Up Plus RNA Kit试剂盒说明书进行,提取的总RNA经2%琼脂糖凝胶电泳检测RNA完整性,Nanodrop 8000紫外分光检测仪检测OD值后直接用于反转录或保存于-80℃备用。

    • 按照TransScript Ⅱ One-Step gDNA Removal and cDNA Synthesis SuperMix试剂盒说明书,各取样品的总RNA 7μL反转录成cDNA第一链。所得cDNA产物直接用于PCR和qRT-PCR或保存于-20℃备用。

    • 使用Primer Premier 5.0软件设计候选内参基因的qRT-PCR引物,引物序列见表 2。引物由昆明硕擎生物技术有限公司合成,采用PAGE纯化。

      表 2  qRT-PCR检测中所用的6个候选内参基因的引物序列

      Table 2.  Primer sequences of six candidate reference genes for qRT-PCR analysis

      基因
      Gene name
      引物序列(5’→3’)
      Primer sequence (5’→3’)
      解链温度
      Tm/℃
      扩增产物长度Length of amplified product/bp
      PP2A F:CAGTCATCCCATGCTCTTGCT 60.6 158
      R:GAAGCTGGCTTACGTGAGATGAA 59.8
      EF-1α F:GCAACAAAACCACGCTTGAGA 59.6 186
      R:GGGTCGTGTTGAGACTGGTG 59.5
      GAPDH F:AGTCCACTGGTGTCTTCACTG 59.6 178
      R:CAATTTGTGGTGCAGCTAGCATT 60.0
      Actin F:TGCCTACATTGCCCTGGACTA 59.6 162
      R:TCCATTCCAATGAAGGATGGCT 59.8
      Tubulin F:GTCGATGCGCGAGAATACCT 59.4 186
      R:GCAGTTGCAACCATCAAGACC 59.6
      TIP-41 F:CAGCGAAATGGAAGTTCAGAAG 59.8 169
      R:TGGTCATCAGTATCCTCCCAAC 59.7
      注:正向引物(F)forwards primer,反向引物(R)reverse primer,解链温度(Tm)melting temperature。
    • 利用表 2中的引物对6个候选内参基因进行普通PCR,体系为2×Taq PCR MasterMix 10 μL,混合cDNA 2.0 μL,正向和反向引物(浓度10 μmol·L-1)各1.0 μL,ddH2O补充体积至20 μL。反应温度程序为:95℃预变性2 min;94℃变性30 s,60℃退火30 s,72℃延伸30 s,45个循环;72℃延伸10 min。反应完成后用2%琼脂糖凝胶电泳检测各PCR产物。

    • qRT-PCR采用20 μL的反应体系:2×TransStart Top Green qPCR SuperMix 10 μL,各样品的cDNA 1.0 μL(阴性对照加ddH2O 1.0 μL),正向和反向引物(浓度10 μmol·L-1)各0.4 μL,用ddH2O补足体积至20 μL。反应程序为:94℃预变性30 s;94℃变性5 s,60℃退火30 s,40个循环;6595℃溶解曲线分析。实验设计3个技术重复。加样完成后迅速放入BioRad CFX 96 real-time PCR仪上进行反应。数据读取由仪器自动完成。

    • qRT-PCR扩增完成后,采用Ct值的方法[12]计算6个内参基因的表达量。利用内参基因评价软件GeNorm和NormFinder分别对不同基因进行稳定性评价,以确定理想的内参基因。

    • 以筛选出的基因为内参,对巨龙竹木质素合成中PTAL(Phenylalanine/tyrosine ammonia-lyase)[13]基因的表达量进行定量分析,验证内参基因的有效性。体系和反应程序同1.2.5,实验设3个生物学重复,每个生物学重复设3个技术重复。根据项目组前期采用相同植物材料进行转录组测序获得的结果(待发表)和实时荧光定量分析之间基因表达量的相关性验证内参基因的有效性。

    • 提取的总RNA经2%琼脂糖凝胶电泳结果显示如图 1,所提总RNA有28S、18S两条清晰的条带,条带无拖尾,RNA完整。经Nanodrop 8000紫外分光检测仪检测,A260/280值在1.72.0范围内,A260/230值大于2.0,RNA的纯度高,表明所制备的总RNA质量较好,满足试验需求。

      图  1  巨龙竹不同秆形各发育时期总RNA电泳图

      Figure 1.  Electrophoresis results of total RNA extracted from different types of D.sinicus at three developmental stages

    • 常规PCR结果显示,各内参基因均可见预期大小的目标条带(如图 2),且条带单一,无非特异性扩增,表明引物特异性良好,可用于后续qRT-PCR实验。

      图  2  6个候选内参基因常规PCR的电泳结果

      Figure 2.  Electrophoresis results of six candidateerence genes

    • 设置5个浓度梯度,3个技术重复,标准曲线相关系数R2均大于0.99,说明了线性关系好。扩增效率E在90%105%。熔解曲线有单一峰值,峰线窄而尖,说明引物特异性好,扩增产物单一。

    • 对于不同样品或组织中基因表达的精确定量分析,选择合适的内参基因是关键,GeNorm软件和NormFinder软件是专业用来筛选内参基因的软件。

    • GeNorm软件程序采用几何平均算法,通过计算出每个内参基因稳定性的M值来筛选出稳定性较好的内参基因,M值越小内参基因的稳定性就越好,反之,稳定性越差。同时,通过归一化因子法计算出最适合的内参基因搭配组合[12, 14],来确定所需内参基因的最适数目,当Vn/(n+1)<1.5时,说明有n个基因可以作为qRT-PCR定量分析的内参基因,当Vn/(n+1)>1.5时,说明有n+1个基因可以作为qRT-PCR定量分析的内参基因。通过GeNorm分析,6个内参基因的稳定性由高到低的顺序为:EF-1α=GAPDH > Actin > TIP-41 > PP2A > Tubulin(图 3)。根据标准差数据分析(图 4),Vn/(n+1)均大于1.5,其中V5/6的M值最小,可以选择6个内参基因做组合。

      图  3  6个候选内参基因的表达稳定性GeNorm分析

      Figure 3.  Expression stabilities of six candidate reference genes analyzed by GeNorm

      图  4  6个候选内参基因标准化因子的配对差异分析

      Figure 4.  Pairwise variance analyses of six candidate reference genes

    • NormFinder软件采用方差估计算法,先获得内参基因表达稳定值,再根据稳定值(M值)大小来筛选最适合的内参基因,M值最小的内参基因为最合适的内参基因[14-15]。通过NormFinder分析,6个内参基因的稳定性由高到低的顺序为:Actin > GAPDH > EF-1α > PP2A > Tubulin > TIP-41

      综上所述:通过使用GeNorm与NormFinder 2种软件对6个候选内参基因的稳定性评估结果得知GAPDHActinEF-1α的稳定性最佳,可作为研究巨龙竹秆形发育过程基因表达分析使用的内参基因。

    • 为了验证筛选的内参基因在巨龙竹茎秆发育过程中应用的有效性,选取木质素合成过程中的酶基因PTAL,分别以ActinEF-1αGAPDH为内参,对两种秆型巨龙竹不同发育时期的材料中的PTAL基因进行定量分析。结果如图 6所示,以Actin定量的PTAL基因表达量最高,GAPDH定量的PTAL基因表达量最低。根据前期转录组测序的结果,在弯曲型巨龙竹中,PTAL基因的FPKM值在发育5天、15天、30天的材料中分别为24.15、73.52和1 029.38;在通直型巨龙竹中分别为245.31、53.82和363.4。以不同内参基因获得的表达量与转录组测序获得的表达量进行相关性分析,相关性从大到小为Actin > EF-1α > GAPDH,相关系数分别为0.99,0.95和0.92,说明筛选出的3个内参基因有效性均较高。

      图  5  6个候选内参基因的表达稳定性NormFinder分析

      Figure 5.  Expression stabilities of six candidate reference genes analyzed by NormFinder

      图  6  不同内参基因对PTAL基因的定量分析

      Figure 6.  Quantitative analysis of different reference genes to PTAL

    • 近年来,随着测序技术的快速突破,“基因组学”的研究越来越便捷,极大地促进了转录组等技术手段的迅猛发展,使之成为探索基因调控网络研究的得力工具。目前,荧光定量PCR(qRT-PCR)技术由于其更高的敏感性和特异性,已经成为基因表达深入研究的基本方法[6]。而合适稳定表达的内参基因的筛选将决定表达量数据的可靠性。本研究对不同变异类型、不同发育时期的巨龙竹竹笋样品进行内参基因筛选的结果对解析“通直型”和“弯曲型”巨龙竹差异发育的分子机理有重要意义。

      在qRT-PCR反应中,合适内参基因的选择是决定目的基因精确定量分析的先决条件。理想的内参基因应在生物体不同的生长发育阶段、组织及各种外界因素条件下均恒定表达。然而,内参基因的适用性并不是普遍的,需要从多个内参基因中筛选出表达相对稳定的一个或一个以上作为研究对象的内参基因[15]。Zeng等[16]在对枸杞属两个种Lycium barbarum L.和L. ruthenicum Murray内参基因的筛选中,确定在不同的组织细胞中,表达最稳定的内参基因并不统一。对牡丹不同发育阶段种子和花瓣进行qRT-PCR内参基因筛选的结果也显示,在不同发育阶段及不同组织,最稳定的内参基因也不同[17]。而相同内参基因在禾本科的大麦、小麦和燕麦间的通用性研究结果表明,GAPDHTUBB、18S rRNA在三者中都表达稳定,而EIF4ATUBA分别在小麦和燕麦中表达稳定,说明并不是所有内参基因在亲缘关系相近的物种中都通用[18]。在竹亚科植物中,相关内参基因筛选的研究十分有限,而且并非所有的竹亚科植物都能通用一种内参基因。对毛竹成花基因PheTFL1表达分析的研究中,选择了PP2A作为毛竹PheTFL1基因定量分析的内参基因[9]。Fan[10]等从14个候选内参基因中选择TIP41和核苷酸通道结合蛋白基因(NTB)作为毛竹开花前后表达量分析的内参基因。

      本研究中,选择EF-1αGAPDHActinTubulinTIP-41PP2A作为巨龙竹的候选内参基因,其中EF-1αGAPDHActinTubilin均为植物学研究中经典的内参基因。EF-1α主要参与快速增殖细胞的生长与增殖[19]GAPDHActinTubilin均与细胞骨架结构的形成有关[20-22]PP2ATIP-41是新挖掘的内参基因,PP2A主要参与细胞功能的调节[23]TIP-41主要参与细胞水分的转运[24]。这些内参基因已在多种植物研究中证明了其有效性[9-10, 18]

    • 通过候选内参基因的稳定性和有效性的验证分析与评估,在不同变异类型的巨龙竹笋期及其不同发育时期,表达稳定性最佳的内参基因是ActinEF-1αGAPDH。下一步将结合巨龙竹秆形发育转录组学的分析,利用筛选出的内参基因对后续筛选的目的基因进行定量分析,研究巨龙竹秆型发育过程中相关基因的表达模式。

参考文献 (24)

目录

    /

    返回文章
    返回