高级搜索

留言板

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

基于转录组测序挖掘甜荞自交不亲和相关功能基因

,  ,  ,  ,  ,  ,  ,  ,  ,  ,  , 

上一篇

下一篇

张建, 黄科慧, 刘洋, 等. 基于转录组测序挖掘甜荞自交不亲和相关功能基因[J]. 西南大学学报（自然科学版）, 2022, 44(11): 99-113. doi: 10.13718/j.cnki.xdzk.2022.11.010
引用本文: 张建, 黄科慧, 刘洋, 等. 基于转录组测序挖掘甜荞自交不亲和相关功能基因[J]. 西南大学学报（自然科学版）, 2022, 44(11): 99-113. doi: 10.13718/j.cnki.xdzk.2022.11.010
shu
ZHANG Jian, HUANG Kehui, LIU Yang, et al. Mining of Common Buckwheat Self-Incompatibility Related Gene Based on RNA-Seq[J]. Journal of Southwest University Natural Science Edition, 2022, 44(11): 99-113. doi: 10.13718/j.cnki.xdzk.2022.11.010
Citation: ZHANG Jian, HUANG Kehui, LIU Yang, et al. Mining of Common Buckwheat Self-Incompatibility Related Gene Based on RNA-Seq[J]. Journal of Southwest University Natural Science Edition, 2022, 44(11): 99-113. doi: 10.13718/j.cnki.xdzk.2022.11.010
shu

基于转录组测序挖掘甜荞自交不亲和相关功能基因

  • 1.

    西南大学 农学与生物科技学院/重庆市荞麦产业创新团队，重庆 400715

  • 2.

    甘肃省古浪县防沙治沙技术推广中心，甘肃 古浪 733100

  • 基金项目: 重庆市自然科学基金项目(CSTB2022NSCQ-MSX1034)；中央高校基本业务费专项(SWU119070，XDJK2019C038，XDJK2019D022)；重庆市农委荞麦产业技术体系项目(CQCJT2021001)
详细信息
    作者简介:

    张建，博士，副教授，主要从事作物分子育种研究 .

    通讯作者: 易泽林，副研究员，硕士研究生导师

  • 中图分类号: S517

