摘要
为研究陆地棉查尔酮异构酶CHI(chalcone isomerase)在彩色棉纤维发育及响应胁迫中的作用,对陆地棉基因组中GhCHI基因家族进行鉴定,开展蛋白理化性质、结构域、启动子顺式作用元件、蛋白质高级结构、系统进化分析,利用转录组数据和实时荧光定量PCR分析其表达特征。结果显示,从陆地棉基因组中鉴定获得12个GhCHI基因家族成员,可分为2个亚家族,具有典型的查尔酮超家族结构域,编码201~452个氨基酸,主要为亲水性蛋白,二级结构主要以α-螺旋和无规则卷曲为主。GhCHI基因启动子顺式作用元件类型主要包含光反应元件、激素响应元件及参与类黄酮生物合成调控的元件等。表达分析结果显示GhCHI1、GhCHI2、GhCHI3、GhCHI4与纤维发育密切相关,尤其在彩色棉纤维发育后期的色素沉积着色时期表达水平较高;这4个基因在叶、花托、雌蕊中也具有较高的表达水平;GhCHI1、GhCHI2、GhCHI3在热胁迫、盐胁迫和干旱胁迫下均受到显著的诱导表达。GhCHI1―GhCHI4蛋白互作分析结果显示,它们可能和参与类黄酮合成的2-氧代戊二酸3-双加氧酶和类黄酮3'-单加氧酶等蛋白质存在相互作用。结果表明,GhCHI1、GhCHI2、GhCHI3和GhCHI4基因在彩色棉的纤维发育和胁迫应答中发挥重要作用。
随着社会发展和生活水平的不断提高,人们的绿色消费观念不断加深,对绿色纺织品的需求量不断增大。作为陆地棉的重要组成部分,天然彩色棉纤维是绿色纺织品的重要原料来源。目前生产上天然彩色棉纤维主要包括棕色棉纤维和绿色棉纤维,与白色棉纤维相比,彩色棉纤维存在质量略差、着色不稳定、色泽不均匀等问题,制约了天然彩色棉产品的应用和推
查尔酮异构酶(CHI)是黄酮类化合物合成途径的关键限速酶,决定了植物细胞内黄酮类化合物的产
目前,关于棉花GhCHI基因已有部分研究,宋成攀
从Pfam数据库(https://pfam-legacy.xfam.org/)下载棉花CHI隐马尔科夫模型Chalcone (PF02431) 作为种子序列,从CottonGen数据库(https://www.cottongen.org/)和CottonFGD数据库(https://cottonfgd.net/)下载棉花基因组数据,对这些基因组数据进行分析筛选GhCHI基因。
1) 陆地棉GhCHI基因家族成员鉴定及其理化性质分析。利用TBtools中的Simple HMM Search工具,确定陆地棉GhCHI基因家族成员的基因ID,通过TBtools中的Fasta Extract工具提取陆地棉GhCHI基因家族成员的蛋白序列。将提取序列提交Clustal Omega(https://www.ebi.ac.uk/Tools/msa/clustala)进行多序列比对;利用ProtParam tool(https://web.expasy.org/)对鉴定蛋白结构域后的基因序列进行理化性质分析。利用MEME(http://alter-nate.meme-suite.org/tools/meme)在线分析棉花GhCHI1―GhCHI12蛋白保守基序。利用NCBI中BatchCD Search工具分析GhCHI1―GhCHI12蛋白结构域并通过TBtools处理得到可视化图形。
2) 陆地棉GhCHI基因启动子分析及系统发育树构建。利用TBtools中的GXF Sequences Extract工具,从基因组文件中提取家族基因起始密码子上游2 000 bp的核酸序列,经Fasta Stats验证后,再用Fasta Extract提取启动子序列。将筛选后的陆地棉GhCHI基因家族的启动子序列提交至PlantCARE(https://bioinformatics.psb.ugent.be/webtools/plantcare/html/)分析顺式作用元件组成,再利用TBtools的Basic Biosequence View进行可视化。利用Mega软件通过邻接法(Neighbor-Joining)构建陆地棉GhCHI与拟南芥AtCHI和水稻OsCHI(蛋白序列来源于JGI,https://phytozome-next.jgi.doe.gov/)系统发育树。
3)GhCHI蛋白质结构预测及蛋白相互作用分析。利用SOMPA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)、SWISS-MODEL(https://swissmodel.expasy.org/)对陆地棉GhCHI蛋白序列进行二级结构、三级结构分析。利用STRING数据库(https://cn.string-db.org/)在线分析基因间的互作关系。
4)陆地棉GhCHI基因家族成员在棉花不同组织、不同胁迫处理以及不同时期多种纤维中的表达分析。利用R包Heatmap,以Log2(FPKM+1)对棉花在开花后不同发育期的胚珠和纤维材料 0 d(胚珠)、5 d(胚珠+纤维)、10 d(纤维)、15 d(纤维)、20 d(纤维)转录组数据进行均一化,分析GhCHI基因在纤维发育不同时期的表达情况。参考NCBI数据库中经过0、1、3、6、12 h冷、热、盐、干旱胁迫处理的不同转录组数据分析GhCHI基因在棉花不同胁迫处理及棉花不同组织中的表达变化。
利用天根生化科技(北京)有限公司的多糖多酚植物总RNA提取试剂盒(DP441,北京)和cDNA第一链合成预混试剂(KR211,北京)提取棕色棉、白色棉和绿色棉不同发育时期的纤维样品(开花后0、5、10、15和20 d)的总RNA,反转录获得cDNA。使用 Primer Premier 5.0软件设计ChCHI基因特异性引物(
基因名称/基因编号 Gene name/Gene ID | 引物名称 Primer name | 引物序列(5′- 3′) Primer sequence (5′- 3′) |
---|---|---|
GhCHI1/Gh_A13G0197 |
GhCHI1-F GhCHI1-R | CCACCAAAACCCTGTTCCTT |
ACTCCACCGCGCTCTTACCC | ||
GhCHI2/Gh_D13G0211 |
GhCHI2-F GhCHI2-R | AGCCTCCTGGTTCCACCAAA |
CACCGCGCTCTTACCCTTCC | ||
GhCHI3/Gh_A05G3491 |
GhCHI3-F GhCHI3-R | CCTTTTCCTCCACAAATCAC |
ATCTTCAGCTAGCACGTTCC | ||
GhCHI4/Gh_D04G1926 |
GhCHI4-F GhCHI4-R | TTCACGACGGAGGGAAAGGA |
GCAAAGCTCAGCGGAGAGGG | ||
GhCHI5/Gh_D13G1602 | GhCHI5-F | CAAATTAACTCAAACCCATT |
GhCHI5-R | CAAAAACTTTCTCCCTGTAC | |
GhCHI6/Gh_A06G0722 | GhCHI6-F | GTGGCAGTCTATCTTTTCTA |
GhCHI6-R | TGATGATCTTCATGGTTCTT | |
GhCHI7/Gh_D08G2001 | GhCHI7-F | GAACCACCCGTTCTTCCAAC |
GhCHI7-R | GCTTTTTTTCCTCAACCCAA | |
GhCHI8/Gh_D06G0839 | GhCHI8-F | CGACACATTCCTTTGCTATC |
GhCHI8-R | GCATCCTTGACCAACTTCAC | |
GhCHI9/Gh_A08G2369 | GhCHI9-F | TTGCTGTCGTCTCGCAGAAT |
GhCHI9-R | CGGGAAAGAAACCCCACTCT | |
GhCHI10/Gh_D07G0160 | GhCHI10-F | CTTCTGAGGTGCTTGTTGGT |
GhCHI10-R | AAAGTCATTGTGTTTGTTCG | |
GhCHI11/Gh_A07G2350 | GhCHI11-F | CTTCTGAGGTGCTTGTTGGT |
GhCHI11-R | AAAGTCATTGTGTTTGTTCG | |
GhCHI12/Gh_D08G2008 | GhCHI12-F | CTACCTTTTTCCAATTTTCT |
GhCHI12-R | CTATATCCCCACTTCACACC | |
GhUBQ7/DQ116441 | GhUBQ7-F | GAAGGCATTCCACCTGACCAAC |
GhUBQ7-R | CTTGACCTTCTTCTTCTTGTGCTTG |
从陆地棉基因组数据库中鉴定获得12个GhCHI基因家族成员(GhCHI1―GhCHI12),陆地棉GhCHI基因的开放阅读框由1 013~7 460个核苷酸组成,编码的蛋白质包含201~451个氨基酸、等电点为4.86~9.38,GhCHI蛋白主要为亲水蛋白,α-螺旋和无规则卷曲是GhCHI蛋白二级结构的主要形式(
基因ID Gene ID | 基因名称 Gene name | 核苷酸数 Nucleotide number | 氨基酸数量 Amino acid number | 等电点 Isoelectric point | 亲水指数 Grand average of hydropathicity | α-螺旋/% α-Helix | 延伸链/% Extended strand | β-折叠/% β-Turn | 无规则 卷曲/% Random coil |
---|---|---|---|---|---|---|---|---|---|
Gh_A13G0197 | GhCHI1 | 2 750 | 227 | 5.02 | -0.043 | 43.61 | 19.38 | 9.25 | 27.75 |
Gh_D13G0211 | GhCHI2 | 3 097 | 227 | 4.92 | -0.058 | 40.53 | 18.50 | 9.25 | 31.72 |
Gh_A05G3491 | GhCHI3 | 1 056 | 216 | 4.86 | -0.090 | 43.52 | 24.54 | 7.41 | 24.54 |
Gh_D04G1926 | GhCHI4 | 1 013 | 203 | 4.98 | -0.172 | 46.80 | 21.18 | 7.39 | 24.63 |
Gh_D13G1602 | GhCHI5 | 1 575 | 278 | 9.13 | -0.063 | 33.45 | 19.42 | 2.52 | 44.60 |
Gh_A06G0722 | GhCHI6 | 2 881 | 432 | 8.52 | -0.045 | 38.43 | 16.67 | 5.09 | 39.81 |
Gh_D08G2001 | GhCHI7 | 2 410 | 297 | 9.13 | -0.154 | 40.07 | 19.87 | 5.72 | 34.34 |
Gh_D06G0839 | GhCHI8 | 2 867 | 432 | 8.71 | -0.080 | 36.8 | 16.90 | 4.86 | 41.44 |
Gh_A08G2369 | GhCHI9 | 7 460 | 451 | 9.38 | -0.116 | 38.58 | 21.06 | 5.54 | 34.81 |
Gh_D07G0160 | GhCHI10 | 2 557 | 439 | 6.95 | -0.031 | 39.64 | 15.26 | 5.01 | 40.09 |
Gh_A07G2350 | GhCHI11 | 6 264 | 243 | 7.00 | -0.299 | 28.40 | 19.34 | 6.58 | 45.68 |
Gh_D08G2008 | GhCHI12 | 1 056 | 201 | 8.90 | 0.131 | 42.79 | 23.38 | 6.47 | 27.36 |

图1 棉花GhCHI基因的基因结构(A)与保守基序(B)分析
Fig.1 Gene structure (A) and conserved motif (B) analysis of cotton GhCHI genes

图2 棉花GhCHI的多序列比对和保守结构域分析
Fig.2 Multiple sequence alignment and conserved domain analysis of GhCHI
为了研究棉花GhCHI基因家族在进化过程中的亲缘关系,选取拟南芥和水稻CHI蛋白与陆地棉GhCHI蛋白构建系统进化树。结果显示,这些CHI可以分为3个组:Group 1(OsCHI2、OsCHI4、GhCHI3、GhCHI4)、Group 2(AtCHI1-AtCHI8、AtCHI10、AtCHI12)、Group 3(OsCHI1、OsCHI3、AtCHI9、GhCHI1、GhCHI2、GhCHI5-GhCHI12)。Group 1含有2个棉花CHI,Group 2 均为拟南芥CHI,都是具有物种特异性的分支,Group 3中的CHI在棉花中具有更紧密的进化关系(

图3 棉花与拟南芥、水稻CHI基因系统进化分析
Fig.3 Phylogenetic tree of CHI gene family members of cotton, rice and Arabidopsis
红色、绿色、蓝色星号分别代表棉花、拟南芥、水稻的CHI蛋白。Red, green and blue asterisks represent CHI proteins in cotton, Arabidopsis and rice, respectively.
选取棉花GhCHI基因起始密码子上游2 000 bp的启动子序列,利用PlantCARE分析基因启动子序列中顺式作用元件的组成,结果显示,棉花GhCHI基因启动子中包含大量的光响应元件,尤其是GhCHI10启动子含有16个光响应元件,暗示了它们与光响应之间的密切联系。此外,大多数GhCHI启动子还含有脱落酸响应元件、茉莉酸响应元件以及厌氧诱导顺式作用元件。GhCHI1、GhCHI6、GhCHI7、GhCHI9的启动子中均含有参与类黄酮生物合成基因调控的MYB结合位点,推测这些基因在类黄酮化合物生物合成中发挥重要作用(

图4 GhCHI基因启动子顺式作用元件分析
Fig.4 Constituent analysis of cis-elements of GhCHI promoters
A:MYB结合位点 MYB binding site;B:参与类黄酮生物合成基因调控的MYB结合位点 MYB binding sites involved in gene regulation of flavonoid biosynthesis;C:玉米醇溶蛋白代谢调节顺式调节元件 Element for regulating corn gliadin metabolism;D:光反应顺式作用元件 Photoreactive response element;E:脱落酸反应顺式作用元件 Abscisic acid response element;F:水杨酸响应元件 Salicylic acid responsive element;G:细胞周期调控顺式作用元件 Elements in cell cycle regulation;H:与光反应相关的MYB结合位点 MYB binding sites related to photoreaction;I:赤霉素反应元件 Gibberellin reaction element;J:分生组织表达调控顺式作用元件 Expression and regulation of element in meristem tissue;K:茉莉酸响应元件 Jasmonic acid responsive element;L:低温响应元件 Low temperature responsive element;M:叶肉细胞分化元件 Mesophyll cell differentiation element;N:无氧诱导顺式作用元件 Anaerobic induced element;O:防御和应激反应响应元件 Defense and stress response element;P:参与干旱诱导相关的MYB结合位点 MYB binding sites involved in drought induction;Q:与胚乳表达相关的顺式作用元件 Element related to endosperm expression;R:昼夜节律调控顺式作用元件 Circadian rhythm regulation element;S:生长素反应 Auxin reaction.
选取彩色棉和白棉纤维不同发育时期(开花后0、5、10、15、20 d)、不同组织、不同胁迫处理的转录组数据进行GhCHI基因的表达特征分析。结果显示,GhCHI1、GhCHI2、GhCHI3、GhCHI4四个基因在棉纤维发育过程中具有较高的表达,在白色棉纤维快速伸长发育时期的10 d累积,在棕色棉和绿色棉纤维中主要在10~20 d高表达,表明这些基因在纤维发育的次生壁发育时期具有重要作用(

图5 GhCHI基因参与纤维发育的表达分析
Fig.5 Expression analysis of GhCHI genes involved in cotton fiber development

图6 GhCHI基因组织特异性表达分析
Fig.6 Expression analysis of GhCHI genes in different cotton tissues

图7 GhCHI基因响应非生物胁迫的表达分析
Fig.7 Expression analysis of GhCHI genes in response to abiotic stress
GhCHI基因在棕色、白色和绿色棉纤维的不同发育时期(开花后0、5、10、15、20 d)的表达情况如

图8 GhCHI基因在棉花纤维发育中的表达特征分析
Fig.8 Expression analysis of GhCHI genes in cotton fiber development
不同小写字母表示组间具有显著差异。Different lowercase letters indicate significant difference between different groups(P<0.05).
棉花GhCHI蛋白三级结构分析显示,α-螺旋、无规则卷曲、延伸链以及β-折叠在GhCHI蛋白中均有分布,且以α-螺旋和无规则卷曲为主,不同的蛋白含有不同数目的螺旋-转角-螺旋结构(

图9 GhCHI蛋白三级结构分析
Fig.9 Tertiairy structure analysis of GhCHI proteins

图10 GhCHI蛋白质相互作用分析(A)及互作蛋白的表达分析(B)
Fig.10 GhCHI protein interaction analysis (A) and interaction protein expression analysis(B)
A0A1U8MSF4:GhCHI1;A0A1U8HHL4:GhCHI2;D6N3G6:GhCHI3;A0A1U8M0E4:GhCHI4;A0A1U8K257:X1亚型2-氧代戊二酸3-双加氧酶 X1 subtype 2-oxoglutarate 3-dioxygenase (Gh_D11G1775);A0A1U8K7C8:2-氧代戊二酸3-双加氧酶 2-Oxoglutarate 3-dioxygenase (Gh_A11G1616);A0A1U8NBV3:X1亚型2-氧代戊二酸3-双加氧酶 X1 subtype 2-oxoglutarate 3-dioxygenase (Gh_D12G0566);A0A1U8NF28:类黄酮3'-单加氧酶 Flavonoids 3'-monooxygenase (Gh_D12G1798);A0A1U8NS60:类黄酮3'-单加氧酶 Flavonoids 3'-Monooxygenase (Gh_A12G2650).
查尔酮异构酶CHI是植物细胞内调控类黄酮化合物合成的关键酶,在植物生长发育和响应环境胁迫过程中发挥着重要作用。目前已在红花(Carthamus tinctorius
CHI基因的表达与植物发育和逆境胁迫响应密切相关。如甘薯中CHI基因的表达在盐和干旱胁迫下显著累
GhCHI基因在天然彩色棉纤维发育后期的次生壁加厚时期具有较高的表达(
参考文献References
LIU H F,LUO C,SONG W,et al.Flavonoid biosynthesis controls fiber color in naturally colored cotton[J/OL] PeerJ,2018,6:e4537 [2023-04-24].https://doi.org/10.7717/peerj.4537. [百度学术]
CHEN M Y,ZHANG T T,HE L,et al.Qualitative analysis of cotton fiber pigment composition[J].Textile research journal,2021,91(3/4):456-463. [百度学术]
MCKHANN H I,HIRSCH A M.Isolation of chalcone synthase and chalcone isomerase cDNAs from alfalfa (Medicago sativa L.):highest transcript levels occur in young roots and root tips[J].Plant molecular biology,1994,24(5):767-777. [百度学术]
GAO J F,SHEN L,YUAN J L,et al.Functional analysis of GhCHS,GhANR and GhLAR in colored fiber formation of Gossypium hirsutum L.[J/OL].BMC plant biology,2019,19(1):455 [2023-04-24].https://doi.org/10.1186/s12870-019-2065-7. [百度学术]
麻新妍,南春利,薛永常.植物查尔酮异构酶的结构与功能研究概述[J].生物学教学,2022,47(1):2-4.MA X Y,NAN C L,XUE Y C. Review on the structure and function of plant chalcone isomerase[J].Biology teaching,2022,47(1):2-4 (in Chinese). [百度学术]
薛永常,麻新妍.黄芩查尔酮异构酶基因的全长cDNA克隆及生物信息学分析[J].中国野生植物资源,2022,41(6):50-53,96.XUE Y C,MA X Y.Cloning and bioinformatic analysis of chalcone isomerase gene full-length cDNA from Scutellaria baicalensis[J].Chinese wild plant resources,2022,41(6):50-53,96(in Chinese with English abstract). [百度学术]
NISHIHARA M,NAKATSUKA T,YAMAMURA S.Flavonoid components and flower color change in transgenic tobacco plants by suppression of chalcone isomerase gene[J].FEBS Letters,2005,579(27):6074-6078. [百度学术]
李国庆.黄酮醇在小麦抗旱节水中的功能分析[D].济南:山东大学,2021.LI G Q.Function analysis of flavonols in drought resistance and water saving of wheat[D].Ji’nan:Shandong University,2021 (in Chinese with English abstract). [百度学术]
宁亚茹,韩民利,张晓东,等.盐胁迫对黄蜀葵不同部位总黄酮含量及抗氧化活性的影响[J].江苏农业科学,2021,49(19):196-200.NING Y R,HAN M L,ZHANG X D,et al.Effects of salt stress on total flavonoids content and antioxidant activity in different parts of Abelmoschus manihot[J].Jiangsu agricultural sciences.,2021,49(19):196-200 (in Chinese). [百度学术]
宋静武,殷德松,赵弟广,等.核桃叶片内多酚黄酮类成分对低温胁迫的响应[J].河北林果研究,2017,32(1):34-41.SONG J W,YIN D S,ZHAO D G,et al.The response of polyphenol flavonoids in Juglans regia leaves to low temperature stress[J].Hebei journal of forestry and orchard research,2017,32(1):34-41 (in Chinese with English abstract). [百度学术]
李军,李洪清,李美茹.大花美人蕉查尔酮异构酶基因的cDNA克隆和序列分析[J].植物生理学通讯,2006,42(3):449-453.LI J,LI H T,LI M R.Cloning and sequence analysis of cDNA of chalcone isomerase gene from Canna generalis Bailey[J].Plant physiology communications,2006,42(3) :449-453 (in Chinese with English abstract). [百度学术]
段怡红.陆地棉开花期相关基因GhMYB44的功能研究[D].保定:河北农业大学,2020.DUAN Y H.Functional identification of flowering-stage relatedgene GhMYB44 from Gossypium hirsutum[D].Baoding:Hebei Agricultural University,2020 (in Chinese with English abstract). [百度学术]
ZHU J H,ZHAO W,LI R S,et al.Identification and characterization of chalcone isomerase genes involved in flavonoid production in Dracaena cambodiana[J/OL].Frontiers in plant science,2021,12:616396[2023-04-24]. https://doi.org/10.3389/fpls.2021.616396. [百度学术]
丁宁,海燕,王晓晖,等.白木香查尔酮异构酶基因的克隆鉴定与表达分析[J].药学学报,2021,56(2):630-638.DING N,HAI Y,WANG X H,et al.Cloning and expression analysis of chalcone isomerase from Aquilaria sinensis[J].Acta pharmaceutica sinica,2021,56(2):630-638 (in Chinese with English abstract). [百度学术]
宋成攀,夏松波,王孝刚,等.嫁接陆地棉查尔酮合成酶与查尔酮异构基因的克隆及表达分析[J].湖北农业科学,2017,56(23):4616-4623.SONG C P,XIA S B,WANG X G,et al.Cloning and expression analysis of chalcone synthase and chalcone isomer genes in grafted upland cotton [J].Hubei agricultural science,2017,56 (23):4616-4623 (in Chinese with English abstract). [百度学术]
任超翔,唐小慧,何雯,等.红花查尔酮异构酶基因的克隆及表达分析[J].天然产物研究与开发,2018,30(9):1521-1525,1574.REN C X,TANG X H,HE W,et al.Cloning and expression analysis of chalcone isomerase gene in safflower[J] Natural product resarch and development,2018,30(9):1521-1525,1574 (in Chinese with English abstract). [百度学术]
PARK S I,PARK H L,BHOO S H,et al.Biochemical and molecular characterization of the rice chalcone isomerase family[J/OL].Plants,2021,10(10):2064[2023-04-24].https://doi.org/10.3390/plants10102064. [百度学术]
CHAO N,WANG R F,HOU C,et al.Functional characterization of two chalcone isomerase (CHI) revealing their responsibility for anthocyanins accumulation in mulberry[J].Plant physiology and biochemistry,2021,161:65-73. [百度学术]
张晓冬,尹彦超,周姗,等.光果甘草查尔酮异构酶基因的克隆及序列分析[J].生物技术通讯,2017,28(4):478-484.ZHANG X D,YIN Y C,ZHOU S,et al.Cloning and sequence analysis of chalcone isomerase gene from Glycyrrhiza glabra L.[J].Letters in biotechnology,2017,28(4):478-484 (in Chinese with English abstract). [百度学术]
史敏莉,唐锐敏,张毅,等.甘薯lncRNA TCONS_00074371及其靶基因CHI的克隆及生物信息学分析[J].山西农业科学,2022,50(12):1599-1607.SHI M L,TANG R M,ZHANG Y,et al.Cloning and bioinformatics analysis of sweet potato lncRNA TCONS_00074371 and its target gene CHI [J].Journal of Shanxi agricultural science,2022,50(12):1599-1607(in Chinese with English abstract). [百度学术]
周平,林志楷,郭瑞,等.低温处理对桃树叶片基因表达及类黄酮合成代谢的影响[J].农业生物技术学报,2021,29(7):1283-1294.ZHOU P,LIN Z K,GUO R,et al.Effects of low temperature treatment on gene expression and flavonoids biosynthesis metabolism in peach (Prunus persica) leaves[J].Journal of agricultural biotechnology,2021,29(7):1283-1294 (in Chinese with English abstract). [百度学术]
WANG Y,SHI Y,LI K,et al.Roles of the 2-oxoglutarate-dependent dioxygenase superfamily in the flavonoid pathway:a review of the functional diversity of F3H,FNS I,FLS,and LDOX/ANS[J/OL].Molecules,2021,26(21):6745 [2023-04-24].https://doi.org/10.3390/molecules26216745. [百度学术]
LV L L,FENG X F,LI W,et al.High temperature reduces peel color in eggplant (Solanum melongena) as revealed by RNA-seq analysis[J].Genome,2019,62(7):503-512. [百度学术]
ZHAO D Q,TAO J,HAN C X,et al.Flower color diversity revealed by differential expression of flavonoid biosynthetic genes and flavonoid accumulation in herbaceous peony (Paeonia lactiflora Pall.)[J].Molecular biology reports,2012,39(12):11263-11275. [百度学术]