留言板

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

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

‘无子瓯柑’CHS基因家族的克隆和表达分析

俞狄虎 张迟 柯甫志 敬露阳 顾雪娇 吴宝玉 张敏

俞狄虎, 张迟, 柯甫志, 敬露阳, 顾雪娇, 吴宝玉, 张敏. ‘无子瓯柑’CHS基因家族的克隆和表达分析[J]. 浙江农林大学学报, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
引用本文: 俞狄虎, 张迟, 柯甫志, 敬露阳, 顾雪娇, 吴宝玉, 张敏. ‘无子瓯柑’CHS基因家族的克隆和表达分析[J]. 浙江农林大学学报, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
YU Dihu, ZHANG Chi, KE Fuzhi, JING Luyang, GU Xuejiao, WU Baoyu, ZHANG Min. Cloning and expression analysis of CHS gene family in Citrus suavissima 'Seedless'[J]. Journal of Zhejiang A&F University, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
Citation: YU Dihu, ZHANG Chi, KE Fuzhi, JING Luyang, GU Xuejiao, WU Baoyu, ZHANG Min. Cloning and expression analysis of CHS gene family in Citrus suavissima 'Seedless'[J]. Journal of Zhejiang A&F University, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013

‘无子瓯柑’CHS基因家族的克隆和表达分析

doi: 10.11833/j.issn.2095-0756.2019.05.013
基金项目: 

浙江省基础公益研究计划项目 LGN18C160005

浙江省基础公益研究计划项目 LGN19C160005

国家自然科学基金资助项目 31000897

国家现代农业产业技术体系专项基金资助项目 CARS-27

详细信息
    作者简介: 俞狄虎, 从事经济林栽培与利用研究。E-mail:1175248409@qq.com
    通信作者: 张敏, 副教授, 博士, 从事经济林栽培与利用研究。E-mail:mzhang@zafu.edu.cn
  • 中图分类号: S666.1;S759.3