Mining of Common Buckwheat Self-Incompatibility Related Gene Based on RNA-Seq

  • 1.

    College of Agronomy and Biotechnology, Southwest University, Chongqing 400715, China

  • 2.

    Gulang County Sand Prevention and Control Technology Extension Center of Gansu, Gulang Gansu 733100, China

  • 摘要: 普通甜荞属于花柱异长型自交不亲和1年生草本作物，营养价值丰富. 为了探明甜荞育性和花型发育的重要功能途径和关键基因，以自交不亲和普通甜荞乌克兰大粒荞和自交亲和甜荞品种甜自21为材料，分别对盛花期的甜自21等柱花，乌克兰大粒荞长柱花和短柱花雌雄蕊进行转录组测序，获得283 170条平均长度为1 351 bp，N50为2 040 bp的Unigenes，其中有68 813条差异表达基因(DEGs). GO功能富集分析显示，自交不亲和甜荞中长、短柱花雌雄蕊中的差异基因主要富集在导管形成和导管发育过程中；KEGG代谢通路分析显示，自交不亲和甜荞(S-型和L-型)与自交亲和甜荞(H-型)雌雄蕊中差异基因主要富集在倍半萜和三萜化合物合成途径中；自交不亲和甜荞长、短柱花雌雄蕊中差异基因(L vs S)主要富集在植物激素信号转导途径中. 倍半萜和三萜化合物合成途径共富集到77个差异基因，其中三萜化合物合成途径基因在自交不亲和甜荞中下调表达，倍半萜化合物合成途径基因以上调表达为主. 植物激素信号转导途径中，与自交亲和型等柱花雌雄蕊比较，生长素、茉莉酸和赤霉素信号通路中的差异表达基因在自交不亲和型长、短柱花雌雄蕊中均高表达. 通过qRT-PCR对15个自交亲和性相关差异表达基因进行分析，结果表明：其表达趋势与转录组结果一致，AGL80和PG在短柱花雌雄蕊中高表达，长柱花雌雄蕊中低表达，可能与自交不亲和甜荞长、短柱花发育有关. 锌指蛋白基因CCCH，转录因子ARF，MYC，SP11，SUI1与甜荞育性发育相关. 该文比较分析了不同花型甜荞雌雄蕊转录组测序结果，获得大量Unigenes数据及差异表达基因的相关信息，发现了导管形成与发育对自交不亲和甜荞花柱发育具有重要影响，明确了萜类化合物、植物激素和转录因子对甜荞育性发育具有重要作用.
    Abstract: Common buckwheat (Fagopyrum esculentum), a heteromorphic self-incompatible (SI) species, is the annual herbaceous crops with rich nutritional value. The objective of this study is to explore the important functional pathways and key genes regulating fertility and style development of self-incompatible common buckwheat, which has great significance to further study the molecular mechanism of self-incompatibilities in common buckwheat. Ukraine daliqiao (UD) (heterostyly and self-incompatible) and Tianzi 21 (TZ) (homostyled and self-compatible) as materials, RNA-seq were carried out with the samples of pistils and stamens from S-morph flower (S) and L-morph flower (L) of UD, and H-morph flower of TZ (H). A total of 283 170 unigenes with an average length of 1 351 bp and N50 of 2 040 bp were obtained, among which 68 813 differentially expressed genes (DEGs) were identified. GO function enrichment analysis showed that the DEGs in L vs S were significantly enriched in "angiogenesis" and "blood vessel development". KEGG metabolic pathway analysis showed that DEGs of L vs H and S vs H were significantly enriched in the pathways of sesquiterpenoid and triterpenoid biosynthesis. The most of DEGs of L vs S were enriched in the pathway of plant hormone signal transduction. There were 77 DEGs in "sesquiterpenoid and triterpenoid biosynthesis", among which the gene expression of triterpenoid synthesis pathway was down-regulated and the gene expression of sesquiterpenoid synthesis pathway was mainly up-regulated in self-incompatible common buckwheat. In the plant hormone signal transduction pathway, compared with the self-compatible common buckwheat, the DEGs in auxin, JA and GA had high expression in the self-incompatible common buckwheat. Fifteen DEGs related to self-incompatibility were analyzed by qRT-PCR, and the expression trend was consistent with the results of RNA-Seq analysis. The expression levels of AGL80 and PG were high in the pistils and stamens of S-morph flower and low in the pistils and stamens of L-morph flower in UD, which may be related to the style development of self-incompatible common buckwheat. CCCH encoding zin finger protein, transcription factors ARF, MYC, SP11, SUI1 are related to the fertility development of common buckwheat. This study acquired the gene expression profiling data of the pistils and stamens of different flower style in common buckwheat, found that the formation and development of vessel have an important effect on the style development of self-incompatible common buckwheat, and the terpenoids, certain plant hormones and transcription factor play important roles in fertility development of common buckwheat.
  • 加载中

  • 图 1  甜荞花表型分析及性状测量

    下载: 全尺寸图片 幻灯片

    图 2  L vs H和S vs H差异表达基因分析

    下载: 全尺寸图片 幻灯片

    图 3  L vs H，S vs H和L vs S中差异表达基因GO分析

    下载: 全尺寸图片 幻灯片

    图 4  倍半萜和三萜化合物合成途径差异表达分析

    下载: 全尺寸图片 幻灯片

    图 5  植物激素信号通路中相关基因差异表达分析

    下载: 全尺寸图片 幻灯片

    图 6  与甜荞育性和花柱发育相关的差异表达基因qRT-PCR验证

    下载: 全尺寸图片 幻灯片

    表 1  不同花型雌雄蕊转录组测序结果

    样品 原始序列 过滤后序列 过滤后碱基数 碱基错误率/% Q30/% GC含量/%
    L_1 56420138 55733278 8.36G 0.02 95.02 45.78
    L_2 51854578 51235976 7.69G 0.02 94.89 45.83
    L_3 54170174 53401442 8.01G 0.03 92.39 45.39
    S_1 63152392 61912778 9.29G 0.02 95.03 45.87
    S_2 66402272 65372118 9.81G 0.02 94.82 45.68
    S_3 59586576 58608248 8.79G 0.03 92.51 45.42
    H_1 59225340 58204792 8.73G 0.02 95.04 45.53
    H_2 53393352 52412146 7.86G 0.02 95.04 45.67
    H_3 56415928 55677058 8.35G 0.03 92.79 45.77
    下载: 导出CSV

    表 2  差异表达基因KEGG通路分析_top20

    条目 节点 DG数目 总基因数目 p值 修正后p值
    L vs H Sesquiterpenoid and triterpenoid biosynthesis ko00909 103 265 1.49E-06 0.000 185 026
    Anthocyanin biosynthesis ko00942 38 76 6.05E-05 0.003 751 799
    Ribosome ko03010 665 2 591 0.000 147 0.006 079 685
    beta-Alanine metabolism ko00410 166 574 0.001 189 0.036 855 917
    Diterpenoid biosynthesis ko00904 55 155 0.001 965 0.044 372 750
    Cutin，suberine and wax biosynthesis ko00073 64 190 0.002 504 0.044 372 750
    Arginine and proline metabolism ko00330 209 761 0.002 505 0.044 372 750
    Tryptophan metabolism ko00380 98 344 0.014 270 0.221 180 526
    Glycine，serine and threonine metabolism ko00260 192 737 0.019 200 0.230 627 731
    Phenylalanine metabolism ko00360 105 378 0.019 818 0.230 627 731
    Limonene and pinene degradation ko00903 75 258 0.020 459 0.230 627 731
    Cysteine and methionine metabolism ko00270 237 932 0.023 183 0.239 562 705
    alpha-Linolenic acid metabolism ko00592 186 723 0.029 347 0.279 923 248
    Other types of O-glycan biosynthesis ko00514 22 61 0.034 566 0.306 159 299
    Histidine metabolism ko00340 85 316 0.055 426 0.458 191 199
    Circadian rhythm-plant ko04712 139 544 0.060 126 0.465 980 375
    C5-Branched dibasic acid metabolism ko00660 28 92 0.082 679 0.572 725 017
    Steroid biosynthesis ko00100 95 370 0.095 350 0.572 725 017
    Phenylalanine，tyrosine and tryptophan biosynthesis ko00400 162 657 0.097 898 0.572 725 017
    Linoleic acid metabolism ko00591 64 241 0.099 799 0.572725017
    S vs H Sesquiterpenoid and triterpenoid biosynthesis ko00909 92 265 0.000 388 0.033 980 646
    Anthocyanin biosynthesis ko00942 35 76 0.000 642 0.033 980 646
    alpha-Linolenic acid metabolism ko00592 210 723 0.001 012 0.033 980 646
    Cutin，suberine and wax biosynthesis ko00073 68 190 0.001 096 0.033 980 646
    Flavonoid biosynthesis ko00941 91 286 0.003 496 0.086 699 505
    Arginine and proline metabolism ko00330 203 761 0.021 725 0.283 804 733
    Circadian rhythm-plant ko04712 149 544 0.023 304 0.283 804 733
    Cysteine and methionine metabolism ko00270 244 932 0.024 288 0.283 804 733
    Diterpenoid biosynthesis ko00904 49 155 0.028 043 0.283 804 733
    Plant hormone signal transduction ko04075 500 2 007 0.028 868 0.283 804 733
    beta-Alanine metabolism ko00410 155 574 0.029 894 0.283 804 733
    Regulation of autophagy ko04140 158 588 0.032 232 0.283 804 733
    Stilbenoid，diarylheptanoid and gingerol biosynthesis ko00945 55 180 0.033 530 0.283 804 733
    Oxidative phosphorylation ko00190 295 1 157 0.037 121 0.283 804 733
    Taurine and hypotaurine metabolism ko00430 33 100 0.041 519 0.283 804 733
    Caffeine metabolism ko00232 25 71 0.041 532 0.283 804 733
    Sulfur metabolism ko00920 70 242 0.043 616 0.283 804 733
    Other types of O-glycan biosynthesis ko00514 22 61 0.044 775 0.283 804 733
    Phenylalanine metabolism ko00360 104 378 0.045 970 0.283 804 733
    Tyrosine metabolism ko00350 131 488 0.047 863 0.283 804 733
    L vs S Carotenoid biosynthesis ko00906 5 385 0.019 206 0.557 989 998
    Caffeine metabolism ko00232 2 71 0.033 378 0.557 989 998
    One carbon pool by folate ko00670 3 183 0.037 015 0.557 989 998
    Monoterpenoid biosynthesis ko00902 2 90 0.050 675 0.557 989 998
    Glyoxylate and dicarboxylate metabolism ko00630 7 842 0.050 843 0.557 989 998
    Plant hormone signal transduction ko04075 13 2 007 0.055 964 0.557 989 998
    DNA replication ko03030 5 536 0.062 567 0.557 989 998
    Phenylalanine metabolism ko00360 4 378 0.063 977 0.557 989 998
    Circadian rhythm-plant ko04712 5 544 0.065 735 0.557 989 998
    SNARE interactions in vesicular transport ko04130 4 383 0.066 427 0.557 989 998
    Regulation of autophagy ko04140 5 588 0.084 727 0.634 355 072
    RNA degradation ko03018 8 1 202 0.104 618 0.634 355 072
    Flavonoid biosynthesis ko00941 3 286 0.104 991 0.634 355 072
    Photosynthesis ko00195 3 303 0.119 043 0.634 355 072
    Phenylalanine，tyrosine and tryptophan biosynthesis ko00400 5 657 0.119 639 0.634 355 072
    Photosynthesis-antenna proteins ko00196 2 155 0.126 108 0.634 355 072
    Tyrosine metabolism ko00350 4 488 0.128 381 0.634 355 072
    Terpenoid backbone biosynthesis ko00900 4 501 0.137 319 0.640 822 483
    Stilbenoid，diarylheptanoid and gingerol biosynthesis ko00945 2 180 0.159 543 0.697 022 324
    Ubiquitin mediated proteolysis ko04120 8 1 374 0.176 527 0.697 022 324
    下载: 导出CSV

    表 3  与甜荞育性和花柱发育相关的差异表达基因qRT-PCR引物

    基因名称 上游引物 下游引物
    AGL80(Cluster-11933.108535) GGCCCTTCGATTCTCGTAACT GATTATGATGACCGGCTTTGCC
    ARF(Cluster-11933.100857) TCCCAGTTTGCAGTTCACATG AGTAAGGCTATTCGGGGTGAC
    ARF(Cluster-11933.94328) CATGTTCGAATCCCCACCTGA ACCGTGCTGATCCTGAAAGG
    CCCH(Cluster-11933.183409) CCTGGAAGCGCCATAAGGAA TTGGCTTTTCAGAAGGACAAGT
    MYC(Cluster-11933.187734) CAAAGAGCAGGGCCACCTAAG GCAAGCCCGCTTACCTTCCT
    PG(Cluster-11933.6467) CTGGAACAGGAGCCACGATT CAACGTCGTTGACTATGGTGCA
    SP11(Cluster-11933.123004) CTCCGCGCTGTTGAATTTGG GGAATCATACAGGATGAACTGGT
    SP11(Cluster-11933.168503) AAGCACGGAATCTCCGGTAA AAAAATGGACATGGCAGATGGT
    SP11(Cluster-11933.66132) GGGTGCAAATTATTCTGAAAAGTCA TCTCGTATGGAGCAATCGGC
    SP11(Cluster-11933.98978) GGAGTGAAGCTTTCGCTCTCT ACCGGAGATGAAGTTCAAGTG
    SUI1(Cluster-11933.145211) AGGGCTAGGGCGGTTCAT AATAGCTGCTACTCGTGGGC
    SUI1(Cluster-11933.126825) AGTTGTGTTCCCCACCTACG GCCACAACGTATGCGAGTTG
    THL(Cluster-11933.24640) CAGCTAATTTCGCGGAAGAAACC ACCTCCAGGACTTCCGATCA
    UN(Cluster-11933.147201) AGTCATTCCGAAAGTTCTTGGGT TGCTTGATCTACCCTTTGGTG
    UN(Cluster-11933.213250) GCTGGCAGTCATTACGGATG TATGGCAATAGCCACCAGCAT
    下载: 导出CSV
  • [1] 柴岩, 张宗文. 荞麦研究进展: 第十届国际荞麦会议论文集[M]. 杨凌: 西北农林科技大学出版社, 2007.
    [2] 蔡齐宗, 王佳蕊, 陈庆富, 等. 苦荞全基因组SSR位点鉴定及分子标记开发[J]. 河南农业大学学报, 2022, 56(3): 392-400. doi: 10.16445/j.cnki.1000-2340.20220506.001
    [3] 李春花, 加英多拉, 田娟, 等. 自交可育红花甜荞种质资源创新利用研究[J]. 南方农业学报, 2021, 52(10): 2751-2757. doi: 10.3969/j.issn.2095-1191.2021.10.015
    [4] KOYAMA M, NAKAMURA C, NAKAMURA K. Changes in Phenols Contents from Buckwheat Sprouts during Growth Stage[J]. Journal of Food Science and Technology, 2013, 50(1): 86-93. doi: 10.1007/s13197-011-0316-1
    [5] 柴岩, 王鹏科, 冯佰利. 中国小杂粮产业发展指南[M]. 杨凌: 西北农林科技大学出版社, 2007.
    [6] GARBER R J, QUISENBERRY K S. Self-Fertilization in Buckwheat[J]. J Agric Res, 1927, 34: 185-190.
    [7] SŁOMKA A, SYCHTA K, DUBERT F, et al. Embryological Background of Low Seed Set in Distylous Common Buckwheat (Fagopyrum esculentum Moench) with Biased Morph Ratios, and Biostimulant-Induced Improvement of it[J]. Crop and Pasture Science, 2017, 68(7): 680-690. doi: 10.1071/CP17009
    [8] OHNISHI O. Search for the Wild Ancestor of Buckwheat. Ⅰ. Description of New Fagopyrum (Polygonaceae) Species and Their Distribution in China and Himalayan Hills[J]. Fagopyrum, 1998, 15: 18-28.
    [9] MATSUI K, TETSUKA T, NISHIO T, et al. Heteromorphic Incompatibility Retained in Self-Compatible Plants Produced by a Cross between Common and Wild Buckwheat[J]. New Phytologist, 2003, 159(3): 701-708. doi: 10.1046/j.1469-8137.2003.00840.x
    [10] WANG Y J, SCARTH R, CAMPBELL G C. Inheritance of Seed Shattering in Interspecific Hybrids between Fagopyrum esculentum and F. homotropicum[J]. Crop Science, 2005, 45(2): 693-697. doi: 10.2135/cropsci2005.0693
    [11] 陈庆富. 荞麦生产状况及新类型栽培荞麦育种研究的最新进展[J]. 贵州师范大学学报(自然科学版), 2018, 36(3): 1-7, 131. doi: 10.16614/j.gznuj.zrb.2018.03.001
    [12] MATSUI K, YASUI Y S. Buckwheat Heteromorphic Self-Incompatibility: Genetics, Genomics and Application to Breeding[J]. Breeding Science, 2020, 70(1): 32-38. doi: 10.1270/jsbbs.19083
    [13] YASUI Y, MORI M, AII J, et al. S-LOCUS EARLY FLOWERING 3 is Exclusively Present in the Genomes of Short-Styled Buckwheat Plants that Exhibit Heteromorphic Self-Incompatibility[J]. PLoS ONE, 2012, 7(2): e31264. doi: 10.1371/journal.pone.0031264
    [14] WOO S H, ADACHI T, JONG S K, et al. Inheritance of Self-Compatibility and Flower Morphology in an Inter-Specific Buckwheat Hybrid[J]. Canadian Journal of Plant Science, 1999, 79(4): 483-490. doi: 10.4141/P98-117
    [15] MATSUI K, NISHIO T, TETSUKA T. Genes Outside the S Supergene Suppress S Functions in Buckwheat (Fagopyrum esculentum)[J]. Annals of Botany, 2004, 94(6): 805-809. doi: 10.1093/aob/mch206
    [16] MIZUNO N, YASUI Y S. Gene Flow Signature in the S-Allele Region of Cultivated Buckwheat[J]. BMC Plant Biology, 2019, 19(1): 125-138. doi: 10.1186/s12870-019-1730-1
    [17] MATSUI K, YASUI Y S. Genetic and Genomic Research for the Development of an Efficient Breeding System in Heterostylous Self-Incompatible Common Buckwheat (Fagopyrum esculentum)[J]. Theoretical and Applied Genetics, 2020, 133(5): 1641-1653. doi: 10.1007/s00122-020-03572-6
    [18] BOWLER C, CHUA N H. Emerging Themes of Plant Signal Transduction[J]. The Plant Cell, 1994, 6(11): 1529-1541.
    [19] FOURQUIN C, FERRÁNDIZ C. Functional Analyses of AGAMOUS Family Members in Nicotiana Benthamiana Clarify the Evolution of Early and Late Roles of C-Function Genes in Eudicots[J]. The Plant Journal, 2012, 71(6): 990-1001. doi: 10.1111/j.1365-313X.2012.05046.x
    [20] LU S W, ZHANG Y, ZHU K J, et al. The Citrus Transcription Factor CsMADS6 Modulates Carotenoid Metabolism by Directly Regulating Carotenogenic Genes[J]. Plant Physiology, 2018, 176(4): 2657-2676. doi: 10.1104/pp.17.01830
    [21] USHIJIMA K, NAKANO R, BANDO M, et al. Isolation of the Floral Morph-Related Genes in Heterostylous Flax (Linum Grandiflorum): The Genetic Polymorphism and the Transcriptional and Post-Transcriptional Regulations of the S Locus[J]. The Plant Journal for Cell and Molecular Biology, 2012, 69(2): 317-331. doi: 10.1111/j.1365-313X.2011.04792.x
    [22] VIERSTRA R D. The Ubiquitin/26S Proteasome Pathway, the Complex last Chapter in the Life of many Plant Proteins[J]. Trends in Plant Science, 2003, 8(3): 135-142. doi: 10.1016/S1360-1385(03)00014-1
    [23] SHI T X, LI R Y, CHEN Q J, et al. De Novo Sequencing of Seed Transcriptome and Development of Genic-SSR Markers in Common Buckwheat (Fagopyrum esculentum)[J]. Molecular Breeding, 2017, 37(12): 1-15.
    [24] FANG X M, ZHANG Y L, ZHANG Y K, et al. De Novo Transcriptome Assembly and Identification of Genes Related to Seed Size in Common Buckwheat (Fagopyrum esculentum M. )[J]. Breeding Sci, 2019, 69(3): 487-497. doi: 10.1270/jsbbs.18194
    [25] XU J M, FAN W, JIN J F, et al. Transcriptome Analysis of Al-Induced Genes in Buckwheat (Fagopyrum esculentum Moench) Root Apex: New Insight into Al Toxicity and Resistance Mechanisms in an Al Accumulating Species[J]. Frontiers in Plant Science, 2017, 8: 1141. doi: 10.3389/fpls.2017.01141
    [26] 黄娟, 邓娇, 陈庆富. 荞麦根的转录组学分析及黄酮合成基因的鉴定[J]. 中国农业科技导报, 2017, 19(2): 9-19. doi: 10.13304/j.nykjdb.2016.206
    [27] YOKOSHO K, YAMAJI N, MA J F. Global Transcriptome Analysis of Al-Induced Genes in an Al-Accumulating Species, Common Buckwheat (Fagopyrum esculentum Moench)[J]. Plant and Cell Physiology, 2014, 55(12): 2077-2091. doi: 10.1093/pcp/pcu135
    [28] LOGACHEVA M D, KASIANOV A S, VINOGRADOV D V, et al. De Novo Sequencing and Characterization of Floral Transcriptome in Two Species of Buckwheat (Fagopyrum)[J]. BMC Genomics, 2011, 12: 30. doi: 10.1186/1471-2164-12-30
    [29] FANG Z W, HOU Z H, WANG S P, et al. Transcriptome Analysis Reveals the Accumulation Mechanism of Anthocyanins in Buckwheat (Fagopyrum esculentum Moench) Cotyledons and Flowers[J]. International Journal of Molecular Sciences, 2019, 20(6): 1493. doi: 10.3390/ijms20061493
    [30] LANGE B M, RUJAN T, MARTIN W, et al. Isoprenoid Biosynthesis: The Evolution of Two Ancient and Distinct Pathways across Genomes[J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(24): 13172-13177. doi: 10.1073/pnas.240454797
    [31] LIU Y, WANG H, YE H C, et al. Advances in the Plant Isoprenoid Biosynthesis Pathway and Its Metabolic Engineering[J]. Journal of Integrative Plant Biology, 2005, 47(7): 769-782. doi: 10.1111/j.1744-7909.2005.00111.x
    [32] TETALI S D. Terpenes and Isoprenoids: a Wealth of Compounds for Global Use[J]. Planta, 2019, 249(1): 1-8. doi: 10.1007/s00425-018-3056-x
    [33] PICHERSKY E, GERSHENZON J. The Formation and Function of Plant Volatiles: Perfumes for Pollinator Attraction and Defense[J]. Current Opinion in Plant Biology, 2002, 5(3): 237-243. doi: 10.1016/S1369-5266(02)00251-0
    [34] BRODERICK S R, WIJERATNE S, WIJERATN A J, et al. RNA-Sequencing Reveals Early, Dynamic Transcriptome Changes in the Corollas of Pollinated Petunias[J]. BMC Plant Biology, 2014, 14: 307. doi: 10.1186/s12870-014-0307-2
    [35] OVERVOORDE P, FUKAKI H, BEECKMAN T. Auxin Control of Root Development[J]. Cold Spring Harbor Perspectives in Biology, 2010, 2(6): a001537.
    [36] CHEN D, ZHAO J. Free IAA in Stigmas and Styles during Pollen Germination and Pollen Tube Growth of Nicotiana Tabacum[J]. Physiologia Plantarum, 2008, 134(1): 202-215. doi: 10.1111/j.1399-3054.2008.01125.x
    [37] HASENSTEIN K H, ZAVADA M S. Auxin Modification of the Incompatibility Response in Theobroma Cacao[J]. Physiologia Plantarum, 2001, 112(1): 113-118. doi: 10.1034/j.1399-3054.2001.1120115.x
    [38] KOVALEVA L, ZAKHAROVA E. Hormonal Status of the Pollen-Pistil System at the Progamic Phase of Fertilization after Compatible and Incompatible Pollination in Petunia Hybrida L[J]. Sexual Plant Reproduction, 2003, 16(4): 191-196. doi: 10.1007/s00497-003-0189-1
    [39] SOLFANELLI C, BARTOLINI S, VITAGLIANO C, et al. Immunolocalization and Quantification of IAA after Self-and Free-Pollination in Olea Europaea L[J]. Scientia Horticulturae, 2006, 110(4): 345-351. doi: 10.1016/j.scienta.2006.06.026
    [40] FERNÁNDEZ-CALVO P, CHINI A, FERNÁNDEZ-BARBERO G, et al. The Arabidopsis BHLH Transcription Factors MYC3 and MYC4 are Targets of JAZ Repressors and Act Additively with MYC2 in the Activation of Jasmonate Responses[J]. The Plant Cell, 2011, 23(2): 701-715. doi: 10.1105/tpc.110.080788
    [41] PARK J H, HALITSCHKE R, KIM H B, et al. A Knock-out Mutation in Allene Oxide Synthase Results in Male Sterility and Defective Wound Signal Transduction in Arabidopsis Due to a Block in Jasmonic Acid Biosynthesis[J]. The Plant Journal, 2002, 31(1): 1-12. doi: 10.1046/j.1365-313X.2002.01328.x
    [42] QI T C, SONG S S, REN Q C, et al. The Jasmonate-ZIM-Domain Proteins Interact with the WD-Repeat/BHLH/MYB Complexes to Regulate Jasmonate-Mediated Anthocyanin Accumulation and Trichome Initiation in Arabidopsis Thaliana[J]. The Plant Cell, 2011, 23(5): 1795-1814. doi: 10.1105/tpc.111.083261
    [43] YOSHIDA Y, SANO R, WADA T, et al. Jasmonic Acid Control of GLABRA3 Links Inducible Defense and Trichome Patterning in Arabidopsis[J]. Development, 2009, 136(6): 1039-1048. doi: 10.1242/dev.030585
    [44] JIANG Y J, LIANG G, YANG S Z, et al. Arabidopsis WRKY57 Functions as a Node of Convergence for Jasmonic Acid-and Auxin-Mediated Signaling in Jasmonic Acid-Induced Leaf Senescence[J]. The Plant Cell, 2014, 26(1): 230-245. doi: 10.1105/tpc.113.117838
    [45] ZHANG X, ZHU Z Q, AN F Y, et al. Jasmonate-Activated MYC2 Represses ETHYLENE INSENSITIVE3 Activity to Antagonize Ethylene-Promoted Apical Hook Formation in Arabidopsis[J]. The Plant Cell, 2014, 26(3): 1105-1117. doi: 10.1105/tpc.113.122002
    [46] SONG S S, QI T C, HUANG H, et al. The Jasmonate-ZIM Domain Proteins Interact with the R2R3-MYB Transcription Factors MYB21 and MYB24 to Affect Jasmonate-Regulated Stamen Development in Arabidopsis[J]. The Plant Cell, 2011, 23(3): 1000-1013. doi: 10.1105/tpc.111.083089
    [47] HUANG H, GONG Y L, LIU B, et al. The DELLA Proteins Interact with MYB21 and MYB24 to Regulate Filament Elongation in Arabidopsis[J]. BMC Plant Biology, 2020, 20(1): 64-72. doi: 10.1186/s12870-020-2274-0
    [48] SHI D Q, TANG C, WANG R Z, et al. Transcriptome and Phytohormone Analysis Reveals a Comprehensive Phytohormone and Pathogen Defence Response in Pear Self-/Cross-Pollination[J]. Plant Cell Reports, 2017, 36(11): 1785-1799. doi: 10.1007/s00299-017-2194-0
    [49] ACHARD P, GENSCHIK P. Releasing the Brakes of Plant Growth: How GAs Shutdown DELLA Proteins[J]. Journal of Experimental Botany, 2009, 60(4): 1085-1092. doi: 10.1093/jxb/ern301
    [50] YAMAGUCHI S. Gibberellin Metabolism and Its Regulation[J]. Annual Review of Plant Biology, 2008, 59: 225-251. doi: 10.1146/annurev.arplant.59.032607.092804
    [51] LEE S, CHENG H, KING K E, et al. Gibberellin Regulates Arabidopsis Seed Germination via RGL2, a GAI/RGA-Like Gene whose Expression is Up-Regulated Following Imbibition[J]. Genes & Development, 2002, 16(5): 646-658.
    [52] TYLER L, THOMAS S G, HU J H, et al. DELLA Proteins and Gibberellin-Regulated Seed Germination and Floral Development in Arabidopsis[J]. Plant Physiology, 2004, 135(2): 1008-1019. doi: 10.1104/pp.104.039578
    [53] DILL A, SUN T P. Synergistic Derepression of Gibberellin Signaling by Removing RGA and GAI Function in Arabidopsis Thaliana[J]. Genetics, 2001, 159(2): 777-785. doi: 10.1093/genetics/159.2.777
  • 加载中
图( 6) 表( 3)
计量
  • 文章访问数:  8578
  • HTML全文浏览数:  8578
  • PDF下载数:  495
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-12-21
  • 刊出日期:  2022-11-20

基于转录组测序挖掘甜荞自交不亲和相关功能基因

    通讯作者: 易泽林，副研究员，硕士研究生导师
    作者简介: 张建，博士，副教授，主要从事作物分子育种研究
  • 1. 西南大学 农学与生物科技学院/重庆市荞麦产业创新团队，重庆 400715
  • 2. 甘肃省古浪县防沙治沙技术推广中心，甘肃 古浪 733100
  • 收稿日期:  2021-12-21
基金项目:  重庆市自然科学基金项目(CSTB2022NSCQ-MSX1034)；中央高校基本业务费专项(SWU119070，XDJK2019C038，XDJK2019D022)；重庆市农委荞麦产业技术体系项目(CQCJT2021001)

摘要: 普通甜荞属于花柱异长型自交不亲和1年生草本作物，营养价值丰富. 为了探明甜荞育性和花型发育的重要功能途径和关键基因，以自交不亲和普通甜荞乌克兰大粒荞和自交亲和甜荞品种甜自21为材料，分别对盛花期的甜自21等柱花，乌克兰大粒荞长柱花和短柱花雌雄蕊进行转录组测序，获得283 170条平均长度为1 351 bp，N50为2 040 bp的Unigenes，其中有68 813条差异表达基因(DEGs). GO功能富集分析显示，自交不亲和甜荞中长、短柱花雌雄蕊中的差异基因主要富集在导管形成和导管发育过程中；KEGG代谢通路分析显示，自交不亲和甜荞(S-型和L-型)与自交亲和甜荞(H-型)雌雄蕊中差异基因主要富集在倍半萜和三萜化合物合成途径中；自交不亲和甜荞长、短柱花雌雄蕊中差异基因(L vs S)主要富集在植物激素信号转导途径中. 倍半萜和三萜化合物合成途径共富集到77个差异基因，其中三萜化合物合成途径基因在自交不亲和甜荞中下调表达，倍半萜化合物合成途径基因以上调表达为主. 植物激素信号转导途径中，与自交亲和型等柱花雌雄蕊比较，生长素、茉莉酸和赤霉素信号通路中的差异表达基因在自交不亲和型长、短柱花雌雄蕊中均高表达. 通过qRT-PCR对15个自交亲和性相关差异表达基因进行分析，结果表明：其表达趋势与转录组结果一致，AGL80和PG在短柱花雌雄蕊中高表达，长柱花雌雄蕊中低表达，可能与自交不亲和甜荞长、短柱花发育有关. 锌指蛋白基因CCCH，转录因子ARF，MYC，SP11，SUI1与甜荞育性发育相关. 该文比较分析了不同花型甜荞雌雄蕊转录组测序结果，获得大量Unigenes数据及差异表达基因的相关信息，发现了导管形成与发育对自交不亲和甜荞花柱发育具有重要影响，明确了萜类化合物、植物激素和转录因子对甜荞育性发育具有重要作用.

English Abstract

    全文HTML

  • 开放科学(资源服务)标志码(OSID):

  • 甜荞Fagopyrum esculentum和苦荞Fagopyrm tataricum都属于蓼科Polygonaceae荞麦属Fagpyrum的1年生草本作物，起源于中国，具有生长期短、抗逆性强的特点[1-3]，富含多种矿质营养元素(铁、钾、锌、硒等)和保健成分(芦丁、槲皮素、不饱和脂肪酸、抗性淀粉、β-葡聚糖等)[4-5]，是重要药食兼用作物. 栽培种甜荞花型为异型花(雌雄蕊不等长)，分为2种类型，雌蕊长于雄蕊(长柱花，pin)和雌蕊短于雄蕊(短柱花，thrin)[6]. 普通栽培种甜荞自花授粉或同类型花植株之间相互授粉不结实，具有自交不亲和现象[7]；同型花雌雄蕊等长，能够自花授粉结实，常见于野生种中. 近年来，国内外多个研究团队利用自交可育野生荞麦F. homotropicum[8]，通过种间杂交的方式(F. esculentum和F. homotropicum)培育出自交可育甜荞，即雌雄蕊等长的同型花甜荞[9-12]. Yasui等[13]指出荞麦花的形态是thrum和pin两种形式，通过一个具有两个等位基因状态的基因(S和s)，控制thrum(Ss)和pin(ss)形式. 与其他二态型作物一样，花的形态和种内不亲和性均由单个S位点决定，且S-型植株为杂合子(S/s)，L-型植株在该位点为纯合隐性(s/s)，未发现S/S基因型[13]. 但随后在中国云南省发现了花柱等长自交亲和野生甜荞. Woo等[14]发现自交亲和性受单基因Sh控制，甜荞花型S基因的优势关系为S＞Sh＞s，自交亲和性等位基因Sh保留了异型性不亲和性，并认为Sh等位基因来源于S超基因复合体[15, 12]. Yasui等[13]在S-Locus基因区注释到32个预测基因，并在这32个预测基因中发现2个与SI相关的候选基因，编码RING/U-box超家族蛋白. Mizuno等[16]利用GBS技术获得甜荞全基因组SNP标记，S等位基因区域内的SNP密度和遗传多样性分析表明该区域具有高度多样性. 由于普通甜荞属于自交不亲和的虫媒传粉作物，花期遇阴雨天气将严重影响结实率，导致每667 m2产量低且十分不稳定(全国平均单产量约为60 kg/667 m2)[11]，而通过种间杂交培育的自交可育甜荞，因携带野生荞麦的特性，表现出结实率低、产量低、极易落粒，自交衰退等特点[11, 17]. 因此，甜荞产量与甜荞育性息息相关，对甜荞育性及花柱发育的分子机理进行深入研究具有重要意义. 本研究以自交亲和甜自21和自交不亲和乌克兰大粒荞为材料，分别对甜自21等柱花、乌克兰大粒荞长柱花和短柱花雌雄蕊进行RNA-seq，通过对转录组数据的分析解析甜荞自交不亲和信号传导网络，挖掘自交不亲和相关基因，阐明自交不亲和的分子机理，揭示甜荞育性及花柱发育的分子机制.

1.   材料与方法

    1.1.   试验材料

  • 2018年秋，自交亲和甜荞品种甜自21和贵甜2号(来源于贵州师范大学)，自交不亲和甜荞品种乌克兰大粒荞和酉荞2号分别隔离种植于重庆市北碚区西南大学校内实验基地，常规田间管理，待开花初期将乌克兰大粒荞和酉荞2号长柱花(L-型)植株和短柱花(S-型)植株移栽后隔离种植. 盛花期分别对甜自21和贵甜2号等柱花，乌克兰大粒荞和酉荞2号长柱花和短柱花的花药直径、花冠直径、雄蕊高、花柱高和花柱到花药的距离进行测量并统计分析. 于清晨7：00~9：00点，分别取甜自21和贵甜2号等柱花，乌克兰大粒荞和酉荞2号长柱花和短柱花雌雄蕊，并进行3次生物学重复，样品经液氮速冻后-80 ℃保存，用于RNA提取.

    • 1.2.   RNA提取及转录组测序

  • 采用博日Bioflux的Biospin多糖多酚植物总RNA提取试剂盒提取各个花型雌雄蕊的总RNA. 通过1%琼脂糖凝胶电泳检测总RNA的完整性，使用Agilent Bioanalyzer 2100 System对RNA精确质检. 将检测合格的甜自21等柱花(H)，乌克兰大粒荞长柱花(L)和短柱花(S)雌雄蕊总RNA送至北京诺禾致源科技股份有限公司测序，构建cDNA文库，在Illumina HiSeq 4000平台上进行测序. 为了保证信息分析的质量，去除测序得到的原始测序序列(原始数据)里面含有带接头的、低质量的读长后，得到有效数据，然后采用Trinity软件对有效数据进行拼接，得到拼接的转录本. 然后利用Corset层次聚类软件对转录本进行层次聚类，Corset层次聚类后得到最长的Cluster被确定为Unigene，用于后续的分析. 得到的Unigene分别与NCBI(National Center for Biotechnology Information)中的NR和Nt数据库、KO(KEGG ortholog database)数据库、Swiss-Prot蛋白质数据库、Pfam蛋白数据库、GO数据库、KOG/COG数据库进行比对，e-value值设置为小于1e-5. 然后，对注释到GO数据库中的Unigene进行GO功能分类，对注释到COG/KOG数据库中的Unigene进行COG/KOG功能分类分析，对注释到KO数据库中的Unigene进行KEGG代谢途径分析.

    • 1.3.   基因表达数据分析

  • 以Trinity拼接得到的转录组作为参考序列，将每个样品的有效数据在参考序列上做比对，采用RSEM软件，设置bowtie2参数为mismatch 0，对bowtie的比对结果进行统计，得到了每个样品比对到每个基因上的读长数目，并对其进行FPKM转换，进而得到了所有基因的表达水平. 再用DEGseq进行差异分析，筛选阈值为qvalue＜0.005且|log2.Fold_change|＞1. 然后，分别对差异表达(上调和下调)基因进行GO功能富集分析及KEGG代谢途径富集分析.

    • 1.4.   实时荧光定量PCR(qRT-PCR)分析

  • 根据转录组分析结果，筛选出可能与甜荞育性和花柱发育相关的差异表达基因，利用Primer-NCBI在线设计引物. 同时，在自交亲和甜自21和贵甜2号等柱花，自交不亲和乌克兰大粒荞和酉荞2号长柱花和短柱花雌雄蕊总RNA中，利用实时荧光定量PCR进一步验证差异基因的表达特性. 将RNA采用第1链cDNA合成试剂盒(Invitrogen)反转录为cDNA，以cDNA为模板，FeActin为内参，采用TaKaRa(日本)公司的SYBR©Premix Ex TaqTMⅡ(Tli RNaseH Plus)试剂盒，在ABI 7500FAST荧光定量PCR仪(ABI公司，美国)进行荧光定量检测. 每个样品设3次技术重复，依照2-ΔΔCT法计算相对表达量.

2.   结果

    2.1.   甜荞花表型分析及性状测量

  • 选取花柱异长型自交不亲和甜荞品种乌克兰大粒荞(UD)和花柱等长型自交亲和甜荞品种甜自21(TZ)作为实验材料(图 1A). 表型分析结果表明，不同花柱类型花的花冠直径差异不明显，短柱花柱头到花药的距离、花药直径、雄蕊高均大于长柱花，花柱高则是长柱花高于短柱花. 等柱花的花药直径、雄蕊高、花柱高均在长柱花和短柱花之间，且短柱花柱头到花药的距离最小(图 1). 方差分析结果显示，除花药直径外，花冠直径、雄蕊高、花柱高、花柱到花药的距离在长柱花、短柱花和等柱花之间差异均具有统计学意义(p＜5%)(图 1).

    • 2.2.   RNA-seq测序及差异基因筛选

  • 甜自21等柱花(H)，乌克兰大粒荞长柱花(L)和短柱花(S)雌雄蕊总RNA样品经转录组Illumina HiSeq 4000测序平台测序后，9个测序样品得到的原始数据，经过严格的质量评估和低质量数据过滤，使其得到高质量的有效数据，Q30(碱基被测错的概率为1‰)均在92%以上，GC含量均在45%左右(表 1)，说明得到的有效数据质量和准确度较高，满足后续的数据分析需求.

    高质量的有效数据进行拼接后共得到326 245个转录本，平均长度为1 209 bp，N50的长度为1 981 bp. 得到的转录本利用Corset软件进行层次聚类，通过层次聚类后得到Clusters，这些Clusters中具有最长序列的被确定为Unigene，最终得到283 170个平均长度为1 351 bp，N50为2 040 bp的Unigenes.

    通过分析长柱花(L-型)乌克兰大粒荞、短柱花(S-型)乌克兰大粒荞和等柱花(H-型)甜自21雌雄蕊中的基因差异表达，以等柱花雌雄蕊(H-型)作为对照，在长柱花(L-型)雌雄蕊中上调表达27 788个，下调表达25 244个；在短柱花(S-型)雌雄蕊中上调表达29 929个，下调表达25 963个；长(L-型)、短柱花(S-型)雌雄蕊同时上调基因20 466个，同时下调基因19 753个. 通过对长柱花(L-型)乌克兰大粒荞和短柱花(S-型)乌克兰大粒荞雌雄蕊中的基因进行差异表达分析，共716个差异基因，其中长柱花(L-型)雌雄蕊相对于短柱花(S-型)雌雄蕊上调表达357个，下调表达359个(图 2).

    • 2.3.   差异表达基因GO功能分析

  • 差异表达基因GO富集分析表明，自交不亲和甜荞雌雄蕊(L-型和S-型)与自交亲和甜荞雌雄蕊(H-型)相比较(L vs H，S vs H)，差异表达基因特异性富集到氧化还原过程(GO：0055114，oxidation-reduction process)和氧化还原酶活性(GO：0016491，oxidoreductase activity；GO：0016705，oxidoreductase activity，acting on paired donors，with incorporation or reduction of molecular oxygen)最为显著，其次在血红素结合(GO：0020037，heme binding)、铁离子结合(GO：0005506，iron ion binding)、单个有机体代谢过程(GO：0044710，single-organism metabolic process)、转移酶活性(GO：0016758，transferase activity，transferring hexosyl groups)和四吡咯结合(GO：0046906，tetrapyrrole binding)等均有显著富集(图 3)，表明这些GO分类中的差异基因与甜荞亲和性有关. 在自交不亲和甜荞中，长柱花(L-型)雌雄蕊与短柱花(S-型)雌雄蕊(L vs S)差异表达基因特异性富集到导管形成过程(GO：0001525，angiogenesis)最为显著，其次为导管发育(GO：0001568，blood vessel development)、导管系统开发(GO：0001944，vasculature development)、导管生成调控(GO：0045765，regulation of angiogenesis)和导管形态发生(GO：0048514，blood vessel morphogenesis)(图 3)，表明导管形成相关基因调控甜荞花柱发育.

    • 2.4.   差异表达基因KEGG通路分析

  • 差异表达基因KEGG富集分析表明，自交不亲和甜荞与自交亲和甜荞雌雄蕊(L vs H，S vs H)差异表达基因均显著富集到倍半萜和三萜生物合成途径(ko00909，Sesquiterpenoid and triterpenoid biosynthesis)，此外还有花色苷合成途径(ko00942，Anthocyanin biosynthesis)，二萜化合物合成过程(ko00904，Diterpenoid biosynthesis)，角质、软木脂和蜡质合成(ko00073，Cutin，suberine and wax biosynthesis)，β-丙氨酸代谢(ko00410，beta-Alanine metabolism)，精氨酸和脯氨酸代谢(ko00330，Arginine and proline metabolism)，半胱氨酸和蛋氨酸代谢(ko00270，Cysteine and methionine metabolism)等(表 2)，表明萜类化合物可能与甜荞亲和性有关. 在自交不亲和甜荞中，未检测到长柱花和短柱花雌雄蕊(L vs S)差异基因显著富集到某个代谢途径，但植物激素信号转导(ko04075，Plant hormone signal transduction)中的差异基因最多.

    • 2.5.   萜类化合物合成途径基因差异表达分析

  • 在倍半萜和三萜化合物合成途径(ko00909)中共检测到77个差异基因，与自交亲和甜荞雌雄蕊(H-型)相比，40个DEGs在长柱花雌雄蕊(L-型)和短柱花雌雄蕊(S-型)中同时上调表达，37个DEGs同时下调表达. 其中，三萜化合物(Triterpenoid)合成途径基因FDFT1(farnesyl-diphosphate farnesyltransferase)和SQLE(squalene monooxygenase)以下调表达基因为主. 倍半萜化合物有多个类型，无环型类倍半萜(Acyclic sesquiterpenoid)合成途径基因AFS1(alpha-farnesene synthase)在自交不亲和甜荞长柱花(L)和短柱花(S)雌雄蕊中以下调表达为主；大根香叶烯(Germacren-type)和蛇麻烯(Humule-type)合成途径相关基因SSTLE1(germacrene C synthase)，GERD((-)-germacrene D synthase)，TPS1(valencene/7-epi-alpha-selinene synthase)，HVS(vetispiradiene synthase)，TPS21(valencene/7-epi-alpha-selinene synthase)多数为上调表达基因(图 4).

    • 2.6.   植物激素相关基因差异表达分析

  • 植物激素信号转导对于激素触发的生化变化非常重要[16]. 本研究发现，长柱花雌雄蕊相对于等柱花雌雄蕊上调表达基因在“植物激素信号转导”途径(Ko04075)中显著富集. 而且，在Ko04075途径中发现，与自交亲和型等柱花雌雄蕊比较，生长素、茉莉酸和赤霉素信号通路中的差异表达基因在自交不亲和型长、短柱花雌雄蕊中均高表达(图 5). 生长素(AUX/IAA)信号通路中，相对于等柱花(H)雌雄蕊，AUX1，TIR1(transport inhibitor response 1)，AUX/IAA，ARF在长柱花(L)和短柱花(S)雌雄蕊中上调表达；茉莉酸相关信号通路与α-亚麻酸代谢密切相关，在茉莉酸信号通路中共发现13个基因，分别为6个GOI1，2个JAR1(jasmonic acid-amino synthetase)和5个MYC2；在赤霉素(GA)信号通路中，与自交亲和型等柱花雌雄蕊比较，自交不亲和型长、短柱花雌雄蕊中表达上调的基因为10个，其中1个GID1(gibberellin receptor GID1)，3个DELLA和6个TF.

    • 2.7.   甜荞花柱异型候选基因的qRT-PCR验证

  • 为消除品种特异性影响并验证转录组测序结果的准确性，本研究分别提取自交不亲和的两个甜荞品种乌克兰大粒荞(UD)和酉荞2号(YQ)L-型雌雄蕊、S-型雌雄蕊和自交亲和的两个甜荞品种甜自21(TZ)和贵甜2号(GT) H-型雌雄蕊总RNA，利用qRT-PCR技术，根据已有参考文献对MADS-box[19-20]、MYB[21]及U-box超基因家族[22]等16个差异表达基因进行了qRT-PCR(表 3)，检测其在不同品种、不同花柱类型中的表达情况，发现与通过转录组测序FPKM方法计算读长数得出的表达水平，在不同花柱类型表现出一致的表达趋势，表明了转录组测序结果的准确性(图 6).

    此外，MADS-box转录因子家族AGAMOUS-like protein (AGL，Cluster-11933.108535)和糖代谢过程中的编码多聚半乳糖醛酸酶基因(Polygalacturonases，PG，Cluster-11933.6467)在短柱花(S-型)雌雄蕊中高表达，长柱花(L-型)雌雄蕊中低表达，等柱花(H-型)中不表达，表明这两个基因可能与甜荞花柱发育有关. 锌指蛋白基因CCCH(Cluster-11933.183409)和U-box超基因家族基因SP11(Cluster-11933.98978)在自交不亲和性甜荞(S-型和L-型)雌雄蕊中表达量较低，甚至不表达，但在自交亲和甜荞(H-型)中均高表达. 相反，转录因子ARF(Cluster-11933.100857，Cluster-11933.100857，Cluster-11933.94328)，MYC(Cluster-11933.187734)，U-box超基因家族基因SP11(Cluster-11933.123004，Cluster-11933.168503，Cluster-11933.66132，Cluster-11933.98978)，转录抑制子SUI1(Cluster-11933.126825)在自交不亲和甜荞雌雄蕊中高表达，在自交亲和甜荞雌雄蕊中低表达，甚至不表达，表明这些基因可能参与调控了甜荞育性.

3.   讨论

  • 甜荞营养价值较高，且耐旱耐瘠薄，是高寒山区重要的经济和药食兼用作物. 甜荞的分子研究起步较晚，研究基础较薄弱. 高通量测序技术的出现，为非模式植物基因组学研究带来了更多的新方法和新方案. 基于高通量测序的从头(donovo)转录组分析可在非模式植物中有效地用于新基因的发现和新分子标记的开发. 目前，已有研究对甜荞籽粒[23-24]、根[25-26]、根和叶[27]、花序[28]、花和子叶[29]进行了转录组测序，并挖掘了相关基因，但有关甜荞自交不亲和性及花柱发育相关的基因鲜有报道. 本研究分别对盛花期自交亲和甜荞品种甜自21等柱花雌雄蕊、自交不亲和甜荞品种乌克兰大粒荞长柱花雌雄蕊和短柱花雌雄蕊进行转录组测序，得到了283，170个Unigenes，其平均长度为1 351 bp，N50为2 040 bp，获得的高质量转录组拼接结果及表达谱数据，将对甜荞育性和花柱发育相关遗传研究提供非常有用的数据资源.

    萜类化合物也称类异戊二烯，是由两个异构的5碳骨架(异戊烯基二磷酸(IPP)和二甲基烯丙基二磷酸(DMAPP))组成的最大类别的天然化合物，萜类化合物在自然界中普遍存在，具有各种结构和功能[30]. 萜烯通过参与植物的主要代谢，例如植物激素(脱落酸，细胞分裂素，赤霉素和油菜素甾体)、光合色素(叶绿素和类胡萝卜素)、电子载体(质体醌和泛醌)和内膜系统，对植物的生长和发育至关重要[31-32]. 当植物开花时，一些低分子量的萜类化合物被释放以吸引昆虫进行授粉[33]. 在本研究中，通过对差异基因进行KEGG通路分析发现，自交不亲和甜荞与自交亲和甜荞雌雄蕊(L vs H，S vs H)差异表达基因显著富集到倍半萜和三萜生物合成途径，该途径共检测到79个差异表达基因，42个DEGs在长柱花(L-型)雌雄蕊和短柱花(S-型)雌雄蕊中同时上调表达，37个DEGs同时下调表达，由此推断萜类化合物可能参与调控了甜荞育性发育.

    植物激素信号转导对于激素触发的生化变化非常重要[18]. 牵牛花授粉诱导多种植物激素，引起生理和分子变异导致其衰老[34]. 植物生长调节剂生长素基本上参与细胞所有发育过程，包括细胞分裂和扩增，对所有植物发育过程都有影响[35]. 生长素可能直接或间接参与烟草花蕊识别和花粉管伸长[36]，在可可[37]、矮牵牛[38]和欧洲橄榄[39]等植物的自交不亲和反应中有重要作用. 在本研究中，植物激素信号转导途径(Ko04075)中的AUX/IAA与色氨酸合成密切相关，该途径中与自交亲和甜荞等柱花雌雄蕊相比，AUX1，TIR1，AUX/IAA，ARF，CH3和SAUR等基因在自交不亲和甜荞长柱花雌雄蕊中均上调表达，推测AUX/IAA相关基因可能参与了甜荞花柱发育.

    茉莉酸(JA)在植物发育过程中起着重要的调节作用，如根的生长[40]、雄蕊的发育[41]、毛的形成[42]、开花[43]、叶片的衰老[44]和顶端细胞的形成[45]，以及控制对非生物和生物胁迫的不同防御反应. 在拟南芥中，JA诱导R2R3 MYB转录因子MYB21以及MYB24[46]，与bHLH转录因子第二亚组(MYC2，MYC3，MYC4和MYC5)[40]来调节后期雄蕊发育[47]. Shi等[48]的研究结果表明，自花授粉梨树花柱中的JA浓度增加，而异花授粉花柱中的JA浓度显著降低. 结合外源JA处理，Shi等[48]认为S-RNase的表达水平降低与JA信号级联有关. 本研究中α-亚麻酸代谢途径在茉莉酸信号作用下，MYC2在自交不亲和甜荞长、短柱花雌雄蕊中上调表达，自交亲和甜荞、雌雄蕊中下调表达. 结果表明JA信号对bHLH转录因子MYC2有影响，从而参与调控甜荞自交亲和性.

    赤霉素(GA)是植物中许多发育过程必不可少的植物激素，包括种子发芽、茎伸长、叶片膨胀、毛状体发育、花粉成熟和开花诱导[49]. 植物中存在的具有生物活性的GA较少，主要包括GA1，GA3，GA4和GA7，是碱性二萜类羧酸骨架的衍生物[50]. DELLA是GA反应的关键细胞内阻遏物，抑制种子发芽、生长和几乎所有已知的GA依赖性过程，而GA则减轻其抑制活性[49]. DELLA是植物特异性GRAS家族中假定的转录调节因子的子集，拟南芥中有5个DELLA，它们在抑制GA反应中起着独特或重叠的功能[51-52]，其中RGA和GAI抑制营养生长和花诱导[53]，RGA，RGL1和RGL2共同调节花的发育[52]. 本研究的Ko04075途径，在GA作用下二萜化合物生物合成通路中，与自交亲和甜荞相比TF，DELLA和GID1基因在自交不亲和甜荞长、短柱花雌雄蕊中均上调表达. 因此，本研究推测赤霉素通过影响萜类化合物的合成，进而参与调控甜荞自交不亲和性.

4.   结论

  • 本研究获得了自交不亲和甜荞长、短柱花和自交亲和甜荞等柱花雌雄蕊的基因表达谱数据，共获得283 170条Unigenes，其中有68 813个差异表达基因(DEGs). 自交不亲和甜荞中长、短柱花雌雄蕊中的差异基因主要富集在导管形成和导管发育过程中；自交不亲和甜荞(S-型和L-型)与自交亲和甜荞(H-型)雌雄蕊中差异基因主要富集在倍半萜和三萜化合物合成途径中；自交不亲和甜荞长、短柱花雌雄蕊中差异基因(L vs S)主要富集在植物激素信号转导途径中，生长素、茉莉酸和赤霉素信号通路中的差异表达基因在自交不亲和型长、短柱花雌雄蕊中均高表达；ARF，MYC，MADS和U-box超基因家族中某些基因参与了甜荞育性发育.

参考文献 (53)

目录

版权所有：西南大学学报（自然科学版）

地址：重庆市北碚区西南大学 电话：023-68254576, 68252540, 68366165,68366917 电子邮箱：lkxb@swu.edu.cn

本系统由北京仁和汇智信息技术有限公司开发

/

返回文章
返回