中国畜牧兽医 ›› 2025, Vol. 52 ›› Issue (2): 512-521.doi: 10.16431/j.cnki.1671-7236.2025.02.002

• 生物技术 • 上一篇    

腹泻猪群肠道微生物宏基因组分析

张洁1, 贾骁晔2, 童泽宇3, 马岩2, 雷晨莹1, 原霖2, 田克恭1,4   

  1. 1. 河南农业大学动物医学院, 郑州 450003;
    2. 北京中科基因技术股份有限公司, 北京 100000;
    3. 扬州大学兽医学院, 扬州 225009;
    4. 国家兽用药品工程技术研究中心, 洛阳 471003
  • 收稿日期:2024-06-11 发布日期:2025-01-24
  • 通讯作者: 原霖, 田克恭 E-mail:yuanlin_1982@163.com;vetvac@126.com
  • 作者简介:张洁,E-mail:zj9308@126.com。
  • 基金资助:
    国家重点研发计划项目(2023YFD1800505、2022YFC2304005)

Metagenomic Analysis of Gut Microbiota in Diarrheal Pig Herds

ZHANG Jie1, JIA Xiaoye2, TONG Zeyu3, MA Yan2, LEI Chenying1, YUAN Lin2, TIAN Kegong1,4   

  1. 1. College of Veterinary Medicine, Henan Agricultural University, Zhengzhou 450003, China;
    2. Beijing Sino-science Gene Technology Co., Ltd., Beijing 100000, China;
    3. College of Veterinary Medicine, Yangzhou University, Yangzhou 225009, China;
    4. National Research Center for Veterinary Medicine, Luoyang 471003, China
  • Received:2024-06-11 Published:2025-01-24

摘要: 【目的】 了解腹泻猪群病毒混合感染情况及病毒感染下肠道微生物群落结构变化。【方法】 收集73份临床猪腹泻样品,使用单端100 bp读长(SE100)进行宏基因组测序。下机数据经FastQC工具质控获得clean data,利用BWA工具去除宿主序列获得unmapped data。使用Kraken 2工具将unmapped data比对到动物病原数据库进行物种分类注释并生成物种分类表。使用主成分分析(PCA)对数据降维处理,并使用每百万读数(RPM)进行数据归一化,利用Z-score统计方法生成热图进行差异分析可视化。使用MEGAHIT工具从头拼接病毒序列,使用BWA将unmapped data比对到contigs评估拼接效果,选择高质量contigs在NCBI上比对找到最相似序列,将ummapped data比对到该序列后提取一致序列,作为该病毒序列,使用IQ-TREE工具进行猪流行性腹泻病毒(PEDV)S基因和猪星状病毒(PAstV)ORF1a基因遗传进化分析。【结果】 73份猪腹泻样品中检出多种病毒,阳性率分别为:PAstV 26.03%、猪嵴病毒(PKV)20.55%、猪轮状病毒(PoRV)19.18%、PEDV 19.18%、猪札幌病毒(PSaV)6.85%、猪德尔塔冠状病毒(PDCoV)4.11%。混合感染情况复杂,主要为PEDV+PAstV组合,在全部样品中占比为17.81%(13/73),在混合感染样品中占比为30.95%(13/42)。不同病毒感染下肠道微生物结构发生变化,PEDV阳性组样品乳杆菌属(Lactobacillus)相对丰度高于阴性组,PAstV阳性组样品梭菌属(Clostridium)和密螺旋体属(Treponema)相对丰度高于阴性组。拼接得到2条PEDV全基因组序列SE18和SE22,均属于GⅡa谱系,以及4条PAstV ORF1a基因组序列SE18、PE31、PE36和PE70,SE18为PAstV2型,其余3条序列与HAstV-HMO-B株亲缘关系较近。【结论】 本研究通过宏基因组学技术揭示了临床腹泻猪群中PAstV和PEDV的高阳性率及混合感染情况,发现PEDV和PAstV感染显著改变了肠道微生物结构,为猪腹泻混合感染的防控提供了重要数据支持。

关键词: 猪; 腹泻; 混合感染; 宏基因组学

Abstract: 【Objective】 This study was aimed to understand the prevalence of viral co-infections in diarrheic pig populations and the alterations in the gut microbiota community structure under viral infection. 【Method】 Seventy-three clinical porcine diarrhea samples were collected and subjected to metagenomic sequencing using single-end 100 bp read length (SE100).The raw data were quality-controlled using FastQC to obtain clean data,and host sequences were removed using BWA to generate unmapped data.Kraken 2 was used to align the unmapped data to an animal pathogen database for species classification annotation and to generate a species classification table.Principal component analysis (PCA) was performed for dimensionality reduction,and the data were normalized using reads per million (RPM).Z-score statistical methods were used to generate heatmaps for differential analysis visualization.MEGAHIT was employed for de novo assembly of viral sequences,and BWA was used to align unmapped data to contigs to assess assembly quality.High-quality contigs were selected and aligned to the NCBI database to find the most similar sequences.The unmapped data were aligned to these sequences to extract consensus sequences,which were used as the viral sequences.Phylogenetic analysis of Porcine epidemic diarrhea virus (PEDV) S gene and Porcine astrovirus (PAstV) ORF1a gene was performed using IQ-TREE. 【Result】 Various viruses were detected in 73 samples of pig diarrhea with the following positivity rates:PAstV 26.03%,Porcine kobuvirus (PKV) 20.55%,Porcine rotavirus (PoRV) 19.18%,PEDV 19.18%,Porcine sapporo virus (PSaV) 6.85% and Porcine deltacoronavirus (PDCoV) 4.11%.The mixed infection situation was complex,mainly consisting of a combination of PEDV and PAstV,accounting for 17.81% (13/73) of all samples and 30.95% (13/42) of mixed infection samples.The gut microbiota structure showed significant changes under different viral infections.The relative abundance of Lactobacillus was higher in PEDV-positive group than that in PEDV-negative group,and the relative abundances of Clostridium and Treponema were higher in PAstV-positive group than that in PAstV-negative group.Two complete PEDV genomes SE18 and SE22 were assembled,belonging to the GⅡa lineage,along with four PAstV ORF1a gene sequences SE18,PE31,PE36 and PE70,of which SE18 belonged to PAstV2 and the other three were closely related to HAstV-HMO-B strain. 【Conclusion】 This study revealed high positive rates and co-infection patterns of PAstV and PEDV in clinically diarrheic pig populations using metagenomic techniques,and also identified significant changes in gut microbiota composition due to PEDV and PAstV infections,providing critical data support for the prevention and control of mixed viral infections in pig diarrhea.

Key words: pig; diarrhea; mixed infection; metagenomics

中图分类号: