Genetic analysis of plant type traits in hybrid populations of Lagerstroemia indica and L. limii
-
摘要:目的
福建紫薇是典型的乔木型树种，为研究福建紫薇改良紫薇乔木性状的潜力。
方法本研究以福建紫薇与灌木型紫薇品种‘千层绯雪’为亲本构建了正、反交群体，采用相关性分析、主基因+多基因混合遗传模型对一年生杂交后代的株高、地径、一级分枝点高度、一级分枝数等七个性状进行了分析。
结果（1）杂交后代的株高在1.50～24.5 cm之间、地径0.30～4.44 mm、一级分枝点高度0.50～6.20 cm、一级分枝数0～4个，多数性状基本呈正态分布，后代中存在潜在的乔木型种质。（2）七个性状的变异系数范围在31.72%～240.42%之间，变异程度最大的是千层绯雪×福建紫薇群体的一级分枝数，最小的是该群体的一级分枝点高度。（3）相关性分析表明，株高与地径、一级分枝数与开张角度、开张角度与侧枝GSA呈极显著正相关。（4）株高、地径、开张角的最适模型为2MG-AD，侧枝GSA、一级分枝点高度的最适模型为2MG-EA。
结论研究结果说明，福建紫薇是培育生长量大、基部不分枝、主干明显的紫薇品种的优良材料，为乔木型紫薇的培育提供了依据。
Abstract:ObjectiveTo investigate the potential of Lagerstroemia limii to improve L. indica plant type traits.
MethodHybrid populations were constructed using L. limii (tree) and L. indica 'Qianceng Feixue' (shrub) as parents, and seven plant type traits (plant height, plant diameter, length of the first branch, et al.) of the one-year old hybrid plants were analyzed using correlation analysis and major gene plus polygene mixed inheritance model.
Result(1) The plant height of hybrid offsprings ranged from 1.50 cm to 24.5 cm, plant diameter from 0.30 mm to 4.44 mm, length of the first branch from 0.50 cm to 6.20 cm, and the number of branches from 0 to 4. The traits were basically normal distributed. Potential arborescent germplasm is present in the offspring. (2) The coefficients of variation of the seven traits ranged from 31.72% to 240.42%. The largest variation was the number of branches in 'Qianceng Feixue'×L. limii population, and the smallest one was the height of the first branch in 'Qianceng Feixue '×L. limii population. (3) Correlation analysis showed that plant height was positively correlated with plant diameter, number of branches was positively correlated with branch angle, branch angle was positively correlated with branch GSA. (4) The optimal model for plant height, plant diameter and branch angle were 2MG-AD, and the optimal model of branch GSA and the length of the first branch was 2MG-EA.
ConclusionIt can be concluded that L. limii is an excellent material for breeding cultivars of L. indica with tall trunk, unbranched and obvious main stem.
-
Keywords:
- Lagerstroemia /
- distant hybridization /
- plant-type traits /
- genetic analysis
-
0. 引言
【研究意义】兔出血症2型（rabbit hemorrhagic disease Type 2, RHD-2）是由兔出血症病毒2型（rabbit hemorrhagic disease virus serum type 2, RHdV-2）引起的一种新型高度接触传染性、急性致死性传染病。新型毒株的RHdV-2的感染性更强，感染范围更广，不同日龄和品种的家兔、野兔均可被感染[1−3]。RHD-2具有较高发病率和致死率[4]，且国内尚未研制出可靠的商品化疫苗供生产使用，一旦发病将严重影响养兔业的健康发展。当前，我国养兔规模化程度仍然十分低，散养户是养兔业的主力。散养户普遍存在技术力量弱、不会操作检测仪器等问题。RHdV-2感染已对我国养兔也造成了巨大经济损失，如何在养殖场特别是散养户进行快速准确的诊断对于科学防控该病具有十分重要的意义。通过建立快速且准确的检测方法，可以有效识别并控制传染源，从而防止疫情的进一步扩散。【前人研究进展】目前，重组酶介导的核酸等温扩增（recombinase aided amplification, RAA）与CRISPR/Cas13a技术联用在病原检测领域已逐渐展现出其巨大的应用潜力。该法先对目标序列进行RAA扩增以获得大量的检测模板，再利用CRISPR系统对其进行特异性识别以激活Cas13a蛋白来剪切反应体系中的荧光探针，从而实现对低载量样品的快速准确检测[5−6]。作为病原检测研究领域的一个热门方向，RAA-CRISPR/Cas13a检测方法通过结合RAA法和CRISPR/Cas13a技术，在灵敏度和特异性方面实现了显著提升。一些学者已经利用这种方法开发出了快速检测动物疫病病原微生物的新方法，如禽腺病毒 4 型[7]、猪流行性腹泻病毒[8]、口蹄疫病毒[9]和禽流感病毒[10]。虽然国内已经建立多种针对RHdV-2的诊断技术，但是均需要一定的专业知识和仪器设备，难以在养殖场一线应用。有研究根据RHdV-2的保守序列VP60，建立了实时荧光定量聚合酶链反应（real-time quantitative polymerase chain reaction, qPCR）检测方法[11]，酶联免疫吸附技术（enzyme Linked Immunosorbent assay, ELISA）也可用于检测RHdV-2抗原或抗体[12]，这些方法虽然其灵敏度和特异性也较好，但它们均需较为昂贵的配备仪器和过长的反应时间。【本研究切入点】本实验室前期建立了基于RAA的侧流层析试纸（lateral flow device, LFD）方法[13]，用于RHdV-2的特异性检测，然而RAA的扩增效率易受引物设计、模板序列及其他因素的制约，限制了其检测灵敏度和稳定性。近年来，CRISPR-Cas核酸酶的精准切割特性为核酸检测提供了新的技术突破，多项研究[14−15]通过联合RAA与CRISPR-Cas系统，显著提升了检测的灵敏度和特异性，并结合LFD实现了检测结果的可视化。基于此，本研究创新性地整合了RAA与CRISPR/Cas13a技术，首次建立了一种兔出血症病毒2型可视化诊断方法RAA-CRISPR/Cas13a(cpf1)-LFD。【拟解决的关键问题】本研究在本实验室前期研究的重组酶等温扩增结合侧流层析试纸条方法[13]的基础上，结合CRISPR/Cas13a技术，建立了新的RHdV-2检测方法，这种方法具备高度的灵敏度和特异性，操作简单，可实现在实验室和基层养殖场现场的快速检测，同时利用LFD实现结果可视化，为RHdV-2的诊断提供高效准确的检测方法。
1. 材料与方法
1.1 试验材料
1.1.1 重组质粒与临床样本
RHdV-2（KU991797.1）VP60重组质粒、RHdV-1（DQ205345.1）VP60重组质粒均由擎科股份有限公司进行合成。31份临床样本由本实验室保存，来自闽侯、武平、大田、连江等多个兔场，病兔主要表现急性死亡，实质器官有出血，采集脾脏淋巴结等实质器官。
1.1.2 主要试验试剂
总RNA提取试剂盒（ER501-01）（北京全式金生物有限公司）；RAA核酸扩增试剂盒，RAA核酸扩增试剂盒（试纸条法）（T00R01)（江苏奇天基因生物科技有限公司）；一次性核酸检测试纸条（JY0301，北京宝盈同汇生物技术有限公司）；1xPBS（兰杰柯科技有限公司）；LwaCas13a蛋白（C2C2）（广州美格生物有限公司）；HiScribeT7快速高效 RNA 合成试剂盒（New England Biolabs）；重组RNase抑制剂（荷瑞生物）；RNA纯化试剂盒（DP412，天根生物有限公司）；LwaCas13a 10×reaction buffer：200 mmol·L−1 HEPES-NaOH（pH6.8），600 mmol·L−1 NaCl，60 mmol·L−1 MgCl2，10 mmol·L−1 DTT。
1.2 试验方法
1.2.1 RAA引物、crRNA设计与合成
RAA特异性引物和RAA-LFD探针按照文献[13]报道的序列合成。在已扩增的特异且保守序列范围内，首先识别符合条件的PAM序列，继而设计长度为20～23 nt的间隔区序列构成LwaCas13a-crRNA，间隔区序列位于PAM序列之后，同时确保LwaCas13a-crRNA目标序列内3'端的PFS序列由A、C、U组成，且目标序列不与RAA引物重叠，最终，将间隔区序列进行Blast分析比对确定其保守特异性。crRNA的单链DNA（即ssDNA）由擎科股份有限公司进行合成。
Lwacas13a-crRNA合成步骤：将crRNA的单链DNA作为模板与相应引物通过PCR扩增获得大量crRNA的双链DNA（即dsDNA），使用琼脂凝胶回收试剂盒回收目的条带进行纯化，纯化产物通过T7 HiScribe T7高效RNA合成试剂盒37 ℃过夜转录，转录产物使用RNA纯化试剂盒纯化后通过nanodrop 2000测定RNA浓度置于−80 ℃保存。CRISPR-Cas13a相关序列见表1。
表 1 CRISPR-Cas13a crRNA及相关引物Table 1. CRISPR-Cas13a crRNA and related primers基因
Genes序列5'-3'
Sequences (5’-3’)Cas13a crRNA1 TAATACGACTCACTATAGGGGATTTAGACTACCCCAAAAACGAAGGGGACTAAAACACTCATAAGCCTGCATGGTCGTGACGTA Cas13a crRNA2 TAATACGACTCACTATAGGGGATTTAGACTACCCCAAAAACGAAGGGGACTAAAACGGTGGTGGTGGGTTGGGGGTTGCTCGGT 斜体为T7启动子序列，下划线碱基为重复序列区。
The italicized sequence represents the T7 promoter, and the underlined bases indicate the repetitive sequence region.PCR 25 μL反应体系：2×San Taq PCR Master Mix （含蓝色染料）12.5 μL，crRNA-dsDNA 1 μL，上下游引物各1 μL（10 μmol·L−1），ddH2O 9.5 μL。PCR反应过程：95 ℃预变性5 min；95 ℃变性30 s，56 ℃退火30 s，72 ℃延伸15 s，30个循环；72 ℃延伸10 min。
转录体系：T7 mix 2 μL，NTP mix 10 μL，模板1 μg，ddH2O补齐至20 μL。
1.2.2 RAA-LFD反应体系
根据RAA核酸扩增试剂盒说明书，选择50 μL体系进行扩增。反应体系为：缓冲液25 μL，纯水15.7 μL，正向与反向引物各取2.1 μL，探针取0.6 μL，乙酸镁（280 mmol·L−1）2.5 μL，模板3 μL。
1.2.3 Cas13a试纸条法反应条件的确立
（1） RAA-CRISPR/Cas13a (Cpf1)-LFD体系优化。1.5 μL crRNA（300 ng·μL−1）与50 nmol·L−1 Cas13a蛋白进行反应，选择反应效果最佳体系为CRISPR-Cas13a反应体系。该体系的总体积为50 µL，包括以下成分：10×反应缓冲液5 µL，crRNA（浓度为300 ng·µL−1）1.5 µL，探针（浓度为10 μmol·L−1）2 µL，模板2.5 µL，LwaCAS13a蛋白（浓度为1 μmol·L−1）3.5 µL，NTP混合物2.5 µL，RNA酶（浓度为300U·µL−1）5 µL，T7混合物0.5 µL，以及双蒸水27.5 µL。
（2）RAA-CRISPR/Cas13a-LFD反应时间优化。基于上述Cas13a体系对RAA-CRISPR/Cas13a-LFD反应时间进行优化，每5 min一个间隔，设置反应时间为10～25 min，Cas13a试纸条法通过试纸条观察结果。
1.2.4 RAA-CRISPR/Cas13a (Cpf1)-LFD试纸条法检测敏感性的评价
RHdV-2 VP60重组质粒经检测其质粒浓度为6.7×103 copies·μL−1，用1×PBS梯度稀释得到6.7×103、6.7×102、6.7×101、6.7×10−1、6.7×10−2 copies·μL−1 5个浓度梯度，以不同浓度的重组质粒为模板进行扩增，扩增产物与Cas13a蛋白以及crRNA等试剂再进行一定时间恒温反应，反应产物经LFD显示，分析敏感性试验结果。
1.2.5 RAA-CRISPR/Cas13a (Cpf1)-LFD检测RHdV-2特异性评价
RHdV-2 VP60重组质粒（6.7×105 ng·μL−1）、RHdV-1 VP60重组质粒（6.7×105 ng·μL−1）为模板进行RAA扩增，产物与Cas13a蛋白以及crRNA再进行一定时间恒温反应，反应产物经LFD显示。
1.2.6 RNA的提取
根据总RNA提取试剂盒说明进行操作提取，具体步骤：将单份临床兔病料75 mg与1 mL Trans Zol混合并研磨裂解样品，后加入200 μL RNA Extraction Agent，振荡混匀5 min后放于4 ℃高速离心机10 000 r·min−1离心15 min；吸取无色水相部分与同体积无水乙醇混合后加入离心柱后离心；重复加入CB9（500 μL）并离心；重复加入WB9（500 μL）并离心；离心去除剩余乙醇，使用75 μL纯水洗脱后，将提取的RNA保存于-80 ℃保存。
1.2.7 临床样品的检测
提取31份临床样本病料的总RNA，应用建立的RAA-CRISPR/Cas13a(Cpf1)-LFD 法分别对31份临床样品进行RAA扩增反应，扩增产物经LFD显示结果。
2. 结果与分析
2.1 Cas13-crRNA筛选结果
分别以crRNA1和crRNA2为探针对同一阳性样本进行检测，结果如图1A和图1B所示，crRNA1显示出更明显的检测线（T线），并且Image J定量分析（图1C）可见crRNA1与crRNA2之间的信号强度差异显著（P＜0.05），说明crRNA1具有更高的信号强度。综上后续选择crRNA1进行试验。
![]()
图 1 Cas13a-crRNA筛选结果A：试纸条；B：Image J定量图；C：样品条带定量柱状图。1：crRNA 1；2: crRNA 2；N：阴性对照。C：质控线；T：检测线；* 表示与对照相比差异显著 P ＜ 0.05。Figure 1. Cas13-crRNA screening resultsA: test strip image; B: image J quantification graph; C: sample band quantification bar chart. 1: crRNA1; 2: crRNA2; N: negative control. C: quality control Line; T:test line; * indicates a statistically significant difference with control (P ＜ 0.05).2.2 RAA-CRISPR/Cas13a (Cpf1)-LFD反应时间优化
从图2中可以看出，在RAA基础反应30 min后，仅再需10 min，RAA-CRISPR/Cas13a(Cpf1)-LFD即可检测出RHdV-2，且显色明显，从量化柱状图可知，反应时间为10 min时，数值最高（P＜0.05），综合分析可知10 min时即可达到较佳检测效果，故后续反应时间均釆用10 min。
![]()
图 2 不同反应时间RAA-CRISPR/Cas13a(Cpf1)-LFD扩增产物A：试纸条图；B：Image J定量图；C：样品条带定量柱状图。10、15、20、25分别为 10、15、20、25 min 扩增产物； N阴性对。C：质控线；T：检测线；*表示与对照相比差异显著（P ＜ 0.05）。Figure 2. Analysis of RAA-CRISPR/Cas13a(Cpf1)-LFD amplification products at different reaction timeA: test strip image; B: image J quantification graph; C: sample band quantification bar chart. 10, 15, 20, 25: amplification products for 10, 15, 20, 25 min with the standard plasmid as template; N: negative control. C: quality control line; T: test line; * indicates significant difference with control (P ＜ 0.05).2.3 RAA-CRISPR/Cas13a (Cpf1)-LFD敏感性
Cas13a试纸条法对RHdV-2 RNA最低可检值可达6.7×101 copies·μL−1，结合量化结果（图3）综合分析在质粒浓度为6.7×101 ng·μL−1时阳性结果最亮且数值最高（P＜0.05），阴性对照结果为阴性。
![]()
图 3 Cas13a敏感性试验结果A：试纸条图；B：Image J定量图;B：样品条带定量柱状图。1：6.7×103 copies·μL−1；2：6.7×102 copies·μL−1；3：6.7×101 copies·μL−1；4：6.7×10−1 copies·μL−1；5：6.7×10−2 copies·μL−1；N：阴性；* 表示与对照相比差异显著（P ＜ 0.05）。Figure 3. Sensitivity test results about Cas13aA: test strip image; B: image J quantification graph; C: sample band quantification bar chart. 1: 6.7×103 copies·μL−1; 2: 6.7×102 copies·μL−1; 3: 6.7×101 copies·μL−1; 4: 6.7×10−1 copies·μL−1; 5: 6.7×10−2 copies·μL−1; N: negative control; * indicates significant difference with control (P ＜ 0.05).2.4 RAA-CRISPR/Cas13a (Cpf1)-LFD特异性试验结果
如图4可知，RHdV-1结果为阴性，RHdV-2在试纸条检测中呈现强阳性，因此所建立的方法特异性优异。
![]()
图 4 Cas13特异性检测结果A：试纸条图；B：Image J定量图； B：样品条带定量柱状图。1：RHdV-2；2：RHdV-1；N：阴性对照；*表示与对照相比差异显著（P ＜ 0.05）。Figure 4. Cas13a specific test resultsA: test strip image; B: image J quantification graph; C: sample band quantification bar chart. 1: RHdV-2; 2: RHdV-1; N: negative control; * indicates significant difference with control (P ＜ 0.05).2.5 RAA-CRISPR/Cas13a (Cpf1)-LFD法检测结果
对31份不同兔场收集的病死兔脾脏组织进行RAA-CRISPR/Cas13a (Cpf1)-LFD试纸条检测，如图5可知，所有样品在C位置均有条带，表明所有样品均合格，样品4、10、13、15、18、19、22、24、27、29在T位置出现了较为明显的条带，并且Image J定量分析（图5C）可见样品4、10、13、15、18、19、22、24、27、29的条带强度与其他样品相比差异显著（P＜0.05），表明样品4、10、13、15、18、19、22、24、27、29均为阳性样品，RAA-CRISPR/Cas13a (Cpf1)-LFD法的阳性检出率为30.3%。
![]()
图 5 RAA-CRISPR/Cas13a (Cpf1)-LFD临床样品检测结果A：试纸条图；B：Image J定量图；C：样品条带定量柱状图。1～31分别为样品编号； N为阴性对照；*表示与对照相比差异显著（P ＜ 0.05）。Figure 5. Clinical samples detected by RAA-CRISPR/Cas13a (Cpf1)-LFDA: test strip image; B: image J quantification graph; C: sample band quantification bar chart. 1—31 are sample numbers; N: negative control; * indicates significant difference with control (P ＜ 0.05).3. 讨论与结论
自2010年在法国首次发现RHdV-2以来[4]，该病毒已在全球范围内流行。2020年，中国四川省首次报告了RHD2病例[16]。RHdV-2对养兔业构成严重威胁，但目前缺乏有效的商业疫苗。因此，开发可靠的病原学检测技术对于控制RHdV-2传播、减轻其对养兔业和生态系统的影响至关重要。
王波[17]在其研究中建立的检测RHdV-2的SYBR Green Ⅰ实时荧光定量RT-PCR方法检测限度达到68个拷贝数，而本研究中我们构建的RHdV-2核酸试纸条检测方法对兔出血症病毒RNA的最低检测限为6.7×101 copies·μL−1，说明本研究建立的方法灵敏度与实时荧光定量RT-PCR方法相当。同时，应用该方法进行检测，结果显示RHdV-1结果为阴性，RHdV-2则呈现强阳性，说明该方法特异性良好，与RHdV-1血清型核酸无交叉反应。本研究使用新建立的RAA-CRISPR/Cas13a (Cpf1)-LFD方法对本实验室前期研究中保存的31份疑似阳性临床样品[13]进行检测后，阳性率从19.2%提高到了30.3%，说明新方法在检测RHdV-2方面具有更高的有效性和准确性。近年来，关于RHdV-2检测方法的研究主要集中在RT-PCR技术上。尽管部分研究[18−19]开发的检测方法能够达到1×101 copies·μL−1的高灵敏度，但这些方法通常需要较长的反应时间，并且依赖于专业设备来判读结果，这限制了它们在临床快速检测RHDV-2中的应用。相比之下，本研究建立的检测方法在40 min内即可完成，并且通过试纸条直接显示结果，使得结果判读变得快速且直观。但本研究建立的方法也存在不足之处，反应过程中需要开盖转移RAA反应产物作为体系中的模板，易产生气溶胶污染。目前已有RAA-CAS一管反应法，将RAA反应试剂与Cas相关试剂一次性加入一个管内，反应时间相应缩短，减少气溶胶污染。然而，一管法灵敏性较低，这可能与RAA中反应蛋白酶接触Cas相关试剂影响了蛋白酶的活性，后期是否可以通过改变体系内试剂浓度进行一管反应还需进一步研究。
本研究前期通过RHdV两种血清型基因组对比，发现RHdV-2的VP60基因最为保守，将其设定为检测靶序列，合成特异性引物，建立了RAA-CRISPR/Cas13a(cpf1)-LFD检测方法。该方法具有敏感性强、特异性高、操作简易等优点，可实现对RHdV-2快速检测，可为RHdV-2的预防与早期诊断、控制提供便捷技术，规避兔养殖业面临巨大疫病风险，助力养兔业蓬勃发展。
-
![]()
图 2 QC×FJ群体杂交后代株型性状的频率分布直方图
Figure 2. Frequency distribution histogram of plant type traits in QC×FJ population
![]()
图 3 FJ×QC群体杂交后代株型性状的频率分布直方图
Figure 3. Frequency distribution histogram of plant type traits in FJ×QC population
表 1 福建紫薇与‘千层绯雪’的杂交结实率
Table 1 seed setting rate of L.limii and L. indica ‘Qianceng Feixue’
母本×父本
Maternal×
Paternal杂交数量
Flower Number
of crosses结实数
Fruit number结实率
Fruiting
rate/%出苗数
Seedling
number‘千层绯雪’×
福建紫薇371 280 75.5 963 福建紫薇×
‘千层绯雪’17 17 100.0 66 合计 388 297 76.5 1029 
下载: 导出CSV
表 2 QC×FJ群体的描述性统计
Table 2 Plant type difference in hybrid population
性状
Traits最大值
Maximum最小值
Minimum均值
Average标准差
SD偏度
Skewness峰度
Kurtosis变异系数
CV/%株高
Plant height/cm24.5 1.50 6.87 357 1.32 2.16 51.97 地径
Diameter/mm4.44 0.30 1.40 0.61 1.09 1.55 43.57 一级分枝点高度
Height at first branching point/cm6.20 0.80 3.31 1.05 0.11 0.19 31.72 一级分枝数
Number of first-order branches/each4 0 0.24 0.577 2.591 6.839 240.42 开张角度
Opening angle/°87.56 8.20 27.90 11.35 1.89 6.43 40.68 侧枝GSA
Lateral branch GSA/°72.25 15.35 35.26 12.42 0.71 0.05 35.22 节间长度
Internode length/cm2.0 0.4 0.73 0.26 1.58 3.23 35.62
下载: 导出CSV
表 3 FJ×QC群体的描述性统计
Table 3 Descriptive statistics of the FJ×QC population
性状
Traits最大值
Maximum最小值
Minimum均值
Average标准差
SD偏度
Skewness峰度
Kurtosis变异系数
CV/%株高
Plant height/cm22.0 2.0 7.54 4.08 1.06 2.19 54.11 地径
Diameter/mm2.98 0.44 1.40 0.60 0.91 0.15 42.86 一级分枝点高度
Height at first branching point/cm5.8 0.5 2.41 1.39 0.99 0.97 57.68 一级分枝数
Number of first-order branches/each2 0 0.26 0.51 1.84 2.69 196.15 开张角度
Opening angle/°59.69 3.03 27.34 19.70 −1.38 1.121 72.06 侧枝GSA
Lateral branch GSA/°48.35 10.9 24.65 11.33 0.83 −0.10 45.96 节间长度
Internode length/cm54.8 17.3 33.40 13.39 0.38 −1.05 40.08
下载: 导出CSV
表 6 不同遗传模型下QC×FJ 群体的 AIC 值
Table 6 AIC values of QC×FJ population under different genetic models
模型 Model implication 株高
Plant height地径
Diameter开张角度
Opening angle侧枝GSA
Lateral branch GSA一级分枝点高度
Height at first branching pointA-0 0MG 5188.884 391.072 1303.792 1334.032 506.3505 A-1 1MG-AD 4899.472 300.2864 1250.446 1313.571 503.0956 A-2 1MG-A 4942.795 306.6786 1248.605 1312.763 501.3233 A-3 1MG-EAD 5032.944 335.7514 1282.295 1318.325 506.3529 A-4 1MG-NCD 4993.272 331.1677 1287.728 1321.534 509.4304 B-1 2MG-ADI 4950.868 325.7444 1278.015 1314.864 515.8107 B-2 2MG-AD 4871.614 294.8911 1245.498 1305.003 497.1486 B-3 2MG-A 4910.321 302.0399 1249.31 1314.201 510.264 B-4 2MG-EA 4888.363 295.9231 1246.738 1303.926 491.0738 B-5 2MG-CD 5192.889 395.0742 1307.791 1338.033 510.3601 B-6 2MG-EAD 5190.889 393.0743 1305.791 1336.033 508.3601 加粗值代表最小 AIC 值；MG 代表主基因模型；A 代表加性效应；D 代表显性效应； I 代表上位性 效应；E 代表相等；N 代表负向；C 代表完全。
The values in bold indicate the minimum AIC value; MG represents the major gene model; A represents the additive effect; D represents the dominance effect; I represents epistatic interaction；E represents equal; N represents negative; C represents complete. The same as below.
下载: 导出CSV
表 7 不同遗传模型下FJ×QC群体的 AIC 值
Table 7 AIC values of FJ×QC population under different genetic models
模型
ModelModel implication 地径
Diameter开张角度
Opening angleA-0 0MG 123.9251 118.3907 A-1 1MG-AD 109.1332 118.0483 A-2 1MG-A 114.9762 118.7385 A-3 1MG-EAD 115.9589 119.0364 A-4 1MG-NCD 110.7568 119.0415 B-1 2MG-ADI 110.3914 122.8118 B-2 2MG-AD 100.7582 106.6611 B-3 2MG-A 111.5877 118.2622 B-4 2MG-EA 107.1847 112.1628 B-5 2MG-CD 127.9204 122.358 B-6 2MG-EAD 125.9204 120.358
下载: 导出CSV
表 8 QC×FJ、FJ×QC 群体入选模型的适合性检验
Table 8 Fitness test for QC×FJ and FJ×QC population inclusion models
群体
Group性状
Traits模型
ModelU12 U22 U32 nW2 Dn QC×FJ 株高
Plant height2MG-AD 0.0251 （0.8741 ）0.0453 （0.8314 ）0.0565 （0.8121 ）0.0393 （0.9369 ）0.0206 （0.7996 ）地径
Diameter2MG-AD 0.0004 （0.9851 ）0.0005 （0.9822 ）0.0003 （0.9866 ）0.012（1） 0.0107 （0.9998 ）开张角度
Opening angle2MG-AD 0.0015 （0.9692 ）0.0004 （0.9927 ）0.0127 （0.9103 ）0.0077 （1.0008 ）0.0184 （1）侧枝GSA
Lateral branch GSA2MG-EA 0.0008 （0.9776 ）0.0013 （0.971）0.0647 （0.7992 ）0.0206 （0.9964 ）0.032（0.993） 一级分枝点高度
Height at first branching point2MG-EA 0.0014 （0.9704 ）0.001（ 0.9748 ）0.0003 （0.986）0.064（ 0.7929 ）0.0621 （0.5055 ）FJ×QC 地径
Diameter2MG-AD 0.0013 （0.9708 ）0.0017 （0.9674 ）0.0005 （0.9828 ）0.0063 （1.0046 ）0.0267 （1）开张角度
Opening angle2MG-AD 0.0217 （0.883）0.0142 （0.9051 ）0.0087 （0.9257 ）0.0152 （0.9996 ）0.1078 （0.9871 ）U12, U22 和 U32 为均匀性检验，nW2 为 Smirnov 检验，Dn 为 Kolmogorov 检验，括号内为概率值。
U12, U22and U32 represent the uniformity tests; nW2 represents the Smirnov test; and Dn represents the Kolmogorov test, values in brackets are probability.
下载: 导出CSV
表 9 QC×FJ、FJ×QC（一年生）群体各性状在最适模型下的遗传参数
Table 9 Genetic parameters of traits in QC×FJ and FJ×QC population under optimal models
群体
Group性状
Traitsda(d) db ha(h) hb i jab jba l h2mg /% QC×FJ 株高
Plant height3.5662 1.5173 −1.452 −0.782 0 地径
Diameter0.3965 0.1399 −0.226 −0.29 0 开张角度
Opening angle8.4046 6.8204 −3.912 −3.831 0 侧枝GSA
Lateral branch GSA8.8543 23.802 一级分枝点高度
Height at first branching point0.6373 63.933 FJ×QC 地径
Diameter0.6005 0.3401 −0.238 −0.257 91.1643 开张角度
Opening angle10.415 7.135 −1.855 −3.439 95.6964 da:第 1 对主基因加性效应；db:第 2 对主基因加性效应；ha:第 1 对主基因显性效应；hb:第 2 对主基因显性效应；i:加×加效应；jab:加×显应；jba:显×加效应；l:显×显效应；h2mg:主基因遗传率。
da: the first major-gene additive effect; db: the second major-gene additive effect; ha: the first major-gene dominant effect; hb: the second major-gene dominant effect; i: additivity × additivity effect; jab: additivity × dominance effect; jba: dominance × additivity effect; l: dominance × dominance` effect; h2mg: indicates major gene heritability.
下载: 导出CSV
-
[1] 中国科学院中国植物志编辑委员会. 中国植物志-第七十三卷,第二分册[M]. 北京:科学出版社,1983. [2] EGOLF D R,ANDRICK A O. The Lagerstroemia handbook/checklist:A guide to crapemyrtle cultivars [J]. USA:American Association of Botanical Gardens and Arboreta,1978.
[3] 王金凤,柳新红,陈卓梅. 紫薇属植物育种研究进展[J]. 园艺学报,2013,40(9) :1795−1804. WANG J F,LIU X H,CHEN Z M. Research progress in breeding of Lagerstroemia plant[J]. Acta Horticulturae Sinica,2013,40(9) :1795−1804. (in Chinese)
[4] 王甜. 紫薇和川黔紫薇远缘杂交亲和性及种子萌发特性研究[D]. 长沙:中南林业科技大学,2022. WANG T. Studies on Compatibility and Seed Germination Characteristics of the Distant Hybridization of Lagerstroemia indica and Lagerstroemia excelsa (Dode) Chun. Changsha:Central South University of Forestry & Technology,2022. (in Chinese)
[5] SEGURA V,CILAS C,COSTES E. Dissecting apple tree architecture into genetic,ontogenetic and environmental effects:Mixed linear modelling of repeated spatial and temporal measures[J]. New Phytologist,2008,178(2) :302−314.
[6] SOCQUET-JUGLARD D,CHRISTEN D,DEVÈNES G,et al. Mapping architectural,phenological,and fruit quality QTLs in apricot[J]. Plant Molecular Biology Reporter,2013,31(2) :387−397.
[7] LIU R Z,AI N J,ZHU X X,et al. Genetic analysis of plant height using two immortalized populations of “CRI12 × J8891” in Gossypium hirsutum L[J]. Euphytica,2014,196(1) :51−61. DOI: 10.1007/s10681-013-1013-0
[8] ZHANG Y G,ZHU J,DAI H Y. Characterization of transcriptional differences between columnar and standard apple trees using RNA-seq[J]. Plant Molecular Biology Reporter,2012,30(4) :957−965.
[9] PETERSEN R,KROST C. Tracing a key player in the regulation of plant architecture:The columnar growth habit of apple trees (Malus × domestica) [J]. Planta,2013,238(1) :1−22. DOI: 10.1007/s00425-013-1898-9
[10] 陈卓梅,沈鸿明,张冬林,等. 从美国专利品种看中国紫薇育种发展趋势[J]. 浙江林业科技,2018,38(1) :49−61. CHEN Z M,SHEN H M,ZHANG D L,et al. Patented cultivars of Lagerstroemia in the US and trends of breeding in China[J]. Journal of Zhejiang Forestry Science and Technology,2018,38(1) :49−61. (in Chinese)
[11] 国家林业局科技发展中心(国家林业局植物新品种保护办公室) . 中国林业植物授权新品种-2016[M]. 北京:中国林业出版社,2017:113–115. [12] 胡琦,何金儒,潘会堂. 紫薇新品种[J]. 中国花卉园艺,2022(6) :32−33. DOI: 10.3969/j.issn.1009-8496.2022.6.zghhyy202206009 HU Q,HE J R,PAN H T. New Cultivar of Lagerstroemia indica[J]. China Flowers & Horticulture,2022(6) :32−33. (in Chinese) DOI: 10.3969/j.issn.1009-8496.2022.6.zghhyy202206009
[13] 赵静媛,陈发棣,滕年军,等. 地被菊匍匐性的遗传分析与RAPD标记研究[J]. 中国农业科学,2009,42(2) :734−741. ZHAO J Y,CHEN F D,TENG N J,et al. Genetic analysis and RAPD marker of creeping habits in ground-cover Chrysanthemum[J]. Scientia Agricultura Sinica,2009,42(2) :734−741. (in Chinese)
[14] 刘阳. 紫薇微卫星标记开发及矮化性状的分子标记[D]. 北京:北京林业大学,2013 LIU Y. Development of Microsatellite of Lagerstroemia and Molecular Markers of Dwarf Trait. Beijing:Beijing Forestry University,2013. (in Chinese)
[15] 石俊,陈之琳,秦波等. 紫薇匍匐平展性状相关表型遗传分析[C]. /见张启翔主编:中国观赏园艺研究进展2016,北京:中国林业出版社,2016:100-109. [16] 焦垚. 大花紫薇与紫薇杂交后代重要观赏性状遗传分析[D]. 北京:北京林业大学,2017 JIAO Y. Phenotypic and Genetic Analysis of Ornamental Traits in Hybrids of Lagerstroemia indica and L. speciosa. Beijing:Beijing Forestry University,2017. (in Chinese)
[17] 黄冬贤,黎可华,黎明星,等. 紫薇的整形修剪及其园林应用[J]. 广东园林,2010,32(4) :65−67. DOI: 10.3969/j.issn.1671-2641.2010.04.016 HUANG D X,LI K H,LI M X,et al. Pruning technology and landscape application of Lagerstroemia indica[J]. Guangdong Landscape Architecture,2010,32(4) :65−67. (in Chinese) DOI: 10.3969/j.issn.1671-2641.2010.04.016
[18] YE Y J,WU J Y,FENG L,et al. Heritability and gene effects for plant architecture traits of crape myrtle using major gene plus polygene inheritance analysis[J]. Scientia Horticulturae,2017,225:335−342.
[19] 姜珊,易星湾,徐庭亮,等. 月季花瓣数量遗传分析[J]. 植物科学学报,2021,39(2) :142−151. DOI: 10.11913/PSJ.2095-0837.2021.20142 JIANG S,YI X W,XU T L,et al. Genetic analysis of petal number in Rosa[J]. Plant Science Journal,2021,39(2) :142−151. (in Chinese) DOI: 10.11913/PSJ.2095-0837.2021.20142
[20] 吴元奇,唐兰,邱贵兰,等. 一份矮生玉米株高主基因+多基因遗传模式分析[J]. 四川农业大学学报,2022,40(3) :353−361. WU Y Q,TANG L,QIU G L,et al. Genetic pattern analysis of the height main gene and multigene of A dwarf maize[J]. Journal of Sichuan Agricultural University,2022,40(3) :353−361. (in Chinese)
[21] 王靖天,张亚雯,杜应雯,等. 数量性状主基因+多基因混合遗传分析R软件包SEA v2.0[J]. 作物学报,2022,48(6) :1416−1424. WANG J T,ZHANG Y W,DU Y W,et al. SEA v2.0:An R software package for mixed major genes plus polygenes inheritance analysis of quantitative traits[J]. Acta Agronomica Sinica,2022,48(6) :1416−1424. (in Chinese)
[22] 张向前,刘景辉,齐冰洁,等. 燕麦种质资源主要农艺性状的遗传多样性分析[J]. 植物遗传资源学报,2010,11(2) :168−174. ZHANG X Q,LIU J H,QI B J,et al. Cluster diversity analysis of the main agronomic traits in oat germplasm[J]. Journal of Plant Genetic Resources,2010,11(2) :168−174. (in Chinese)
[23] 秦波,李素珍,胡玲,等. 基于转录组的紫薇SSR分子标记开发[C]//中国观赏园艺研究进展2018. 哈尔滨,2018:526–532. [24] 郭长花,康向阳. 树木发育中的阶段转变研究进展[J]. 生物技术通讯,2008,19(5) :784−786. GUO C H,KANG X Y. Progress in phase change during tree development[J]. Letters in Biotechnology,2008,19(5) :784−786. (in Chinese)
[25] GREENWOOD M S. Phase change in loblolly pine:Shoot development as a function of age[J]. Physiologia Plantarum,1984,61(3) :518−522. DOI: 10.1111/j.1399-3054.1984.tb06366.x
[26] 戴思兰. 园林植物育种学[M]. 北京:中国林业出版社,2007. [27] 叶福民,马策,张晓波,等. 多头切花秋菊杂交F1群体观赏性状主基因+多基因混合遗传分析[J]. 辽宁农业科学,2022(6) :7−12. YE F M,MA C,ZHANG X B,et al. Heterosis of ornamental characters and mixed inheritance of major gene and polygene in hybrid progenies of cut Chrysanthemum[J]. Liaoning Agricultural Sciences,2022(6) :7−12. (in Chinese)
[28] 张琳,郭丽丽,郭大龙,等. 牡丹杂交F1代性状分离规律及混合遗传分析[J]. 南京林业大学学报(自然科学版) ,2018,42(6) :51−60. ZHANG L,GUO L L,GUO D L,et al. Separation analysis and mixed genetic analysis of phenotypic traits in F1 progenies of tree peony[J]. Journal of Nanjing Forestry University (Natural Sciences Edition) ,2018,42(6) :51−60. (in Chinese)
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 1
