留言板

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

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

菜豆黄单胞菌侵染下木薯内参基因筛选与评价

李慈云 杨静 卢东莹 杨健飞 杨琦 王菲 陈银华 牛晓磊

李慈云,杨静,卢东莹,等. 菜豆黄单胞菌侵染下木薯内参基因筛选与评价[J]. 热带生物学报,2025, 16(0):1−8. doi:10.15886/j.cnki.rdswxb.20240201 doi:  10.15886/j.cnki.rdswxb.20240201
引用本文: 李慈云,杨静,卢东莹,等. 菜豆黄单胞菌侵染下木薯内参基因筛选与评价[J]. 热带生物学报,2025, 16(0):1−8. doi:10.15886/j.cnki.rdswxb.20240201 doi:  10.15886/j.cnki.rdswxb.20240201
LI Ciyun, YANG Jing, LU Dongying, YANG Jianfei, YANG Qi, WANG Fei, Chen Yinhua, NIU Xiaolei. Screening and Evaluation of Reference Genes in Cassava unpon Infection by Xanthomonas phaseoli pv. manihotis[J]. Journal of Tropical Biology. doi: 10.15886/j.cnki.rdswxb.20240201
Citation: LI Ciyun, YANG Jing, LU Dongying, YANG Jianfei, YANG Qi, WANG Fei, Chen Yinhua, NIU Xiaolei. Screening and Evaluation of Reference Genes in Cassava unpon Infection by Xanthomonas phaseoli pv. manihotis[J]. Journal of Tropical Biology. doi: 10.15886/j.cnki.rdswxb.20240201

菜豆黄单胞菌侵染下木薯内参基因筛选与评价

doi: 10.15886/j.cnki.rdswxb.20240201
基金项目: 海南省国际科技合作研发项目(GHYF2022005);海南省高等学校教育教学改革研究重点项目(Hnjg2019ZD-2)
详细信息
    第一作者:

    李慈云(2000—),女,海南大学热带农林学院2022级硕士研究生。E-mail:1695841062@qq.com

    通信作者:

    牛晓磊(1979—),男,副教授,博士生导师。研究方向:作物遗传育种。E-mail:ninterxll@hainanu.edu.cn

  • 中图分类号: S476

Screening and Evaluation of Reference Genes in Cassava unpon Infection by Xanthomonas phaseoli pv. manihotis

  • 摘要: 为筛选出菜豆黄单胞菌(Xanthomonas phaseoli pv. manihotis, Xpm)侵染木薯(Manihot esculenta)叶片后稳定表达的内参基因,确保后续基因功能研究的可靠性,以木薯抗病品种‘JG1301’和感病品种‘SC9’叶片为研究材料,在病原菌侵染后0、6、12、24、72和120 h时采样,对常用内参基因NascentEF1aACTGTPbTUB在接种不同时间点后的稳定性进行分析。结果表明:在Xpm侵染抗病品种‘JG1301’和感病品种‘SC9’木薯叶片后,5个内参基因表达量变化均有差异;使用Delta Ct、GeNorm、NormFinder和 BestKeeper软件对这些内参基因稳定性进行排序,并结合RefFinder软件 进行各内参基因稳定性的综合排名,结果发现,Xpm侵染木薯叶片后表达最稳定的内参基因为 EF1α。本研究筛选出了Xpm侵染木薯叶片后表达稳定的内参基因,为木薯响应Xpm基因表达和抗病基因的进一步挖掘奠定了基础。
  • 图  2  5个候选内参基因Ct值分布

    横跨该框的线表示中位数。箱部上方–上四分位数,箱部下方–下四分位数,上下误差线表示最大值和最小值。

    Fig.  2  Ct value distribution of five candidate reference genes

    The line across the box represents the median. Upper part of the box-upper quartile, lower part of the box-lower quartile, The upper and lower error bars indicate the maximum and minimum values.

    图  3  Genorm分析

    A. 内参基因稳定性;B. 配对变异分析。

    Fig.  3  Genorm analysis

    A. Reference gene stability; B. Pairwise variation analysis.

    表  1  内参基因及其对应的引物

    Table  1  Reference genes and their corresponding primers

    基因符号
    Gene
    引物序列(5′−3′)
    Primer sequences(5′−3′)
    基因登录号
    Accession No.
    扩增长度/bp
    Length/bp
    引物扩增效率/%
    Efficiency/%
    决定系数R2
    Correlation R2
    参考文献
    References
    TUB F: gtggaggaactggttctgga
    R: tgcactcatctgcattctcc
    Manes.08G061700 180 105.5 0.997 [15]
    Nascent F: ttttgaagctgctgcagagg
    R: tgaagcggatctgaacaacc
    Manes.09G005100 192 100.7 0.997 [18]
    GTPb F: ggtgggtaacaaggcagaca
    R: tgagtcggtatctgcaagcc
    Manes.09G086600 169 102.4 0.993 [19]
    ACT F: ttcgtgtcaaggtgtcgtga
    R: gccctctcatttgctgcaat
    Manes.13G086400 149 106.7 0.999 [20]
    EF1α F: tgaaccaccctggtcagattggaa
    R: aacttgggctccttctcaagctct
    Manes.15G054800 139 109.4 0.998 [21, 22]
    下载: 导出CSV

    表  2  Delta Ct软件分析结果

    Table  2  Results of Delta Ct Software

    基因
    Gene
    平均标准差
    m SD
    排名
    Ranking
    Nascent 0.75 1
    GTPb 0.77 2
    ACT 0.92 3
    EF1α 1.00 4
    TUB 1.18 5
    下载: 导出CSV

    表  3  Normfinder 软件分析结果

    Table  3  Results of Normfinder Software

    基因
    Gene
    稳定值
    Stable value
    排名
    Ranking
    EF1α 0.348 1
    GTPb 0.381 2
    Nascent 0.444 3
    ACT 0.520 4
    TUB 0.714 5
    下载: 导出CSV

    表  4  Bestkeeper 软件分析结果

    Table  4  Results of BesKeeper Software

    参数
    Parameter
    基因 Gene
    TUBNascentGTPbACTEF1α
    几何平均值[CP]
    geo Mean [CP]
    24.65 22.64 24.044 23.36 20.37
    算术平均值[CP]
    ar Mean [CP]
    24.68 22.65 24.05 23.38 20.39
    最小值[CP]
    min [CP]
    22.30 21.02 22.59 21.1 18.46
    最大值[CP]
    max [CP]
    26.89 24.15 25.36 24.84 22.48
    标准差[± CP]
    SD [± CP]
    0.97 0.61 0.66 0.75 0.78
    变异系数[% CP]
    CV [% CP]
    3.91 2.68 2.73 3.21 3.84
    调节系数标准差[± x-fold]
    std dev [± x-fold]
    1.95 1.52 1.58 1.68 1.72
    下载: 导出CSV

    表  5  RefFinder分析结果

    Table  5  Results of RefFinder Software

    基因
    Gene
    稳定值
    Stable value
    排名
    Ranking
    EF1α 1.41 1
    Nascent 1.73 2
    GTPb 2.21 3
    ACT 3.72 4
    TUB 5.00 5
    下载: 导出CSV
  • [1] NOLAN T, HANDS R E, BUSTIN S A. Quantification of mRNA using real-time RT-PCR[J]. Nature Protocols, 2006, 1(3): 1559 − 1582. doi:  10.1038/nprot.2006.236
    [2] SOWA S, SOZONIUK M, TOPOROWSKA J, et al. Validation of reference genes as an internal control for studying Avena sativa–Puccinia coronata interaction by RT-qPCR[J]. Scientific Reports, 2022, 12(1): 14601. doi:  10.1038/s41598-022-18746-z
    [3] SU B, LI Z, LIU H, et al. Identification and validation of reference genes for RT-qPCR analysis in Iris domestica under Cd stress[J]. Heliyon, 2024, 10(17): e36923. doi:  10.1016/j.heliyon.2024.e36923
    [4] YANG Y, CAO G, TANG L. Selection and validation of reference genes for qRT-PCR normalization in dayflower (Commelina communis) based on the transcriptome profiling[J]. BMC Plant Biology, 2024, 24(1): 1131 − 1143. doi:  10.1186/s12870-024-05853-4
    [5] VANGUILDER H D, VRANA K E, FREEMAN W M. Twenty-Five Years of Quantitative PCR for Gene Expression Analysis[J]. BioTechniques, 2008, 44(5): 619 − 626. doi:  10.2144/000112776
    [6] TANG X, ZHANG N, SI H, et al. Selection and validation of reference genes for RT-qPCR analysis in potato under abiotic stress[J]. Plant Methods, 2017, 13(1): 85. doi:  10.1186/s13007-017-0238-7
    [7] ZHANG J, HE X, ZHOU J, et al. Selection and Verification of Standardized Reference Genes of Angelica dahurica under Various Abiotic Stresses by Real-Time Quantitative PCR[J]. Genes, 2024, 15(1): 79 − 95. doi:  10.3390/genes15010079
    [8] JEONG Y M, MUN J H, LEE I, et al. Distinct Roles of the First Introns on the Expression of Arabidopsis Profilin Gene Family Members[J]. Plant Physiology, 2006, 140(1): 196 − 209. doi:  10.1104/pp.105.071316
    [9] 方佳, 濮文辉, 张慧坚. 国内外木薯产业发展近况[J]. 中国农学通报, 2010, 26(16): 353 − 361.
    [10] 梁露锋, 玉琼广, 刘洁. 木薯产业发展动态及展望前景[J]. 大众科技, 2011, 128(6): 128 − 130. doi:  10.3969/j.issn.1008-1151.2011.06.057
    [11] 严华兵, 叶剑秋, 李开绵. 中国木薯育种研究进展[J]. 中国农学通报, 2015, 31(15): 63 − 70. doi:  10.11924/j.issn.1000-6850.casb14110159
    [12] 刘佳, 张箭. 美洲木薯在非洲的引种和推广[J]. 世界农业, 2019, 14(2): 83 − 88.
    [13] 陈松笔, 蔡杰, 安飞飞, 等. 木薯育种现状及发展趋势[J]. 中国科学: 生命科学, 2024, 54(10): 1833 − 1842.
    [14] AMELEWORK A B, BAIRU M W. Advances in Genetic Analysis and Breeding of Cassava (Manihot esculenta Crantz): A Review[J]. Plants, 2022, 11(12): 1617 − 1635. doi:  10.3390/plants11121617
    [15] CARLOS A. Z -C, CRUZ G D, VERDIER V, et al. Cassava diseases caused by Xanthomonas phaseoli pv. manihotis cassavae[J]. Molecular Plant Pathology, 2021, 22(12): 1520 − 1537. doi:  10.1111/mpp.13094
    [16] 李超萍, 时涛, 刘先宝, 等. 国内木薯病害普查及细菌性萎蔫病安全性评估[J]. 热带作物学报, 2011, 32(1): 116 − 121.
    [17] CARLOSA Z C, MOUFID Y, LóPEZ C E, et al. First Report of Cassava Bacterial Blight Caused by Xanthomonas phaseoli pv. manihotis in the Amazonian Forest of Ecuador[J]. Plant Disease, 2024, 108(6): 1879.
    [18] HU M, HU W, XIA Z, et al. Validation of Reference Genes for Relative Quantitative Gene Expression Studies in Cassava (Manihot esculenta Crantz) by Using Quantitative Real-Time PCR[J]. Frontiers in Plant Science, 2016, 7: 680.
    [19] MORENO I, GRUISSEM W, VANDERSCHUREN H. Reference genes for reliable potyvirus quantitation in cassava and analysis of Cassava brown streak virus load in host varieties[J]. Journal of Virological Methods, 2011, 177(1): 49 − 54. doi:  10.1016/j.jviromet.2011.06.013
    [20] CZECHOWSKI T, STITT M, ALTMANN T, et al. Genome-Wide Identification and Testing of Superior Reference Genes for Transcript Normalization in Arabidopsis[J]. Plant Physiology, 2005, 139(1): 5 − 17. doi:  10.1104/pp.105.063743
    [21] CHANG Y, BAI Y, WEI Y, et al. CAMTA3 negatively regulates disease resistance through modulating immune response and extensive transcriptional reprogramming in cassava[J]. Tree Physiology, 2020, 40(11): 1520 − 1533. doi:  10.1093/treephys/tpaa093
    [22] ZHENG L, GAO S, BAI Y, et al. NF‐YC15 transcription factor activates ethylene biosynthesis and improves cassava disease resistance[J]. Plant Biotechnology Journal, 2024, 22(9): 2424 − 2434. doi:  10.1111/pbi.14355
    [23] ZHU S, PAN Y, LI K, et al. Complete Genome Sequence of Xanthomonas phaseoli pv. manihotis Strain CHN01, the Causal Agent of Cassava Bacterial Blight[J]. Plant Disease, 2022, 106(3): 1039 − 1041. doi:  10.1094/PDIS-09-21-2016-A
    [24] WEI Y, SHI H, XIA Z, et al. Genome-Wide Identification and Expression Analysis of the WRKY Gene Family in Cassava[J]. Frontiers in Plant Science, 2016, 7: 25.
    [25] TONG Z, GAO Z, WANG F, et al. Selection of reliable reference genes for gene expression studies in peach using real-time PCR[J]. BMC Molecular Biology, 2009, 10: 71. doi:  10.1186/1471-2199-10-71
    [26] ZHONG H Y, CHEN J W, LI C Q, et al. Selection of reliable reference genes for expression studies by reverse transcription quantitative real-time PCR in litchi under different experimental conditions[J]. Plant Cell Reports, 2011, 30(4): 641 − 653. doi:  10.1007/s00299-010-0992-8
    [27] NICOT N, HAUSMAN J F, HOFFMANN L, et al. Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress[J]. Journal of Experimental Botany, 2005, 56(421): 2907 − 2914. doi:  10.1093/jxb/eri285
    [28] GAO K, KHAN WU, LI J, et al. Identification and Validation of Reliable Reference Genes for Gene Expression Studies in Koelreuteria paniculata[J]. Genes, 2022, 13(5): 714 − 725. doi:  10.3390/genes13050714
    [29] LI X, CHENG J, ZHANG J, et al. Validation of Reference Genes for Accurate Normalization of Gene Expression in Lilium davidii var. unicolor for Real Time Quantitative PCR[J]. Plos One, 2015, 10(10): 28 − 43.
    [30] LI G, SUN X, ZHU X, et al. Selection and Validation of Reference Genes in Virus-Infected Sweet Potato Plants[J]. Genes, 2023, 14(7): 1477 − 1488. doi:  10.3390/genes14071477
    [31] MASCIA T, SANTOVITO E, GALLITELLI D, et al. Evaluation of reference genes for quantitative reverse‐transcription polymerase chain reaction normalization in infected tomato plants[J]. Molecular Plant Pathology, 2010, 11(6): 805 − 816. doi:  10.1111/j.1364-3703.2010.00646.x
    [32] 唐枝娟, 刘秦, 肖晓蓉, 等. 白叶枯病菌侵染下的水稻内参基因稳定性[J]. 分子植物育种, 2017, 15(1): 300 − 306.
  • [1] 张博文, 徐浩然, 曾鸿秋.  MeJAZ2.2基因的克隆及互作蛋白筛选 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20240169
    [2] 唐璐芝, 吴竞远, 陈元来, 朱彬彬, 赵惠萍.  木薯RXam2的原核蛋白表达及抗病功能分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20240147
    [3] 胡文成, 朱寿松, 王艺璇, 阳达, 陈银华.  木薯茉莉酸氧化酶JOXs基因家族的鉴定及表达模式分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20240199
    [4] 缪蝶, 张舒心, 王鹏.  番茄不同组织及不同激素处理下内参基因的筛选 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20240140
    [5] 陈雨, 程雨果, 蔡静, 陶均.  GNAT类乙酰转移酶调控水稻细菌性条斑病菌的生长和致病力 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230050
    [6] 贾新蕾, 黄增朝, 杨林狄, 吕静, 李妍萍, 简纪常, 黄郁葱.  不同培养温度的鱼源海豚链球菌转录组分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230027
    [7] 赵金铃, 潘晋龙, 吴宇垚, 韩昊祥, 张怀方, 陈新, 曾长英.  不同氮源处理对4个木薯品种幼苗生长发育的影响 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230018
    [8] 张逸杰, 王思琦, 陆小静, 宋记明, 王睿, 曹敏, 张瑞, 王红刚, 吴金山.  木薯块根灰分和水分近红外光谱预测模型的构建与优化 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230031
    [9] 冯亚亭, 张逸杰, 林南方, 陈银华, 骆凯.  木薯不同组织内生细菌多样性的比较分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230077
    [10] 贾素行, 朱寿松, 符仁稳, 李春霞, 陈银华.  木薯MePPD3基因的克隆及功能分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230035
    [11] 李昉泽, 詹剑.  木薯渣和甘蔗渣基生物炭对砖红壤的改良效果 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230071
    [12] 徐春华, 李超萍, 时涛, 王国芬, 蔡吉苗, 李博勋, 黄贵修.  木薯细菌性萎蔫病菌2个候选抗铜基因簇的克隆和分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20220009
    [13] 戴云鹏, 潘晋龙, 吴宇垚, 韩昊祥, 张怀方, 陈新, 曾长英.  2种氮素的分根处理对木薯幼苗生长发育和氮效率的影响 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.20230008
    [14] 刘满意, 王禹童, 孙铭泽, 李荣, 王蓓蓓.  套作白三叶草对香蕉枯萎病发病率及土壤微生物群落的影响 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2021.02.011
    [15] 邓雅洁, 刘云豪, 孙琪, 邱越, 宋家明, 赖杭桂.  在干旱胁迫下有性四倍体木薯叶片的蛋白质组学研究 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2021.02.010
    [16] 黄文枫, 胡艳平, 朱白婢.  海南薄皮甜瓜枯萎病病原菌鉴定及抗病砧木的筛选 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2020.03.008
    [17] 赵克学, 许淑娟, 马瑞仙, 耿金静, 刘杨, 李琼毅, 冯若飞.  TRIM21基因SYBR GreenⅠ实时荧光定量PCR检测方法的建立 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2020.01.017
    [18] 梁峻玮, 曾凡云, 漆艳香, 彭军, 张欣, 谢培兰, 谢艺贤.  香蕉枯萎病菌4号生理小种中性海藻糖酶NTH1基因敲除和功能分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2020.02.010
    [19] 赵锐, 樊若晨, 丁凯旋, 唐枝娟, 王红刚, 黄思源, 蒋凌雁, 陈银华, 耿梦婷.  木薯MeNOSIP与MeRbohD蛋白的互作验证及其对外源激素的响应 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2020.02.007
    [20] 宋娜, 张刘宁颖, 曹敏, 郭文雅, 洪雨慧, 吴金山, 陈银华, 于晓惠.  木薯ERF转录因子调控的靶基因筛选与表达分析 . 热带生物学报, doi: 10.15886/j.cnki.rdswxb.2020.02.008
  • 加载中
