全州禾花鲤体质量和体长遗传参数估算及其选育进展

杜雪松, 陈忠, 文露婷, 潘贤辉, 黄姻, 周康奇, 林勇, 李哲, 罗辉, 覃俊奇

杜雪松, 陈忠, 文露婷, 潘贤辉, 黄姻, 周康奇, 林勇, 李哲, 罗辉, 覃俊奇. 全州禾花鲤体质量和体长遗传参数估算及其选育进展[J]. 南方水产科学, 2022, 18(6): 69-77. DOI: 10.12131/20210353
引用本文: 杜雪松, 陈忠, 文露婷, 潘贤辉, 黄姻, 周康奇, 林勇, 李哲, 罗辉, 覃俊奇. 全州禾花鲤体质量和体长遗传参数估算及其选育进展[J]. 南方水产科学, 2022, 18(6): 69-77. DOI: 10.12131/20210353
DU Xuesong, CHEN Zhong, WEN Luting, PAN Xianhui, HUANG Yin, ZHOU Kangqi, LIN Yong, LI Zhe, LUO Hui, QIN Junqi. Estimation of genetic parameters and its breeding progress for body mass and body length of Cyprinus carpio var. Quanzhounensis[J]. South China Fisheries Science, 2022, 18(6): 69-77. DOI: 10.12131/20210353
Citation: DU Xuesong, CHEN Zhong, WEN Luting, PAN Xianhui, HUANG Yin, ZHOU Kangqi, LIN Yong, LI Zhe, LUO Hui, QIN Junqi. Estimation of genetic parameters and its breeding progress for body mass and body length of Cyprinus carpio var. Quanzhounensis[J]. South China Fisheries Science, 2022, 18(6): 69-77. DOI: 10.12131/20210353

全州禾花鲤体质量和体长遗传参数估算及其选育进展

基金项目: 广西特色淡水鱼产业创新团队 (nycytxgxcxtd-21-08);国家自然科学基金项目 (31960730)
详细信息
    作者简介:

    杜雪松 (1988—),男,助理研究员,硕士,研究方向为鱼类育种。E-mail: 627380756@qq.com

    通讯作者:

    覃俊奇 (1980—),男,高级工程师,研究方向为水产养殖。E-mail: 28434483@qq.com

  • 中图分类号: S 917.4