Cloning and expression analysis of CHS gene family in Citrus suavissima 'Seedless'

  • 摘要: 为研究查尔酮合成酶基因(CHS)家族在‘无子瓯柑’Citrus suavissima‘Seedless’雄性不育发生过程中的作用,以克里曼丁橘Citrus clementina基因组数据库为参照,在‘无子瓯柑’和瓯柑Citrus suavissima的转录组和蛋白质组测序结果中筛选出4个CHS同源差异表达基因,并对其克隆和表达量分析,对克里曼丁橘CHS基因家族成员进行生物信息学分析。结果表明:‘无子瓯柑’和瓯柑CHS基因编码区核苷酸序列相似度达98%以上。小孢子母细胞发育各时期CHS基因表达量在瓯柑和‘无子瓯柑’间差异显著。与瓯柑相比,‘无子瓯柑’小孢子母细胞时期Ciclev10005133mCiclev10030093mCiclev10015535m的表达量显著下调;减数分裂时期Ciclev10005133mCiclev10001405mCiclev10015535m的表达量显著下调;四分体时期Ciclev10001405mCiclev10030093m的表达量显著减少。在花粉粒成熟时期的花蕾中,Ciclev10001405mCiclev10005133m在花药中特异性表达。CHS同源基因在花药不同时期中的差异表达可能是‘无子瓯柑’花药发育异常的重要原因,并最终导致‘无子瓯柑’雄性不育。
  • 图  1  ‘无子瓯柑’和瓯柑CHS基因在花药发育过程中的基因表达分析

    Figure  1  Expression level of anther in different development stages in C. suavissima 'Seedless' and C. suavissima

    图  2  ‘无子瓯柑’和瓯柑花蕾不同部位的CHS基因的相对表达量

    Figure  2  Expression of CHS genes among different floral organs in C. suavissima 'seedless' and C. suavissima

    图  3  4个物种CHS蛋白系统进化树

    Figure  3  Phyogenetic trees of CHS proteins from four species

    表  1  ‘无子瓯柑’ CHS基因家族成员及其小孢子母细胞时期相对表达量

    Table  1.   CHS gene family and expression in C. suavissima 'Seedless' at microsporcyte

    基因名称 差异表达倍数 错误发现率/假设概率
    Ciclev10015535m 0.76 0.004 6b
    Ciclev10005133m 0.65 0.007 5a
    Ciclev10001405m 1.32 0.006 3a
    Ciclev10001413m 5.50 0.595 4a
    Ciclev10001395m 0.94 0.717 6a
    Ciclev10028604m 5.50 0.161 3a
    Ciclev10030093m 0.71 0.007 9a
    Ciclev10028605m 0.91 0.144 1a
    Ciclev10030398m 5.50 0.399 2a
    Ciclev10001905m 0.90 0.064 2a
    说明: a表示错误发现率; b表示假设概率
    下载: 导出CSV

    表  2  基因克隆引物

    Table  2.   Sequences of primer for gene cloning

    基因名称 上游引物(5′→3′) 下游引物(5′→3′)
    Ciclev10005133m ATGGTGACCGTCGATGAAG CAGTGTTGCCGCTGCTTAA
    Ciclev10001405m ATGGAGAAAGTTAAAGATG TCCCTACTGTTACAACCTAG
    Ciclev10030093m ATGACGACAGTGAAAAGTAA ACAAATTTCACCAACTACTGA
    Ciclev10015535m ATGGCAACCGTTCAAGAGAT GTGTCCCCATCAAAGCTTGA
    下载: 导出CSV

    表  3  qRT-PCR引物

    Table  3.   Sequences of primer for real-time PCR

    基因ID 上游引物(5′→3′) 下游引物(5′→3′)
    GU911361 ATCTGCTGGAAGGTGCTGAG CCAAGCAGCATGAAGATCAA
    Ciclev10015535m AAGAGCGAGCATATGACGGA CAGCTTCTTTCCCGAGCTTC
    Ciclev10001405m AATTGTGTCGGGTGCACAAA ATTCCGAACGGACTAAACGC
    Ciclev10005133m AGCCGAGAACAACAAAGG ATGGGCTTCTCGATCTCAGG
    Ciclev10030093m CAAGGACCAACAGCAACGAT CTTTCAGCTCGGTCTTGTGG
    下载: 导出CSV

    表  4  瓯柑和‘无子瓯柑’的CHS同源基因

    Table  4.   Alignments of CsCHS nucleotide sequences between Citrus suavissima 'Seedless' and C. suavissima

    基因名称 瓯柑 ‘无子瓯柑’ 开放阅读框长度/bp 同源基因在瓯柑与‘无子瓯柑’间的相似度/% 差异碱基位/bp 氨基酸的变化
    命名 GenBank
    登录号
    命名 GenBank
    登录号
    Ciclev10015535m CsCHS1 MK070533 CsCHS5 MK070537 1 176 99 501
    Ciclev10001405m CsCHS6 MK070536 CsCHS6 MK070536 1 188 100
    Ciclev10005133m CsCHS3 MK070534 CsCHS7 MK070538 1 173 98 549~954 922 bp处苏氨酸突变为丙氨酸
    Ciclev10030093m CsCHS4 MK070535 CsCHS8 MK070539 1 179 98 342,495
    下载: 导出CSV

    表  5  CHS基因及其表达

    Table  5.   Basic information of CHS genes and proteins

    基因名称 氨基酸数量/个 分子量/Da 等电点
    Ciclev10015535m 391 42 592.17 6.47
    Ciclev10005133m 390 42 676.30 5.84
    Ciclev10015900m 328 35 477.11 6.11
    Ciclev10001405m 395 42 918.26 6.95
    Ciclev10001413m 394 43 364.92 6.56
    Ciclev10001395m 396 43 377.16 6.77
    Ciclev10028604m 395 43 383.97 5.48
    Ciclev10030093m 392 43 406.82 5.82
    Ciclev10028605m 395 43 436.08 6.17
    Ciclev10030398m 390 429 51.35 6.10
    Ciclev10001905m 312 33 958.16 5.44
    Ciclev10001806m 330 36 069.43 6.33
    Ciclev10003127m 309 33 988.20 5.59
    下载: 导出CSV

    表  6  CHS蛋白二级结构分析

    Table  6.   Secondary structure analysis of CHS proteins

    蛋白名称 α螺旋 延伸链 β转角 无规则卷曲
    数量/个 比例/% 数量/个 比例/% 数量/个 比例/% 数量/个 比例/%
    Ciclev10015535m 161 41.18 65 16.62 32 8.18 133 34.02
    Ciclev10005133m 170 43.59 60 15.38 31 7.95 129 33.08
    Ciclev10015900m 142 43.29 59 17.99 29 8.84 98 29.88
    Ciclev10001405m 160 40.51 63 15.95 28 7.09 144 36.46
    Ciclev10001413m 170 43.15 63 15.99 27 6.85 134 34.01
    Ciclev10001395m 175 44.19 62 15.66 26 6.57 133 33.59
    Ciclev10028604m 169 42.78 65 16.46 25 6.33 136 34.43
    Ciclev10030093m 156 39.80 64 16.33 31 7.91 141 35.97
    Ciclev10028605m 161 40.76 66 16.71 32 8.10 136 34.43
    Ciclev10030398m 165 42.31 61 15.64 28 7.18 136 34.87
    Ciclev10001905m 131 41.99 56 17.95 26 8.33 99 31.73
    Ciclev10001806m 133 40.30 53 16.06 22 6.67 122 36.97
    Ciclev10003127m 119 38.51 57 18.45 25 8.09 108 34.95
    下载: 导出CSV
  • [1] 曹庆芹, 伊华林, 邓秀新.果树雄性不育研究进展[J].果树学报, 2005, 22(6):678-681.

    CAO Qingqin, YI Hualin, DENG Xiuxin. Advances in research on male sterility in fruit crops[J]. J Fruit Sci, 2005, 22(6):678-681.
    [2] HU Zhiyong, ZHANG Min, WEN Qigen, et al. Abnormal microspore development leads to pollen abortion in a seedless mutant of 'Ougan' mandarin (Citrus suavissima Hort. ex Tanaka)[J]. J Amer Soc Hort Sci, 2007, 132(6):777-782.
    [3] 张迟, 张敏, 朱铨, 等. '瓯柑'及其无子突变体花粉发育的细胞学观察[J].果树学报, 2014, 31(2):265-269.

    ZHANG Chi, ZHANG Min, ZHU Quan, et al. Cytological observation of pollen development in 'Ougan' (Citrus suavissima Hort. ex Tanaka) and its seedless mutant[J]. J Fruit Sci, 2014, 31(2):265-269.
    [4] NAPOLI C A, FAHY D, WANG Huayu, et al. White anther:a petunia mutant that abolishes pollen flavonol accumulation, induces male sterility, and is complemented by a chalcone synthase transgene[J]. Plant Physiol, 1999, 120(2):615-622.
    [5] ZHENG Rui, YUE Sijun, XU Xiaoyan, et al. Proteome analysis of the wild and YX-1 male sterile mutant anthers of wolfberry (Lycium barbarum L.)[J]. PLoS One, 2012, 7(7):e41861.
    [6] 陈遥, 彭彦, 刘峰, 等.辣椒CHS基因表达的定量分析及与雄性不育的相关性[J].湖南农业大学学报(自然科学版), 2013, 39(3):259-264.

    CHEN Yao, PENG Yan, LIU Feng, et al. The quantitative analysis of CHS in Capsicum annuum L. and its correlation to male sterility[J]. J Hunan Agric Univ Nat Sci, 2013, 39(3):259-264.
    [7] XIE Yang, ZHANG Wei, WANG Yan, et al. Comprehensive transcriptome-based characterization of differentially expressed genes involved in microsporogenesis of radish CMS line and its maintainer[J]. Funct Integr Genomic, 2016, 16(5):529-543.
    [8] 王志彬, 申晚霞, 朱世平, 等.柑橘CHS基因序列多态性及表达水平对类黄酮生物合成的影响[J].园艺学报, 2015, 42(3):435-444.

    WANG Zhibin, SHEN Wanxia, ZHU Shiping, et al. Polymorphism and expression of chalcone synthase gene in citrus related to the flavonoids content[J]. Acta Hortic sin, 2015, 42(3):435-444.
    [9] 朱咪咪, 张迟, 常爱玲, 等. '无籽'瓯柑小孢子母细胞减数分裂特性基因RAD51和MS1的表达差异分析[J].浙江农林大学学报, 2016, 33(6):921-927.

    ZHU Mimi, ZHANG Chi, CHANG Ailing, et al. Expression of microsporocyte meiosis with special genes RAD51 and MS1 in Citrus suavissima 'seedless'[J]. J Zhejiang A&F Univ, 2016, 33(6):921-927.
    [10] QIU Lingling, JIANG Bo, FANG Jia, et al. Analysis of transcriptome in hickory (Carya cathayensis), and uncover the dynamics in the hormonal signaling pathway during graft process[J]. BMC Genomics, 2017, 17(1):935.
    [11] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method[J]. Methods, 2001, 25(4):402-408.
    [12] YU Chinsheng, LIN Chinjen, HWANG Jennkang. Predicting subcellular localization of proteins for gram-negative bacteria by support vector machines based on n-peptide compositions[J]. Protein Sci, 2004, 13(5):1402-1406.
    [13] ZHANG Yi, QU Lijia, LIU Meihua, et al. Structural and expressional analysis of a cDNA that expresses predominantly in rice stamens[J]. Chin Sci Bull, 1998, 43(9):765-769.
    [14] AKADA S, DUBE S K. Organization of soybean chalcone synthase gene clusters and characterization of a new member of the family[J]. Plant Mol Biol, 1995, 29(2):189-199.
    [15] ZHENG Honghong, QU Lijia, LIU Meijia, et al. An anther-specific chalcone synthase-like gene D5 related to rice pollen development[J]. Chin Sci Bull, 2000, 45(21):1921-1926.
  • [1] 张耀, 王家璇, 蔡璇, 曾祥玲, 杨洁, 陈洪国, 邹晶晶.  桂花OfACOs基因家族鉴定及表达分析 . 浙江农林大学学报, 2023, 40(3): 492-501. doi: 10.11833/j.issn.2095-0756.20220783
    [2] 叶潇铃, 赵宇虹, 姜楠, 张敏, 张迟.  ‘无子瓯柑’E3泛素连接酶基因CsRNF217对转基因烟草育性的影响 . 浙江农林大学学报, 2023, 40(6): 1181-1187. doi: 10.11833/j.issn.2095-0756.20220710
    [3] 赖梦霞, 杜长霞, 樊怀福.  黄瓜AQP基因家族的鉴定与生物信息学分析 . 浙江农林大学学报, 2022, 39(2): 318-328. doi: 10.11833/j.issn.2095-0756.20210361
    [4] 火艳, 招雪晴, 黄厚毅, 黄贤斌, 许云方, 祝遵凌, 苑兆和.  观赏石榴表型遗传多样性分析 . 浙江农林大学学报, 2020, 37(5): 939-949. doi: 10.11833/j.issn.2095-0756.20190619
    [5] 洪森荣, 曾清华, 谭鑫, 陈永华, 郑亚娇, 徐迎昕, 邱梦琴.  上饶早梨‘六月雪’和‘黄皮消’全基因组重测序分析 . 浙江农林大学学报, 2019, 36(2): 227-235. doi: 10.11833/j.issn.2095-0756.2019.02.003
    [6] 钱武兵, 李建设, 高艳明, 周文波.  营养液中添加不同盐类对水培番茄果实糖组分和相关酶活性的影响 . 浙江农林大学学报, 2018, 35(6): 1120-1127. doi: 10.11833/j.issn.2095-0756.2018.06.016
    [7] 龚小强, 李素艳, 李燕, 孙向阳.  绿化废弃物好氧堆肥和蚯蚓堆肥作为蔬菜育苗基质研究 . 浙江农林大学学报, 2016, 33(2): 280-287. doi: 10.11833/j.issn.2095-0756.2016.02.013
    [8] 朱咪咪, 张迟, 常爱玲, 党婉誉, 周彩红, 俞狄虎, 吴莹莹, 张敏.  ‘无籽’瓯柑小孢子母细胞减数分裂特性基因RAD51和MS1的表达差异分析 . 浙江农林大学学报, 2016, 33(6): 921-927. doi: 10.11833/j.issn.2095-0756.2016.06.001
    [9] 王文鹏, 周莉花, 赵宏波, 包志毅.  蜡梅与光叶红蜡梅和夏蜡梅属间杂交亲和性初步研究 . 浙江农林大学学报, 2015, 32(5): 756-762. doi: 10.11833/j.issn.2095-0756.2015.05.015
    [10] 雷燕, 李庆卫, 李文广, 景珊, 陈俊愉.  2个地被菊品种对不同遮光处理的生理适应性 . 浙江农林大学学报, 2015, 32(5): 708-715. doi: 10.11833/j.issn.2095-0756.2015.05.008
    [11] 张敏, 刘志辉, 宋雪恩, 朱咪咪, 党婉誉, 颜福花, 张迟.  ‘无籽’瓯柑CsLTPCsLOX基因的克隆与表达分析 . 浙江农林大学学报, 2014, 31(6): 823-830. doi: 10.11833/j.issn.2095-0756.2014.06.001
    [12] 王国立, 安华明, 秦巧平, 李孟娇, 刘真真, 陈佳莹, 周倩, 张岚岚.  柑橘果实成熟特异基因CsPMEI/InvI的克隆与序列分析 . 浙江农林大学学报, 2013, 30(3): 336-342. doi: 10.11833/j.issn.2095-0756.2013.03.005
    [13] 秦巧平, 林飞凡, 张岚岚.  枇杷果实糖酸积累的分子生理机制 . 浙江农林大学学报, 2012, 29(3): 453-457. doi: 10.11833/j.issn.2095-0756.2012.03.021
    [14] 孟赐福, 曹志洪, 姜培坤, 徐秋芳, 周国模.  叶面施硼矫正杨梅缺硼的研究 . 浙江农林大学学报, 2008, 25(5): 543-547.
    [15] 余学军, 窦可, 章兆福, 柳德坚.  不同保鲜预处理对绿竹笋呼吸速率的影响 . 浙江农林大学学报, 2007, 24(4): 424-427.
    [16] 周建, 袁德义, 张琳, 李佑武.  黄金梨叶片营养诊断及施肥标准 . 浙江农林大学学报, 2007, 24(1): 39-43.
    [17] 曾艳玲, 谭晓风, 张党权, 曾晓峰, 李秀根, 刘先雄.  5个中国砂梨品种S基因型的确定 . 浙江农林大学学报, 2007, 24(6): 654-660.
    [18] 胡滨, 陈观平, 施农农, 卢江杰, 王慧中.  植物细胞质雄性不育及其在育种中的应用 . 浙江农林大学学报, 2006, 23(6): 689-693.
    [19] 王齐瑞, 谭晓风, 张琳.  覆草栽培对甜樱桃生长及光合速率的影响 . 浙江农林大学学报, 2006, 23(1): 24-28.
    [20] 何福基, 余象煌, 李平.  杉木雄性不育性遗传机理初步分析* . 浙江农林大学学报, 1995, 12(2): 219-220.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2019.05.013

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2019/5/943