计量
  • 文章访问数:  26
  • HTML全文浏览量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-12-30
  • 录用日期:  2025-02-23
  • 修回日期:  2025-02-16

菜豆黄单胞菌侵染下木薯内参基因筛选与评价

doi: 10.15886/j.cnki.rdswxb.20240201
    基金项目:  海南省国际科技合作研发项目(GHYF2022005);海南省高等学校教育教学改革研究重点项目(Hnjg2019ZD-2)
    作者简介:

    李慈云(2000—),女,海南大学热带农林学院2022级硕士研究生。E-mail:1695841062@qq.com

    通讯作者: 牛晓磊(1979—),男,副教授,博士生导师。研究方向:作物遗传育种。E-mail:ninterxll@hainanu.edu.cn
  • 中图分类号: S476

摘要: 为筛选出菜豆黄单胞菌(Xanthomonas phaseoli pv. manihotis, Xpm)侵染木薯(Manihot esculenta)叶片后稳定表达的内参基因,确保后续基因功能研究的可靠性,以木薯抗病品种‘JG1301’和感病品种‘SC9’叶片为研究材料,在病原菌侵染后0、6、12、24、72和120 h时采样,对常用内参基因NascentEF1aACTGTPbTUB在接种不同时间点后的稳定性进行分析。结果表明:在Xpm侵染抗病品种‘JG1301’和感病品种‘SC9’木薯叶片后,5个内参基因表达量变化均有差异;使用Delta Ct、GeNorm、NormFinder和 BestKeeper软件对这些内参基因稳定性进行排序,并结合RefFinder软件 进行各内参基因稳定性的综合排名,结果发现,Xpm侵染木薯叶片后表达最稳定的内参基因为 EF1α。本研究筛选出了Xpm侵染木薯叶片后表达稳定的内参基因,为木薯响应Xpm基因表达和抗病基因的进一步挖掘奠定了基础。

