-
竹亚科牡竹属的巨龙竹(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]优点被广泛用于基因表达分析。在植物学研究中,常选用稳定表达的内参基因转录延伸因子-1α(EF-1α)、甘油醛-3-磷酸-脱氢酶(GAPDH)、肌动蛋白基因(Actin)、微管蛋白基因(Tubilin)为内参基因[7-8]。随着基因芯片和基因组学的发展,又挖掘了新的内参基因蛋白磷酸酶2A(PP2A)和液泡膜内嵌蛋白基因(TIP-41)等[8]。各内参基因在不同组织、不同发育时期、不同处理中的稳定性也存在差异[8]。近期有关竹亚科植物内参基因的筛选已有一些相关报道,主要集中在毛竹等温带散生竹中,筛选出PP2A、TIP41等作为毛竹不同组织基因表达定量分析的内参基因[9-11],但这些内参基因在不同发育时期的表达稳定性是未知的。针对巨龙竹不同秆形发育特点,本研究用qRT-PCR评价EF-1α、GAPDH、Actin、Tubulin、TIP-41、PP2A 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个技术重复。根据项目组前期采用相同植物材料进行转录组测序获得的结果(待发表)和实时荧光定量分析之间基因表达量的相关性验证内参基因的有效性。
巨龙竹秆形发育过程实时荧光定量PCR内参基因的筛选
Reference Gene Selection for Quantitative Real-Time PCR in Studying Culm Shape Development of Dendrocalamus sinicus
-
摘要:
目的 利用荧光定量PCR(qRT-PCR)和表达分析软件筛选适用于巨龙竹秆形发育研究的qRT-PCR的稳定内参基因。 方法 应用常规PCR和qRT-PCR, 分析EF-1α、GAPDH、Actin、Tubulin、TIP-41、PP2A共6个内参基因在不同秆形巨龙竹("通直型"和"弯曲型")竹笋3个不同发育时期的表达情况。利用GeNorm、NormFinder软件对各内参基因的表达稳定性进行评估; 以筛选的内参基因对巨龙竹木质素合成中的PTAL基因的表达量进行定量分析, 验证所筛选内参基因的有效性。 结果 PCR和qRT-PCR的结果表明, 6个候选内参基因PCR目的片段电泳条带单一, 熔解曲线具有明显的单一峰。内参基因Actin、EF-1α、GAPDH的稳定性和相关性表现最佳; 以它们为内参, 对茎秆发育过程中PTAL基因进行定量分析, 结果显示它们均具有高效性。 结论 Actin、EF-1α、GAPDH可以作为巨龙竹秆形发育研究中目的基因定量分析的内参基因。为进一步分析巨龙竹秆形发育过程中关键基因表达量的变化提供研究基础。 Abstract:Objective To select suitable quantitative real-time PCR (qRT-PCR) reference genes for studying culm shape development of Dendrocalamus sinicus Chia et J. L. Sun. Method PCR and qRT-PCR were used to analyze the mRNA expression stability of six candidate reference genes (EF-1α, GAPDH, Actin, TIP-41, Tubulin and PP2A) in a set of six bamboo shoot samples at three development stages of the straight-culmed and bending-culmed D. sinicus respectively. The software GeNorm and NomFinder were employed to evaluate the data for reference genes. Quantitative analysis of PTAL was used to test the effectiveness of the candidate reference genes. Result The results showed that six candidate reference genes were of expected size and single peaking melting curve. Of them, Actin, EF-1α and GAPDH had the best stability and correlativity. Meanwhile, these three candidate reference genes had the highest efficiency based on quantitative analysis results of the PTAL gene during the culm development. Conclusion The Actin, EF-1α and GAPDH were optional inference genes for normalizing purpose genes in studying culm shape development of D. sinicus. The study can provide references for gene expression analysis in the culm development of D. sinicus. -
Key words:
- Dendrocalamus sinicus
- / culm shape development
- / reference genes
- / qRT-PCR
-
表 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 表 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。 -
[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.