图(3) / 表(6)
计量
  • 文章访问数:  3063
  • HTML全文浏览量:  658
  • PDF下载量:  62
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-10
  • 修回日期:  2019-01-23
  • 刊出日期:  2019-10-20

‘无子瓯柑’CHS基因家族的克隆和表达分析

doi: 10.11833/j.issn.2095-0756.2019.05.013
    基金项目:

    浙江省基础公益研究计划项目 LGN18C160005

    浙江省基础公益研究计划项目 LGN19C160005

    国家自然科学基金资助项目 31000897

    国家现代农业产业技术体系专项基金资助项目 CARS-27

    作者简介:

    俞狄虎, 从事经济林栽培与利用研究。E-mail:1175248409@qq.com

    通信作者: 张敏, 副教授, 博士, 从事经济林栽培与利用研究。E-mail:mzhang@zafu.edu.cn
  • 中图分类号: S666.1;S759.3

摘要: 为研究查尔酮合成酶基因(CHS)家族在‘无子瓯柑’Citrus suavissima‘Seedless’雄性不育发生过程中的作用,以克里曼丁橘Citrus clementina基因组数据库为参照,在‘无子瓯柑’和瓯柑Citrus suavissima的转录组和蛋白质组测序结果中筛选出4个CHS同源差异表达基因,并对其克隆和表达量分析,对克里曼丁橘CHS基因家族成员进行生物信息学分析。结果表明:‘无子瓯柑’和瓯柑CHS基因编码区核苷酸序列相似度达98%以上。小孢子母细胞发育各时期CHS基因表达量在瓯柑和‘无子瓯柑’间差异显著。与瓯柑相比,‘无子瓯柑’小孢子母细胞时期Ciclev10005133mCiclev10030093mCiclev10015535m的表达量显著下调;减数分裂时期Ciclev10005133mCiclev10001405mCiclev10015535m的表达量显著下调;四分体时期Ciclev10001405mCiclev10030093m的表达量显著减少。在花粉粒成熟时期的花蕾中,Ciclev10001405mCiclev10005133m在花药中特异性表达。CHS同源基因在花药不同时期中的差异表达可能是‘无子瓯柑’花药发育异常的重要原因,并最终导致‘无子瓯柑’雄性不育。