English Abstract

李慈云,杨静,卢东莹,等. 菜豆黄单胞菌侵染下木薯内参基因筛选与评价[J]. 热带生物学报,2025, 16(0):1−8. doi:10.15886/j.cnki.rdswxb.20240201 doi:  10.15886/j.cnki.rdswxb.20240201
引用本文: 李慈云,杨静,卢东莹,等. 菜豆黄单胞菌侵染下木薯内参基因筛选与评价[J]. 热带生物学报,2025, 16(0):1−8. doi:10.15886/j.cnki.rdswxb.20240201 doi:  10.15886/j.cnki.rdswxb.20240201
LI Ciyun, YANG Jing, LU Dongying, YANG Jianfei, YANG Qi, WANG Fei, Chen Yinhua, NIU Xiaolei. Screening and Evaluation of Reference Genes in Cassava unpon Infection by Xanthomonas phaseoli pv. manihotis[J]. Journal of Tropical Biology. doi: 10.15886/j.cnki.rdswxb.20240201
Citation: LI Ciyun, YANG Jing, LU Dongying, YANG Jianfei, YANG Qi, WANG Fei, Chen Yinhua, NIU Xiaolei. Screening and Evaluation of Reference Genes in Cassava unpon Infection by Xanthomonas phaseoli pv. manihotis[J]. Journal of Tropical Biology. doi: 10.15886/j.cnki.rdswxb.20240201
  • 实时荧光定量PCR(Reverse transcription quantitative real-time PCR, qRT-PCR)是基因表达研究中最准确、最可靠的技术之一[12]。该方法通过荧光信号测定PCR 循环扩增产物总量,并利用内参基因作为参考对特定cDNA序列进行定量分析[3, 4]。qRT-PCR常用于检测基因在不同处理、不同生长发育时期的表达量变化,而内参基因的稳定性决定了检测基因表达量变化的准确性。因此,使用稳定表达的内参基因对目的基因表达水平的标准化至关重要[5]。好的内参基因参与细胞的基础代谢功能,并在整个过程中保持稳定表达,不受外部条件改变的影响[6, 7]。然而,基因在生物体中不同组织和不同环境下行使功能的差异会导致没有任何一个内参基因能够在所有条件下都恒定表达[8]。因此,应根据具体的实验条件从而筛选出表达量相对稳定的内参基因来对靶标基因的表达进行分析,从而获得可靠的试验结果。

    木薯(Manihot esculenta),是大戟科(Euphorbiaceae)木薯属植物,与马铃薯(Solanum tuberosum )、红薯(Ipomoea batatas )并列为世界重要的薯类作物,被誉为“地下粮仓”和“淀粉之王”[911]。木薯具有多种优良特性,如耐瘠、耐旱、适应性强、管理粗放、容易栽培和淀粉含量高等特点[12];块根淀粉含量丰富且种植易成活的优势,成为了热带地区低收入国家的主要人类粮食[13]。木薯还能作为酒精的主要生产原料,也可以用以制作沼气、纺织、医疗药品等[1314]。但病害是制约木薯产业健康、持续发展的重要因素之一。世界范围内,危害木薯生产的病害高达30多种,其中木薯细菌性枯萎病(Cassava Bacterial Blight,CBB)的危害尤其严重[15]。CBB是由菜豆黄单胞菌木薯萎蔫致病变种(Xanthomonas phaseoli pv. manihotis, Xpm)引起的检疫性病害,导致木薯严重减产甚至绝收[1617]。尽管对CBB的研究已取得一定进展,但大多停留在生理水平,其在分子层面上的抗病机制仍有待进一步阐明。目前,尚未有木薯在Xpm侵染条件下稳定表达内参基因的研究报道,因此,筛选在Xpm侵染木薯过程中稳定表达的内参基因是研究木薯响应细菌性枯萎病相关功能基因定量表达准确性的关键因素,可为木薯抗病分子育种基因功能研究提供衡量基因表达水平的标尺。

    • 木薯抗病品种‘JG1301’、感病品种‘SC9’和菜豆黄单胞菌Xpm均由海南大学热带农林学院热带作物与微生物互作团队提供。

    • 随机选择生长良好、高度和长势一致的3株生长45 d的木薯‘JG1301’和SC9’,采用叶片注射法接种病原菌Xpm,每个品种接种3株,每株接种3个掌叶,每片掌叶接种2个复叶,接种后分别于0、6、12、72和120 h 采取叶片,每个时间点从3株木薯各采1复叶(3次重复),采集后放入离心管中,液氮速冻后,−80 ℃保存备用。

    • 使用天根RNAprep Pure多糖多酚植物总RNA提取试剂盒(DP432)按照说明书方法进行木薯叶片总RNA的提取。用微型分光光度计(Nano-300 Micro-Spectrophotometer)检测RNA浓度和质量。通过w=1.0% 琼脂凝胶电泳检测RNA完整度。采用天根FastKing一步法除基因组cDNA第一链合成预混试剂(KR118),按照说明书方法进行反转录,使用1 μg总RNA建立20 μL cDNA逆转录反应体系。反转录后的cDNA于−80 ℃冰箱保存备用。

    • 根据文献[15, 1822]表1)查找木薯研究中qRT-PCR实验常用的内参基因,并最终筛选出实验所需要的 5个内参基因,分别为TUBNascentGTPbACTEF1a,通过在Phytozome数据库中查询基因序列,确定各候选基因的CDS序列,根据基因的CDS序列使用在线工具Primer3 web(https://primer3.ut.ee/)进行引物设计,通过Phytozome Blast鉴定引物特异性后,送北京擎科生物科技股份有限公司进行引物合成,引物序列见表1

      表 1  内参基因及其对应的引物

      Table 1.  Reference genes and their corresponding primers

      基因符号
      Gene
      引物序列(5′−3′)
      Primer sequences(5′−3′)
      基因登录号
      Accession No.
      扩增长度/bp
      Length/bp
      引物扩增效率/%
      Efficiency/%
      决定系数R2
      Correlation R2
      参考文献
      References
      TUB F: gtggaggaactggttctgga
      R: tgcactcatctgcattctcc
      Manes.08G061700 180 105.5 0.997 [15]
      Nascent F: ttttgaagctgctgcagagg
      R: tgaagcggatctgaacaacc
      Manes.09G005100 192 100.7 0.997 [18]
      GTPb F: ggtgggtaacaaggcagaca
      R: tgagtcggtatctgcaagcc
      Manes.09G086600 169 102.4 0.993 [19]
      ACT F: ttcgtgtcaaggtgtcgtga
      R: gccctctcatttgctgcaat
      Manes.13G086400 149 106.7 0.999 [20]
      EF1α F: tgaaccaccctggtcagattggaa
      R: aacttgggctccttctcaagctct
      Manes.15G054800 139 109.4 0.998 [21, 22]
    • 对每一对内参基因引物制作对应的标准曲线,计算相应引物的扩增效率。以ddH2O作为阴性对照,用于检测实验过程中的试剂或者人为污染。所有样品重复3次。采用耶拿荧光定量PCR qTOWER384G(Analytik Jena)进行qRT-PCR实验,根据测定的Ct值和对应的不同稀释梯度模板浓度经线性回归(其斜率值即为 slope),计算得各候选基因的引物扩增效率(E=10−1/slope)及决定系数R2

    • 采用耶拿荧光定量PCR qTOWER384G(Analytik Jena)进行qRT-PCR实验,每个反应均包含正向和反向引物各2 μL(10 μM),不同稀释倍数的cDNA模板5 μL(20 ng/μL),10 μL 2× SYBR Greenq PCR Mix(2×SYBR Green qPCR Mix,思科捷),加 ddH2O至20 μL反应体系。 每个反应体系进行3次重复。qRT-PCR 反应条件为:94 ℃预变性 3 min,随后进行 94 ℃变性10 s,60 ℃延伸34 s,循环40次的二步法程序。在每个实验结束时,使用默认参数在60 ~94 °C以0.3 ℃的增量持续60 s的周期进行熔解曲线分析。

    • Xpm侵染木薯‘JG1301’和‘SC9’叶片不同时间点cDNA等量混合做为模板,进行普通PCR扩增。利用PCR对内参基因引物特异性进行检测,在2%琼脂糖凝胶电泳中,内参基因扩增产物在100~250 bp间,且目标条带单一,无非特异性条带和引物二聚体,扩增片段条带大小符合预期 (图1−A)。以Xpm侵染木薯‘JG1301’和‘SC9’叶片后不同时间点的cDNA为模板,经过qRT-PCR扩增后,对候选扩增片段的溶解曲线进行分析,结果表明,所选择的5个木薯候选内参基因引物特异性好,无杂峰,可用于后续内参基因稳定性研究(图1−B)。

      图  1  内参基因引物检测

      Figure 1.  Reference gene primer detection

      通过5个木薯候选内参基因引物的标准曲线得到引物的扩增效率和决定系数,结果表明,所选择的5个木薯候选内参基因荧光定量引物扩增效率在100.7%~109.4%,扩增效率最高的基因是EF1α,引物扩增效率为109.4%,引物扩增效率最低的是Nascent基因,扩增效率为100.7%。5个引物的决定系数R2 在0.993~0.999间,决定系数最高的基因是ACT基因,决定系数为0.999,决定系数最低的是GTPb基因,决定系数为0.993(表1)。

    • 对 5 对木薯内参基因引物以Xpm侵染抗病品种‘JG1301’和感病品种‘SC9’木薯叶片后不同时间点的cDNA为模板进行qRT-PCR实验,通过检测其表达量得到对应的 Ct 值。Ct 值与基因的表达丰度成反比,Ct 值越大,基因的表达丰度越低,反之Ct 值越小,则基因的表达丰度越高。从图2可知,5个木薯候选内参基因中表达丰度最高的基因是EF1α,Ct值在样本检测中为18.46~22.48;表达丰度最低的基因是TUB,Ct值在样本检测中为22.30~26.89。

      图  2  5个候选内参基因Ct值分布

      Figure 2.  Ct value distribution of five candidate reference genes

    • Delta Ct算法是用平均标准差(m SD)来评价内参基因的稳定性。m SD值的大小能够反应内参基因的稳定性,数值越小代表内参基因越稳定,越大则越不稳定。通过Delta Ct方法计算5个木薯候选内参基因的稳定性可知,5个候选内参基因的m SD排序为TUB > EF1α > ACT > GTPb > Nascent, 其中m SD最大的是TUB,其m SD为1.18;m SD最小的是GTPb,其m SD为0.75,表明稳定性最好的是Nascent基因,稳定性最差的是TUB基因(表2)。

      表 2  Delta Ct软件分析结果

      Table 2.  Results of Delta Ct Software

      基因
      Gene
      平均标准差
      m SD
      排名
      Ranking
      Nascent 0.75 1
      GTPb 0.77 2
      ACT 0.92 3
      EF1α 1.00 4
      TUB 1.18 5
    • GeNorm 程序是根据平均变异度M值来衡量基因的稳定性,程序默认的取舍值为1.5,当M值<1.5 时,基因被认为相对稳定,M值越低则表示基因越稳定。该程序采用成对比较法(Pairwise comparison approach)计算出基因表达稳定性的M值,根据Genorm默认变异系数(Vn/Vn+1)的值来确认最优的内参基因数目。当Vn/Vn+1<0.15时,则认为最佳的内参基因数为n个,若Vn/Vn+1>0.15还可引入1个新的内参基因后计算标准化因子的平均配对变异V值,并根据Vn/Vn+1 值来确定所需最适内参基因数目。

      图3所示,5个木薯内参候选基因的M值都小于1.5(图3−A),但是表达稳定性有较大差异,M值排序为TUB > ACT > GTPb >EF1α= Nascent,其中M值最大的是ACTM值为1.00,M值最小的是EF1α和Nascent,M值为0.66,说明稳定性最好的是EF1αNascent,稳定性最差的是TUB。其中V4/5配对变异系数为0.141<0.15(图3−B),说明最适的内参基因数4个,不需要引入第5个基因进行校准。因此,本实验中最优基因数组合为4个,即EF1αNascentGTPbACT

      图  3  Genorm分析

      Figure 3.  Genorm analysis

    • 该程序基于模型分析,获得内参基因的表达稳定值,再根据稳定值大小来筛选出最合适的内参基因。NormFinder 程序根据候选内参基因的稳定值进行评价,稳定值越高表示稳定性越差,稳定值最小的基因,则为最稳定的基因。结果表明5个木薯候选内参基因的稳定值排序为TUB >ACT > Nascent > GTPb > EF1α,其中稳定值最大的是TUB,稳定值为0.714;稳定值最小的是EF1α,稳定值为0.348,说明稳定性最好的是EF1α基因,稳定性最差的是TUB基因(表3)。

      表 3  Normfinder 软件分析结果

      Table 3.  Results of Normfinder Software

      基因
      Gene
      稳定值
      Stable value
      排名
      Ranking
      EF1α 0.348 1
      GTPb 0.381 2
      Nascent 0.444 3
      ACT 0.520 4
      TUB 0.714 5
    • BestKeeper 程序以内参基因 Ct 值的标准差(SD)来衡量基因的表达稳定性,SD值越小,表达越稳定,SD 临界值默认为1,当 SD值>1 时,说明该基因表达不稳定。采用重复配对相关分析和回归分析,通过计算Ct值的标准偏差(SD)来确定稳定性较好的内参基因。结果表明,5个木薯候选内参的SD值均小于1,SD值排序为TUB > EF1α > ACT > GTPb > Nascent,其中SD值最大的是TUBSD值为0.97;SD值最小的是NascentSD值为0.61,说明稳定性最好的是Nascent,稳定性最差的是TUB表4)。

      表 4  Bestkeeper 软件分析结果

      Table 4.  Results of BesKeeper Software

      参数
      Parameter
      基因 Gene
      TUBNascentGTPbACTEF1α
      几何平均值[CP]
      geo Mean [CP]
      24.65 22.64 24.044 23.36 20.37
      算术平均值[CP]
      ar Mean [CP]
      24.68 22.65 24.05 23.38 20.39
      最小值[CP]
      min [CP]
      22.30 21.02 22.59 21.1 18.46
      最大值[CP]
      max [CP]
      26.89 24.15 25.36 24.84 22.48
      标准差[± CP]
      SD [± CP]
      0.97 0.61 0.66 0.75 0.78
      变异系数[% CP]
      CV [% CP]
      3.91 2.68 2.73 3.21 3.84
      调节系数标准差[± x-fold]
      std dev [± x-fold]
      1.95 1.52 1.58 1.68 1.72
    • RefFinder 在线软件可以对上述 4 种算法分析所得到的结果进行综合排序,赋予每个候选内参基因一个适当的权重,然后计算出其几何平均值,通过赋分综合排名,得分越低越稳定,筛选出表现最稳定的内参基因。结果显示RefFinder综合以上四种算法计算出来的稳定值排序为TUB > ACT > GTPb > Nascent > EF1α,其中稳定值最大的是TUB,稳定值为5.00;稳定值最小的是EF1α,稳定值为1.41,说明稳定性最好的是EF1α基因,稳定性最差的是TUB基因(表5)。

      表 5  RefFinder分析结果

      Table 5.  Results of RefFinder Software

      基因
      Gene
      稳定值
      Stable value
      排名
      Ranking
      EF1α 1.41 1
      Nascent 1.73 2
      GTPb 2.21 3
      ACT 3.72 4
      TUB 5.00 5
    • 木薯种植和生产领域中,CBB是目前最为严重的木薯病害之一,根据种植环境条件的不同,引起的损失至高可达100%。因此,充分挖掘木薯抗病基因已成为木薯产业亟待解决的问题之一,而选择稳定的内参基因是木薯抗病基因挖掘的基础。木薯中稳定的内参基因筛选已有部分报道,如Moreno 等[19]研究了5对内参基因在木薯褐条病毒侵染木薯品种‘TMS60444’、‘Ebwanateraka’和‘AR34’叶片中的稳定性,发现GTPb基因在感染木薯褐条病的条件下表达最稳定。Hu等[18]研究了26对内参基因在木薯‘Ku50’、‘SC124’、‘Rongyong9’、‘SC5’和‘Arg7’等品种中不同发育时期的叶片、块根、叶柄及花以及干旱胁迫下的不同时间点叶片合计21个样本中的内参稳定性,发现Nascent基因在绝大多数情况下都稳定表达[18]。在哥伦比亚Xpm生理小种UA681和CFBP1851侵染感病木薯品种60444后受诱导表达的木薯基因功能研究中,Carlos[15]等使用TUB基因作为内参基因,我们所使用的Xpm生理小种 CHN 01 是在海南省分离得到的,将TUB作为候选内参基因研究Xpm 生理小种CHN 01侵染木薯下TUB在木薯中表达是否稳定[23]ACT基因被认为是在拟南芥抗病研究中最稳定的内参基因[20],本研究选择了序列相似度最高的木薯直系同源基因(Manes.13G086400,ACT)作为候选内参基因。EF1α基因也被用于CBB抗病基因筛选的内参基因使用[2122]

      目前,木薯黄单胞菌侵染后的木薯内参基因稳定性研究较少,因此选择了上述5对文献报道的木薯基因GTPb、Nascent、TUB、ACT、EF1α作为后续内参基因研究。本研究采用了Delta Ct、Genorm、NormFinder、BestKeeper以及在线软件Reffinder对所选用的5对内参基因在木薯抗感品种(‘JG1301’和‘SC9’)不同接种时间点进行了稳定性打分,进一步出筛选在Xpm侵染木薯的条件下仍然能稳定表达的内参基因。结果表明,qRT-PCR数据显示5个候选内参基因都有不同幅度的变化,研究发现,Nascent基因在木薯不受外界胁迫的情况下可以稳定表达,但是在干旱胁迫下表达不稳定[18] ,结合Nascent基因在Xpm侵染木薯叶片后在木薯叶片中稳定性排名第2,说明外界生物与非生物胁迫影响了Nascent基因的表达量。GTPb基因在Xpm侵染木薯叶片后表达稳定性综合排名第3,而在木薯褐条病毒侵染木薯条件下较为稳定[19],可能原因为褐条病毒侵染机理与Xpm作为维管束病害的侵染机制不同,从而造成了GTPb基因在不同生物胁迫条件下的表达量变化差异。

      Carlos在木薯抗病基因研究中使用了TUB作为内参基因,同时Wei等[24]在进行木薯耐盐研究实验中也使用TUB作为内参基因 。Tong等[25]研究发现,TUB基因在桃的不同组织中表达稳定 ,Zhong等[26]研究发现TUB基因在荔枝的不同组织和遮光处理下都能稳定表达,但是Nicot等[27]的研究表明,TUB在马铃薯侵染晚疫病和冷胁迫实验中表达量不稳定 。本研究结果表明,TUBXpm侵染木薯后稳定性最差,说明TUB在木薯受到病原菌胁迫时表达不稳定。

      ACT基因能在栾树的不同组织中稳定[28]ACT也能在百合花不同发育时期和不同组织中稳定表达[29],但Nicot等[27]研究发现,在马铃薯侵染晚疫病和盐处理下ACT基因在马铃薯叶片中表达最不稳定,Li等[30]在甘薯内参基因筛选实验中发现,在病毒感染和非病毒感染的条件下甘薯叶片中的ACT基因表达量不稳定。本研究发现,Xpm侵染木薯叶片后ACT基因在木薯叶片中表达不稳定,说明ACT基因虽然在多种植物中被广泛用作内参基因,但是不适合在部分薯类作物中当作内参基因。

      Zheng和Chang等[2122]在进行CBB研究时均使用EF1α作为内参基因,Nicot等[27]研究发现,在马铃薯晚疫病侵染下和冷胁迫下马铃薯叶片中表达最稳定的内参基因是EF1α。Mascia等[31]研究发现, EF1α基因在菊苣叶片和根组织中的表达稳定性最高。唐枝娟等[32]研究发现,强弱致病力的白叶枯病菌PXO99和PXO61侵染下的水稻EF1α基因表达最稳定。本研究结果也表明,在Xpm侵染下,木薯EF1α基因在抗感木薯不同时间点表达量最为稳定。综合上述结果可知,在多种植物和多种外界胁迫下,EF1α基因都能在植物叶片中稳定表达,可作为一个稳定的内参基因使用。我们的实验结果为木薯抗细菌性枯萎病相关抗病基因挖掘研究中稳定内参基因的使用奠定了理论基础

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回