瓦氏黄颡鱼全基因组微卫星的分布特征及其定位的初步研究

彭冶, 李杰, 王涛, 张凯, 宁先会, 暨杰, 尹绍武

彭冶, 李杰, 王涛, 张凯, 宁先会, 暨杰, 尹绍武. 瓦氏黄颡鱼全基因组微卫星的分布特征及其定位的初步研究[J]. 南方水产科学, 2022, 18(1): 90-98. DOI: 10.12131/20210168
引用本文: 彭冶, 李杰, 王涛, 张凯, 宁先会, 暨杰, 尹绍武. 瓦氏黄颡鱼全基因组微卫星的分布特征及其定位的初步研究[J]. 南方水产科学, 2022, 18(1): 90-98. DOI: 10.12131/20210168
PENG Ye, LI Jie, WANG Tao, ZHANG Kai, NING Xianhui, JI Jie, YIN Shaowu. Preliminary study on distribution characteristics and positioning of microsatellites in whole genome of Pelteobagrus vachelli[J]. South China Fisheries Science, 2022, 18(1): 90-98. DOI: 10.12131/20210168
Citation: PENG Ye, LI Jie, WANG Tao, ZHANG Kai, NING Xianhui, JI Jie, YIN Shaowu. Preliminary study on distribution characteristics and positioning of microsatellites in whole genome of Pelteobagrus vachelli[J]. South China Fisheries Science, 2022, 18(1): 90-98. DOI: 10.12131/20210168

瓦氏黄颡鱼全基因组微卫星的分布特征及其定位的初步研究

基金项目: 国家自然科学基金青年基金(32102760); 江苏省农业重大新品种创制项目(PZCZ201742); 江苏省农业科技自主创新资金项目 [CX(19)2034]; 江苏省高校自然科学研究面上项目(18KJB240001, 20KJD240001); 浙江省淡水水产研究所农业部淡水渔业健康养殖重点实验室开放课题(重大)(ZJK202003)
详细信息
    作者简介:

    彭 冶 (1999—),男,硕士研究生,研究方向为鱼类种质资源与遗传育种。E-mail: 17839730970@163.com

    通讯作者:

    尹绍武 (1969—),男,教授,博士,从事鱼类种质资源与遗传育种研究。E-mail: yinshaowu@163.com

  • 中图分类号: S 917.4

Preliminary study on distribution characteristics and positioning of microsatellites in whole genome of Pelteobagrus vachelli

  • 摘要: 文章探究了瓦氏黄颡鱼 (Pelteobagrus vachelli) 全基因组微卫星的分布特征及其规律,旨在为相关功能性微卫星分子标记的筛选提供依据。利用MISA (MIcroSAtellite identification tool) 软件对其全基因组微卫星进行筛查和分析,并对外显子中含有微卫星的基因进行了GO (Gene Ontology) 注释和KEGG (Kyoto Encyclopedia of Genes and Genomes) 富集分析。在瓦氏黄颡鱼全基因组 (约663.53 Mb) 中筛查到6种完整型微卫星共417 724个,相对丰度为630 个·Mb−1,在全基因组总长度中占比1.48%。其中二碱基重复类型的微卫星个数最多,占微卫星总数的43.36%,其他依次分别是单碱基 (39.02%)、四碱基 (9.05%)、三碱基 (7.34%)、五碱基 (1.12%)和六碱基 (0.12%)。对筛选得到的完整型微卫星进行初步定位,发现其中10 924个微卫星分布在外显子区,共定位到5 788个基因上。GO分析表明注释到生物过程的基因数量最多,GO富集较为显著的条目为结合活性和细胞大分子代谢过程。KEGG富集表明这些基因共富集到273条通路,其中黄酮与黄酮醇生物合成等通路最为显著 (P=0)。联合分析预测瓦氏黄颡鱼定位到外显子上的微卫星和其体内的生物代谢过程密切相关。
    Abstract: In this study, MISA (MIcroSAtellite identification tool) was used to screen and analyze the distribution characteristics of microsatellites in the whole genome of Pelteobagrus vachelli, aiming to provide a basis for the selection of functional microsatellite markers. The genes containing microsatellites in the exon regions were subjected to GO annotation and KEGG enrichment. In the whole genome of P. vachelli (Approximately 663.53 Mb), 417 724 perfect microsatellites were identified, accounting for 1.48% of the total length. The relative abundance of microsatellites in P. vachelli was 630 pcs·Mb−1. Among the six repeat types, dinucleotides were the most frequent, accounting for 43.36% of the total microsatellites, followed by mononucleotides (39.02%), tetranucleotides (9.05%), trinucleotides (7.34%), pentanucleotides (1.12%) and hexanucleotides (0.12%). By investigating the locations of microsatellites in the genome, we found that 10 924 microsatellites which belonged to 5 788 genes were located in the exons. The GO annotation shows that the number of genes annotated to biological process was the largest, mainly associated with binding activity and cellular macromolecular metabolism. KEGG enrichment analysis shows that these genes were enriched in 273 pathways, among which, flavonoids and flavonol biosynthesis (P=0) were the most significantly enriched pathways. Integrated analysis indicates that the microsatellites located in the exons of genes in P. vachelli were closely related to the biological metabolism processes.
  • 鳗败血假单胞菌 (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.   Distribution of different copy numbers of various types of microsatellites in P. vachelli

    图  2   瓦氏黄颡鱼微卫星分布于外显子的基因GO功能注释

    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  2.   GO function annotation of genes with microsatellites in exons in P. vachelli

    1. Binding; 2. Catalytic activity; 3. Transporter activity; 4. Nucleic acid binding transcription factor activity; 5. Molecular transducer activity; 6. Signal transducer activity; 7. Molecular function regulator; 8. Structural molecule activity; 9. Transcription factor activity, protein binding; 10. Translation regulator activity; 11. Cell; 12. Cell part; 13. Membrane; 14. Organelle; 15. Membrane part; 16. Macromolecular complex;17. Organelle part; 18. Extracellular region; 19. Membrane-enclosed lumen; 20. Extracellular region part; 21. Cellular process;22. Metabolic process; 23. single-organism process; 24. Biological regulation; 25. Regulation of biological process; 26. Response to stimulus; 27. Signaling; 28. Localization; 29. Cellular component organizationor biogenesis; 30. Multicellular organismal process; 31. Developmental process; 32. Negative regulation of biological process; 33. Positive regulation of biological process; 34. Biological adhesion; 35. Locomotion.

    表  1   瓦氏黄颡鱼基因组中不同类型微卫星统计

    Table  1   Summary of different types of microsatellite in genome of P. vachelli

    重复类型
    Repeat type
    数量
    Total
    number/个
    占比
    Proportion/%
    总长度
    Total
    length/bp
    占比
    Proportion/%
    平均长度
    Average
    length/bp
    频率
    Frequency/
    (个·Mb−1)
    密度
    Density/
    (bp·Mb−1)
    单核苷酸 Mononucleotide 162 987 39.02 2 137 066 21.75 13.11 245.63 3 220.73
    二核苷酸 Dinucleotide 181 107 43.36 5 291 784 53.85 29.22 272.94 7 975.15
    三核苷酸 Trinucleotide 30 661 7.34 956 385 9.73 31.19 46.21 1 441.35
    四核苷酸 Tetranucleotide 37 787 9.05 1 228 528 12.50 32.51 56.95 1 851.49
    五核苷酸 Pentanucleotide 4 689 1.12 188 890 1.92 40.28 7.07 284.67
    六核苷酸 Hexanucleotide 493 0.12 23 472 0.24 47.61 0.74 35.37
    下载: 导出CSV

    表  2   次数排名前10的重复碱基类别

    Table  2   Top ten types of SSRs with most repeated copy numbers

    SSR类别
    SSR type
    SSR出现次数
    SSR occurrence
    占各碱基类型比例
    Proportion of each base type/%
    A 158 915 97.50
    AC 130 224 71.90
    AG 31 118 17.18
    AT 19 593 10.82
    AAT 17 101 55.77
    AAAT 9 490 25.11
    AATG 9 210 24.37
    ATC 4 106 13.39
    AGAT 3 755 9.94
    AAAG 2 865 7.58
    下载: 导出CSV

    表  3   瓦氏黄颡鱼基因组中排名前3的微卫星重复拷贝类别

    Table  3   Top three dominant base classes in each base repeat type in P. vachelli genome

    单碱基
    Mononucleotide
    二碱基
    Dinucleotide
    三碱基
    Trinucleotide
    四碱基
    Tetranucleotide
    五碱基
    Pentanucleotide
    六碱基
    Hexanucleotide
    A
    158 915a
    97.50%b
    AC
    13 002a
    71.90%b
    AAT
    17 101a
    55.77%b
    AAAT
    9 490a
    25.11%b
    AATCT
    365a
    7.78%b
    GGGTTA
    111a
    22.5%b
    C
    4 072a
    2.50%b
    AG
    31 118a
    17.18%b
    ATC
    4 106a
    13.39%b
    AATG
    9 210a
    24.37%b
    AAAGA
    305a
    6.50%b
    CTAACC
    73a
    14.8%b
    AT
    19 593a
    10.82%b
    AAG
    3 060a
    9.98%b
    AGAT
    3 755a
    9.94%b
    AAAAT
    210a
    4.48%b
    TGTAAA
    54a
    10.95%b
    注:a. 微卫星的个数;b. 该种微卫星占其碱基类型的比例。 Note: a. Number of microsatellites; b. Proportion of the microsatellites in the base type.
    下载: 导出CSV

    表  4   瓦氏黄颡鱼微卫星分布于外显子的基因GO富集

    Table  4   GO enrichment of genes with microsatellites located in exons in P. vachelli

    条目
    Term
    GO码
    GO ID
    输入数
    Input number/个
    背景数
    Background number/个
    P
    结合 Binding 0005488 2 660 8 375 7.20×10−17
    细胞大分子代谢过程 Cellular macromolecule Metabolic process 0044260 976 2 782 1.49×10−14
    核 Nucleus 0005634 371 920 3.38×10−14
    大分子代谢过程 Macromolecule metabolic process 0043170 1 095 3 196 3.70×10−13
    细胞膜结合细胞器 Intracellular membrane-bounded organelle 0043231 475 1 263 5.81×10−12
    细胞代谢过程 Cellular metabolic process 0044237 1 134 3 353 6.36×10−12
    膜结合细胞器 Membrane-bounded organelle 0043227 486 1 306 2.13×10−11
    基因表达调控 Regulation of gene expression 0010468 364 939 3.10×10−11
    RNA生物合成过程 RNA biosynthetic process 0032774 375 972 3.13×10−11
    转录、DNA模板化 Transcription, DNA-templated 0006351 374 970 3.65×10−11
    下载: 导出CSV

    表  5   瓦氏黄颡鱼微卫星分布于外显子的基因KEGG富集

    Table  5   KEGG enrichment of genes with microsatellites in exons in P. vachelli

    功能分类   
    Functional classification  
    通路
    Pathway
    输入数
    Input number/个
    背景数
    Background number/个
    P
    机体系统
    Organismal system
    甲状腺激素信号通路
    Thyroid hormone signaling pathway
    66 147 1.45×10−5
    胰岛素信号通路
    Insulin signaling pathway
    73 173 7.57×10−5
    神经营养因子信号通路
    Neurotrophin signaling pathway
    64 156 4.83×10−4
    昼夜节律
    Circadian rhythm
    22 46 2.18×10−3
    背腹轴形成
    Dorso-ventral axis formation
    16 32 3.61×10−3
    细胞过程
    Cellular process
    黏合连接
    Adherens junction
    59 115 1.51×10−7
    内吞作用
    Endocytosis
    135 345 1.86×10−5
    黏着
    Focal adhesion
    104 280 1.29×10−3
    凋亡-果蝇
    Apoptosis-fly
    29 63 1.37×10−3
    自噬-酵母
    Autophagy-yeast
    32 75 4.08×10−3
    环境信息处理
    Environmental information processing
    ErbB信号通路
    ErbB signaling pathway
    54 115 1.52×10−5
    Hedgehog 信号通路-果蝇
    Hedgehog signaling pathway-fly
    21 34 1.70×10−5
    Hippo信号通路
    Hippo signaling pathway
    84 200 3.07×10−5
    FoxO信号通路
    FoxO signaling pathway
    76 179 4.34×10−5
    Notch信号通路
    Notch signaling pathway
    35 69 4.48×10−5
    代谢
    Metabolism
    黄酮与黄酮醇生物合成
    Flavone and flavonol biosynthesis
    1 1 0.00
    单巴坦生物合成
    Monobactam biosynthesis
    2 2 0.00
    角质的生物合成
    Cutin,suberine and wax biosynthesis
    1 1 0.00
    安莎霉素的生物合成
    Biosynthesis of ansamycins
    2 2 0.00
    香叶醇降解
    Geraniol degradation
    2 2 0.00
    遗传信息处理
    Genetic information processing
    mRNA监控通路
    mRNA surveillance pathway
    42 96 7.33×10−4
    基础转录因子
    Basal transcription factors
    19 38 1.90×10−3
    剪接体
    Spliceosome
    49 124 4.54×10−3
    RNA转运
    RNA transport
    63 168 6.96×10−3
    真核生物核糖体的生物合成
    Ribosome biogenesisin eukaryotes
    30 74 1.22×10−2
    下载: 导出CSV
  • [1]

    TAUTZ D, RENZ M. Simple sequences are ubiquitous repetitive components of eukaryotic genomes[J]. Nucleic Acids Res, 1984, 12(10): 4127-4138. doi: 10.1093/nar/12.10.4127

    [2] 何平. 真核生物中的微卫星及其应用[J]. 遗传, 1998, 20(4): 44-49.
    [3] 罗文永, 胡骏, 李晓方. 微卫星序列及其应用[J]. 遗传, 2003, 25(5): 615-619. doi: 10.3321/j.issn:0253-9772.2003.05.027
    [4]

    ZHAO X Y, TAN Z Y, FENG H P, et al. Microsatellites in different Potyvirus genomes: survey and analysis[J]. Gene, 2011, 488(1/2): 52-56.

    [5] 张增翠, 侯喜林. SSR分子标记开发策略及评价[J]. 遗传, 2004(5): 763-768. doi: 10.3321/j.issn:0253-9772.2004.05.037
    [6]

    WU X P, FENG Y W, JIANG H L, et al. A preliminary genetic linkage map of Sinonovacula constricta (Lamarck, 1818) based on microsatellites derived from RAD sequencing[J]. J Ocean U China, 2018, 17(4): 947-956. doi: 10.1007/s11802-018-3543-5

    [7] 赵娜, 常剑波, 陶江平, 等. 基于微卫星标记的中华鲟亲子关系判别及案例分析[J/OL]. 水生态学杂志. [2021-01-15]. http://kns.cnki.net/kcms/detail/42.1785.X.20200915.1323.002.html.
    [8] 苏胜彦, 张林兵, 李海洋, 等. 基于微卫星标记的大口黑鲈 (Micropterus salmoides) 原种和养殖群体遗传多样性和结构分析[J]. 浙江大学学报 (农业与生命科学版), 2020, 46(6): 687-698.
    [9] 蔡焰值, 蔡烨强, 何长仁. 瓦氏黄颡鱼生物学的初步研究[J]. 北京水产, 2003(6): 24-29.
    [10] 李明锋. 瓦氏黄颡鱼研究进展及前景展望[J]. 现代渔业信息, 2011, 26(1): 5-12.
    [11] 巩高瑞, 张晋, 丹成, 等. 应用DNA分子标记鉴定黄颡鱼、瓦氏黄颡鱼及其杂交种的研究[J]. 水生生物学报, 2017, 41(2): 321-325. doi: 10.7541/2017.39
    [12] 葛学亮, 尹洪滨, 毕冰, 等. 黄颡鱼遗传图谱构建及生长相关性状的QTL定位[J]. 水产学报, 2010, 34(2): 185-193.
    [13] 李大宇, 殷倩茜, 侯宁, 等. 黄颡鱼 (Pelteobagrus fulvidraco) 不同生态地理分布群体遗传多样性的微卫星分析[J]. 海洋与湖沼, 2009, 40(4): 460-469. doi: 10.3321/j.issn:0029-814X.2009.04.011
    [14] 李林, 梁宏伟, 李忠, 等. 瓦氏黄颡鱼线粒体全基因组序列分析及系统进化[J]. 遗传, 2011, 33(6): 627-635.
    [15] 郑翔, 徐杰杰, 张佳佳, 等. 4个瓦氏黄颡鱼群体遗传多样性的微卫星分析[J]. 水产科学, 2020, 39(5): 657-668.
    [16] 何晶晶, 黄建华, 牛红艳, 等. 猪蛔虫全基因组微卫星分子标记开发与特征分析[J]. 中国兽医杂志, 2020, 56(9): 11-14, 19.
    [17] 张雪莲, 王红梅, 王磊, 等. 草地贪夜蛾基因组微卫星的分布规律[J]. 应用昆虫学报, 2020, 57(6): 1287-1298.
    [18] 徐杰杰, 毕宜慧, 程景颢, 等. 中华绒螯蟹 (Eriocheir sinensis) 全基因组微卫星分布特征研究[J/OL]. 基因组学与应用生物学. [2021-03-21]. http://61.175.198.136:8083/rwt/CNKI/http/NNYHGLUDN3WXTLUPMW4A/kcms/detail/45.1369.Q.20210118.1225.006.html.
    [19] 段永楠, 刘奕, 胡隐昌, 等. 美丽硬仆骨舌鱼全基因组微卫星分布规律特征[J]. 中国农学通报, 2019, 35(23): 152-158. doi: 10.11924/j.issn.1000-6850.casb18030101
    [20] 王耀嵘, 杨尉, 任席林, 等. 金钱鱼基因组微卫星分布特征分析及多态性标记开发[J]. 广东海洋大学学报, 2020, 40(4): 7-14. doi: 10.3969/j.issn.1673-9159.2020.04.002
    [21] 梁霞, 王慧琪, 马宇璇, 等. 鲤鱼 (Cyprinus carpio) 全基因组微卫星分布特征研究[J]. 南京师范大学学报(自然科学版), 2021, 44(3): 103-111.
    [22]

    CONESA A, GÖTZ S, GARCÍA-GÓMEZ J M, et al. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research[J]. Bioinformatics, 2005, 21(18): 3674-3676. doi: 10.1093/bioinformatics/bti610

    [23]

    XIE C, MAO X, HUANG J, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases[J]. Nucleic Acids Res, 2011, 39(Suppl2): W316-W322. doi: 10.1093/nar/gkr483

    [24]

    SUBRAMANIAN S, MISHRA R K, SINGH L. Genome-wide analysis of microsatellite repeats in humans: their abundance and density in specific genomic regions[J]. Genome Biol, 2003, 4(2): 1-10. doi: 10.1186/gb-2003-4-2-p1

    [25] TONG X L, DAI F Y, LI B, et al. Microsatellite repeats in mouse: abundance, distribution and density[J]. 动物学报, 2006, 52(1): 138-152.
    [26] 戚文华, 蒋雪梅, 肖国生, 等. 牛和绵羊全基因组微卫星序列的搜索及其生物信息学分析[J]. 畜牧兽医学报, 2013, 44(11): 1724-1733.
    [27] 徐杰杰, 郑翔, 张鑫宇, 等. 4种河鲀全基因组微卫星分布特征分析研究[J/OL]. 基因组学与应用生物学. [2021-01-25]. http://kns.cnki.net/kcms/detail/45.1369.q.20191126.1019.002.html.
    [28] 徐杰杰, 郑翔, 李杰, 等. 黄颡鱼 (Pelteobagrus fulvidraco) 全基因组微卫星分布特征分析[J]. 基因组学与应用生物学, 2020, 39(12): 5488-5498.
    [29] 王月月, 刘雪雪, 董坤哲, 等. 7种家养动物全基因组微卫星分布的差异研究[J]. 中国畜牧兽医, 2015, 42(9): 2418-2426.
    [30] 阮晓红. 大菱鲆(Turbot)微卫星标记的筛选与应用[D]. 青岛: 中国海洋大学, 2009: 53-54.
    [31] 黄杰, 杜联明, 李玉芝, 等. 红原鸡全基因组中微卫星分布规律研究[J]. 四川动物, 2012, 31(3): 358-363.
    [32] 崔凯, 岳碧松. 绿尾虹雉全基因组微卫星分布规律研究[J]. 四川动物, 2018, 37(5): 533-540. doi: 10.11984/j.issn.1000-7083.20180096
    [33] 黄杰, 刘磊, 杨波, 等. 普通鸬鹚基因组微卫星分布规律研究[J]. 野生动物学报, 2020, 41(1): 108-114. doi: 10.3969/j.issn.1000-0127.2020.01.015
    [34]

    TÓTH G, GÁSPÁRI Z, JURKA J. Microsatellites in different eukaryotic genomes: survey and analysis[J]. Genome Res, 2000, 10(7): 967-981. doi: 10.1101/gr.10.7.967

    [35] 上官清, 陈昆慈, 刘海洋, 等. 斑鳢基因组中微卫星分布特征及野生种群遗传结构分析[J]. 南方水产科学, 2020, 16(3): 47-60.
    [36]

    SCHORDERET D F, GARTLER S M. Analysis of CpG suppression in methylated and nonmethylated species[J]. P Natl Acad Sci USA, 1992, 89(3): 957-961. doi: 10.1073/pnas.89.3.957

    [37]

    PEARSON C E, SINDEN R R. Alternative structures in duplex DNA formed within the trinucleotide repeats of the myotonic dystrophy and fragile X loci[J]. Biochemistry-US, 1996, 35(15): 5041-5053. doi: 10.1021/bi9601013

    [38]

    XU Y, LI W, HU Z, et al. Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species[J]. Gene, 2018, 643: 124-132. doi: 10.1016/j.gene.2017.12.008

    [39]

    XIAO T W, YU J Z, LIAO Q, et al. Comparative analyses of simple sequence repeats (SSRs) in 23 mosquito species genomes: identification, characterization and distribution (Diptera: Culicidae)[J]. Insect Sci, 2019, 26: 607-619. doi: 10.1111/1744-7917.12577

    [40]

    ELLEGREN H. Heterogeneous mutation processes in human microsatellite DNA sequences[J]. Nature Genet, 2000, 24(4): 400-402. doi: 10.1038/74249

    [41]

    WIERDL M, DOMINSKA M, PETES T D. Microsatellite instability in yeast: dependence on the length of the microsatellite[J]. Genetics, 1997, 146(3): 769-779. doi: 10.1093/genetics/146.3.769

推荐阅读
面粉类型对油炸裹糊罗非鱼肉品质的影响
云万里 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)  /  表(5)
计量
  • 文章访问数:  869
  • HTML全文浏览量:  216
  • PDF下载量:  52
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-06-03
  • 修回日期:  2021-06-29
  • 录用日期:  2021-07-22
  • 网络出版日期:  2021-08-03
  • 刊出日期:  2022-02-04

目录

/

返回文章
返回