Variation and correlation analysis of florescence and main ornamental traits of Chrysanthemum × morifolium Ramat. varieties
-
摘要:
以菊花（Chrysanthemum × morifolium Ramat.）大菊品种自然群体及7个杂交群体为实验材料，分析比较了不同品种开花期和7个主要观赏性状的变异特点及其相关性，为筛选杂交亲本和选育早花优良品种提供依据。结果显示：大菊自然群体的8个性状中，除开花期外，其余性状变异系数均高于15%；相关性分析结果表明，开花期和封顶叶片数之间呈显著负相关；开花持续期与开花期呈极显著负相关，与舌状花轮数呈极显著正相关。杂交群体的8个性状中，除开花期外，变异系数均高于15%，但各杂交群体的开花期性状变异系数明显高于自然群体；杂交群体Ⅲ和Ⅴ的瓣型与花色更为丰富，且花期也较早。相关性分析结果表明，开花期和瓣型、舌状花轮数、封顶叶片数之间均呈极显著负相关，和花色之间呈显著负相关；开花持续期和舌状花轮数、舌状花数量之间呈极显著正相关，和瓣型之间呈极显著负相关，和花色之间呈显著负相关；开花期和开花持续期之间呈极显著正相关。研究结果表明，通过有性杂交可以获得开花期变异更广泛的子代植株，从而筛选开花期更早的大菊新品种。
Abstract:Using natural and hybrid populations of chrysanthemum as experimental materials, we analyzed the variation characteristics and correlations of florescence and ornamental characters among different varieties, then compared and analyzed the results to provide a basis for screening hybrid parents and breeding good varieties of early-flowering chrysanthemums. Results showed that the coefficients of variation of the eight traits in the natural chrysanthemum population were higher than 15%, except for flowering stage traits. Analysis showed a significant negative correlation between flowering stage and number of stem leaves. There was a highly significant negative correlation between flowering duration and flowering stage, and a highly significant positive correlation between flowering duration and number of ray flower turns. Among the eight traits in the hybrid populations, the coefficients of variation were higher than 15%, except for flowering traits, but the coefficients of variation for florescence traits of each hybrid population were also significantly higher than that of the natural population. Compared with other hybrid populations, the petal types and flower colors of hybrid populations Ⅲ and Ⅴ were more abundant, with earlier flowering stages. Analysis showed a highly significant negative correlation between flowering stage and petal type, number of ray flower turns, and number of stem leaves, and a significant negative correlation between flowering stage and flower color. There was a highly significant positive correlation between flowering duration and number of tubular flower turns and number of tubular flowers, a highly significant negative correlation between flowering duration and petal type, and a significant negative correlation between flowering duration and flower color. There was a highly significant positive correlation between flowering stage and flowering duration. Therefore, progeny plants with wider variation in flowering stage can be obtained through sexual hybridization, with new chrysanthemum varieties showing earlier flowering. Early flowering parents can be selected according to the characteristics of fewer leaves, more ray flower turns, petals close to flat petals, and flower color close to orange and pink lines. Sexual hybridization with varieties exhibiting better flower type and color can be carried out to obtain better early flowering varieties of Chrysanthemum with stable flowering stage from late September to early October.
-
Keywords:
- Chrysanthemum /
- Florescence traits /
- Main ornamental traits /
- Variation /
- Correlation
-
地黄（Rehmannia glutinosa Libosch.）为玄参科地黄属植物，多年生草本，以块根入药，为著名的“四大怀药”之一。地黄最早记载见《神农本草经》，被列为上品，迄今已有2000余年的应用历史。根据炮制方法的不同，地黄药材分为鲜地黄、生地黄和熟地黄[1]。2020版《中国药典》记载，鲜地黄具有清热生津、凉血、止血的功效，生地黄具有清热凉血、养阴生津的功效，熟地黄具有补血滋阴，益精填髓的功效[2]。地黄富含环烯醚萜类、苯乙醇苷类、紫罗兰酮类、三萜类、黄酮类和糖类等，对人体心脑血管、血液、中枢神经和免疫系统等均有显著作用[3]。
毛蕊花糖苷（Acteoside）是地黄中含量较高的苯乙醇苷类化合物[4]，具有抗氧化、免疫调节、抗炎、保肝、抗肿瘤、增强记忆力等生物活性[5]，是2010版、2015版《中国药典》规定的地黄药材质量控制的指标性成分之一。地黄毛蕊花糖苷的含量易受品种[6]、产地[7]、收获时期[5]、种植密度[8]和光照条件[9]等因素的影响，造成某些年份部分地黄药材的毛蕊花糖苷含量达不到《中国药典》规定的要求。地黄毛状根中也含有丰富的毛蕊花糖苷，本课题组前期研究表明，在毛状根诱导的特定时期添加水杨酸（SA）可显著促进毛蕊花糖苷的含量[10]。然而，地黄生长发育过程中叶面喷施SA对毛蕊花糖苷的含量是否有影响还未见报道。
本研究以大田栽培的地黄为材料，采用叶面喷施方法分析SA对地黄叶片和块根毛蕊花糖苷含量的影响，并利用转录组测序技术分析地黄块根中的基因表达特征，研究结果旨在为生产中应用外源激素提高毛蕊花糖苷的含量提供理论依据。
1. 材料与方法
1.1 实验材料
供试材料为地黄‘温85-5’，种植在河南省武陟县蔡庄村（35°2′51″N，113°18′34″E），经河南农业大学王丰青教授鉴定为Rehmannia glutinosa Libosch.。将浓度为100 μmol/L的SA水溶液均匀喷施在生长180 d处于膨大后期长势一致的地黄叶片上，以叶片完全湿润且无水滴落下为准。分别于处理后1、3和6 h进行取样，取样时选择位置相同的叶片。对照组喷施蒸馏水，与处理组同时取样。样品清洗干净后，将块根切成1 cm3左右的小块，叶片剪碎，55℃烘干后打粉，过3号药典筛。样品粉末存放于干净的8号自封袋中备用。
1.2 毛蕊花糖苷含量测定
对照品毛蕊花糖苷（批号MUST-18032725）购于成都曼斯特科技有限公司，纯度均 ≥ 98%。采用 Agilent1260高效液相色谱仪（美国安捷伦科技有限公司）进行检测。
色谱条件：采用的色谱柱型号为 Dikma Diamonsil C18（4.6 × 250 mm，5 μm），柱温30℃，流速1 min/mL。毛蕊花糖苷的流动相为乙腈−0.1%醋酸水（16 : 84)，检测波长334 nm，进样量为20 μL。
供试样品溶液制备：精密称取地黄叶和块根样品粉末0.8 g，放入锥形瓶中，精密吸取50 mL甲醇加入锥形瓶中，称重，并在65℃加热回流提取1.5 h，放凉至室温称重，用甲醇补足失重后，摇匀过滤。于蒸发皿中精密吸取滤液20 mL进行毛蕊花糖苷分析，在电热恒温水浴锅上浓缩至近干，残渣用流动相溶解，转移至5 mL容量瓶中，用流动相稀释至刻度，摇晃均匀，用0.22 μm微孔的滤膜过滤，滤液装入2 mL的进样瓶待测。
含量计算：以本实验室建立的标准曲线Y = 30024X − 110.9来计算毛蕊花糖苷含量，Y为峰面积积分值，X为样品的质量浓度。
1.3 样品RNA提取及高通量测序
用TRIzol试剂提取样品的总RNA，用核酸测定仪检测RNA的浓度和质量。使用带有Oligo dT的磁珠富集具有polyA尾巴的mRNA，然后将RNA片段化，反转录后再合成cDNA第2链，形成双链cDNA。双链cDNA经过末端修复、3′末端加A、添加测序接头、多轮扩增、热变性成单链及单链环化等一系列步骤后，完成测序文库的制备。测序委托华大基因科技有限公司进行，测序平台为BGISEQ-500。
对测序得到的原始序列（Raw reads）进行质控处理，去除低质量、接头及污染序列，获取过滤后的测序序列（Clean reads）。使用Bowtie2将Clean reads比对到课题组前期获得的地黄叶和根的参考基因序列集[11]，统计不同样品的片段序列比对率以及分布，之后再使用RSEM 计算基因的表达水平，表达量用FPKM表示。对比分析水杨酸喷施前后地黄块根中基因的表达水平，获取水杨酸喷施处理的特异响应基因，并对其进行GO注释和KEGG注释，根据注释结果进行功能分类和KEGG pathway分类，并使用R软件中的phyper函数进行富集分析，获得水杨酸喷施处理下地黄块根内的关键分子响应进程。
1.4 实时荧光定量（qRT-PCR）分析
用TaKaRa反转录试剂盒对RNA进行反转录，合成cDNA，反应体系包括1 μL oligo dT primer，1 μL dNTP Mixture，2 μg 模板RNA，加水补足体积到10 μL，65℃保温5 min后，冰上迅速冷却。再加0.5 μL的RNase抑制剂，1 μL PrimeScript Ⅱ RTase，4 μL 5 × PrimeScript Ⅱ Buffer，加水补足体积到20 μL。反应程序为42℃ 60 min，95℃ 5 min。以RgTIP41为内参基因，用实时荧光定量PCR检测基因表达水平。所用试剂盒为SYBR® Premix Ex Taq™ Ⅱ（Tli RNaseH Plus) （Takara，大连），使用仪器为Bio-Rad IQ5（上海伯乐公司）。定量反应体系为25 μL，包含2 μL 上述反转录cDNA产物，上、下游引物各1 μL，12.5 μL SYBR® Premix Ex Taq，8.5 μL ddH2O。反应程序为：95℃变性30 s；然后95℃ 5 s， 60℃ 30 s，40个循环。结束反应后获得不同样品的扩增循环数Ct，使用2－ΔΔCt法计算不同基因的相对表达量。
2. 结果与分析
2.1 SA处理对地黄毛蕊花糖苷含量的影响
对SA处理的地黄叶和块根中的毛蕊花糖苷含量进行测定，结果表明，SA能够显著提高毛蕊花糖苷的含量（图1）。在叶中， SA处理1 ~ 6 h后，毛蕊花糖苷含量分别提高了11.2% ~ 19.3%。块根中毛蕊花糖苷的提升幅度远高于叶片，分别较对照提高了0.9 ~ 1.4倍，处理 6 h后的毛蕊花糖苷含量最高，达0.5 mg/g，远超2015版《中国药典》规定的0.02%，说明地黄叶面喷施SA可显著提高叶和块根中毛蕊花糖苷的含量。
![]()
图 1 SA处理下地黄叶和块根中的毛蕊花糖苷含量A：叶中毛蕊花糖苷含量；B：块根中毛蕊花糖苷含量。* 表示处理与对照差异显著，P < 0.05。Figure 1. Acteoside content in SA-treated leaves and tuberous roots of Rehmannia glutinosaA: Acteoside content in leaves; B: Acteoside content in tuberous roots. * indicates significant difference at P < 0.05 level.2.2 地黄块根RNA测序分析
利用Agilent 2100 Bioanalyzer和Fragment Analyzer分别对提取的各样品总RNA质量进行检测，结果显示，12个样品的总RNA浓度在430 ~ 1260 ng/μL，总RNA质量在8.6 ~ 25.2 μg，RNA的浓度和总量满足建库需求。RIN值在8.2 ~ 9.9，28S/18S > 1.6，说明RNA较为完整，符合建库要求。
为了分析SA处理后地黄块根相关基因的表达特性，对SA处理1、3和6 h后的地黄块根进行RNA-seq分析，结果表明（表1），每个测序样本获得的总原始读段量均为21.94 M，去除低质量、接头污染及未知碱基N含量过高的reads，获得的高质量reads在20.98 ~ 21.27 M，碱基数均在1.05 ~ 1.06 Gb，测序数据量基本一致。将每个样品的测序数据匹配地黄参考转录组，发现匹配率在85.04% ~ 87.70%，特异匹配率在47.35% ~ 51.31%，测序数据能够较好地反映细胞中基因表达的真实情况，说明测序质量良好，可以进行后续基因表达分析。
表 1 测序数据统计结果Table 1. Statistics of sequenced data样本
Sample总原始序列
Total raw reads / M总测序序列
Total clean reads / M总测序碱基数
Total clean bases / Gb测序序列比率
Clean read ratio / %总匹配率
Total mapped / %特异匹配率
Uniquely mapped / %Control_1 21.94 21.13 1.06 96.3 87.70 51.31 Control_2 21.94 21.07 1.05 96.04 86.92 51.18 Control_3 21.94 21.08 1.05 96.08 85.72 51.09 SA 1h_1 21.94 21.27 1.06 96.93 85.44 47.35 SA 1h_2 21.94 21.08 1.05 96.05 87.07 48.03 SA 1h_3 21.94 21.15 1.06 96.37 87.28 49.26 SA 3h_1 21.94 20.98 1.05 95.59 86.85 50.36 SA 3h_2 21.94 21.05 1.05 95.93 85.74 49.92 SA 3h_3 21.94 21.05 1.05 95.95 86.81 50.47 SA 6h_1 21.94 21.09 1.05 96.11 86.50 50.25 SA 6h_2 21.94 21.06 1.05 96.00 85.04 49.42 SA6h_3 21.94 21.03 1.05 95.86 86.68 49.97 2.3 SA处理前后差异表达基因分析及筛选
为了分析SA处理后块根中基因的表达特性，采用FPKM值比较基因丰富度的差异（图2）。使用以下标准对差异表达基因（DEGs）进行识别和筛选：校正P值 < 0.001且log2值 ≥ 2。分析SA处理不同时间后与对照样品中的差异表达基因（图2：A），发现SA处理1 h后834个基因上调，2017个基因下调；处理3 h后1587个基因上调，2726个基因下调；处理6 h后1602个基因上调，2063个基因下调。进一步分析SA处理后的3个时间点与CK相比的共同差异表达基因，发现共有603个基因是共同差异表达的（图2：B），其中上调表达和下调表达的基因数分别为195（图2：C）和366个（图2：D）。
![]()
图 2 SA处理后基因显著差异表达A：SA处理过程中上调和下调的基因数目； B ~ D：SA处理后不同时间点鉴别出的总DGEs（B）、上调DGEs（C）和下调DGEs的维恩图。Figure 2. Significant DEGs in response to SA treatmentA: Up-regulated and down-regulated gene numbers during SA treatment; B − D: Venn diagram of total DEGs (B), up-regulated DEGs (C), and down-regulated DEGs (D) identified at different time points after SA treatment.对603个DGEs进行GO功能分类（图3：A），结果显示其共分为分子功能（Molecular function）、细胞功能（Cellar function）和生物功能（Biological function）3个大类。其中分子功能分类中的催化活性（Catalytic activity）占比最多，其次是ATP结合；在细胞功能分类中，细胞（Cell）、细胞膜（Membrane）、细胞膜构件（Membrane part）较多；生物功能分类中，细胞过程（Cellular process）和代谢过程（Metabolic process）所占比重最多。进一步对603个DEGs进行KEGG代谢通路富集分析，图3：B展示了最显著的前20个代谢通路，这些通路涉及各项生命活动。其中首先被富集的是苯乙醇苷生物合成通路（Phenylpropanoid biosynthesis），其次是淀粉和蔗糖代谢合成通路（Starch and sucrose metabolism）、植物MARK信号通路（MAPK signaling pathway plant）及RNA聚合酶通路（RNA polymerase）。这说明SA喷施对地黄块根中次生代谢物的积累产生了较大影响，且调控了碳水化合物和MAPK等多个代谢通路。
![]()
图 3 SA处理后3种比对均差异表达基因的GO分类及KEGG通路富集A：603个基因的GO分类；B：KEGG富集的前20个代谢通路。圆点大小和颜色分别表示通路中DEGs的数量和Q值范围。1：分子转导活性；2：分子载体活性；3：信号转导活性；4：结构分子活性；5：转运活性；6：转录调节活性；7：抗氧化活性；8：结合；9：分子功能调节；10：催化活性；11：膜封闭腔；12：超分子复合物；13：共质体；14：细胞连接；15：细胞组分；16：细胞外区域；17：膜组分；18：膜；19：细胞器组分；20：细胞器；21：大分子复合物；22：细胞；23：生殖过程；24：繁殖；25：解毒作用；26：多细胞生物过程；27：发育过程；28：多生物体过程；29：信号；30：定域化；31：生物调节；32：对刺激的反应；33：代谢过程；34：细胞过程；35：细胞成分组织或生物合成。Figure 3. GO classification and KEGG pathway enrichment of co-DEGs in three comparisonsA: GO classification of 603 genes; B: Top 20 enriched KEGG pathways among 603 genes. Size and color of dot represent number and scope of DEGs in pathway, respectively. 1: Molecular transducer activity; 2: Molecular carrier activity; 3: Signal transducer activity; 4: Structural molecule activity; 5: Transporter activity; 6: Transcription regulator activity; 7: Antioxidant activity; 8: Binding; 9: Molecular function regulator; 10: Catalytic activity; 11: Membrane-enclosed lumen; 12: Supramolecular complex; 13: Symplast; 14: Cell junction; 15: Cell part; 16: Extracellular region; 17: Membrane part; 18: Membrane; 19: Organelle part; 20: Organelle; 21: Macromolecular complex; 22: Cell; 23: Reproductive process; 24: Reproduction; 25: Detoxification; 26: Multicellular organismal process; 27: Developmental process; 28: Multi-organism process; 29: Signaling; 30: Localization; 31: Biological regulation; 32: Response to stimulus; 33: Metabolic process; 34: Cellular process; 35: Cellular component organization or biogenesis.2.4 毛蕊花糖苷合成相关催化酶基因表达分析
研究表明，在植物体内，毛蕊花糖苷由苯丙氨酸途径的咖啡酰辅酶A（Caffeoyl CoA）和酪氨酸途径的羟基酪醇苷（Hydroxtrosol glucoside）经缩合糖苷化后生成[8]。在地黄转录组中鉴定出可能参与毛蕊花糖苷合成的编码催化酶的基因215个，但SA处理后，地黄块根中仅有少数催化酶基因上调表达（图4）。在酪氨酸途径中，编码乙醛脱氢酶（ALDH）的基因CL7331.Contig2，在SA处理1 h和3 h后的地黄块根中表达量显著升高，另一个ALDH编码基因Unigene9172在SA处理后的表达量也有一定程度上升。编码糖苷转移酶（UGT）的基因CL4677.Contig1随着SA处理时间的延长其表达量逐渐升高，在SA处理6 h后表达量最高。编码多酚氧化酶（PPO）的两个基因Unigene12832和Unigene12832在SA处理3 h和6 h后的表达量增加较为明显。而苯丙氨酸途径的催化酶基因表达变化不明显。这说明SA处理后毛蕊花糖苷含量的增加可能主要与酪氨酸途径部分催化酶基因的表达量增加有关。
![]()
图 4 SA处理后毛蕊花糖苷合成相关的催化酶基因的表达Figure 4. Expression characteristics of enzyme genes involved in acteoside biosynthesis after SA treatment2.5 差异表达转录因子筛选
对SA处理不同时间后地黄块根中差异表达的转录因子进行分析，结果表明（表2），25种转录因子的编码基因在处理前后呈差异表达。在处理1和3 h后，下调表达的转录因子较多，处理6 h后则上调表达的转录因子较多。WRKY、MYB、bHLH、AP2-EREBP、NAC和GRAS转录因子的差异表达基因较多，其中AP2-EREBP、WRKY和MYB的差异表达基因最多，且均表现为SA处理1 h后下调的基因较多， 3 h和6 h后上调表达的基因较多。进一步分析发现，共有20个转录因子编码基因在SA处理后的3个时间点均上调表达（表3），其中NAC和AP2-EREBP基因均为4个，WRKY基因有3个，MYB、GRAS、PLATZ基因各2个，bHLH、MADS和C2C2-CO-lik基因各1个。具有调控毛蕊花糖苷合成功能的RgWRKY37编码基因CL394.Contig2在SA处理1、3、6 h后的Log2(SA处理/CK)的值分别为0.79、0.82和0.57，Log2(SA处理/CK)的值虽然小于1，但其Q-value和P-value均达到了显著水平，说明在SA处理后CL394.Contig2上调表达。
表 2 SA处理后差异表达的转录因子数Table 2. Number of differentially expressed transcription factors (TFs) after SA treatment转录因子
Transcription factorCK-vs-1 h CK-vs-3 h CK-vs-6 h 共同差异表达的基因数
Number of common DEGs下调 Down 上调 Up 下调 Down 上调 Up 下调 Down 上调 Up zf-HD 2 1 1 1 0 0 0 WRKY 11 5 8 16 5 15 5 TUB 1 0 1 1 0 0 0 Trihelix 1 0 1 0 0 0 0 Tify 1 0 3 0 3 0 1 PLATZ 0 4 0 4 1 3 2 SRS 0 0 2 0 2 0 0 OFP 2 0 3 0 1 0 0 NAC 2 8 2 6 0 8 4 MYB 15 10 5 7 6 12 3 mTERF 1 0 2 0 1 1 0 MADS 2 2 3 4 2 2 1 LOB 1 0 4 2 2 0 0 HSF 5 2 4 6 3 6 2 GRAS 1 8 0 4 0 7 2 G2-like 2 0 2 1 0 3 0 CPP 2 0 2 0 4 0 2 C2H2 3 1 8 1 3 1 2 C2C2-GATA 2 0 1 0 5 0 0 C2C2-Dof 1 1 3 0 2 1 0 C2C2-CO-like 0 1 0 2 4 3 1 bZIP 1 2 0 0 1 0 0 bHLH 12 3 15 4 5 4 2 AP2-EREBP 15 12 12 18 9 21 6 ABI3VP1 5 2 7 1 1 1 0 总数 88 62 89 78 60 88 33 表 3 SA处理后上调表达的转录因子基因Table 3. Up-regulated transcription factor genes after SA treatment转录因子
Transcription factor基因
GeneLog2(SA/CK) 功能
Function1 h 3 h 6 h AP2-EREBP CL1637.Contig3 2.82 1.29 2.42 Ethylene-responsive transcription factor ERF071 AP2-EREBP CL4501.Contig2 1.66 1.75 2.16 Pathogenesis-related genes transcriptional activator PTI6 AP2-EREBP CL7827.Contig1 1.77 3.51 4.24 Ethylene-responsive transcription factor ERF106-like AP2-EREBP Unigene2558 1.62 1.06 1.74 Ethylene-responsive transcription factor 2-like MYB CL1983.Contig1 4.55 4.19 3.12 Transcription factor TFIIIB component B''-like MYB CL4303.Contig1 3.76 2.61 3.47 Single MYB histone protein NAC CL2945.Contig2 6.82 7.15 7.51 NAC domain-containing protein 82-like isoform X1 NAC CL4851.Contig1 3.15 1.99 2.44 NAC transcription factor 29 NAC CL4851.Contig2 3.15 2.27 2.82 NAC transcription factor 29 NAC Unigene5193 2.09 1.38 2.37 NAC domain-containing protein 72 C2C2-CO-like CL5505.Contig3 2.94 2.66 4.31 Zinc finger protein CONSTANS-LIKE 4-like PLATZ CL5569.Contig1 2.31 1.49 2.53 Interleukin-1 receptor-associated kinase 4 PLATZ CL5569.Contig3 2.47 1.71 2.38 Interleukin-1 receptor-associated kinase 4 GRAS CL645.Contig2 1.93 1.65 2.91 Scarecrow-like protein 14 GRAS Unigene10453 1.03 1.03 1.08 Scarecrow-like protein 15 WRKY CL6521.Contig1 2.84 1.66 2.35 Probable WRKY transcription factor 25 WRKY CL7324.Contig3 1.13 1.59 1.21 Probable WRKY transcription factor 35 WRKY CL791.Contig6 2.59 2.28 2.29 Probable WRKY transcription factor 40 bHLH Unigene12420 2.20 2.20 2.29 Phytochrome-interacting factor 3 MADS Unigene23315 2.17 3.01 2.27 MADS-box transcription factor 2.6 qRT-PCR验证基因表达差异
为了验证转录组测序对基因表达量分析的可靠性，随机选取CL7331.Contig2、Unigene1886、CL985.Contig1、CL5931.Contig2和CL379.Contig3等5个基因，利用qRT-PCR方法检测其在不同样本中的相对表达量。结果表明（图5），5个基因的定量结果与转录组获得的FPKM值变化趋势基本一致，其中CL7331.Contig2和CL5931.Contig2为SA处理后上调表达的基因，Unigene1886、CL985.Contig1和CL379.Contig3为SA处理后下调表达的基因。qRT-PCR与转录组测序的Pearson相关系数为0.602，相关性达显著水平，说明利用转录组测序分析SA处理后的基因表达量结果可靠。
![]()
图 5 差异表达基因的qRT-PCR验证A ~ E分别为CL7331.Contig2、Unigene1886、CL985.Contig1、CL5931.Contig2和CL379.Contig3的FPKM值与相对表达量；F：FPKM值与qRT-PCR的相关性分析。Figure 5. Validation of DEGs by qRT-PCRA–E: Represent expression and FPKM values of CL7331.Contig2, Unigene1886, CL985.Contig1, CL5931.Contig2, and CL379.Contig3; F: Correlation analysis between FPKM and qRT-PCR data.3. 讨论
植物次生代谢产物的积累既受自身遗传因素的控制，也受生长过程中生物与非生物环境条件的影响。一般而言，若药用植物的次生代谢产物在细胞中的含量相对较低，则会制约次生代谢产物的最终含量。近年来，人们常利用生物和非生物诱导子对植物进行处理，以提高植物特定次生代谢产物的生产[12]。如花生四烯酸（AA）、SA、茉莉酸甲酯（MeJA）和AgNO3均能够提高紫杉醇的含量[13]。50 µmol/L的乙烯利可以显著促进人参（Panax ginseng C. A. Meyer）根的生长和人参皂苷的积累[14]。诱导子提高苯乙醇苷含量的研究亦有报道，如外源添加Ag + 和腐胺均可以提高肉苁蓉（Cistanche deserticola Y. C. Ma）细胞培养物中松果菊苷和毛蕊花糖苷的含量[15]。本课题组前期研究发现，在地黄毛状根培养的培养基中添加25 μmol/L的SA可将毛蕊花糖苷的含量提高2.28倍[10]。本研究发现，叶面喷施100 μmol/L的SA可显著提高地黄叶片和块根中毛蕊花糖苷的含量，在块根中最高可提高1.4倍，说明在大田中叶面喷施诱导子可有效提高地黄块根中的次生代谢产物，有助于改善中药材的品质。
转录组测序分析不但可以高通量地获得基因表达的有关信息，还能够揭示基因表达与生命现象之间的内在联系，从而表征生命体的生理活动规律并确定其代谢特征[16]。目前，转录组测序不仅用于模式植物和大田作物生长发育及逆境胁迫响应关键基因的筛选[17-19]，在药用植物次生代谢产物合成调控的结构基因和转录因子基因的挖掘中也有广泛应用[20-22]。由于地黄为同源四倍体物种，其基因组测序虽有报道[23]，但作为参考基因组仍存在一些问题。因此，本研究利用课题组前期获得的地黄根、叶转录组为参考基因集进行分析，发现测序数据的特异匹配率偏低（50%左右），与地黄毛状根转录组测序的结果[10]类似，可能与其为同源四倍体物种有关。地黄叶片表面喷施SA后，上调表达的基因数少于下调表达，与SA处理的地黄毛状根结果[10]不同，这可能与本研究以大田地黄材料进行SA处理有关。本研究还发现，利用RNA-seq分析基因表达与qRT-PCR分析的结果相关系数仅为0.602，虽然达到显著相关，但未达到极显著相关水平，可能与地黄的基因组较大（约2.6 Gb），而采用RNA-seq测序获得的数据量较小有关。因此，对于基因组较大的物种，建议提高RNA-seq测序的深度，以获得更多的基因表达信息，提高基因表达量分析的准确性。
毛蕊花糖苷的生物合成途径目前已经比较清楚，其羟基酪醇基团来源于酪氨酸途径，咖啡酰基团来源于苯丙氨酸途径[24]。本课题组进一步推导、优化了毛蕊花糖苷的生物合成途径，认为毛蕊花糖苷是由羟基酪醇苷和咖啡酰辅酶A在莽草酸邻羟基肉桂酰转移酶（HCT）/毛蕊花糖苷合酶（AcS）和UGT的催化下合成[10]。周延清等[25] 基于地黄代谢组学分析获得了KEGG途径中的香豆酸-3-羟化酶（C3H），并克隆了其全长编码序列。李欣容等[26]根据SA处理下毛状根中催化酶基因的表达特性，鉴定并克隆了响应SA诱导的毛蕊花糖苷合酶基因RgAcS1。Yang等[27] 鉴定了4个酪氨酸脱羧酶（TyDC）基因，遗传转化发现过量表达RgTyDC2和RgTyDC4的地黄块根、纤维根、茎、嫩叶和成熟叶中的毛蕊花糖苷含量均显著高于野生型。Wang等[28]筛选了1个响应SA和H2O2诱导的WRKY转录因子基因RgWRKY37，功能研究发现RgWRKY37过量表达的毛状根转化体中毛蕊花糖苷和总苯乙醇苷的含量均显著高于对照。本研究发现，在SA处理的地黄块根中2个ALDH基因、1个UGT基因和2个PPO基因均上调表达，可能与块根中毛蕊花糖苷的含量增加有关。同时，在SA处理后的3个时间点，编码WKRY、NAC和AP2-EREBP等转录因子的20个基因均显著上调表达，其中RgWRKY37的表达量均明显增加。本研究结果为进一步探讨SA诱导毛蕊花糖苷合成的分子机理奠定了基础。
-
![]()
图 1 大菊杂交亲本遗传距离聚类分析
Figure 1. Cluster analysis of genetic distance of chrysanthemum hybrid parents
![]()
图 2 早花大菊杂交亲本
A：‘玉环公主’；B：‘东篱知秋’；C：‘绿宝石’；D：‘滦风飒爽’；E：‘滦凤洒爽’；F：‘红魁’；G：‘唐宇金秋’。
Figure 2. Hybrid parents of early flowering chrysanthemums
A: ‘Yuhuangongzhu’; B: ‘Donglizhiqiu’ ; C: ‘Lvbaoshi’ ; D: ‘Luanfengsashuang’ ; E: ‘Luanfengsashuang’ ; F: ‘Hongkui’ ; G: ‘Tangyujinqiu’.
表 1 大菊自然群体花期和主要观赏性状的变异分析
Table 1 Analysis of major variation in chrysanthemum varieties and ornamental characters
性状
Trait最小值
Minimum最大值
Maximum标准差
Standard deviation变异系数
Variation coefficient / %花径 / cm 8.00 29.30 3.00 19.00 瓣型 1.00 8.00 1.15 57.33 舌状花轮数 / 轮 0.00 25.00 4.08 51.50 舌状花数量 / 朵 0.00 2531.00 269.01 73.01 封顶叶片数 / 片 10.00 47.00 5.73 25.02 花色 2.00 9.00 1.83 33.50 开花期 / d 145.00 210.00 9.44 5.10 开花持续期 / d 13.00 49.00 6.90 22.30 
下载: 导出CSV
表 2 大菊自然群体主要观赏性状和花期性状间的相关系数
Table 2 Correlation coefficients between main ornamental and flowering characters of chrysanthemum varieties
性状
Trait观赏性状和花期性状间的相关系数 Correlation coefficients between ornamental and flowering characters 花径
Flower diameter瓣型
Petal type舌状花轮数
No. of ray flower rounds舌状花数量
No. of ray flowers封顶叶片数
No. of leaves花色
Color开花期
Flowering period开花持续期
Flowering duration花径 1.00 瓣型 0.11* 1.00 舌状花轮数 −0.08 −0.06 1.00 舌状花数量 −0.04 −0.09 0.12* 1.00 封顶叶片数 −0.01 −0.14* 0.03 −0.01 1.00 花色 0.02 0.03 −0.02 −0.07 −0.02 1.00 开花期 0.06 0.04 −0.07 0.06 −0.12* 0.05 1.00 开花持续期 −0.12 −0.06 0.40** 0.01 0.13 0.04 −0.38** 1.00 注：*、**分别表示在0.05和0.01水平上差异显著。下同。 Note: * and ** indicate significant difference at 0.05 and 0.01 levels, respectively. Same below.
下载: 导出CSV
表 3 初步筛选大菊杂交亲本的性状表现
Table 3 Preliminary screening of characters of chrysanthemum hybrid parents
杂交亲本
Hybrid parents花径
Flower diameter / cm瓣型
Petal type舌状花轮数
No. of ray flower rounds舌状花数量
No. of ray flowers封顶叶片数
No. of leaves花色
Color开花期
Flowering period / d开花持续期
Flowering duration / d定植日期
Planting date开花日期
Flowering date得分
Score等级
Grade‘东篱知秋’ 15.3 2.0 11.0 434.0 21.0 3.0 152.0 31.0 4月23日 9月22日 86.78 1 ‘唐宇金秋’ 17.4 5.0 12.0 457.0 27.0 7.0 166.0 33.0 4月23日 10月06日 86.37 1 ‘绿宝石’ 16.2 7.0 12.0 459.0 19.0 8.0 161.0 27.0 4月23日 10月01日 86.13 1 ‘滦风飒爽’ 17.6 2.0 25.0 516.0 21.0 6.0 171.0 39.0 4月23日 10月11日 84.80 1 ‘飞流直下’ 16.0 8.0 15.0 480.0 24.0 6.0 166.0 39.0 4月23日 10月06日 84.15 1 ‘红 魁’ 18.2 1.0 25.0 625.0 28.0 5.0 173.0 33.0 4月23日 10月13日 81.18 1 ‘滦凤洒爽’ 15.0 4.0 10.0 420.0 27.0 4.0 172.0 32.0 4月23日 10月12日 80.82 1 ‘唐宇秋荷’ 15.7 1.0 10.0 418.0 33.0 6.0 168.0 40.0 4月23日 10月08日 79.60 1 ‘粉旋球’ 17.8 5.0 8.0 228.0 21.0 6.0 171.0 34.0 4月23日 10月11日 77.14 1 ‘粉十八’ 16.5 7.0 2.0 50.0 17.0 6.0 172.0 36.0 4月23日 10月12日 76.21 1 ‘玉环公主’ 10.9 1.0 2.0 107.0 15.0 3.0 145.0 19.0 4月23日 9月15日 72.31 2 ‘晚月莲香’ 11.5 7.0 11.0 460.0 22.0 6.0 173.0 25.0 4月23日 10月13日 72.30 2 ‘天成秀色’ 16.4 5.0 12.0 458.0 27.0 3.0 165.0 32.0 4月23日 10月05日 72.11 2
下载: 导出CSV
表 4 早花大菊杂交组合
Table 4 Early flowering chrysanthemum hybrid combinations
杂交亲本
Hybrid parents杂交组合编号
Hybrid combination number母本 × 父本 ‘玉环公主’ Ⅰ ‘玉环公主’ × ‘东篱知秋’ ‘东篱知秋’ Ⅱ ‘东篱知秋’ × ‘玉环公主’ ‘绿宝石’ Ⅲ ‘绿宝石’ × ‘玉环公主’ ‘滦风飒爽’ Ⅳ ‘滦风飒爽’ × ‘玉环公主’ ‘滦凤洒爽’ Ⅴ ‘滦凤洒爽’ × ‘玉环公主’ ‘唐宇金秋’ Ⅵ ‘唐宇金秋’ × ‘东篱知秋’ ‘红魁’ Ⅶ ‘红魁’ × ‘东篱知秋’
下载: 导出CSV
表 5 大菊杂交群体花期和主要观赏性状的变异分析
Table 5 Variation analysis of flowering period and ornamental characters of chrysanthemum hybrid populations
性状
Trait杂交群体
Hybrid combination最小值
Minimum最大值
Maximum标准差
Standard deviation变异系数
Variation coefficient / %花径 / cm Ⅰ 5.20 16.00 2.38 23.00 Ⅱ 8.70 17.70 2.21 18.18 Ⅲ 7.00 17.80 1.95 17.25 Ⅳ 7.00 14.90 1.52 14.09 Ⅴ 6.90 17.20 2.06 18.39 Ⅵ 9.90 18.07 2.54 19.35 Ⅶ 6.00 17.90 2.02 18.50 瓣型 Ⅰ 1.00 5.00 0.53 25.88 Ⅱ 1.00 4.00 0.61 30.87 Ⅲ 1.00 7.00 1.21 67.22 Ⅳ 1.00 5.00 1.07 68.98 Ⅴ 1.00 8.00 1.25 71.05 Ⅵ 1.00 8.00 2.17 74.45 Ⅶ 1.00 5.00 0.78 50.09 舌状花轮数 Ⅰ 1.00 11.00 2.28 55.94 Ⅱ 3.00 24.00 4.58 49.41 Ⅲ 1.00 14.00 3.38 57.41 Ⅳ 1.00 15.00 2.89 48.69 Ⅴ 1.00 12.00 2.25 40.58 Ⅵ 2.00 4.00 3.67 67.29 Ⅶ 1.00 11.00 2.26 44.54 舌状花数量 Ⅰ 32.00 540.00 64.58 73.27 Ⅱ 36.00 664.00 156.56 66.79 Ⅲ 14.00 526.00 117.67 65.66 Ⅳ 29.00 344.00 72.94 54.48 Ⅴ 14.00 399.00 76.49 59.27 Ⅵ 42.00 524.00 128.07 80.55 Ⅶ 29.00 398.00 79.55 66.97 封顶叶片数 Ⅰ 8.00 42.00 5.92 26.93 Ⅱ 16.00 34.00 3.01 13.60 Ⅲ 11.00 33.00 6.63 28.97 Ⅳ 14.00 30.00 4.19 20.46 Ⅴ 12.00 31.00 4.71 23.23 Ⅵ 12.00 35.00 7.02 31.50 Ⅶ 11.00 33.00 4.61 23.69 花色 Ⅰ 2.00 7.00 1.57 39.83 Ⅱ 3.00 7.00 1.44 27.46 Ⅲ 2.00 7.00 1.57 36.52 Ⅳ 2.00 9.00 1.24 31.09 Ⅴ 2.00 8.00 2.04 41.15 Ⅵ 2.00 8.00 1.85 43.27 Ⅶ 2.00 7.00 1.46 36.38 开花期 / d Ⅰ 140.00 192.00 17.12 11.80 Ⅱ 131.00 172.00 16.17 10.70 Ⅲ 131.00 161.00 14.01 10.50 Ⅳ 132.00 168.00 15.00 9.85 Ⅴ 131.00 152.00 12.00 9.10 Ⅵ 140.00 192.00 16.15 11.30 Ⅶ 133.00 190.00 23.15 9.26 开花持续期 / d Ⅰ 14.00 30.00 5.61 25.70 Ⅱ 14.00 45.00 5.72 26.00 Ⅲ 14.00 45.00 6.73 26.10 Ⅳ 15.00 43.00 5.60 22.30 Ⅴ 16.00 45.00 6.32 21.23 Ⅵ 14.00 33.00 6.10 28.33 Ⅶ 15.00 42.00 6.71 35.35
下载: 导出CSV
表 6 大菊杂交群体主要观赏性状和花期性状间的相关系数
Table 6 Correlation coefficients between main ornamental and flowering period characters of hybrid populations of chrysanthemum
性状
Trait主要观赏性状和花期性状间的相关系数 Correlation coefficients between main ornamental and flowering characters 花径
Flower diameter瓣型
Petal
type舌状花轮数
No. of ray flower rounds舌状花数量
No. of ray flowers封顶叶片数
No. of leaves花色
Color开花期
Flowering period开花持续期
Flowering duration花径 1.00 瓣型 0.02 1.00 舌状花轮数 0.17** −0.12** 1.00 舌状花数量 0.16** −0.06 0.79** 1.00 封顶叶片数 0.08 −0.05 0.02 0.09* 1.00 花色 −0.23 0.24 0.28* 0.25 0.09 开花期 0.03 −0.15** −0.15** −0.03 −0.16** −0.35* 1.00 开花持续期 −0.03 −0.44** 0.20** 0.26** −0.07 −0.29* 0.18** 1.00
下载: 导出CSV
-
[1] 陈俊愉. 中国菊花过去和今后对世界的贡献[J]. 中国园林,2005,21(9):73−75. Chen JY. Contributions of Chinese chrysanthemum to the world in the past and future[J]. Chinese Landscape Architecture,2005,21 (9):73−75. Chen JY. Contributions of Chinese chrysanthemum to the world in the past and future[J]. Chinese Landscape Architecture, 2005, 21（9）: 73-75.
[2] 吴丽娜. 《红楼梦》菊花诗的拟题、分配与意象[J]. 明清小说研究,2021(4):121−135. doi: 10.3969/j.issn.1004-3330.2021.04.008 [3] 戴思兰,王文奎,黄家平. 菊属系统学及菊花起源的研究进展[J]. 北京林业大学学报,2002,24(5-6):234−238. Dai SL,Wang WK,Huang JP. Advances of researches on phylogeny of Dendranthema and origin of chrysanthemum[J]. Journal of Beijing Forestry University,2002,24 (5-6):234−238. Dai SL, Wang WK, Huang JP. Advances of researches on phylogeny of Dendranthema and origin of chrysanthemum[J]. Journal of Beijing Forestry University, 2002, 24（5-6）: 234-238.
[4] 戴思兰. 中国菊花的魅力[J]. 中国园林,2012,28(8):46−48. Dai SL. The fascinations of Chinese chrysanthemum[J]. Chinese Landscape Architecture,2012,28 (8):46−48. Dai SL. The fascinations of Chinese chrysanthemum[J]. Chinese Landscape Architecture, 2012, 28（8）: 46-48.
[5] 秦忠文. 中国传统菊花栽培起源与花文化发展[D]. 武汉: 华中农业大学, 2006: 20−33. [6] 李学松, 季玉山, 戴思兰. 百菊图话[M]. 北京: 化学工业出版社, 2013: 10−40. [7] 张树林, 戴思兰. 中国菊花全书[M]. 北京: 中国林业出版社, 2013: 120−133. [8] 赵小刚. 日中性小菊新品种选育及小菊开花期遗传分析[D]. 北京: 北京林业大学, 2019: 1−10. [9] 秦贺兰. 夏花型小菊新品[J]. 中国花卉园艺,2006(20):46−49. [10] 宋彬,刘景安. 宿根花卉花期调控的研究进展[J]. 北方园艺,2020(1):122−127. Song B,Liu JA. Research advance of flowering time regulator for perennial flower[J]. Northern Horticulture,2020 (1):122−127. Song B, Liu JA. Research advance of flowering time regulator for perennial flower[J]. Northern Horticulture, 2020（1）: 122-127.
[11] 张秋玲,刘海鹏,高康,孔德元,戴思兰. 盆栽菊花反季节开花调控技术研究[J]. 黑龙江农业科学,2021(9):62−67. Zhang QL,Liu HP,Gao K,Kong DY,Dai SL. Research on anti-seasonal flowering control technology of potted chrysanthemum[J]. Heilongjiang Agricultural Sciences,2021 (9):62−67. Zhang QL, Liu HP, Gao K, Kong DY, Dai SL. Research on anti-seasonal flowering control technology of potted chrysanthemum[J]. Heilongjiang Agricultural Sciences, 2021（9）: 62-67.
[12] 张莉俊,戴思兰. 菊花种质资源研究进展[J]. 植物学报,2009,44(5):526−535. Zhang LJ,Dai SL. Research advance on germplasm resources of Chrysanthemum × morifolium[J]. Bulletin of Botany,2009,44 (5):526−535. Zhang LJ, Dai SL. Research advance on germplasm resources of Chrysanthemum × morifolium[J]. Bulletin of Botany, 2009, 44（5）: 526-535.
[13] 雒新艳,宋雪彬,戴思兰. 中国传统大菊品种数量性状变异及其概率分级[J]. 北京林业大学学报,2016,38(1):101−111. Luo XY,Song XB,Dai SL. Variation and probability grading of quantitative characters of traditional chrysanthemum cultivars[J]. Journal of Beijing Forestry University,2016,38 (1):101−111. Luo XY, Song XB, Dai SL. Variation and probability grading of quantitative characters of traditional chrysanthemum cultivars[J]. Journal of Beijing Forestry University, 2016, 38（1）: 101-111.
[14] 雒新艳. 大菊品种资源遗传多样性研究[D]. 北京: 北京林业大学, 2009: 32−42. [15] 牛雅静,张蒙蒙,孟昕. 观赏植物托桂花型的研究进展[J]. 安徽农业科学,2021,49(24):29−32. Niu YJ,Zhang MM,Meng X. Research progress on anemone type of ornamental plants[J]. Journal of Anhui Agricultural Sciences,2021,49 (24):29−32. Niu YJ, Zhang MM, Meng X. Research progress on anemone type of ornamental plants[J]. Journal of Anhui Agricultural Sciences, 2021, 49（24）: 29-32.
[16] 宋雪彬,高康,戴思兰,季玉山,王朔,刘昊. 菊花新品种‘东篱知秋’[J]. 园艺学报,2020,47(S2):3024−3025. Song XB,Gao K,Dai SL,Ji YS,Wang S,Liu H. A new large-flowered chrysanthemum cultivar ‘Dongli Zhiqiu’[J]. Acta Horticulturae Sinica,2020,47 (S2):3024−3025. Song XB, Gao K, Dai SL, Ji YS, Wang S, Liu H. A new large-flowered chrysanthemum cultivar ‘Dongli Zhiqiu’[J]. Acta Horticulturae Sinica, 2020, 47（S2）: 3024-3025.
[17] 张玉鸽, 高康, 刘芷兰, 宋雪彬, 戴思兰. 不同开花期大菊品种观赏性状综合分析[C]//中国观赏园艺研究进展2018. 哈尔滨: 中国园艺学会, 2018: 50−58. Zhang YG, Gao K, Liu ZL, Song XB, Dai SL. Analysis of florescence ornamental characters on chrysanthemum cultivars[C]//Advances in Ornamental Horticulture of China 2018. Harbin: Chinese Society for Horticultural Science, 2018: 50-58.
[18] 雒新艳,张俊丽,张二海. 大菊主要数量性状分析及其应用探讨[J]. 江苏农业科学,2015,43(7):161−163. [19] 张蒙蒙, 王青, 戴思兰, 季玉山, 王朔. 盆栽小菊表型性状筛选与品种分类研究[C]//中国观赏园艺研究进展2014. 青岛: 中国园艺学会, 2014: 111−117. Zhang MM, Wang Q, Dai SL, Ji YS, Wang S. The screening of the phenotypic traits and taxonomic study for potted chrysanthemum[C]//Advances in Ornamental Horticulture of China 2014. Qingdao: Chinese Society for Horticultural Science, 2014: 111−117.
[20] 许兰杰,梁慧珍,余永亮,谭政委,杨青,李磊. 菊花种质综合评价体系构建及优异种质筛选[J]. 北方园艺,2022(12):55−63. Xu LJ,Liang HZ,Yu YL,Tan ZW,Yang Q,Li L. Establishment of comprehensive evaluation of chrysanthemum germplasm and selection of elite resource[J]. Northern Horticulture,2022 (12):55−63. Xu LJ, Liang HZ, Yu YL, Tan ZW, Yang Q, Li L. Establishment of comprehensive evaluation of chrysanthemum germplasm and selection of elite resource[J]. Northern Horticulture, 2022（12）: 55-63.
[21] 高康, 宋雪彬, 戴思兰, 季玉山, 王朔. 基于AHP法的大菊杂种F1代新品种筛选[C]//中国观赏园艺研究进展2016. 长沙: 中国园艺学会, 2016: 127−133. Gao K, Song XB, Dai SL, Ji YS, Wang S. AHP-based screening of traditional chrysanthemum hybrid F1 for new variety[C]//Research Progress of Ornamental Horticulture in China 2016. Changsha: Chinese Society for Horticultural Science, 2016: 127-133.
[22] 宋雪彬,高康,黄河,刘芷兰,戴思兰,嵇彧. 中国传统大菊叶片形态的数量化定义与分类[J]. 植物学报,2021,56(1):10−24. Song XB,Gao K,Huang H,Liu ZL,Dai SL,Ji Y. Quantitative definition and classification of leaves in large-flowered Chinese chrysanthemum based on the morphological traits[J]. Bulletin of Botany,2021,56 (1):10−24. Song XB, Gao K, Huang H, Liu ZL, Dai SL, Ji Y. Quantitative definition and classification of leaves in large-flowered Chinese chrysanthemum based on the morphological traits[J]. Bulletin of Botany, 2021, 56（1）: 10-24.
[23] 盖钧益, 章元明, 王建康. 植物数量性状遗传体系[M]. 北京: 科学出版社, 2003: 1−20. [24] 冯俊杰,赵文达,张新全,刘英杰,袁帅,等. 引种日本多花黑麦草标准品种DUS性状变异分析及应用[J]. 中国农业科学,2022,55(12):2447−2460. Feng JJ,Zhao WD,Zhang XQ,Liu YJ,Yuan S,et al. DUS traits variation analysis and application of standard varieties of Lolium multiflorum introduced from Japan[J]. Scientia Agricultura Sinica,2022,55 (12):2447−2460. Feng JJ, Zhao WD, Zhang XQ, Liu YJ, Yuan S, et al. DUS traits variation analysis and application of standard varieties of Lolium multiflorum introduced from Japan[J]. Scientia Agricultura Sinica, 2022, 55（12）: 2447-2460.
[25] 中华人民共和国农业农村部. GB/T 19557.19-2018 植物品种特异性、一致性和稳定性测试指南 菊花[S]. 北京: 中国标准出版社, 2018. Ministry of Agriculture and Rural Affairs of the People’s Republic of China. GB/T 19557.19-2018 Guidelines for the conduct of tests for distinctness, uniformity and stability-Chrysanthemum × morifolium Ramat. （Chrysanthemum × grandiflorum Ramat.）, Chrysanthemum pacificum Nakai （Ajania pacifica Bremer and Humphries） and hybrids between them[S]. Beijing: Standards Press of China, 2018.
[26] 杜世坤,赵宝勰,赵振宁,高玉芳,李雨阳,等. 间作用大豆亲本的相关性、主成分及聚类分析[J]. 分子植物育种,2022,20(1):266−275. Du SK,Zhao BX,Zhao ZN,Gao YF,Li YY,et al. Correlation analysis,principal component analysis and cluster analysis in parent selection of intercropping soybean[J]. Molecular Plant Breeding,2022,20 (1):266−275. Du SK, Zhao BX, Zhao ZN, Gao YF, Li YY, et al. Correlation analysis, principal component analysis and cluster analysis in parent selection of intercropping soybean[J]. Molecular Plant Breeding, 2022, 20（1）: 266-275.
[27] 杨群力,李思锋. 大丽花三种园艺性状之间相关性的研究[J]. 中国农学通报,2009,25(23):295−302. Yang QL,Li SF. Study on relativity among three gardening characters of dahlia pinnata[J]. Chinese Agricultural Science Bulletin,2009,25 (23):295−302. Yang QL, Li SF. Study on relativity among three gardening characters of dahlia pinnata[J]. Chinese Agricultural Science Bulletin, 2009, 25（23）: 295-302.
[28] 徐肇友,肖德卿,沈斌,肖纪军,陈杏林,周志春. 秀丽四照花叶色参数和叶形性状的变异及相关性分析[J]. 植物资源与环境学报,2021,30(1):61−68. Xu ZY,Xiao DQ,Shen B,Xiao JJ,Chen XL,Zhou ZC. Variation and correlation analyses on leaf color parameters and leaf shape traits of Cornus hongkongensis subsp. elegans[J]. Journal of Plant Resources and Environment,2021,30 (1):61−68. Xu ZY, Xiao DQ, Shen B, Xiao JJ, Chen XL, Zhou ZC. Variation and correlation analyses on leaf color parameters and leaf shape traits of Cornus hongkongensis subsp. elegans[J]. Journal of Plant Resources and Environment, 2021, 30（1）: 61-68.
[29] 赵芳,魏玮,张晓磊,宋国亮,王晓明,等. 224个谷子品种农艺性状聚类和相关性分析[J]. 种子,2022,41(1):74−83. Zhao F,Wei W,Zhang XL,Song GL,Wang XM,et al. Clustering and correlation analysis of agronomic traits of 224 millet varieties[J]. Seed,2022,41 (1):74−83. Zhao F, Wei W, Zhang XL, Song GL, Wang XM, et al. Clustering and correlation analysis of agronomic traits of 224 millet varieties[J]. Seed, 2022, 41（1）: 74-83.
[30] 吴瑞,刘文辉,张永超,秦燕,魏小星,刘敏洁. 青藏高原老芒麦落粒性及农艺性状相关性研究[J]. 草业学报,2021,30(4):130−139. Wu R,Liu WH,Zhang YC,Qin Y,Wei XX,Liu MJ. A study of the correlation between seed shattering and agronomic traits of Elymus sibiricus on the Qinghai-Tibetan Plateau[J]. Acta Prataculturae Sinica,2021,30 (4):130−139. Wu R, Liu WH, Zhang YC, Qin Y, Wei XX, Liu MJ. A study of the correlation between seed shattering and agronomic traits of Elymus sibiricus on the Qinghai-Tibetan Plateau[J]. Acta Prataculturae Sinica, 2021, 30（4）: 130-139.
[31] 王旭,高致明,张红瑞,张景景,马彦秋,元玉碧. 6个药用菊花栽培类型生长势及抗性综合评价[J]. 河南农业大学学报,2012,46(2):131−135. Wang X,Gao ZM,Zhang HR,Zhang JJ,Ma YQ,Yuan YB. Comprehensive evaluation of the growth potential and resistance in six cultivars of medicinal Chrysanthemum morifolium[J]. Journal of Henan Agricultural University,2012,46 (2):131−135. doi: 10.3969/j.issn.1000-2340.2012.02.004 Wang X, Gao ZM, Zhang HR, Zhang JJ, Ma YQ, Yuan YB. Comprehensive evaluation of the growth potential and resistance in six cultivars of medicinal Chrysanthemum morifolium[J]. Journal of Henan Agricultural University, 2012, 46（2）: 131-135, 151. doi: 10.3969/j.issn.1000-2340.2012.02.004
[32] 吴盼婷,王江民,沈佳逾,杨友翠,管志勇,等. 不同菊花品种根系、地上部和叶片相关指标分析及抗逆性评价[J]. 植物资源与环境学报,2017,26(2):46−54. Wu PT,Wang JM,Shen JY,Yang YC,Guan ZY,et al. Analyses on related indexes of root,above-ground part and leaf of different cultivars of Chrysanthemum morifolium and stress resistance evaluation[J]. Journal of Plant Resources and Environment,2017,26 (2):46−54. Wu PT, Wang JM, Shen JY, Yang YC, Guan ZY, et al. Analyses on related indexes of root, above-ground part and leaf of different cultivars of Chrysanthemum morifolium and stress resistance evaluation[J]. Journal of Plant Resources and Environment, 2017, 26（2）: 46-54.
[33] 马杰,徐婷婷,苏江硕,杨信程,房伟民,等. 菊花F1代舌状花耐寒性遗传变异与QTL定位[J]. 园艺学报,2018,45(4):717−724. Ma J,Xu TT,Su JS,Yang XC,Fang WM,et al. Genetic variation and QTL mapping for cold tolerance of ray florets in an F1 population of Chrysanthemum morifolium[J]. Acta Horticulturae Sinica,2018,45 (4):717−724. Ma J, Xu TT, Su JS, Yang XC, Fang WM, et al. Genetic variation and QTL mapping for cold tolerance of ray florets in an F1 population of Chrysanthemum morifolium[J]. Acta Horticulturae Sinica, 2018, 45（4）: 717-724.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 183
- HTML全文浏览量: 32
- PDF下载量: 43
- 被引次数: 1