Estimation of genetic parameters and its breeding progress for body mass and body length of Cyprinus carpio var. Quanzhounensis

  • 摘要: 为提高全州禾花鲤 (Cyprinus carpio var. Quanzhounensis) 养殖群体的生长速度,采用家系选育方法对其开展了连续2世代的选育,测量了G0、G1和G2代121个家系共3 699尾个体的体质量和体长指标,利用混合线性模型估计体质量和体长性状的遗传力、表型相关和遗传相关,并通过基于目标性状育种值和最小二乘值估计的方法评估了G0—G2代的选择进展。结果显示,禾花鲤体质量变异系数为0.50~0.72,体长性状变异系数为0.18~0.22。体长和体质量的遗传力估计值分别为0.132 (P<0.01) 和0.122 (P<0.01),遗传相关和表型相关分别为0.921 (P<0.01)和0.995 (P<0.01)。采用目标性状最小二乘值和目标性状育种值估计的体质量性状选择反应为12.82%和15.45%,体长性状选择反应为2.76%和6.60%。表明禾花鲤体质量均匀度较低,不利于养殖生产,但具有较高的选育改良潜力;该群体在完成2个世代的家系选育后获得了有效遗传进展,但总体进展不显著。在今后选育过程中可采用分子辅助选育 (Molecular-assisted selection MAS)、大规模混合家系选育等选育策略,以提高选育效率。
    Abstract: In order to improve the growth rate of the cultured population of Cyprinus carpio var. Quanzhounensis in Quanzhou, two generations of breeding were carried out by family breeding method. The body mass and body length indexes of 3 699 individuals in 121 families of G0, G1 and G2 generations of C. carpio were measured. The heritability, phenotypic correlation and genetic correlation of body mass and body length traits were estimated by mixed linear model. The selection progress of G0−G2 generations was evaluated by the method of breeding value and least square value estimation based on target traits. The results show that the coefficient of variation of body mass and body length was 0.50−0.72 and 0.18−0.22 respectively. The heritability estimates were 0.132 (P<0.01) and 0.122, respectively (P<0.01). The genetic correlation and phenotypic correlation were 0.921 (P<0.01) and 0.995 (P<0.01), respectively. The value of selective reaction in body mass estimated by breeding value based on target traits and least square estimation was 12.82% and 15.45%, respectively, while 2.76% and 6.60% in body length, respectively, which indicates that the evenness of body mass was low, not conductive to breeding production, but had high breeding and improvement potential. In general, this population ofC. carpio var. Quanzhounensis obtained effective genetic progress after two successive generations by the family breeding, but with a modest progress. In is suggested that molecular-assisted selection (MAS) and marker-assisted large-scale mixed families' selection should be adopted in the future breeding so as to improve the breeding efficiency.
  • 鳗败血假单胞菌 (Pseudomonas anguilliseptica, PA) 于1972年由 Wakabayashi和Egusa[1]在日本鳗鲡 (Anguilla japonica) 上首次分离并报道,其后发现该菌主要影响欧洲地区的欧洲鳗鲡 (A. anguilla)[2]、大西洋鲷 (Sparus aurata)[3]、大西洋鳕 (Gadus morhua)[4]、圆鳍鱼 (Cyclopterus lumpus)[5]等海洋鱼类。因其感染的流行季节主要在水温低于12 ℃的冬季,且被感染宿主呈现体表出血性红斑的特征,因此被称为“冬季病”[6]和“红斑病”[7]。近年来,笔者研究发现PA是大口黑鲈 (Micropterus salmoides) 的一种新病原菌 (另文报道),对大口黑鲈的健康养殖造成威胁,但目前其在我国大口黑鲈中的流行情况尚不完全清楚。

    PA的生长条件苛刻、培养时间长[8],同时其感染的临床症状与嗜水气单胞菌 (Aeromonas hydrophila)[9]、维氏气单胞菌 (A. veronii)[10]等相似,通过病原菌分离鉴定及临床症状的传统方法难以进行快速诊断,往往出现漏检、错检、误诊的情况,影响了PA感染的有效防控。建立快速、特异的检测方法不仅对PA的早期临床诊断和防控具有重要意义,还能为开展流行病学调查提供条件。目前针对水生动物细菌性病原的检测方法主要有环介导等温扩增技术 (LAMP)[11]、核酸序列依赖性扩增检测技术 (NASBA)[12]、免疫磁珠分离技术 (IMS)[13]、酶联免疫吸附测定技术 (ELISA)[14]、重组酶介导等温扩增技术 (RAA)[15]、实时荧光定量PCR (qPCR)[16]等,综合考虑经济性、操作简便性、时效性等因素后,本研究采用SYBR Green I 实时荧光定量PCR (SYBR Green I real-time quantitative PCR) 和重组酶介导等温扩增技术结合侧流层析试纸条(Recombinase-mediated isothermal amplification combined with lateral flow dipstick, RAA-LFD),以PA的管家基因recA为靶基因建立2种检测方法,以期为PA的早期、快速检测提供方法参考。

    健康大口黑鲈 [外部无异常、PA特异性检测结果为阴性,体质量为 (18.4±1.5) g)]购自四川某苗种生产企业;各菌株 (表1) 由本实验室保存。主要试剂为:pUC18-T载体 (宝生物工程有限公司,大连);细菌基因组DNA提取试剂盒、胶回收试剂盒及质粒提取试剂盒 (天根生化科技有限公司,北京);Taq Pro Universal SYBR qPCR Master Mix、2×Rapid Taq Master Mix (诺唯赞生物科技有限公司,南京);RAA核酸扩增试剂盒及侧流层析试纸条 (乐尚生物科技有限公司,无锡)。

    表  1  研究所用的细菌菌株
    Table  1  Bacterial strains used in this study
    菌株名称
    Bacterial strain
    菌株编号
    Strain No.
    鳗败血假单胞菌 Pseudomonas anguilliseptica SICAU PA001
    无乳链球菌 Streptococcus agalactiae GY104
    拟态弧菌 Vibrio mimicus SCCF01
    嗜水气单胞菌 Aeromonas hydrophila YYL
    鮰爱德华菌 Edwardsiella ictaluri CNHYC01
    鰤诺卡菌 Nocardia seriolae DY01
    荧光假单胞菌 P. fluorescens ATCC®13525
    柱状黄杆菌 Flavobacterium cloumnare ATCC®23463
    维氏气单胞菌 A. veronii ATCC®35624
    下载: 导出CSV 
    | 显示表格

    根据Mjølnerød等[17]对53株PA管家基因的分析结果,选择其中具有完全同源性的recA基因 (登录号:MW685241.1—MW685293.1) 作为建立2种检测方法的靶基因,并采用验证recA基因的特异性引物 (recA-PCR-F/R) 作为本研究PCR和构建标准品质粒的引物。使用Primer Premier 5.0软件,根据SYBR Green I real-time quantitative PCR和RAA引物设计原则设计多对引物,并设计RAA探针。成功筛选出引物recA-qPCR-F/R、recA-RAA-F/R和探针recA-RAA-Probe,并在recA-RAA-R的5' 端标记生物素。本研究所用引物和探针 (表2) 均由生工生物工程 (上海) 有限公司合成。将recA-RAA-F/R配合recA-RAA-Probe的RAA扩增产物和recA-qPCR-F/R的PCR 扩增产物送至生工生物工程 (上海) 有限公司测序,在GenBank中获取相应序列并与测序结果进行比对,以验证引物的特异性。

    表  2  研究所用引物及探针
    Table  2  Primers used in this study
    检测方法
    Detection method
    引物名称
    Primer name
    序列 (5'—3')
    Sequence (5'−3')
    产物大小
    Length of product/bp
    PCRrecA-PCR-F/RF: CTCGCTTGGTCTGGACATCG
    R: CCTTGCCTTGGCCGATCTT
    740
    qPCRrecA-qPCR-F/R F: GGCAGCTCTTACACCCAAAG
    R: GGCTCCCAAACATAACACCTAT
    179
    RAA-LFDrecA-RAA-F/RF: CATGCACTTGACCCAATCTACGCAGCAAAG
    R: Biotin CCCGTAATCTTACGCATTGCCTGACTCATC
    251
    recA-RAA-ProbeFAM AAGCTGAAATCGAAGGTGAAATTGGGGATA/idSp/
    GCACGTTGGGTTGGC C3 Spacer
    下载: 导出CSV 
    | 显示表格

    以PA的基因组DNA为模板,用recA-PCR-F/R引物扩增recA基因目的片段,与pUC18-T载体连接后转化于DH5α感受态细胞,获得重组标准品质粒pUC18-recA。用微量核酸蛋白检测仪测定质粒浓度和纯度,并根据摩尔定律计算出单位体积质粒所含的DNA拷贝数。

    2种方法均采用2.816×106 拷贝·μL−1的标准品质粒进行优化。qPCR引物终浓度梯度设置为0.2~1.0 μmol·L−1,退火温度梯度为 58~64 ℃,以获得qPCR最小的平均Ct值和标准差;RAA引物终浓度梯度设置为 6.25~200 μmol·L−1,RAA扩增反应温度梯度为29~54 ℃,并通过观察不同条件下侧流层析试纸条第30分钟时的显色情况确定RAA扩增的最佳反应体系和反应条件。

    qPCR方法的建立:将标准品质粒10倍稀释成9组浓度梯度 (2.816×109~2.816×101 拷贝·μL−1),采用优化后的反应条件 (反应体系总体积10 μL:recA-qPCR-F/R各0.2 μL,2×Taq Pro Universal SYBR qPCR Mix 5 μL,ddH2O 3.6 μL,模板DNA 1 μL;反应条件为:95 ℃预变性30 s;95 ℃变性10 s,60 ℃退火20 s,72 ℃延伸10 s,40个循环) 进行qPCR检测,建立以质粒拷贝数对数值为x轴和以Ct值为y轴的qPCR标准曲线;每个浓度梯度的标准品质粒设置3次qPCR检测,以进行组内重复性分析,并以相同的标准品质粒为模板,在不同时间点开展3次独立重复实验,进行组间重复性分析;通过Ct值变异系数评价方法的重复性。

    RAA-LFD方法的建立:采用优化后的反应条件 (反应体系总体积50 μL:recA-RAA-Probe 0.6 μL,recA-RAA-LF/LR各2.1 μL,基础缓冲液25 μL,模板DNA 2 μL,ddH2O 15.2 μL;各组分混匀后,加入含冻干酶粉的反应管中溶解,然后加入280 mmol·L−1的醋酸镁溶液3 μL作为启动剂混匀,39 ℃水浴20 min,立即置于冰上5 min) 对标准品质粒 (2.816×109 拷贝·μL−1) 进行RAA扩增,将侧流层析试纸条样品垫末端垂直插入含有50 μL稀释5倍的RAA核酸扩增产物的EP管中,观察试纸条C线和T线在30 min内的颜色变化,并参考Saxena等[18]的标准判定结果。

    各取101~108 拷贝·μL−1的标准品质粒1 μL作为模板,分别进行qPCR、RAA和普通PCR检测,分析3种方法的灵敏性。

    取质量浓度为10 ng·μL−1的各株细菌DNA 1 μL作为模板,同时进行qPCR、RAA和普通PCR扩增,分析3种方法的特异性。

    通过建立的qPCR和RAA-LFD方法以及普通PCR方法,对10份疑似感染PA的大口黑鲈肾组织样本与30份腹腔注射感染PA的大口黑鲈肾组织样本进行检测并比较检出率的差异。

    选取15尾健康大口黑鲈 (共3组,每组5尾),腹腔注射PA,于感染后的第7天分别采集每组病鱼的肝、脾、肾、心、肌肉、脑、眼和鳃组织并提取DNA,采用qPCR方法和标准曲线方程测定每组病鱼组织中的平均PA载量。

    PA提取的DNA通过qPCR引物recA-qPCR-F/R扩增后得到179 bp的单一目的条带 (图1-a),通过RAA引物recA-RAA-F/R配合recA-RAA-Probe扩增后得到251 bp的单一目的条带 (图1-b)。将目的条带的测序结果与GenBank中获取的相应序列进行比对,结果显示2段产物序列与其他序列的同源性均为100% (图2)。

    图  1  基于鳗败血假单胞菌recA基因的引物特异性与标准品质粒鉴定结果
    注:a. qPCR引物特异性结果;b. RAA引物及探针特异性结果;c. 标准品质粒鉴定结果。
    Fig. 1  Results based on primers specificity and standard plasmid identification of recA gene of P. anguilliseptica
    Note: a. Primer specificity results for qPCR; b. Primer and probe specificity results for RAA; c. Results of standard plasmid identification.

    标准品质粒扩增获得的目的片段大小为740 bp (图1-c),与预期大小相符。重组质粒pUC18-recA的质量浓度为1.057×102 ng·μL−1,经质粒拷贝数计算公式换算后为2.816×1010 拷贝·μL−1

    qPCR扩增条件的优化结果显示,qPCR引物浓度为0.6 μmol·L−1、退火温度为60 ℃时平均Ct值 (13.259) 和标准差 (0.045) 最小 (表3),故选择该引物浓度和退火温度为最佳反应条件。

    表  3  鳗败血假单胞菌qPCR扩增条件的优化
    Table  3  Optimization of qPCR amplification conditions for P. anguilliseptica
    退火温度
    Annealing temperature/℃
    引物浓度
    Primer concentration/(μmol·L−1)
    Ct值
    Ct value
    平均Ct值
    Mean Ct value
    标准差
    Standard deviation
    58 0.2 14.583 15.418 14.845 14.949 0.349
    0.4 13.876 14.356 14.230 14.154 0.203
    0.6 13.336 13.576 13.429 13.447 0.099
    0.8 13.745 13.401 13.313 13.486 0.187
    1.0 13.756 13.444 13.895 13.698 0.189
    60 0.2 14.419 15.074 14.355 14.616 0.325
    0.4 14.079 13.832 13.786 13.899 0.129
    0.6 13.204 13.259 13.314 13.259 0.045
    0.8 13.349 13.608 13.573 13.510 0.115
    1.0 13.749 14.023 13.988 13.920 0.122
    62 0.2 16.812 16.033 16.267 16.370 0.326
    0.4 15.305 15.466 15.077 15.283 0.160
    0.6 14.275 14.119 14.439 14.277 0.131
    0.8 14.206 14.401 15.046 14.551 0.359
    1.0 14.748 14.207 14.805 14.587 0.269
    64 0.2 15.745 16.784 15.695 16.075 0.502
    0.4 14.139 14.916 14.878 14.644 0.358
    0.6 15.061 15.328 15.457 15.282 0.165
    0.8 15.740 15.112 15.557 15.470 0.264
    1.0 15.963 15.369 15.844 15.725 0.257
    下载: 导出CSV 
    | 显示表格

    RAA扩增条件的优化结果如图3所示。试纸条显色第30分钟时的T线灰度值与引物浓度成正比 (图3-a—3-b),由于引物浓度过高可能会产生大量引物二聚体,且在25~100 μmol·L−1时T线灰度值差异不显著 (p>0.05),故选择50 μmol·L−1为最佳引物浓度。T线灰度值在反应温度介于29~54 ℃时呈正态分布,在39~44 ℃时达到最高,差异不显著 (p>0.05,图3-c—3-d),故选择其中获得最高T线灰度值的39 ℃为最适反应温度。

    图  2  鳗败血假单胞菌qPCR及RAA扩增序列同源性对比
    Fig. 2  Comparison of homology between qPCR and RAA amplified sequences of P. anguilliseptica
    图  3  鳗败血假单胞菌RAA扩增条件的优化
    注:a. 不同引物浓度的RAA-LFD显色结果;b. 灰度分析结果 (对应3-a);c. 不同反应温度的RAA-LFD显色结果;d. 灰度分析结果 (对应3-c);不同字母表示差异显著 (p<0.05)。
    Fig. 3  Optimization of RAA amplification conditions for P. anguilliseptica
    Note: a. Results of RAA-LFD color development with different primer concentrations; b. Grayscale analysis results (Corresponding to 3-a); c. Results of RAA-LFD color development at different reaction temperatures; d. Grayscale analysis results (Corresponding to 3-c); different letters represent significant differences (p<0.05).

    采用优化后的反应条件对9组10倍梯度稀释的标准品质粒 (2.816×109~2.816×101 拷贝·μL−1) 进行qPCR检测,建立质粒拷贝数对数值与Ct值对应关系的定量标准曲线 (图4-a),标准曲线介于2.816×107~2.816×102 拷贝·μL−1之间Ct值线性良好。标准曲线方程为y= –3.158x+32.213,r2=0.999 2。熔解曲线显示,6个浓度的标准品质粒在熔解温度为 (84±0.5) ℃时均出现特异性细窄单峰 (图4-b),表明扩增产物具有良好的特异性和均一性。标准曲线重复性检测结果表明,组内变异系数≤0.730%,组间变异系数≤1.443% (表4),该方法重复性和稳定性较好。

    图  4  鳗败血假单胞菌qPCR方法建立结果
    注:a. qPCR标准曲线;b. qPCR熔解曲线,1—6. 质粒标准品浓度为2.816×107~2.816×102 拷贝·μL−1
    Fig. 4  Results established by qPCR for P. anguilliseptica
    Note: a. Standard curve for qPCR; b. Melting curves for qPCR, 1−6. Concentration of standard plasmid was 2.816×107−2.816×102 copies·μL−1.

    采用优化后的反应条件对标准品质粒 (2.816×109 拷贝·μL−1) 进行RAA扩增,并通过侧流层析试纸条对扩增产物进行检测。结果显示,T线灰度值与显色时长成正比,在15 min后试纸条T线显色情况较好 (T线灰度值≥40.0),且各时间点无显著性差异 (p>0.05,图5-a—5-b),表明该方法最快15 min可读取检测结果,且显色较为稳定。

    表  4  鳗败血假单胞菌qPCR方法重复性试验结果
    Table  4  Results of qPCR reproducibility test for P. anguilliseptica
    标准品质粒浓度
    Standard quality plasmid
    concentration/
    (拷贝·μL−1)
    组内重复性试验
    Repeatability test within group
    组间重复性试验
    Repeatability test between group
    平均Ct值±标准差
    Mean Ct value±standard deviation
    变异系数
    Variable coefficient/%
    平均Ct值±标准差
    Mean Ct value±standard deviation
    变异系数
    Variable coefficient/%
    2.816×107 10.300±0.032 0.311 10.328±0.149 1.443
    2.816×106 13.213±0.042 0.318 13.257±0.122 0.920
    2.816×105 16.710±0.122 0.730 16.606±0.212 1.277
    2.816×104 19.511±0.102 0.523 19.523±0.111 0.569
    2.816×103 22.609±0.081 0.358 22.618±0.168 0.743
    2.816×102 26.083±0.130 0.498 26.120±0.210 0.804
    下载: 导出CSV 
    | 显示表格
    图  5  鳗败血假单胞菌RAA-LFD反应结果
    注:a. 不同反应时间的RAA-LFD显色;b. 灰度分析 (对应5-a),不同字母表示差异显著 (p<0.05)。
    Fig. 5  Results of RAA-LFD reaction of P. anguilliseptica
    Note: a. Results of RAA-LFD color development at different reaction time; b. Grayscale analysis results (Corresponding to 5-a), different letters represent significant differences (p<0.05).

    由普通PCR灵敏性试验结果可知,模板浓度低于2.816×105 拷贝·μL−1时,无法产生扩增条带 (图6-a),故普通PCR的最低检测浓度为2.816×105 拷贝·μL−1。由2.816×107 拷贝·μL−1浓度起始的6个10倍梯度稀释的标准品质粒扩增曲线可知,模板浓度在低于2.816×102 拷贝·μL−1时,无法产生扩增曲线 (图6-b),故qPCR的最低检测浓度为2.816×102 拷贝·μL−1。由RAA-LFD灵敏性试验结果可知,模板浓度低于2.816×104 拷贝·μL−1时扩增产物在试纸条上显色极不明显 (T线灰度值≤5.0,图6-c—6-d),故RAA-LFD的最低检测浓度为2.816×104 拷贝·μL−1。综上,构建的2种方法灵敏性均高于普通PCR,且qPCR方法的灵敏性最高。

    图  6  鳗败血假单胞菌灵敏性试验结果
    注:a. 普通PCR灵敏性电泳图,M. 标准分子量DL 2000,NC. 阴性对照,1. 2.816×108 拷贝·μL−1,2. 2.816×107 拷贝·μL−1,3. 2.816×106拷贝·μL−1,4. 2.816×105 拷贝·μL−1,5. 2.816×104 拷贝·μL−1;b. qPCR标准质粒扩增曲线,1—6. 质粒标准品浓度为2.816×107~2.816×102 拷贝·μL−1;c. 不同浓度标准质粒的RAA-LFD显色;d. 灰度分析 (对应6-c),不同字母表示差异显著 (p<0.05)。
    Fig. 6  Results of sensitivity test for P. anguilliseptica
    Note: a. Sensitive electrophoresis pattern of PCR, M. Standard molecular weight is DL 2000, NC. Negative control, 1. 2.816×108 copies·μL−1, 2. 2.816×107 copies·μL−1, 3. 2.816×106copies·μL−1, 4. 2.816×105 copies·μL−1, 5. 2.816×104 copies·μL−1; b. Amplification curves of qPCR standard plasmids, 1−6. Concentration of standard plasmid was 2.816×107−2.816×102 copies·μL−1; c. Results of RAA-LFD color development with different concentrations of standard plasmids; d. Grayscale analysis results (Corresponding to 6-c), different letters represent significant differences (p<0.05).

    qPCR特异性扩增曲线结果显示,PA的DNA扩增结果为阳性,而其他8株菌的结果为阴性 (图7-a)。RAA-LFD特异性试验结果显示,PA的DNA扩增产物显色较为明显 (T线灰度值为27.0),而其他8株菌显色极不明显 (T线灰度值≤3.0,图7-b—7-c)。普通PCR特异性试验结果显示,PA的DNA扩增后均出现单一条带(740 bp),而其他8株菌DNA均未出现扩增条带 (图7-d)。综上,3种检测方法均出现明显的特异性扩增结果。

    图  7  鳗败血假单胞菌特异性试验结果
    注:a. qPCR特异性扩增曲线;b. 不同菌株的RAA-LFD显色结果;c. 灰度分析结果 (对应7-b),不同字母表示差异显著 (p<0.05);d. 普通PCR特异性电泳图。
    Fig. 7  Results of specificity test for P. anguilliseptica
    Note: a. Specific amplification curves of qPCR; b. RAA-LFD chromogenic results of different strains; c. Grayscale analysis results (Corresponding to 7-b), different letters represent significant differences (p<0.05); d. Specific electrophoresis pattern of PCR.

    通过建立的2种方法及普通PCR方法对采集的10份疑似感染PA的大口黑鲈肾组织样本与30份腹腔注射感染PA的大口黑鲈肾组织样本进行检测。普通PCR方法的阳性样本检出率为62.50%,而qPCR和RAA-LFD方法的阳性样本检出率分别为87.50%和85.00% (表5),较普通PCR方法有明显提高。且qPCR与RAA-LFD方法具有高一致性,较普通PCR方法更可靠。

    表  5  鳗败血假单胞菌3种检测方法的应用结果对比
    Table  5  Comparison of application results of three detection methods for P. anguilliseptica
    检测方法
    Detection method
    样本总数
    Total sample number
    阳性样本检出率
    Detection rate of positive samples
    共同阳性/阴性样本数
    Number of common positive/negative samples
    Kappa值
    Kappa value
    qPCR/RAA-LFD 40/40 87.50%/85.00% 34/5 0.895
    qPCR/PCR 40/40 87.50%/62.50% 25/5 0.385
    RAA-LFD/PCR 40/40 85.00%/62.50% 25/6 0.455
    注:Kappa值的取值范围为0~1.00,其中1.00表示完全一致,0表示完全不一致。常见标准为:Kappa值在0~0.20表示一致性极低,0.21~0.40表示一致性较低,0.41~0.60表示一致性中等,0.61~0.80表示一致性较高,0.81~1.00表示一致性极高。 Note: The Kappa values ranged from 0 to 1.00, with 1.00 indicating complete agreement and 0 indicating complete disagreement. Common criteria are as follows: Kappa of 0−0.20 represents very low consistency; 0.21−0.40 represents low consistency; 0.41−0.60 represents medium consistency; 0.61−0.80 represents high consistency; 0.81−1.00 represents extremely high consistency.
    下载: 导出CSV 
    | 显示表格

    使用qPCR方法测得PA载量 (图8),PA在大口黑鲈肾中载量最高,达3.533×107 拷贝·ng−1,在PA总载量中占比达70.84%;在肝、脾、眼和脑等组织器官中的载量均较高,但在鳃中载量极低。

    图  8  受感染大口黑鲈各部位的鳗败血假单胞菌载量结果
    注:组1、组2、组3表示3个平行的PA感染组,感染量为1×107 CFU·尾−1
    Fig. 8  Results of P. anguilliseptica load in each part of infected M. salmoides
    Note: Group 1, 2 and 3 represent three parallel PA infected groups at 1×107 CFU·ind−1.

    PA曾被认为是海水鱼类的重要病原菌,近年来其感染谱有向淡水鱼类扩展的趋势。Ajengaekanurmaningdyah和Kurniasih[19]于2018年首次报道了PA感染淡水鱼类金鱼 (Cyprinus carpio) 引起死亡;本实验室近年发现PA已成为危害大口黑鲈的新细菌性病原,给大口黑鲈养殖产业造成严重威胁 (另文报道)。PA感染后病鱼处于无明显症状的亚临床期较长,发展到明显症状的临床期往往导致高死亡率,由于PA仅在TSA和BHI等营养丰富的平板上生长[8],且出现菌落的时间超过1周,采用传统病原菌分离鉴定及临床症状进行诊断的难度大、效率低,建立快速、特异、灵敏的PA检测方法在该病的防控中就显得尤为重要。目前,国内外对于PA检测方法的研究较少,除了传统的分离鉴定外,有采用普通PCR进行检测的报道[20-21],但在笔者实验室的应用中发现普通PCR存在假阴性率高的缺点。

    qPCR方法具有重复性好、特异性强和灵敏性高等优点[22]。相较于TaqMan探针法,SYBR Green I法只需在反应混合液中加入引物和待测样品的DNA,无需设计合成探针,简单高效且成本较低[23],目前已应用于鼠伤寒沙门氏菌 (Salmonella typhimurium)[24]、金黄色葡萄球菌 (Staphylococcus aureus)[25]等病原菌的检测。RAA方法对比NASBA、LAMP等温扩增方法有扩增温度较低、仅需1对引物和探针、反应时间较短的优势[26-27],已在铜绿假单胞菌 (Pseudomonas aeruginosa)[28]、鰤诺卡菌 (Nocardia seriolae)[29]等病原菌的检测上得到了应用。RAA-LFD 方法则更利于脱离实验室环境,应用到现场检测中,目前该技术已在副结核分枝杆菌 (Mycobacterium avium subsp. paratuberculosis)[30]、霍乱弧菌 (Vibrio cholerae)[31]、停乳链球菌 (Streptococcus dysgalactiae)[32]等病原菌上开展了研究,为这些病原菌的现场快速检测提供了重要的支持。本研究基于PA管家基因recA成功建立了SYBR Green I real-time quantitative PCR和RAA-LFD 2种快速、特异的检测方法,检出率相较PCR方法提高了至少22.50%。qPCR方法的最低检测浓度为2.816×102 拷贝·μL−1,RAA-LFD方法的最低检测浓度为2.816×104 拷贝·μL−1,尽管建立的qPCR方法的灵敏性高于RAA-LFD方法,但2种方法检测结果的一致性较高,且RAA-LFD方法更具有开展现场检测的优势。这2种方法的建立为PA的早期快速检测与流行病学调查提供了重要的方法依据。

    据报道,PA感染金鱼致肾、肝、脾、肌肉、心、脑等组织器官病变,且肾的病变最为严重,表现为严重坏死[19];PA感染比目鱼 (Psetta maxima) 可致肾、肝、脾、肌肉、脑、眼等组织器官出现病变[33]。由此可见PA感染宿主致全身多组织器官损伤,但对于这些组织器官出现的病变与病原菌的分布及载量之间的关系尚不清楚。本研究建立的SYBR Green I real-time quantitative PCR方法可对PA 进行定量检测,使用该方法对PA感染大口黑鲈后病原菌在各组织的分布与载量进行分析,结果发现PA感染后菌体在肾、肝、脾、眼、脑、心等全身多组织器官分布,其中在肾中的载量最高,其感染致宿主组织器官的病理损伤程度与病原菌的载量有直接关系,表明肾可能是PA感染最重要的侵袭靶器官。

  • 图  1   体质量和体长性状的数据分布

    Figure  1.   Normal distribution of body mass and body length

    图  2   体质量和体长的幂相关图

    Figure  2.   Power correlation plot of body mass and body length

    表  1   各世代家系构建和培育

    Table  1   Family construction and breeding for each generation

    世代
    Generation
    群体
    Population
    家系数
    Number of
    families
    个体总数
    Number of
    individuals
    各家系平均个体数
    Number of individuals
    per family
    养殖池塘
    Aquaculture
    pond
    养殖密度
    Stocking density/
    (尾·m−2)
    养殖天数
    Culturing
    days/d
    G0 选育群 31 907 29.30 P5# 2.3 168
    G1 选育群 37 1 087 29.40 P8#、P12# 2.1 (P8#) 152
    对照群 6 216 36.00 2.3 (P12#) 150
    G2 选育群 41 1 251 30.50 P4#、P7# 2.2 (P4#) 140
    对照群 6 238 39.70 2.0 (P7#) 142
    合计 Total 121 3 699 32.98
    下载: 导出CSV

    表  2   体质量和体长的描述性统计量

    Table  2   Descriptive statistics of body mass and body length

    世代
    Generation
    群体
    Population
    性别
    Sex
    样本数
    Number of samples/尾
    体质量 Body mass 体长 Body length
    平均值
    Mean/g
    标准差
    SD
    变异系数
    CV
    平均值
    Mean/mm
    标准差
    SD
    变异系数
    CV
    G0 选育群 全部 906 37.44 27.01 0.72 97.79 20.65 0.21
    449 40.88 21.35 0.52 100.52 22.04 0.22
    457 33.94 24.60 0.72 95.01 18.74 0.20
    G1 选育群 全部 1087 39.60 25.94 0.66 99.81 21.47 0.22
    534 40.02 26.99 0.67 99.82 21.88 0.22
    553 39.20 24.90 0.64 99.80 21.08 0.21
    对照群 全部 216 33.59 17.15 0.51 95.23 17.62 0.19
    99 33.76 17.56 0.52 95.74 17.68 0.18
    117 33.44 16.75 0.50 94.80 17.63 0.19
    G2 选育群 全部 1281 39.61 27.13 0.68 99.43 21.72 0.22
    653 39.71 27.49 0.69 99.58 22.16 0.22
    628 39.52 26.77 0.68 99.27 21.30 0.21
    对照群 全部 209 34.54 20.60 0.60 96.03 20.09 0.21
    116 35.06 20.36 0.60 95.50 20.38 0.21
    93 34.12 20.98 0.60 96.68 19.81 0.20
    下载: 导出CSV

    表  3   体质量与体长性状的正态性检验

    Table  3   Normality test of body mass and body length

    性状
    Trait
    世代
    Generation
    W检验 W test 偏度
    Skewness
    峰度
    Kurtosis
    Statistic df Sig.
    体质量
    Body mass/g
    G0 0.107 906 0.000 2.417±0.081 13.44±0.162
    G1 0.113 1303 0.000 1.869±0.680 5.589±0.135
    G2 0.124 1490 0.000 2.297±0.063 9.439±0.127
    G0—G2 0.110 3699 0.000 2.186±0.040 9.089±0.080
    体长
    Body length/mm
    G0 0.047 906 0.000 0.388±0.081 0.300±0.162
    G1 0.035 1303 0.001 0.529±0.068 0.355±0.135
    G2 0.040 1490 0.000 0.581±0.063 0.378±0.127
    G0—G2 0.037 3699 0.000 0.520±0.040 0.361±0.080
    下载: 导出CSV

    表  4   方差分量和遗传力

    Table  4   Variance components and heritability

    性状
    Trait
    世代
    Generation
    加性遗传方差
    σa2
    环境方差
    σc2
    残差方差
    σe2
    表型方差
    σp2
    遗传力
    h2±SE
    体质量
    Body mass/g
    G0 59.12 11.4 539.92 610.44 0.097±0.047**
    G1 70.11 2.89 522.33 595.33 0.118±0.066*
    G2 48.55 6.88 611.46 666.89 0.074±0.071
    G0—G2 82.35 17.15 540.11 639.61 0.132±0.037**
    体长
    Body length/mm
    G0 70.21 23.15 328.57 421.93 0.176±0.325
    G1 37.13 35.56 367.62 440.31 0.091±0.069
    G2 37.52 15.56 341.33 394.41 0.099±0.068
    G0—G2 47.93 23.67 346.51 418.11 0.122±0.043**
    注:**. 极显著性水平(P<0.01);*. 显著性水平(P<0.05);下表同此。 Note: **. Extremely significant difference (P<0.01); *. Significance difference (P<0.05); the same case in the following table.
    下载: 导出CSV

    表  5   遗传相关和表型相关

    Table  5   Genetic correlation and phenotypic correlation

    世代
    Generation
    遗传相关
    Genetic correlation
    表型相关
    Phenotypic correlation
    G0 0.764±0.00** 0.875±0.039**
    G1 0.781±0.00** 0.837±0.056**
    G2 0.757±0.00** 0.794±0.067**
    G0—G2 0.921±0.002** 0.995±0.004**
    下载: 导出CSV

    表  6   体质量和体长性状的遗传进展估计

    Table  6   Estimation of predicted genetic gain for body mass and body length

    性状
    Trait
    世代
    Generation
    群体
    Population
    方法A
    Method A
    方法B
    Method B
    最小二乘值
    Least squares
    mean
    选择反应
    Selection response
    平均育种值
    Average
    breeding
    value
    选择反应
    Selection response
    遗传获得
    Genetic gain
    百分比
    Percentage/%
    遗传获得
    Genetic gain
    百分比
    Percentage/%
    体质量
    Body mass/g
    G0 选育群 39.95 −0.423
    G1 选育群 41.57 3.23 8.42 1.467 1.89 4.92
    对照群 38.34
    G2 选育群 41.21 2.51 6.48 4.382 2.91 7.52
    对照群 38.70
    合计 Total 5.74 15.45 4.80 12.82
    体长
    Body length/
    mm
    G0 选育群 99.6 −0.329
    G1 选育群 101.7 4.2 4.27 0.723 1.05 1.07
    对照群 97.5
    G2 选育群 100.7 2.2 2.23 2.362 1.64 1.66
    对照群 98.5
    合计 Total 6.4 6.60 2.69 2.76
    下载: 导出CSV
  • [1] 史君洁. 鲤种质分子鉴定技术构建[D]. 上海: 上海海洋大学, 2017: 5-6.
    [2]

    REN W, HU L, GUO L, et al. Preservation of the genetic diversity of a local common carp in the agricultural heritage rice-fish system[J]. PNAS, 2018: 546-554.

    [3] 支宇, 刘其根, 吴嘉敏, 等. 人工诱导青田田鱼雌核发育子代形态和遗传分析[J]. 上海海洋大学学报, 2022, 31(4): 839-848.
    [4]

    PENG J X, ZENG D G, HE P P, et al. mRNA and microRNA transcriptomics analyses in intermuscular bones of two carp species, rice flower carp (Cyprinus carpio var. Quanzhounensis) and Jian carp (Cyprinus carpio var. Jian)[J]. Comp Biochem Physiol D, 2019, 30: 71-80.

    [5] 鲁翠云, 杜雪松, 郑先虎, 等. 金边鲤群体的遗传结构及个体间遗传差异的QTL标记分析[J]. 淡水渔业, 2021, 51(4): 58-64. doi: 10.3969/j.issn.1000-6907.2021.04.008
    [6] 李妹娟, 钟旭华, 梁开明, 等. 广东省稻鱼共生生态种养发展现状与对策建议[J]. 广东农业科学, 2021, 48(10): 111-120.
    [7]

    CHRISTOS P, MARTIN K, MARTIN P, et al. Accuracy of genomic evaluations of juvenile growth rate in common carp (Cyprinus carpio) using genotyping by sequencing[J]. Front Genet, 2018, 9: 82. doi: 10.3389/fgene.2018.00082

    [8] 赵新春. 微卫星标记对鲤的两个育成品种遗传背景分析和遗传评价体系建立[D]. 上海: 上海海洋大学, 2015: 6-17.
    [9] 樊佳佳, 马冬梅, 朱华平, 等. F5  代禾花鲤形态性状对体质量的影响[J]. 广东农业科学, 2021, 48(8): 124-130.
    [10] 禾花鲤“乳源1号”[J]. 中国水产, 2021(10): 96-101.
    [11] 金万昆. 淡水养殖鱼类种质资源库[M]. 北京: 中国农业科学技术出版社, 2011: 42-45.
    [12] 吴志强, 李坚明. 推进稻渔生态种养, 促进广西农业绿色发展[J]. 广西农学报, 2019, 34(2): 1-5. doi: 10.3969/j.issn.1003-4374.2019.02.002
    [13]

    SUI J, LUAN S, LUO K, et al. Genetic parameters and response to selection of harvest body weight of the Chinese shrimp Fenneropenaeus chinensis after five generations of multi-trait selection[J]. Aquaculture, 2016, 452: 134-141. doi: 10.1016/j.aquaculture.2015.08.011

    [14] 于爱清. “新吉富”罗非鱼选育后期世代生长和遗传变异分析[D]. 上海: 上海海洋大学, 2011: 9-13.
    [15] 唐首杰, 何安元, 李思发, 等. “新吉富”罗非鱼选育后期世代F13—F15的生长性能比较研究[J]. 上海海洋大学学报, 2013, 22(1): 1-6.
    [16] 卢薛, 孙际佳, 王海芳, 等. 鳜鱼生长性状遗传参数的估计[J]. 中国水产科学, 2016, 23(6): 1268-1278.
    [17] 李祥孔, 田永胜, 李洪, 等. 牙鲆(Paralichthys olivaceus)家系生长性状遗传效应分析[J]. 渔业科学进展, 2017, 38(3): 39-50. doi: 10.11758/yykxjz.20160302003
    [18]

    NEIRA R, LHORENTE J, ARANEDA C, et al. Studies on carcass quality traits in two populations of coho salmon (Oncorhynchus kisutch): phenotypic and genetic parameters[J]. Aquaculture, 2004, 241(14): 117-131.

    [19]

    WANG C H, LI S F, XIANG S P, et al. Genetic parameter estimates for growth-related traits in Oujiang color common carp (Cyprinus carpio var. color)[J]. Aquaculture, 2006, 259(14): 103-107.

    [20] 王俊, 匡友谊, 佟广香, 等. 不同温度下哲罗鲑幼鱼生长性状的遗传参数估计[J]. 中国水产科学, 2011, 18(1): 75-82.
    [21] 姜鹏, 韩林强, 白俊杰, 等. 草鱼生长性状的遗传参数和育种值估计[J]. 中国水产科学, 2018, 25(1): 18-25.
    [22] 金武, 李家乐, 付龙龙, 等. 三角帆蚌早期阶段生长性状遗传参数估计[J]. 水产学报, 2012, 36(8): 1209-1214.
    [23] 袁志发, 常智杰, 郭满才, 等. 数量性状遗传分析[M]. 北京: 科学出版社, 2015: 109-123.
    [24] 唐瞻杨, 陈文治, 罗永巨, 等. 尼罗罗非鱼不同月龄数量性状遗传力估计[J]. 海洋与湖沼, 2015, 46(5): 1180-1185.
    [25] 罗坤, 夏永涛, 王斌, 等. 俄罗斯鲟早期生长性状遗传参数的估计[J]. 中国水产科学, 2015, 22(3): 426-432.
    [26]

    HANNE M N, JØRGEN Ø, INGRID O, et al. Genetic analysis of common carp (Cyprinus carpio) strains: I: genetic parameters and heterosis for growth traits and survival[J]. Aquaculture, 2010, 304: 14-21. doi: 10.1016/j.aquaculture.2010.03.016

    [27] 韦信键, 刘贤德, 王志勇. 32个大黄鱼家系早期阶段生长性状比较及遗传参数估计[J]. 集美大学学报 (自然科学版), 2013, 18(5): 321-328.
    [28] 曾聪, 曹小娟, 高泽霞, 等. 团头鲂生长性状的遗传力和育种值估计[J]. 华中农业大学学报, 2014, 33(2): 89-95. doi: 10.3969/j.issn.1000-2421.2014.02.017
    [29] 苟盼盼, 王秀利, 窦冬雨, 等. 红鳍东方鲀不同家系群体的形态性状差异与相关性分析[J]. 大连海洋大学学报, 2019, 34(5): 674-679.
    [30]

    GJEDREM T. Genetic improvement of cold-water fish species[J]. Aquac Res, 2000, 31(1): 25-33. doi: 10.1046/j.1365-2109.2000.00389.x

    [31]

    CAMPOS E C, OLIVEIRA C, ARAÚJO F C T, et al. Genetic parameters and response to selection for growth in tambaqui[J]. Animal, 2020, 14(9): 1-9.

    [32]

    DONG Z, NGUYEN N H, ZHU W. Genetic evaluation of a selective breeding program for common carp Cyprinus carpio conducted from 2004 to 2014[J]. BMC Genet, 2015, 16(1): 94. doi: 10.1186/s12863-015-0256-2

    [33] 马爱军, 王新安, 黄智慧, 等. 大菱鲆家系选育F2早期选择反应和现实遗传力估计[J]. 海洋与湖沼, 2012, 43(1): 57-61. doi: 10.11693/hyhz201201009009
    [34] 莫日馆, 肖述, 秦艳平, 等. 深凹壳型香港牡蛎选育群体生长性状的遗传参数估计[J]. 水产学报, 2020, 44(1): 33-42.
    [35]

    BENTSEN H B, GJERDE B, EKNATH A E, et al. Genetic improvement of farmed tilapias: response to five generations of selection for increased body weight at harvest in Oreochromis niloticus and the further impact of the project[J]. Aquaculture, 2017, 468(1): 206-217.

    [36]

    THODESEN J, RYE M, WANG Y X, et al. Genetic improvement of tilapias in China: genetic parameters and selection responses in fillet traits of Nile tilapia (Oreochromis niloticus) after six generations of multi-trait selection for growth and fillet yield[J]. Aquaculture, 2012, 366/367: 67-75. doi: 10.1016/j.aquaculture.2012.08.028

    [37] 鲁翠云, 李超, 曹顶臣, 等. 标记指导大规模混合家系选育技术应用实例[J]. 水产学杂志, 2019, 032(5): 8-12. doi: 10.3969/j.issn.1005-3832.2019.05.002
    [38] 程磊, 何苹萍, 韦嫔媛, 等. 基于线粒体D-loop区和Cyt b基因分析广西禾花鲤三个群体遗传结构[J]. 水生生物学报, 2021, 45(1): 54-59. doi: 10.7541/2021.2020.041
推荐阅读
面粉类型对油炸裹糊罗非鱼肉品质的影响
云万里 et al., 南方水产科学, 2025
Lh原油和0# 柴油乳化液对凡纳滨对虾肝胰腺抗氧化酶活性及相关功能基因表达的影响
沈楚焰 et al., 南方水产科学, 2025
基于声学技术定点监测海洋牧场鱼类资源昼夜变化研究
程高 et al., 南方水产科学, 2024
草鱼subfatin分子鉴定及表达特性分析
杨博雅 et al., 南方水产科学, 2024
低盐水体na+/k+对凡纳滨对虾生长、体成分与肝胰腺、鳃组织结构的影响
祝华萍 et al., 水产学报, 2024
发酵豆粕部分替代鱼粉后添加晶体赖氨酸对凡纳滨对虾肝胰腺代谢的转录调控
郭旭 et al., 水产学报, 2023
Subtypes of avoidant/restrictive food intake disorder in children and adolescents: a latent class analysis
Sanchez-Cerezo, Javier et al., ECLINICALMEDICINE, 2024
Personalizing eating disorder treatment using idiographic models: an open series trial
Levinson, Cheri A., JOURNAL OF CONSULTING AND CLINICAL PSYCHOLOGY, 2023
Mc-lr aggravates liver lipid metabolism disorders in obese mice fed a high-fat diet via pi3k/akt/mtor/srebp1 signaling pathway
TOXINS
Deep learning models in classifying primary bone tumors and bone infections based on radiographs
NPJ PRECISION ONCOLOGY, 2025
Powered by
图(2)  /  表(6)
计量
  • 文章访问数:  766
  • HTML全文浏览量:  159
  • PDF下载量:  71
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-25
  • 修回日期:  2022-01-08
  • 录用日期:  2022-02-06
  • 网络出版日期:  2022-04-11
  • 刊出日期:  2022-12-04

目录

/

返回文章
返回