Early growth model fitting and growth rhythm of Eucalyptus pellita
-
摘要:目的
利用来自巴布亚新几内亚南部和澳大利亚东北部9个粗皮桉Eucalyptus pellita种源的试验林材料开展生长节律研究，为粗皮桉良种选育及经营管理提供科学依据。
方法通过随机区组设计试验，利用巢氏方差分析法分析粗皮桉9个种源共98个家系的生长节律变异规律。获取各种源地气候环境因子信息，运用Pearson法分析生长性状与地理及气候因子的相关性。对比Logistic、Gompertz和Von Bertalanffy 3种模型的拟合结果，选择最佳模型计算粗皮桉的生长节律。采用系统聚类方法开展优良家系筛选。
结果2.5年生粗皮桉树高和胸径在种源、家系水平差异显著(P<0.05)。种源树高与等温性呈显著正相关，呈典型的纬向变异模式。Logistic模型对粗皮桉树高和胸径的拟合效果最佳，利用Logistic模型计算生长节律得出树高的最大加速时期(T1)、最大减速时期(T2)和线性生长时期(L)平均值分别为96、636和540 d，胸径的T1、T2和L平均值分别为165、681和516 d。利用2.5年生粗皮桉的生长性状进行系统聚类分析，从98个家系中筛选出一般家系29个、中等家系30个、优秀家系39个；结合生长节律T1的聚类分析，将3个等级的家系均再划分为Ⅰ、Ⅱ、Ⅲ 3类。
结论2.5年生粗皮桉在种源、家系水平上生长性状与生长节律变异丰富，利用Logistic模型计算的生长节律可以为粗皮桉早期抚育管理及良种选育提供理论依据。
Abstract:ObjectiveTo investigate growth rhythm using nine Eucalyptus pellita provenances in southern Papua New Guinea and northeastern Australia, provide a scientific basis for the breeding and management of E. pellita.
MethodA randomized block design experiment was conducted, and nested analysis of variance was utilized to analyze the growth rhythm variation patterns of 98 families from nine E. pellita provenances. The information of climate and environment factors in various sources were obtained, and Pearson method was used to reveal the correlation between growth traits and geography/climate factors. The fitting results of Logistic, Gompertz and Von Bertalanffy models were compared to select the optimal model, and calculate the growth rhythm of E. pellita. The excellent family selection was carried out by systematic clustering method.
ResultSignificant differences were observed in the tree height and DBH of 2.5-year-old E. pellita at provenance and family levels (P<0.05). Furthermore, a significant and positive correlation was found between isothermality and tree height of the provenances, exhibiting typical zonal variation patterns. The Logistic model had the best fitting effect on tree height and DBH of E. pellita. The Logistic model calculated the growth rhythms, with the average values of the maximum acceleration period (T1), maximum deceleration period (T2), and linear growth period (L) for tree height being 96, 636 and 540 d respectively, as well as for DBH being 165, 681 and 516 d respectively. A systematic cluster analysis was conducted using the growth traits of 2.5-year-old E. pellita, and 29 general families, 30 medium families and 39 outstanding families were selected from 98 families. All the families of the three grades were further classified into three categories of I, II and III according to the cluster of growth rhythm T1.
ConclusionThe growth traits and growth rhythm of 2.5-year-old E. pellita exhibited considerable variation at provenance and family levels. The growth rhythm calculated using Logistic model can provide a theoretical basis for optimizing seedling rearing and management, as well as improved variety breeding of E. pellita.
-
Keywords:
- Eucalyptus pellita /
- Provenance /
- Family /
- Growth model /
- Growth rhythm
-
猪繁殖与呼吸综合征(Porcine reproductive and respiratory syndrome，PRRS)是一种以母猪繁殖障碍和各种生长阶段猪的呼吸系统疾病为特征的病毒性传染病，俗称“猪蓝耳病”，其病原为猪繁殖与呼吸综合征病毒(Porcine reproductive and respiratory syndrome virus，PRRSV)[1-2]。PRRS于20世纪80年代末首先在美国被发现，临床上主要表现为母猪晚期繁殖障碍和新生仔猪的呼吸道疾病[3]。随后该病相继在北美、欧洲、亚洲的养猪国家广泛流行，并引发了空前的“流产风暴”[4]，给世界养猪业造成了巨大的经济损失。
2006年，在我国南方部分省份出现的“猪无名高热”病，经鉴定其病原为以JXA1株为代表的高致病性PRRSV，引发了严重的PRRS疫情[5-6]。在PRRSV的遗传进化树中，HP-PRRSV和之前流行的经典毒株同属于Lineage 8，其重要特征是Nsp2编码区存在30个不连续氨基酸的缺失[7-8]，这种新出现的毒株具有毒力强、致死率高、传播速度快等特点，至今依然是国内流行的优势毒株[9]，给我国养猪业造成非常严重的损失。
反向遗传学技术是研究RNA病毒的有力工具，该技术通过对相关病毒基因进行直接改造后构建出病毒的感染性克隆，并经RNA-launched和DNA-launched 2种途径对病毒进行拯救[10-11]。RNA-launched途径是RNA病毒基因组先反转录成cDNA，接着克隆到合适的转录载体上，通过克隆到T7或者SP6启动子的下游，先在体外转录出病毒全长的基因组RNA，随后将这些体外的转录本转染宿主细胞，在宿主细胞内包装出活的病毒粒子[12-13]。DNA-launched途径则是将经反转录的病毒cDNA克隆到真核启动子如CMV启动子的下游，并将构建好的质粒不经体外转录，直接转染宿主细胞，利用宿主细胞的核功能转录出病毒基因组，进而完成病毒粒子的包装[14]。相比RNA-launched途径，DNA-launched途径避开了体外转录步骤，不仅操作简便，还降低了转染过程中RNA降解的风险，使转染的效率大大提高。目前，反向遗传技术在研究PRRSV的结构与功能、致病机制以及研发新型疫苗等方面有着广泛应用[15-17]。
本研究运用反向遗传操作技术将HP-PRRSV JXA1株的全基因组分段克隆至改造过的低拷贝载体pOKq上，成功构建了JXA1株的全长感染性克隆。并采取基于DNA-launched途径成功获得拯救的病毒，并对拯救病毒进行生物学特性的分析。HP-PRRSV代表株JXA1株反向遗传平台的构建，为HP-PRRSV后续的结构功能探索、致病机理研究以及新型疫苗研发奠定了基础。
1. 材料与方法
1.1 材料
1.1.1 载体、毒株、细胞、抗体
克隆载体pJET1.2购买自Thermo Fisher Scientific。由pOK12改造具有多克隆位点的低拷贝载体pOKq、XH-GD全长感染性克隆质粒、PRRSV易感细胞Marc-145细胞、转染用BHK-21细胞均由华南农业大学兽医学院传染病教研室保存。HP-PRRSV JXA1毒株由河南农业大学田克恭教授惠赠。PRRSV N蛋白单克隆抗体购买自韩国金诺公司。FITC标记的山羊抗小鼠IgG二抗购买自北京中杉金桥生物技术有限公司。
1.1.2 主要试剂
RNA抽提试剂盒购自上海飞捷生物技术有限公司，反转录试剂盒、Taq酶、高保真DNA聚合酶购自宝日医生物技术(北京)有限公司。限制性核酸内切酶购自NEB公司。T4 DNA连接酶购自Thermo Fisher Scientific。DL2000 DNA marker、250 bp DNA ladder (Dye plus)、克隆菌感受态细胞Escherichia coli DH5α、克隆载体pMD18-T vector购自宝日医生物技术(北京)有限公司；DNA凝胶回收试剂盒、质粒抽提试剂盒购自OMEGA公司。DMEM高糖培养基、胰蛋白酶购自Biological Industries。胎牛血清购自Gibco公司。二甲基亚砜(DMSO)购自Amresco公司。
1.2 方法
1.2.1 全长cDNA构建方案与引物设计
构建方案：JXA1毒株全长cDNA的构建见图1。先将CMV、A1、A2片段融合，融合产物再与A3片段连接，构成A片段。将D1片段和终止信号肽序列BGH片段融合构成D片段，然后将D、C、B、A共4个片段依次亚克隆到低拷贝载体pOKq上。
![]()
图 1 JXA1毒株全长cDNA克隆的构建策略Figure 1. The strategy for construction of the full-length cDNA clone of JXA1 strain根据GenBank中收录的PRRSV JXA1毒株全基因组序列，利用Oligo7.0软件 设计6对引物对病毒的全长基因进行分段克隆，获得A1、A2、A3、B、C和D1段。同时以PRRSV XH-GD全长感染性克隆质粒为模板，设计引物扩增CMV和BGH片段。同时在扩增A1片段时引入1个沉默突变，将全基因组第510位碱基由A突变为G，引入遗传标记FseI，以便与亲代病毒相区别。引物序列由上海英潍捷基生物技术有限公司合成。各段引物的序列、扩增片断大小及引物名称见表1。
表 1 构建JXA1株全长感染性克隆设计的引物序列Table 1. Primer sequences used for construction of the full-length infectious clone of JXA1 strain引物名称 Primer name 引物序列(5′→3′)1) Primer sequence 产物长度/bp Product length CMV-F AGTC GCGGCCGCCGATGTACGGGCCAGATATAC 757 CMV-R ATAGAGCCAACACCTATACGTCA A1-F TGACGTATAGGTGTTGGCTCTAT 528 A1-R TAGGGGTGAGTTGGCCGGCCCTGTA A2-F ATCTACAGGGCCGGCCAACTCAC 2 193 A2-R CGGCGGTGTCTCGAGAATCATCT A3-F CAAAGATGATTCTCGAGACACC 3 096 A3-R AAGGCATAGGTGCTTAAGT B-F CATCCGTTACAGGTGGTAATG 2 951 B-R AGCACAGCTGTCCATGAGTATTGC C-F GAAGACAACCATCACAGACTCAC 2 713 C-R AAACCTCATGCTGGTGGCATT D1-F AGGACTGGGAGGATTACAATGAT 3 341 D1-R GTGCCAAAGAATGCCAGCCCATCATG BGH-F CATGATGGGCTGGCATTCTTTGGCAC 533 BGH-R GGGATCGA GGTACCCAGAAGC 1)下划线部分为酶切位点
1) The underlined portion is the restricted enzyme digestion site1.2.2 JXA1毒株全基因组各片段的扩增
以XH-GD全长感染性克隆质粒为模板，分别用引物(CMV-F、CMV-R)和(BGH-F、BGH-R)扩增CMV启动子序列和BGH终止信号肽序列；按上海飞捷生物技术有限公司的RNA抽提试剂盒操作说明书提取JXA1细胞病毒的总RNA，用适量RNase free水溶解，获得总RNA，用SuperScriptTM III逆转录酶和OligodT(20)(Invitrogen，USA)，按照操作说明书进行cDNA合成，然后将获得的cDNA作为模板，用表1中扩增引物进行各片段的扩增，经过RT-PCR，获得A1、A2、A3、B、C和D1片段，参照DNA purification kit说明书纯化扩增的片段，利用融合PCR技术，用引物(A1-F、A2-R)融合A1和A2片段，构建点突变以引入全长感染性克隆的遗传标记FseI酶切位点，以回收的CMV片段和A1、A2融合片段为模板，用引物(CMV-F、A2-R)融合CMV、A1、A2片段并回收；以回收的BGH片段和D1片段为模板，用引物(D1-F、BGH-R)融合D1和BGH片段构成D片段并回收，PCR按常规方法进行，并设立空白对照。PCR产物用10 g/L琼脂糖凝胶(含0.5 μg/mL EB)电泳检测并回收纯化。另外，为了提高目的片段的保真性，构建感染性克隆过程中的所有PCR反应均使用宝日医生物技术(北京)有限公司的PrimeSTAR® Max高保真DNA聚合酶。
1.2.3 质粒pOK-A、pJET-B、pJET-C和pJET-D的构建
将CMV、A1、A2的融合回收产物连接到pOKq空载体上形成载体pOK-CMV-A1-A2，再将扩增的A3片段连接到pOK-CMV-A1-A2载体上，构建出pOK-A中间载体。将D、C和B这3个片段的纯化产物分别与pJET1.2载体连接。取5 μL的连接产物转化克隆菌感受态细胞E. coli DH5α后，用灭菌枪头挑取可疑菌落于含氨苄青霉素的LB液体培养基中，于37 ℃条件下振荡培养4~6 h，取适量菌液进行PCR鉴定。将PCR鉴定为阳性的菌液送上海英潍捷基生物技术有限公司进行测序验证。测序正确的重组质粒分别命名为pJET-B、pJET-C和pJET-D，并按照OMEGA公司的质粒提取试剂盒说明书进行质粒抽提。
1.2.4 JXA1毒株全长cDNA的构建
将pJET-D、pJET-C、pJET-B和pOK-A按照图1所示的顺序依次亚克隆到pOKq载体，构建的质粒经PCR鉴定为阳性的重组菌液送上海英潍捷基生物技术有限公司进行测序，将测序正确的克隆质粒命名为pOK-JXA1，并参照OMEGA公司去内毒素质粒提取试剂盒说明书进行去内毒素质粒抽提。
1.2.5 JXA1株病毒的拯救
将JXA1株病毒全长cDNA的重组质粒pOK-JXA1转染BHK-21细胞，转染使用Lipofectamine 3000(后文简称Lipo 3000)脂质体，具体操作步骤参照说明书。转染72 h后，将6孔板反复冻融3次，收集细胞液，4 ℃、7 000 r/min离心7 min，取上清接种于单层Marc-145细胞，吸附1 h后补加含2%(φ)FBS的DMEM，置于37 ℃、CO2体积分数为5%的培养箱培养，观察细胞病变情况(Cytopathic effect，CPE)，出现病变后冻融收毒，连续传12代。将拯救的病毒命名为rJXA1。
1.2.6 拯救病毒的鉴定
为验证所看到的细胞病变是拯救病毒所致，排除亲本病毒污染的可能，利用引入的沉默突变FseI酶切位点进行鉴定。取第3代的病毒细胞培养悬液进行总RNA的提取，再取2 μg RNA进行反转录。以得到的cDNA为模板，以A1-F和A2-R为引物进行PCR扩增。PCR产物用10 g/L琼脂糖凝胶电泳检测并回收纯化，将纯化后的PCR产物送上海英潍捷基生物技术有限公司进行测序鉴定。
同时将拯救病毒的第3代rJXA1-P3和亲本病毒的第3代JXA1-P3接种至长满单层的Marc-145细胞6孔板中并设1孔阴性对照，置于37 ℃、CO2体积分数为5%培养箱中培养72 h后对拯救的病毒进行间接免疫荧光(Indirect immunofluoresence assay，IFA)试验鉴定。
1.2.7 拯救病毒的稳定性
将拯救病毒rJXA1在Marc-145细胞上连续传至12代，观察每一代的CPE情况。同时从拯救获得的病毒中，选择第1、3、6、9和12代病毒进行TCID50测定。将Marc-145细胞接种到96孔细胞培养板上，待细胞汇合度达到90%以上时，将获得的相应代次的拯救病毒进行10倍梯度倍比稀释，即用DMEM进行10−1~10−8倍比稀释，分别接种于Marc-145细胞中，每个稀释度接种8孔，每孔100 μL。同时设未接毒的空白对照，将细胞板放置37 ℃、CO2体积分数为5%细胞培养箱作用90 min。弃去病毒液，加入含2%(φ)FBS的DMEM，放置于37 ℃、CO2体积分数为5%细胞培养箱。连续观察2~7 d，记录每个稀释度细胞病变的孔数，本试验重复3次求其平均值，并按照Reed-Muench法计算各个样品的TCID50。
1.2.8 拯救病毒的生长曲线
测定rJXA1第3代病毒和亲本病毒JXA1第3代病毒在Marc-145细胞中的生长曲线。将rJXA1第3代病毒和亲本病毒JXA1第3代病毒用DMEM培养基稀释至感染复数(Multiplicity of infection，MOI)为0.1，接种于长成单层的Marc-145细胞，置于37 ℃、CO2体积分数为5%培养箱中孵育1 h，每个病毒各做3次重复。孵育结束后弃去上清，用PBS缓冲液洗涤细胞2次，加入5 mL含2%(φ)FBS的DMEM细胞维持液，置于37 ℃、CO2体积分数为5%的培养箱中培养。分别于接种病毒12、24、36、48、60、72、84和96 h后，从每瓶细胞中吸取200 μL上清，然后测定各时相上清中病毒的TCID50，绘制病毒的生长曲线。
2. 结果与分析
2.1 JXA1毒株全基因组各片段的扩增
以XH-GD全长感染性克隆质粒为模板扩增CMV启动子序列和BGH终止信号肽序列，以JXA1反转录后的cDNA为模板扩增JAX1的A1、A2、A3、B、C和D1段。PCR产物经琼脂糖凝胶电泳后，可见相应的目的条带，其中CMV、BGH和A1片段的扩增产物大小分别约为750、500和500 bp，与预期大小相符(图2A)。A2、A3、B、C和D1片段的扩增产物大小分别约为2 200、3 000、3 000、2 700和3 300 bp，与预期大小相符(图2B)。
![]()
图 2 JXA1株各基因片段的PCR扩增与构建M1：DL2000 DNA marker；M2：DL5000 DNA marker；M3：DL10000 DNA marker；1：CMV片段；2：BGH片段；3：A1片段；4：B片段；5：C片段；6:A3片段；7：A3片段酶切产物；8：D1片段；9：D1和BGH融合片段；10：A2片段；11~12：CMV、A1和A2融合片段；13：pOK-CMV-A1-A2的酶切产物；14~18：菌液PCR鉴定；N：阴性对照Figure 2. PCR amplification and construction of each gene fragment of JXA1 strainM1: DL2000 DNA marker; M2: DL5000 DNA marker；M3：DL10000 DNA marker; 1: CMV fragment; 2: BGH fragment; 3 A1 fragment; 4: B fragment; 5: C fragment; 6: A3 fragment; 7: Enzyme digestion of A3 fragment; 8: D1 fragment; 9: D1 and BGH fusion fragment; 10: A2 fragment; 11−12: CMV, A1 and A2 fusion fragment; 13: Enzyme digestion of pOK-CMV-A1-A2; 14−18: PCR identification of the bacterial liquid; N: Negative control利用融合PCR技术融合A1和A2片段，构建点突变来引入全长感染性克隆的遗传标记FseI酶切位点，以回收的CMV片段和A1、A2融合片段为模板，对CMV、A1、A2片段进行融合PCR扩增并回收，PCR产物经琼脂糖凝胶电泳后，可见约为3 800 bp扩增片段，与预期大小相符(图2C)。将CMV、A1、A2的融合回收产物连接到pOKq空载体上形成载体pOK-CMV-A1-A2，再将扩增的A3片段连接到pOK-CMV-A1-A2载体上，并进行菌液PCR鉴定(图2C)，构建出pOK-A中间载体；以回收的BGH和D1段为模板，对D1和BGH片段进行融合PCR扩增并回收，PCR产物经琼脂糖凝胶电泳后，可见约为3 800 bp扩增片段，与预期大小相符(图2B)，然后将D、C、B 片段依次克隆到pJET1.2载体上。
2.2 JXA1毒株全长cDNA的构建与鉴定
将pJET-D、pJET-C、pJET-B和pOK-A按照图1所示的顺序依次连接到pOKq载体，构建全长cDNA的重组质粒pOK-JXA1，用引物CMV-F、A1-R对重组质粒pOK-JXA1进行PCR鉴定，PCR产物经琼脂糖凝胶电泳后，可见约为1 300 bp的扩增条带，与预期大小相符(图3)，鉴定阳性的重组菌液送往上海英潍捷基生物技术有限公司进行测序验证，测序正确的重组质粒命名为pOK-JXA1。
![]()
图 3 pOK-JXA1菌液PCR鉴定M：250 bp DNA ladder (Dye plus)；1~8：pOK-JXA1菌液PCR扩增产物；N：阴性对照Figure 3. PCR identification of pOK-JXA1 bacterial liquidM: 250 bp DNA ladder (Dye plus); 1−8: PCR amplification of pOK-JXA1 bacterial liquid; N: Negative control2.3 JXA1株病毒的拯救
将质粒pOK-JXA1转染BHK-21细胞72 h后，冻融收获病毒，离心取上清后接种至Marc-145细胞，置于37 ℃、CO2体积分数为5%的培养箱中培养72 h，显微镜下可见与亲本病毒JXA1类似的CPE(图4)，将拯救病毒命名为rJXA1。
![]()
图 4 拯救病毒产生的细胞病变A：拯救病毒rJXA1产生的CPE；B：亲本病毒JXA1产生的CPE；C：对照细胞Figure 4. Cytopathic effect (CPE) produced in the rescued virusesA: CPE produced in the rescued virus rJXA1; B: CPE produced in the parental virus JXA1; C: Control cell2.4 拯救病毒的鉴定
将拯救获得的rJXA1在Marc-145细胞中连续传代。取第3代病毒(rJXA1-P3)的培养液抽提RNA进行RT-PCR鉴定，PCR产物经琼脂糖凝胶电泳后可检测到强阳性条带，大小与预期相符(图5)。将纯化的PCR产物送测序，测序结果显示拯救病毒rJXA1全基因组的第510位核苷酸由A突变为G，对应的阅读框由AGA变为AGG，推导的氨基酸均为精氨酸，产生了新的酶切位点FseI(GGCCGGCC)，与预期结果一致，表明拯救病毒的遗传标记被成功引入。
![]()
图 5 拯救病毒的遗传标记位点RT-PCR鉴定M：250 bp DNA ladder (Dye plus)；1~3：rJXA1-P3的3个重复的扩增产物；N：阴性对照Figure 5. Identification of genetic marker of the rescued virus by RT-PCRM: 250 bp DNA ladder (Dye plus); 1−3: Amplification products of three replicates of rJXA1-P3; N: Negative control用PRRSV N蛋白单克隆抗体对亲本病毒以及拯救病毒进行间接免疫荧光鉴定，结果显示：接种亲本病毒和拯救病毒的细胞在荧光显微镜下可观察到特异性的绿色荧光，而阴性对照细胞没有荧光产生(图6)，证明病毒拯救成功。
![]()
图 6 拯救病毒的间接免疫荧光鉴定A：亲本病毒；B：拯救病毒；C：空白对照Figure 6. Identification of the rescued virus through indirect immunofluoresence assayA: Parental virus; B: Rescued virus; C:Blank control2.5 拯救病毒的稳定性
将拯救病毒rJXA1在Marc-145细胞上连续传至12代，每一代均可以产生明显CPE。图7为第1、3、6、9和12代的细胞病变情况。同时对第1、3、6、9和12代病毒进行TCID50的测定，结果表明，拯救病毒可在Marc-145细胞上稳定传代，病毒滴度介于106.0~107.3 TCID50/mL。
![]()
图 7 rJXA1在Marc-145细胞中引起的细胞病变A：第1代；B：第3代；C：第6代；D：第9代；E：第12代；F：空白对照Figure 7. Cytopathic effect (CPE) induced by rJXA1 in Marc-145 cellsA: The first generation; B: The third generation; C: The sixth generation; D: The ninth generation; E: The 12th generation; F: Blank control2.6 拯救病毒的生长曲线
对拯救病毒rJXA1的第3代病毒和亲本病毒JXA1的第3代病毒进行生长曲线的测定，结果表明拯救病毒rJXA1第3代和亲本病毒JXA1第3代的生长曲线相似，二者达到最高滴度的时间相同，均为72 h(图8)。
![]()
图 8 拯救病毒和亲本病毒的生长曲线Figure 8. The growth curves of the rescued virus and parental virus.3. 讨论与结论
2006年以来，以感染猪高热、高发病率和高死亡率为特征的高致病性猪繁殖与呼吸综合征在我国爆发，并在全国范围流行。截至目前，该病依然在国内广泛流行，给我国养猪业带来巨大的经济损失[18]。
PRRSV基因组为不分节段的单股正链RNA，全长15 kb左右。构建其全长感染性克隆时需要避免基因突变、缺失等破坏保真性的情况出现。为克服这一难点，本研究在构建过程中全程使用PrimeSTAR® Max高保真聚合酶进行基因的扩增或融合PCR扩增。此外，为了避免长片段扩增导致的保真性降低，本研究还将整个病毒的基因组分为4个节段分段克隆，逐步连接到严谨性较高的低拷贝载体pOKq上，成功构建了JXA1株的全长感染性克隆。并采取基于DNA-launched途径进行病毒拯救，成功构建了HP-PRRSV JXA1株的反向遗传平台，获得了感染性拯救病毒rJXA1。该拯救病毒与亲本毒株JXA1具有相同的体外增殖能力。
为了区分亲本病毒和拯救病毒，需要在拯救病毒的基因中加入遗传标记。本研究选择了20株有代表性的基因II型PRRSV毒株，使用MegAlign软件对其进行全基因组序列比对。最终选择第510位碱基作为遗传标记位点。20株基因II型代表株在该位点的碱基均为A，对应的氨基酸密码子为AGA，编码精氨酸，将其突变为G后不影响编码的氨基酸，同时可以形成一个新的单酶切位点FseI(GGCCGGCC)，便于对拯救病毒和亲本病毒的鉴定。
准确且完整的全长感染性克隆的构建是病毒拯救的关键[19]。本研究采取基于DNA-launched途径进行病毒拯救，避开了RNA体外转录的步骤，不仅减少了试验操作的影响，也简化了整个试验流程。另外，本研究还在JXA1病毒的基因组两端添加了核酶序列，全长质粒转染至BHK-21细胞后在细胞内转录出RNA，在核酶的作用下直接剪切，形成不含外源基因的病毒基因组。PRRSV基因组3′端存在1个ploy A尾，ploy A的个数对病毒转录的稳定性有重要作用，如果ploy A尾上A的个数过少，有可能导致转录失败，无法产生具有感染性的病毒粒子。基于广东省动物源性人兽共患病预防与控制重点实验室以前建立RNA-launched反向遗传平台时的经验，本研究选择了42个A构建病毒基因组的ploy A尾。
虽然国内已使用高致病性PRRSV弱毒活疫苗来控制HP-PRRSV的流行，但HP-PRRSV毒株依然在国内广泛流行，未呈现明显的减弱趋势[20]，其他亚群(如NADC30-like亚群和GM2-like亚群)呈现上升趋势[21-22]，所以很有必要对其原因进行研究。本研究成功构建的HP-PRRSV JXA1株的反向遗传平台将为进一步研究HP-PRRSV的致病机理、基因功能以及新型疫苗研发奠定了基础，同时可为PRRSV的防控提供技术参考。
致谢：感谢河南农业大学田克恭教授惠赠HP-PRRSV JXA1毒株，感谢华南农业大学兽医学院陈耀博士在本研究的选题上给予的指导与建议！
-
![]()
图 1 粗皮桉树高和胸径表型变异
各小图中，黑色圆点表示离群值，灰色钟形曲线表示数据分布。
Figure 1. Phenotypic variation of tree height and DBH in Eucalyptus pellita
In each figure, the black dots represent the outliers, and the gray bell curves represent the data distribution.
![]()
图 2 粗皮桉树高和胸径Logistic拟合曲线
各小图中，黑点表示实际观测值，灰线表示684株个体的拟合曲线。
Figure 2. Logistic fitting curve of tree height and DBH in Eucalyptus pellita
In each figure, the black dots represent the actual observed values, and the gray lines represent the fitted curves of 684 samples.
![]()
图 3 粗皮桉生长性状与生长节律的相关性
“*”和“**”分别表示树高或胸径与生长节律在0.05和0.01水平显著相关(Pearson法)。
Figure 3. Correlation between growth trait and growth rhythm of Eucalyptus pellita
“*” and “**” indicate tree height or DBH is significantly correlated with growth rhythm at 0.05 and 0.01 levels respectively (Pearson method).
![]()
图 4 粗皮桉基于生长性状(A)与生长节律(B、C、D)的聚类分析
B、C和D分别表示一般、中等和优秀家系，黑色虚线表示分类阈值。
Figure 4. Cluster analysis according to growth traits (A) and growth rhythm (B, C, D) in Eucalyptus pellita
B, C and D represent general, medium and outstanding families respectively; Black dotted lines represent classification thresholds.
表 1 本研究粗皮桉种源信息
Table 1 The information of Eucalyptus pellita provenances in this study
种源
Provenance地区1)
District纬度
Latitude经度
Longitude家系数量
No. of families种批号
Seed lot code海拔/m
AltitudeConn L.A./Cardwell (Con) AUS 18º26′S 146º07′E 7 1111 29 Mt Ray (MR) AUS 15º12′S 144º58′E 11 1109 252 CREB Track (CT) AUS 16º07′S 145º18′E 13 1162 223 Black Mountain Road (BMR) AUS 16º41′S 145º31′E 17 1163 502 Helenvale (Hel) AUS 15º44′S 145º12′E 14 1161 214 SW Tully (ST) AUS 17º54′S 145º41′E 12 1110 41 WONGA-DAINTREE (WD) AUS 16º16′S 145º22′E 5 18750 15 MERU WP (MW) PNG 08º27′S 141º28′E 7 18750 40 PNG Seed Production Area (PSP) PNG 08º20′S 141º26′E 12 688 30 1) AUS：澳大利亚，PNG：巴布亚新几内亚。
1) AUS: Australia, PNG: Papua New Guinea.
下载: 导出CSV
表 2 t7时期粗皮桉树高和胸径的方差分析
Table 2 Variance analysis of tree height and DBH in Eucalyptus pellita at t7 period
生长性状
Growth trait变异来源
Source of variation平方和
Sum of square自由度
df均方
Mean squareF P 树高 Tree height 种源 25.192 8 3.149 3.482 0.001 家系(种源) 117.906 89 1.216 1.344 0.022 胸径 DBH 种源 45.652 8 5.706 2.966 0.003 家系(种源) 349.122 89 3.599 1.871 0.000
下载: 导出CSV
表 3 种源地气候因子与t7时期粗皮桉生长性状的相关性分析
Table 3 Correlation analysis between provenance climate factors and growth traits of Eucalyptus pellita at t7 period
生物气候变量
Bioclimate variable皮尔逊相关系数1) Pearson correlation coefficient 树高 Tree height 胸径 DBH 等温性 Isothermality 0.785* 0.587 季节性气温 Seasonal temperature −0.738* −0.659 年温差范围 Annual temperature range −0.671* −0.677* 纬度 Latitude −0.759* −0.543 经度 Longitude −0.767* −0.503 树高 Tree height 1.000 0.711* 胸径 DBH 0.711* 1.000 1) “*”表示t7时期树高或胸径与生物气候变量在0.05水平显著相关(Pearson法)。
1) “*” indicates that tree height or DBH at t7 period is significantly correlated with bioclimatic variables at 0.05 level (Pearson method).
下载: 导出CSV
表 4 粗皮桉树高和胸径拟合曲线模型比较
Table 4 Comparison of tree height and DBH fitting curve models in Eucalyptus pellita
模型
Model生长性状
Growth trait生长参数1) Growth parameter 决定系数
R2赤池信息量准则
AIC均方根误差
RMSEa b r Logistic 树高 7.65 5.22 0.004 68 0.987 3.137 0.171 胸径 8.09 7.42 0.004 90 0.999 −10.895 0.063 Gompertz 树高 8.22 2.20 0.003 02 0.982 5.555 0.203 胸径 8.77 2.71 0.003 10 0.999 −10.422 0.065 Von Bertalanffy 树高 8.57 0.56 0.002 48 0.980 6.313 0.214 胸径 9.22 0.66 0.002 50 0.998 −8.661 0.074 1) a：极值，b：待定系数，r：生长速率。
1)a: Extreme value, b: Undetermined coefficient, r: Growth rate.
下载: 导出CSV
表 5 粗皮桉树高生长参数以及生长节律1)
Table 5 Growth parameters and growth rhythm of tree height in Eucalyptus pellita
参数
Argumenta/m r T1/d T2/d L/d MGR LGR TLG 平均值 Mean 7.8 0.005 13 96 636 540 0.009 8 0.008 7 4.481 6 最小值 Min 4.2 0.001 74 2 301 240 0.003 5 0.003 1 2.437 9 最大值 Max 15.9 0.010 99 349 1 556 1 510 0.016 9 0.015 0 9.168 5 四分位数 Quartile 1.6 0.001 23 75 145 124 0.002 0 0.001 8 0.897 8 变异系数/% CV 17 20 58 24 27 18 18 17 1) a：极值，r：生长速率，T1：最大加速时期，T2：最大减速时期，L：线性生长时期，MGR：最大线性生长速率，LGR：线性生长速率，TLG：线性生长量。
1)a: Extreme value, r: Growth rate, T1: Timing of the maximum acceleration, T2: Timing of the maximum deceleration, L: Duration of linear growth, MGR: The maximum linear growth rate, LGR: Linear growth rate, TLG: Total linear growth.
下载: 导出CSV
表 6 粗皮桉胸径生长参数以及生长节律1)
Table 6 Growth parameters and growth rhythm of DBH in Eucalyptus pellita
参数
Argumenta/m r T1/d T2/d L/d MGR LGR TLG 平均值 Mean 8.4 0.005 53 165 681 516 0.011 0 0.009 8 4.789 5 最小值 Min 4.0 0.001 77 6 320 183 0.004 1 0.003 7 2.314 1 最大值 Max 44.5 0.014 40 1 016 2 502 1486 0.021 5 0.019 1 13.884 2 四分位差 Quartile 2.3 0.001 48 93 147 133 0.002 9 0.002 6 1.317 2 变异系数/% CV 31 27 52 28 33 21 21 26 1) a：极值，r：生长速率，T1：最大加速时期，T2：最大减速时期，L：线性生长时期，MGR：最大线性生长速率，LGR：线性生长速率，TLG：线性生长量。
1)a: Extreme value, r: Growth rate, T1: Timing of the maximum acceleration, T2: Timing of the maximum deceleration, L: Duration of linear growth, MGR: The maximum linear growth rate, LGR: Linear growth rate, TLG: Total linear growth.
下载: 导出CSV
表 7 基于T1时期树高和胸径聚类的粗皮桉家系分类
Table 7 Classification of Eucalyptus pellita families based on tree height and DBH cluster at T1 period
类别
Category一般家系(29)
General family中等家系(30)
Medium family优秀家系(39)
Outstanding familyⅠ ST12、Hel10、Con6、BMR16、BMR14 MW1、BMR2、BMR3、CT2、PSP4、Hel2、MR1、MR3 WD3、CT6、CT10、PSP5、PSP10 Ⅱ WD5、ST9、ST11、ST10、PSP12、MR9、MR8、MR10、Hel14、Hel13、Hel12、CT12、CT11、Con5、Con3、BMR17、BMR15、BMR12 WD1、MW2、MW3、MW5、BMR1、BMR5、BMR6、CT1、CT3、ST1、ST2、ST3、ST4、PSP1、PSP2、PSP3、Hel3、Hel4、MR2 WD2、WD4、MW6、MW7、BMR7、BMR10、BMR11、CT5、CT8、ST6、ST7、ST8、PSP6、PSP8、PSP9、PSP11、Con1、Hel5、Hel7、Hel8、Hel9、MR4、MR5、MR6、MR7 Ⅲ MR11、Hel11、CT13、Con7、Con4、BMR13 MW4、BMR4、Hel1 BMR8、BMR9、CT4、CT7、CT9、ST5、PSP7、Con2、Hel6
下载: 导出CSV
-
[1] 刘涛, 谢耀坚. 中国桉树人工林快速发展动因分析与展望[J]. 桉树科技, 2020, 37(4): 38-47. [2] 梁庭辉. 桉树大径材培育的探索与思考[J]. 农业与技术, 2018, 38(15): 99-100. [3] HARWOOD C E. Eucalyptus pellita: An annotated bibliography[M]. Victoria, Australia: CSIRO Forestry and Forest Products, 1998.
[4] 候宽昭, 陈焕镛, 钟观光, 等. 广州植物志[M]. 北京: 科学出版社, 1956. [5] GUIMARÃES L M D S, MIRANDA T, DOUGLAS L, et al. Eucalyptus pellita as a source of resistance to rust, ceratocystis wilt and leaf blight[J]. Crop Breeding and Applied Biotechnology, 2010, 10(2): 124-131. doi: 10.12702/1984-7033.v10n02a04
[6] ZANATA M, FREITAS M L M, SILVA M T, et al. Genetic parameters and gains with selection in open pollinated progeny test of Eucalypyus pellita, in Batatais-SP[J]. Revista do lnstituto Florestal, 2010, 22(2): 233-242. doi: 10.24278/2178-5031.2010222264
[7] THAVAMANIKUMAR S, ARNOLD R J, LUO J Z, et al. Genomic studies reveal substantial dominant effects and improved genomic predictions in an open-pollinated breeding population of Eucalyptus pellita[J]. G3-Genes Genomes Genetics, 2020, 10(10): 3751-3763. doi: 10.1534/g3.120.401601
[8] LAPAMMU M, WARBURTON P M, JAPARUDIN Y D, et al. Verification of tolerance to infection by ceratocystis manginecans in clones of Acacia mangium[J]. Journal of Tropical Forest Science, 2023, 35(Special issue): 42-50.
[9] 廖柏勇, 刘丽婷, 莫晓勇, 等. 10年生粗皮桉种源家系选择分析[J]. 华南农业大学学报, 2011, 32(4): 72-77. [10] 邓冬丽, 韦吉伟, 覃若飞, 等. 粗皮桉不同无性系的生长对比试验[J]. 桉树科技, 2023, 40(1): 21-25. [11] 门中华, 李生秀. 植物生物节律性研究进展[J]. 生物学杂志, 2009, 26(5): 53-55. [12] 宋淑媛, 顾宸瑞, 李春旭, 等. 应用曲线模型解析施肥对白桦苗期年高生长节律的影响[J]. 东北林业大学学报, 2021, 49(4): 17-23. doi: 10.3969/j.issn.1000-5382.2021.04.004 [13] 李峰卿, 陈焕伟, 周志春, 等. 红豆树优树种子和幼苗性状的变异分析及优良家系的初选[J]. 植物资源与环境学报, 2018, 27(2): 57-65. [14] ZHANG M M, LU N, JIANG L B, et al. Multiple dynamic models reveal the genetic architecture for growth in height of Catalpa bungei in the field[J]. Tree Physiology, 2022, 42(6): 1239-1255. doi: 10.1093/treephys/tpab171
[15] 麻文俊, 王军辉, 张守攻, 等. 楸树无性系苗期年生长参数的分析[J]. 东北林业大学学报, 2010, 38(1): 4-7. [16] 杨志玲, 杨旭, 谭梓峰, 等. 厚朴不同种源苗期生长模型的拟合[J]. 西北农林科技大学学报(自然科学版), 2011, 39(4): 60-68. [17] 秦莉, 赵有科, 黄荣凤, 等. 8年生粗皮桉生长应变及生长遗传变异[J]. 中南林业科技大学学报, 2008, 28(1): 58-63. [18] 刘晓华, 罗建中, 卢万鸿, 等. 两个连续世代粗皮桉生长与抗风能力遗传特征[J]. 分子植物育种, 2017, 15(12): 5103-5111. [19] 王楚彪, 杨艳, 白卫国, 等. 粗皮桉近红外光谱差异与其遗传差异间的关系[J]. 光谱学与光谱分析, 2021, 41(11): 3399-3404. [20] 王楚彪, 王建忠, 李华强, 等. 16个粗皮桉种源遗传特征分析及选优[J]. 热带亚热带植物学报, 2022, 30(1): 54-62. [21] 胡兴峰, 吴帆, 孙晓波, 等. 38年生马尾松种源生长及材性联合分析[J]. 南京林业大学学报(自然科学版), 2022, 46(3): 203-212. [22] 何霞, 邓成, 杨嘉麒, 等. 苦楝种源间生长性状的早期地理变异分析[J]. 北京林业大学学报, 2018, 40(7): 45-54. [23] 张沛健. 雷琼地区尾细桉纸浆林立地质量评价及生长规律研究[D]. 北京: 中国林业科学研究院, 2023. [24] 凡美玲, 方水元, 冯俊娇, 等. 4个竹种高生长模型的比较[J]. 竹子学报, 2018, 37(2): 64-70. [25] ZHENG Z, FANG K, CHEN Y, et al. Is the Pinus massoniana lamb. tree-ring latewood formation influenced by the diurnal temperature range in humid subtropical China[J]. Forests, 2022, 13(9): 1439. doi: 10.3390/f13091439
[26] 杨斌. 柳树苗期年生长模型的研究[J]. 西北林学院学报, 2006, 21(6): 97-99. [27] 万志兵, 冯刚, 朱成磊, 等. 不同柳树无性系一年生生长差异分析[J]. 分子植物育种, 2018, 16(7): 2358-2363. [28] 童洁, 石玉立. 加格达奇3种森林类型树高−胸径的曲线拟合[J]. 东北林业大学学报, 2017, 45(2): 6-11. [29] 和滢埝, 唐军荣, 李亚麒, 等. 氮磷添加对云南松苗木生长节律的影响[J]. 云南农业大学学报(自然科学), 2023, 38(3): 465-475. [30] 朱雅静, 王雪, 王丹等. 不同干型云南松子代苗木生长性状与异速生长分析[J]. 云南农业大学学报(自然科学), 2021, 36(6): 1044-1050. [31] XIAO L, FANG Y, ZHANG H, et al. Natural variation in the prolyl 4-hydroxylase gene PtoP4H9 contributes to perennial stem growth in Populus[J] The Plant Cell, 2023, 35(11): 4046-4065.
[32] 王玉虓, 史雯茜, 李倩中, 等. 元宝枫1年生苗木生长模型的拟合[J/OL]. 分子植物育种, (2023-12-06) [2024-04-26]. https://link.cnki.net/urlid/46.1068.S.20231205.1107.004. [33] 张力斌, 何明珠, 张珂. 柠条锦鸡儿生物量分配规律与异速生长对氮、磷添加的响应[J]. 生态学报, 2023, 43(16): 6627-6636. [34] 乔栋, 刘勇, 田书勇, 等. 不同土壤水势对毛白杨苗木生长节律和苗木质量的影响[J]. 北京林业大学学报, 2022, 44(4): 12-23. [35] 付威, 韦素云, 陈赢男. 植物生长发育动态QTL解析研究进展[J]. 生物技术通报, 2024, 40(2): 9-19. [36] 邬荣领, 王明庥, 黄敏仁, 等. 黑杨派新无性系研究: Ⅵ: 苗期年生长的动态分析[J]. 南京林业大学学报(自然科学版), 1988(4): 1-12. -
期刊类型引用(1)
1. 吴迪，栗云鹏，巴洪宇，杜鹃，吴惠明，张启龙，冯小宇，周德刚. 美洲型猪繁殖与呼吸综合征病毒实时荧光RT-PCR定量检测试剂盒测量不确定度研究. 中国动物检疫. 2020(05): 94-100 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 88
- HTML全文浏览量: 10
- PDF下载量: 39
- 被引次数: 3