English Abstract

俞狄虎, 张迟, 柯甫志, 敬露阳, 顾雪娇, 吴宝玉, 张敏. ‘无子瓯柑’CHS基因家族的克隆和表达分析[J]. 浙江农林大学学报, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
引用本文: 俞狄虎, 张迟, 柯甫志, 敬露阳, 顾雪娇, 吴宝玉, 张敏. ‘无子瓯柑’CHS基因家族的克隆和表达分析[J]. 浙江农林大学学报, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
YU Dihu, ZHANG Chi, KE Fuzhi, JING Luyang, GU Xuejiao, WU Baoyu, ZHANG Min. Cloning and expression analysis of CHS gene family in Citrus suavissima 'Seedless'[J]. Journal of Zhejiang A&F University, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
Citation: YU Dihu, ZHANG Chi, KE Fuzhi, JING Luyang, GU Xuejiao, WU Baoyu, ZHANG Min. Cloning and expression analysis of CHS gene family in Citrus suavissima 'Seedless'[J]. Journal of Zhejiang A&F University, 2019, 36(5): 943-949. doi: 10.11833/j.issn.2095-0756.2019.05.013
  • 雄性不育在开花植物中普遍存在,主要表现为雄蕊发育不正常,不能产生具有正常功能的花粉[1]。‘无子瓯柑’Citrus suavissima ‘Seedless’是瓯柑Citrus suavissima的芽变品种,保留了瓯柑肉质饱满、香味特异、耐储藏的优良品质,因果实无核被人们青睐。研究认为,雄性不育是‘无子瓯柑’无核的重要原因之一[2]。张迟等[3]发现‘无子瓯柑’花粉败育始于小孢子母细胞时期,推测其雄性不育与能量代谢异常和营养物质缺乏有关。对矮牵牛Petunia hybrid[4],枸杞Lycium barbarum[5],辣椒Capscum annuum[6]和萝卜Raphanus sativus[7]等植物的研究发现,查尔酮合成酶基因(CHS)的异常表达会引起雄性不育;CHS是黄酮类化合物代谢通路中的关键基因,通过影响查尔酮的产生从而影响黄酮类化合物的生物合成[4],而黄酮类化合物的缺乏/过量是植物雄性不育的重要原因之一。有研究将克隆自10个柑橘种质的CHS基因与温州蜜柑‘国庆4号’Citrus unshiu ‘Guoqing 4’,柚‘冯威’Citrus maxima ‘Fengwei’和甜橙‘红宝石’Citrus sinensis ‘Ruby’的CHS比对,发现不同品种柑橘的CHS基因编码区核苷酸序列相似度极高,达98%以上[8]。本研究以前期工作为基础,以克里曼丁橘Citrus clementina数据库为参照,对‘无子瓯柑’和瓯柑小孢子母细胞时期的花药进行转录组和蛋白质组测序,筛选CHS同源差异表达基因进行克隆和表达量分析,并对CHS基因家族进行生物信息学分析,以期为‘无子瓯柑’雄性不育机理的深入研究打下基础。

    • ‘无子瓯柑’和瓯柑小孢子母细胞时期的花药转录组测序数据源自课题组上传的NCBI数据(NCBI登录号PRJNA430695)。

    • ‘无子瓯柑’和瓯柑采自浙江省丽水市富岭街道南寨自然村,采样时间参照前期研究[9]。采集‘无子瓯柑’和瓯柑成熟时期的花蕾进行基因克隆。采集两者小孢子母细胞时期(Ⅰ,花蕾直径2.0~2.4 mm)、减数分裂时期(Ⅱ,花蕾直径2.4~2.8 mm)和四分体时期(Ⅲ,花蕾直径2.8~3.1 mm)的花药,用于分析不同发育时期的基因表达量。采集两者成熟花粉粒时期花蕾(花蕾直径6.5~6.9 mm)的4个不同部位(花药、花丝、雌蕊和花瓣),用于分析不同部位的基因表达量。所有材料均保存于-80 ℃。

    • 利用RNA提取试剂盒(TaKaRa MiniBEST Plant RNA Extraction Kit,TaKaRa)提取瓯柑、‘无子瓯柑’不同发育时期的花药和成熟花粉粒时期花蕾不同部位的RNA,并对RNA样品质量进行检测。按反转录试剂盒(PrimeScript RT reagent Kit with gDNA Eraser,TaKaRa)说明书,将得到的RNA反转录成cDNA,并于-20 ℃保存。

      以克里曼丁橘基因组数据库为参考,对小孢子母细胞时期‘无子瓯柑’和瓯柑的花药进行转录组和蛋白质组分析,鉴定得到CHS基因家族成员;以差异表达倍数(FC)>1.2,错误发现率(FDR)<0.01为标准[10],筛选到CHS同源差异表达基因;以差异表达倍数(FC)>1.2,假设概率(P)<0.05为标准,筛选CHS同源差异表达蛋白。

      对得到的差异表达CHS基因进行克隆及定量引物设计。对‘无子瓯柑’和瓯柑的CHS基因的蛋白质编码区(CDS, Coding Sequence)进行克隆。聚合酶链式反应(PCR, Polymerase Chain Reaction)程序为95 ℃,5 min;95 ℃,30 s,56 ℃,40 s,38个循环;72 ℃,10 min。PCR产物在质量分数为1.5%的琼脂糖凝胶上电泳。利用胶回收试剂盒(TaKaRa)纯化胶回收产物并连接到pMD-18质粒(TaKaRa),转化感受态大肠埃希菌Escherichia coli DH5α(TaKaRa)中,37 ℃震荡培养后涂板挑菌。通过菌液PCR及电泳检测后,选出条带大小正确的菌液送生工生物工程(上海)有限公司测序,获得基因的核苷酸序列。

      以克里曼丁橘序列为模板,在线设计qRT-PCR特异性引物(https://www.genscript.com),以甜橙Citrus sinensisActin基因(GU911361)[9]为内参基因,引物序列由生工生物工程(上海)有限公司合成。利用CFX96 real-time system实时荧光定量PCR仪(Bio-Rad)及荧光染料SYBR® Premix Ex TaqTMⅡ(TaKaRa)说明书进行实时定量PCR,反应程序为95 ℃,30 s,95 ℃,5 s,57 ℃,30 s,39个循环,65~95 ℃升温检测扩增产物的溶解曲线[9]。设置3个重复,根据2-ΔΔCt方法,以瓯柑小孢子母细胞时期花药的表达量为基准计算基因相对表达量[11],比较‘无子瓯柑’和瓯柑在花药不同发育时期和花粉粒成熟期花蕾不同部位的基因表达量。

    • 柑橘CHS基因编码区核苷酸序列相似度极高[8],因此以克里曼丁橘为例开展CHS基因及蛋白质的生物信息学分析。登录克里曼丁橘数据库,下载CHS基因全长序列、基因CDS序列、家族蛋白序列等信息。所有蛋白质序列通过在线数据库ExPASY(https://web.expasy.org/protparam)进行蛋白生理生化分析。基于已有报道[12],利用CELLOv2.5进行亚细胞定位。通过在线数据库Prabi(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)进行蛋白二级结构分析。利用MEGA6.0软件采用邻接法(Neighbor-Joining)构建克里曼丁橘(https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Cclemen-tina)、甜橙(https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Csinensis)、拟南芥Arabidopsis thalianahttps://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Athaliana)和水稻Oryza satavahttps://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Osatava)的CHS同源蛋白系统进化树,设置Bootstrap值为1 000,去除支持率低于50%的节点,并显示各分支长度。

    • 对小孢子母细胞时期的‘无子瓯柑’和瓯柑花药转录组和蛋白质组数据分析,共得到10个CHS同源基因(表 1);转录组测序结果满足差异表达倍数>1.2且错误发现率<0.01的有3个基因(Ciclev10005133mCiclev10001405mCiclev10030093m),蛋白质组测序结果满足差异表达倍数>1.2且假设概率<0.05标准的有1个基因(Ciclev10015535m)。得到的4个差异表达CHS基因引物序列如表 2。qRT-PCR特异性引物序列如表 3

      表 1  ‘无子瓯柑’ CHS基因家族成员及其小孢子母细胞时期相对表达量

      Table 1.  CHS gene family and expression in C. suavissima 'Seedless' at microsporcyte

      基因名称 差异表达倍数 错误发现率/假设概率
      Ciclev10015535m 0.76 0.004 6b
      Ciclev10005133m 0.65 0.007 5a
      Ciclev10001405m 1.32 0.006 3a
      Ciclev10001413m 5.50 0.595 4a
      Ciclev10001395m 0.94 0.717 6a
      Ciclev10028604m 5.50 0.161 3a
      Ciclev10030093m 0.71 0.007 9a
      Ciclev10028605m 0.91 0.144 1a
      Ciclev10030398m 5.50 0.399 2a
      Ciclev10001905m 0.90 0.064 2a
      说明: a表示错误发现率; b表示假设概率

      表 2  基因克隆引物

      Table 2.  Sequences of primer for gene cloning

      基因名称 上游引物(5′→3′) 下游引物(5′→3′)
      Ciclev10005133m ATGGTGACCGTCGATGAAG CAGTGTTGCCGCTGCTTAA
      Ciclev10001405m ATGGAGAAAGTTAAAGATG TCCCTACTGTTACAACCTAG
      Ciclev10030093m ATGACGACAGTGAAAAGTAA ACAAATTTCACCAACTACTGA
      Ciclev10015535m ATGGCAACCGTTCAAGAGAT GTGTCCCCATCAAAGCTTGA

      表 3  qRT-PCR引物

      Table 3.  Sequences of primer for real-time PCR

      基因ID 上游引物(5′→3′) 下游引物(5′→3′)
      GU911361 ATCTGCTGGAAGGTGCTGAG CCAAGCAGCATGAAGATCAA
      Ciclev10015535m AAGAGCGAGCATATGACGGA CAGCTTCTTTCCCGAGCTTC
      Ciclev10001405m AATTGTGTCGGGTGCACAAA ATTCCGAACGGACTAAACGC
      Ciclev10005133m AGCCGAGAACAACAAAGG ATGGGCTTCTCGATCTCAGG
      Ciclev10030093m CAAGGACCAACAGCAACGAT CTTTCAGCTCGGTCTTGTGG

      对克隆得到的‘无子瓯柑’和瓯柑的CHS同源基因序列(CsCHS)进行比对,结果发现:Ciclev10015535mCiclev10005133mCiclev10030093m在瓯柑和‘无子瓯柑’间存在核苷酸序列变异,但仅Ciclev10030093m在编码氨基酸水平发生了改变(表 4)。

      表 4  瓯柑和‘无子瓯柑’的CHS同源基因

      Table 4.  Alignments of CsCHS nucleotide sequences between Citrus suavissima 'Seedless' and C. suavissima

      基因名称 瓯柑 ‘无子瓯柑’ 开放阅读框长度/bp 同源基因在瓯柑与‘无子瓯柑’间的相似度/% 差异碱基位/bp 氨基酸的变化
      命名 GenBank
      登录号
      命名 GenBank
      登录号
      Ciclev10015535m CsCHS1 MK070533 CsCHS5 MK070537 1 176 99 501
      Ciclev10001405m CsCHS6 MK070536 CsCHS6 MK070536 1 188 100
      Ciclev10005133m CsCHS3 MK070534 CsCHS7 MK070538 1 173 98 549~954 922 bp处苏氨酸突变为丙氨酸
      Ciclev10030093m CsCHS4 MK070535 CsCHS8 MK070539 1 179 98 342,495
    • 图 1表明:小孢子母细胞时期,‘无子瓯柑’与瓯柑相比Ciclev10005133mCiclev10030093mCiclev10015535m显著下调,Ciclev10001405m显著上调;减数分裂时期,Ciclev10005133mCiclev10001405mCiclev10015535m显著下调;四分体时期,Ciclev10001405mCiclev10030093m显著下调,Ciclev10005133m显著上调。

      图  1  ‘无子瓯柑’和瓯柑CHS基因在花药发育过程中的基因表达分析

      Figure 1.  Expression level of anther in different development stages in C. suavissima 'Seedless' and C. suavissima

      图 2可知:成熟花粉粒时期,‘无子瓯柑’和瓯柑的花丝中CHS基因表达量相仿。花药中Ciclev10001405mCiclev10005133m表达量显著高于其他花器官;Ciclev10030093mCiclev10005133mCiclev10015535m的表达在‘无子瓯柑’和瓯柑的花药中存在显著差异,可能会引起该时期‘无子瓯柑’和瓯柑的花药中黄酮含量的差异。总的来说,花药中Ciclev10001405m表达量远高于花瓣、雌蕊和花丝,推测Ciclev10001405m的表达在花药中具有特异性。

      图  2  ‘无子瓯柑’和瓯柑花蕾不同部位的CHS基因的相对表达量

      Figure 2.  Expression of CHS genes among different floral organs in C. suavissima 'seedless' and C. suavissima

    • 克里曼丁橘基因组数据库共鉴定到13个CHS基因家族成员。对CHS蛋白质理化性质的分析表明:编码的氨基酸数量最少的是Ciclev10003127m(309个),最多的是Ciclev10001395m(396个);平均等电点为6.12,所有蛋白质均为酸性蛋白。相比之下,‘无子瓯柑’和瓯柑中,Ciclev10005133m(390个)编码的氨基酸数量最少,Ciclev10001405m(395个)编码的氨基酸数量最多;平均等电点为6.27(表 5)。亚细胞定位结果表明:这些CHS蛋白都分布在细胞质中,说明CHS蛋白在细胞中的分布具有特异性。蛋白质二级结构分析表明:α-螺旋比例和数量均最高,β-转角则相对较少;表现为α-螺旋数量>无规则卷曲数量>延伸链数量>β-转角数量。‘无子瓯柑’与瓯柑中筛选的4个CHS蛋白也表现出相似趋势(表 6),推测α-螺旋在蛋白质二级结构中起主导作用。

      表 5  CHS基因及其表达

      Table 5.  Basic information of CHS genes and proteins

      基因名称 氨基酸数量/个 分子量/Da 等电点
      Ciclev10015535m 391 42 592.17 6.47
      Ciclev10005133m 390 42 676.30 5.84
      Ciclev10015900m 328 35 477.11 6.11
      Ciclev10001405m 395 42 918.26 6.95
      Ciclev10001413m 394 43 364.92 6.56
      Ciclev10001395m 396 43 377.16 6.77
      Ciclev10028604m 395 43 383.97 5.48
      Ciclev10030093m 392 43 406.82 5.82
      Ciclev10028605m 395 43 436.08 6.17
      Ciclev10030398m 390 429 51.35 6.10
      Ciclev10001905m 312 33 958.16 5.44
      Ciclev10001806m 330 36 069.43 6.33
      Ciclev10003127m 309 33 988.20 5.59

      表 6  CHS蛋白二级结构分析

      Table 6.  Secondary structure analysis of CHS proteins

      蛋白名称 α螺旋 延伸链 β转角 无规则卷曲
      数量/个 比例/% 数量/个 比例/% 数量/个 比例/% 数量/个 比例/%
      Ciclev10015535m 161 41.18 65 16.62 32 8.18 133 34.02
      Ciclev10005133m 170 43.59 60 15.38 31 7.95 129 33.08
      Ciclev10015900m 142 43.29 59 17.99 29 8.84 98 29.88
      Ciclev10001405m 160 40.51 63 15.95 28 7.09 144 36.46
      Ciclev10001413m 170 43.15 63 15.99 27 6.85 134 34.01
      Ciclev10001395m 175 44.19 62 15.66 26 6.57 133 33.59
      Ciclev10028604m 169 42.78 65 16.46 25 6.33 136 34.43
      Ciclev10030093m 156 39.80 64 16.33 31 7.91 141 35.97
      Ciclev10028605m 161 40.76 66 16.71 32 8.10 136 34.43
      Ciclev10030398m 165 42.31 61 15.64 28 7.18 136 34.87
      Ciclev10001905m 131 41.99 56 17.95 26 8.33 99 31.73
      Ciclev10001806m 133 40.30 53 16.06 22 6.67 122 36.97
      Ciclev10003127m 119 38.51 57 18.45 25 8.09 108 34.95
    • 现有研究筛选到克里曼丁橘13个CHS蛋白、甜橙7个CHS蛋白、拟南芥4个CHS蛋白和水稻3个CHS蛋白。按照遗传距离可以将CHS蛋白聚类为5个亚家族(图 3)。Group Ⅰ包含3个拟南芥CHS蛋白和2个水稻CHS蛋白。Group Ⅱ包含5个克里曼丁橘CHS蛋白和3个甜橙CHS蛋白。Group Ⅲ包含2个克里曼丁橘CHS蛋白。Group Ⅳ包含3个克里曼丁橘CHS蛋白和3个甜橙CHS蛋白。Group Ⅴ包含上述4个物种的CHS蛋白成员。其中Ciclev10005133m蛋白和Ciclev10015535m蛋白均聚类于Group Ⅴ,Ciclev10030093m蛋白和Ciclev10001405m蛋白则分别聚类于Group Ⅱ和Group Ⅳ。

      图  3  4个物种CHS蛋白系统进化树

      Figure 3.  Phyogenetic trees of CHS proteins from four species

    • 查尔酮合成酶(CHS)催化丙二酰-CoA和香豆酰-CoA结合形成查尔酮,是黄酮类化合物代谢通路(ko00941)中的第1个限速酶;查尔酮是各种黄酮类化合物的基本骨架,其数量异动会影响黄酮类化合物的形成,进而造成花药颜色、形态异常[4-5]。一般认为黄酮类化合物是柑橘具有高抗氧化性和良好药用功效的原因[8];矮牵牛CHS突变后花药颜色由黄色变成白色,花药功能异常,引起雄性不育[4];枸杞CHS的下调表达会影响花药功能,造成枸杞雄性不育[5];水稻中CHS的过量表达会引起花药表面色素大量沉淀,造成水稻雄性不育[13]

      CHS基因家族中成员不多,不同植物中数量也不同。如:拟南芥基因组数据库中仅鉴定到4个同源基因,大豆Glycine max[14]中略多(8个)。对已筛选的CHS同源基因的研究发现[13]X89859.1基因与水稻绒毡层的发育极其紧密,而后者发育正常与否直接关系着花粉的育性[15]。本研究在CHS蛋白系统进化树分析中发现,柑橘CHS同源蛋白Ciclev10005133m和水稻CHS同源蛋白X89895.1聚类于Group Ⅴ,推测两者具有相似的功能。在‘无子瓯柑’中,Ciclev10005133m同源基因在小孢子母细胞至减数分裂时期的表达显著低于瓯柑,因此,Ciclev10005133m同源基因的异常表达可能对‘无子瓯柑’小孢子母细胞发育过程产生不利影响。此外,Ciclev10015535mCiclev10030093m同源基因的表达也在小孢子母细胞时期出现显著下调,因此,CHS同源基因起始于小孢子母细胞时期的下调表达可能成为后期‘无子瓯柑’小孢子败育的重要原因。

      本研究以克里曼丁橘数据库为基础,从‘无子瓯柑’和瓯柑的转录组和蛋白质组数据中筛选出4个CHS同源差异表达基因,其核苷酸序列在‘无子瓯柑’和瓯柑中的相似度达到98%,并且3个同源基因在‘无子瓯柑’中表达的显著下调起始于小孢子母细胞时期。CHS的异常表达可能是‘无子瓯柑’小孢子发育异常从而造成雄性不育的原因。

参考文献 (15)

目录

    /

    返回文章
    返回