ACTA THERIOLOGICA SINICA ›› 2026, Vol. 46 ›› Issue (3): 337-348.DOI: 10.16829/j.slxb.151032
• ORIGINAL PAPERS • Previous Articles Next Articles
Ying MA1, Youwen WEI1, Feng JIANG2, Aiping ZHANG1, Hongmei GAO2, Jide LI1, Wenlong WANG1, Penghui ZHANG1, Lingwen LI1, Xinyue GU1, Tongzuo ZHANG2(
)
Received:2024-11-28
Accepted:2025-02-24
Online:2026-05-30
Published:2026-06-04
Contact:
Tongzuo ZHANG
马英1, 魏有文1, 江峰2, 张爱萍1, 高红梅2, 李积德1, 王文龙1, 张鹏辉1, 李玲雯1, 谷昕玥1, 张同作2(
)
通讯作者:
张同作
作者简介:马英(1973- ),女,硕士,主要从事鼠疫防控和动物遗传学研究;基金资助:CLC Number:
Ying MA, Youwen WEI, Feng JIANG, Aiping ZHANG, Hongmei GAO, Jide LI, Wenlong WANG, Penghui ZHANG, Lingwen LI, Xinyue GU, Tongzuo ZHANG. Metagenomic analysis of intestinal virome in Marmota himalayana in eastern Qinghai’s agricultural area[J]. ACTA THERIOLOGICA SINICA, 2026, 46(3): 337-348.
马英, 魏有文, 江峰, 张爱萍, 高红梅, 李积德, 王文龙, 张鹏辉, 李玲雯, 谷昕玥, 张同作. 青海东部农业区喜马拉雅旱獭肠道病毒种群宏基因组学分析[J]. 兽类学报, 2026, 46(3): 337-348.
Add to citation manager EndNote|Ris|BibTeX
种群 Population | 采集地点 Region | 样本量(性别) Number (Gender) | 经度 Longitude | 纬度 Latitude | 海拔 Elevation/m |
|---|---|---|---|---|---|
| HDHZ | 海东市互助县 Huzhu County,Haidong City | 3(2♂ 1♀) | 102.2763°E | 36.7299°N | 2 983 |
| HDMH | 海东市民和县 Minhe County,Haidong City | 3(3♀) | 102.6256°E | 36.2035°N | 2 592 |
| HDPA | 海东市平安区 Pingan District,Haidong City | 3(1♂ 2♀) | 101.8478°E | 36.3422°N | 2 803 |
| XNHY | 西宁市湟源县 Huangyuan County,Xining City | 3(3♀) | 101.2291°E | 36.4249°N | 3 334 |
| XNHZ | 西宁市湟中县 Huangzhong County,Xining City | 3(1♂ 2♀) | 101.7416°E | 36.3320°N | 2 934 |
| XNDT | 西宁市大通县 Datong County,Xining City | 3(3♀) | 101.6220°E | 37.1058°N | 2 930 |
Table 1 Sampling information of Marmota himalayana used in the present study
种群 Population | 采集地点 Region | 样本量(性别) Number (Gender) | 经度 Longitude | 纬度 Latitude | 海拔 Elevation/m |
|---|---|---|---|---|---|
| HDHZ | 海东市互助县 Huzhu County,Haidong City | 3(2♂ 1♀) | 102.2763°E | 36.7299°N | 2 983 |
| HDMH | 海东市民和县 Minhe County,Haidong City | 3(3♀) | 102.6256°E | 36.2035°N | 2 592 |
| HDPA | 海东市平安区 Pingan District,Haidong City | 3(1♂ 2♀) | 101.8478°E | 36.3422°N | 2 803 |
| XNHY | 西宁市湟源县 Huangyuan County,Xining City | 3(3♀) | 101.2291°E | 36.4249°N | 3 334 |
| XNHZ | 西宁市湟中县 Huangzhong County,Xining City | 3(1♂ 2♀) | 101.7416°E | 36.3320°N | 2 934 |
| XNDT | 西宁市大通县 Datong County,Xining City | 3(3♀) | 101.6220°E | 37.1058°N | 2 930 |
样品 Samples | 原始数据 Raw data(G) | 原始reads数 Number of raw reads(PE) | 高质量数据 Clean data(G) | 高质量reads数 Number of clean reads(PE) | 高质量占比 Proportion of high quality/% |
|---|---|---|---|---|---|
| HDHZ | 12.5 | 41826533 | 8.9 | 29913139 | 71.52 |
| HDMH | 12.3 | 41085030 | 9.9 | 33035768 | 80.41 |
| HDPA | 11.4 | 37919759 | 8.1 | 27187912 | 71.70 |
| XNHY | 10.0 | 33475840 | 8.5 | 28559824 | 85.31 |
| XNHZ | 9.8 | 32780458 | 5.9 | 19772197 | 60.32 |
| XNDT | 11.3 | 37633105 | 9.1 | 30345487 | 80.64 |
Table 2 Summary of macro virus sequencing reads data
样品 Samples | 原始数据 Raw data(G) | 原始reads数 Number of raw reads(PE) | 高质量数据 Clean data(G) | 高质量reads数 Number of clean reads(PE) | 高质量占比 Proportion of high quality/% |
|---|---|---|---|---|---|
| HDHZ | 12.5 | 41826533 | 8.9 | 29913139 | 71.52 |
| HDMH | 12.3 | 41085030 | 9.9 | 33035768 | 80.41 |
| HDPA | 11.4 | 37919759 | 8.1 | 27187912 | 71.70 |
| XNHY | 10.0 | 33475840 | 8.5 | 28559824 | 85.31 |
| XNHZ | 9.8 | 32780458 | 5.9 | 19772197 | 60.32 |
| XNDT | 11.3 | 37633105 | 9.1 | 30345487 | 80.64 |
样品 Samples | 质控后的reads Clean reads | 比对到宿主基因组的read数 Number of reads mapped to the host genome | 最后剩余reads数 Remaining reads | 最后剩余reads数占比 Proportion of remaining reads/% |
|---|---|---|---|---|
| HDHZ | 29913139 | 487552 | 29425587 | 98.37 |
| HDMH | 33035768 | 484506 | 32551262 | 98.53 |
| HDPA | 27187912 | 1115333 | 26072579 | 95.90 |
| XNHY | 28559824 | 55236 | 28504588 | 99.81 |
| XNHZ | 19772197 | 56003 | 19716194 | 99.72 |
| XNDT | 30345487 | 124262 | 30221225 | 99.59 |
Table 3 Statistics of reads data after removing host sequences
样品 Samples | 质控后的reads Clean reads | 比对到宿主基因组的read数 Number of reads mapped to the host genome | 最后剩余reads数 Remaining reads | 最后剩余reads数占比 Proportion of remaining reads/% |
|---|---|---|---|---|
| HDHZ | 29913139 | 487552 | 29425587 | 98.37 |
| HDMH | 33035768 | 484506 | 32551262 | 98.53 |
| HDPA | 27187912 | 1115333 | 26072579 | 95.90 |
| XNHY | 28559824 | 55236 | 28504588 | 99.81 |
| XNHZ | 19772197 | 56003 | 19716194 | 99.72 |
| XNDT | 30345487 | 124262 | 30221225 | 99.59 |
样品 Samples | 组装获得的序列 碱基总数 Assembled base pairs | 组装contigs条数 Number of assembled contigs | contig长度最大值 Maximum contig length | contig长度最小值 Minimum contig length | N50值 Contig N50 | contigs GC平均含量 Average GC content of contigs |
|---|---|---|---|---|---|---|
| HDHZ | 138.51 | 150 585 | 132 879 | 300 | 1 009 | 42.17 |
| HDMH | 316.15 | 382 618 | 208 313 | 300 | 833 | 45.87 |
| HDPA | 84.25 | 92 097 | 232 212 | 300 | 984 | 43.12 |
| XNHY | 355.78 | 415 327 | 196 346 | 300 | 888 | 47.58 |
| XNHZ | 161.14 | 183 488 | 172 565 | 300 | 916 | 45.00 |
| XNDT | 212.15 | 234 266 | 304 333 | 300 | 976 | 41.98 |
Table 4 Statistics of assembly results after removing host sequences
样品 Samples | 组装获得的序列 碱基总数 Assembled base pairs | 组装contigs条数 Number of assembled contigs | contig长度最大值 Maximum contig length | contig长度最小值 Minimum contig length | N50值 Contig N50 | contigs GC平均含量 Average GC content of contigs |
|---|---|---|---|---|---|---|
| HDHZ | 138.51 | 150 585 | 132 879 | 300 | 1 009 | 42.17 |
| HDMH | 316.15 | 382 618 | 208 313 | 300 | 833 | 45.87 |
| HDPA | 84.25 | 92 097 | 232 212 | 300 | 984 | 43.12 |
| XNHY | 355.78 | 415 327 | 196 346 | 300 | 888 | 47.58 |
| XNHZ | 161.14 | 183 488 | 172 565 | 300 | 916 | 45.00 |
| XNDT | 212.15 | 234 266 | 304 333 | 300 | 976 | 41.98 |
基因名称 Name of gene | 样品Samples | |||||
|---|---|---|---|---|---|---|
| HDHZ | HDMH | HDPA | XNHY | XNHZ | XNDT | |
| SCD|contig_3238 | 0.53 | 0.02 | 0.86 | 0.03 | 0.07 | 0.01 |
| HDPA|contig_89993 | 0.42 | 0 | 0.63 | 0.01 | 0.03 | 0 |
| HNZK|contig_7503 | 0.02 | 0.16 | 33 327.41 | 0.03 | 0.07 | 0.06 |
| XNHY|contig_402915 | 293.08 | 1.95 | 0 | 319.12 | 39 730.32 | 46.19 |
| SMD|contig_147452 | 0 | 0.01 | 0 | 0 | 0 | 0 |
| HNZK|contig_65526 | 0 | 0 | 1.92 | 0.01 | 0.08 | 0.02 |
| HDPA|contig_81697 | 0 | 0.04 | 1.82 | 0 | 0.06 | 0 |
| SMD|contig_147636 | 0 | 1.62 | 0 | 172.60 | 34 127.26 | 0 |
| HDMH|contig_380818 | 0.02 | 1.97 | 0 | 176.59 | 33 716.08 | 0 |
| HNZK|contig_4700 | 0 | 0.06 | 21 873.11 | 0.31 | 0.13 | 12.75 |
Table 5 Abundance statistics of the top 10 contigs by RPKM
基因名称 Name of gene | 样品Samples | |||||
|---|---|---|---|---|---|---|
| HDHZ | HDMH | HDPA | XNHY | XNHZ | XNDT | |
| SCD|contig_3238 | 0.53 | 0.02 | 0.86 | 0.03 | 0.07 | 0.01 |
| HDPA|contig_89993 | 0.42 | 0 | 0.63 | 0.01 | 0.03 | 0 |
| HNZK|contig_7503 | 0.02 | 0.16 | 33 327.41 | 0.03 | 0.07 | 0.06 |
| XNHY|contig_402915 | 293.08 | 1.95 | 0 | 319.12 | 39 730.32 | 46.19 |
| SMD|contig_147452 | 0 | 0.01 | 0 | 0 | 0 | 0 |
| HNZK|contig_65526 | 0 | 0 | 1.92 | 0.01 | 0.08 | 0.02 |
| HDPA|contig_81697 | 0 | 0.04 | 1.82 | 0 | 0.06 | 0 |
| SMD|contig_147636 | 0 | 1.62 | 0 | 172.60 | 34 127.26 | 0 |
| HDMH|contig_380818 | 0.02 | 1.97 | 0 | 176.59 | 33 716.08 | 0 |
| HNZK|contig_4700 | 0 | 0.06 | 21 873.11 | 0.31 | 0.13 | 12.75 |
Fig. 1 Relative abundance of virus composition at the phylum, order, and family levels of Marmota himalayana from the agricultural region of eastern in Qinghai Province
Fig. 2 Shannon index for intestinal viruses of Marmota himalayana from the agricultural region of eastern in Qinghai Province. ns, no significant difference
Fig. 3 PCoA analysis among intestinal microbiomes of Marmota himalayana from the agricultural region of eastern in Qinghai Province. HD, Haidong City; XN, Xining City
Fig. 5 The top 10 host information at the family level for intestinal viruses in Marmota himalayana from the agricultural region of eastern in Qinghai Province.A: HDHZ; B: HDMH; C: HDPA; D: XNHY; E: XNHZ; F: XNDT
Fig. 6 The top 10 host information at the genus level for intestinal viruses in Marmota himalayana from the agricultural region of eastern in Qinghai Province. A: HDHZ; B: HDMH; C: HDPA; D: XNHY; E: XNHZ; F: XNDT
Fig. 7 The top 10 host information at the species level for intestinal viruses in Marmota himalayana from the agricultural region of eastern in Qinghai Province. A: HDHZ; B: HDMH; C: HDPA; D: XNHY; E: XNHZ; F: XNDT
| Bolger A M, Lohse M, Usadel B,2014. Trimmomatic:a flexible trimmer for Illumina Sequence Data[J]. Bioinformatics,30(15):2114‑2120. DOI:10.1093/bioinformatics/btu170 . | |
| Dai X, Shang G, Lu S, Yang Y, Xu J G,2018. A new subtype of eastern tick‑borne encephalitis virus is covered in Qinghai‑Xizang Plateau,China[J]. Emerg Microbes Infect,7(1):74. | |
| Dion M B, Oechslin F, Moineau S,2020. Phage diversity,genomics and phylogeny[J]. Nat Rev Microbiol,18(3):125‑138. DOI:10.1038/s41579-019-0311-5 . | |
| Drewes S, Straková P, Drexler J F, Jacob J, Ulrich R G,2017. Assessing the diversity of rodent borne viruse:exploring of high throughput sequencing and classical amplification /sequencing approaches[J]. Adv Virus Res,99:61‑108. | |
| Guo J, Bolduc B, Zayed A A, Varsani A, Dominguez‑Huerta G, Delmont T O, Pratama A A, Gazitúa M C, Dean Vik D, Sullivan M B,2021. VirSorter2:a multi‑classifier,expert‑guided approach to detect diverse DNA and RNA viruses[J]. Microbiome,9(1):37. DOI:10.1186/s40168-020-00990-y . | |
| He W T, Hou X, Zhao J, Sun J, He H, Si W, Wang J, Jiang Z W, Yan Z Q, Xing G, Lu M, Suchard M A, Ji X, Gong W J, He B, Li J, Lemey P, Guo D Y, Tu C C, Holmes E C, Shi M, Su S,2022. Virome characterization of game animals in China reveals a spectrum of emerging pathogens[J]. Cell,185(7):1117‑1129. DOI:10.1016/j.cell.2022.02.014 . | |
| Hess M, Sczyrba A, Egan R, Kim T W, Chokhawala H, Schroth G, Luo S, Clark D S, Chen F, Zhang T, Mackie R I, Pennacchio L A, Tringe S G, Visel A, Woyke T, Wang Z, Rubin E M,2011. Metagenomic discovery of biomass‑degrading genes and genomes from cow rumen[J]. Science,331(6016):463‑467. DOI:10.1126/science.1200387 . | |
| Hulo C, E de Castro, Masson P, Bougueleret L, Bairoch A, Xenarios I, Mercier P L,2011. Viral Zone:a knowledge resource to understand virus diversity[J]. Nucleic Acids Res,39:576‑582. DOI:10.1093/nar/gkq901 . | |
| Karesh W B, Dobson A, Lloyd‑Smith J O, Lubroth J, Dixon M A,2012. Ecology of zoonoses:natural and unnatural histories[J]. Lancet,380(9857):1936‑1945. DOI:10.1016/s0140-6736(12)61678-x . | |
| Li H, Durbin R,2009. Fast and accurate short read alignment with burrows‑wheeler transform[J]. Bioinformatics,25(14):1754‑1760. DOI:10.1093/bioinformatics/btp324 . | |
| Li D, Luo R, Liu C M, Leung C M, Ting H F, Sadakane K, Yamashita H, Lam T W,2016. MEGAHIT v1.0:A fast and scalable metagenome assembler driven by advanced methodologies and community practices[J]. Methods,102 :3‑11. DOI:10.1016/j.ymeth.2016.02.020 . | |
| McCormick J B, Webb P A, Krebs J W, Johnson K M, Smith E S,1987. A prospective study of the epidemiology and ecology of Lassa fever[J]. J Infect Dis,155(3):437‑444. | |
| Nayfach S, Camargo A P, Schulz F, Eloe‑Fadrosh E, Roux S, Kyrpides N C,2021. CheckV assesses the quality and completeness of metagenome‑assembled viral genomes[J]. Nature Biotechnology,39(5):578‑585. DOI:10.1038/S41587- 020-00774-7 . | |
| Nemeth N M, Root J J, Hartwig A E, Bowen R A, Bosco‑Lauth A M,2021. Powassan virus experimental infections in three wild mammal species[J]. Am J Trop Med Hyg,104(3):1048‑1054. DOI:10.4269/ajtmh.20-0105 . | |
| Pons J C, Paez‑Espino D, Riera G, Ivanova N, Kyrpides N C, Llabrés M,2021. VPF‑Class:Taxonomic assignment and host prediction of uncultivated viruses based on viral protein families[J]. Bioinformatics,37(13):1805‑1808. DOI:10.1093/bioinformatics/btab026 . | |
| Reyes A, Haynes M, Hanson N, Angly F E, Heath A C, Rohwer F, Gordon J I,2010. Viruses in the faecal microbiota of monozygotic twins and their mothers[J]. Nature,466(7304):334‑338. DOI:10.1038/nature09199 . | |
| Shang J, Sun Y,2022. CHERRY:a Computational metHod for accuratE pRediction of virus‑pRokarYotic interactions using a graph encoder‑decoder mode[J]. Briefings in Bioinformatics,23(5):1‑16. DOI:10.48550/arXiv.2201.01018 . | |
| Shi L, Yang G, Zhang Z, Xia L, Liang Y, Tan H, He J, Xu J, Song Z, Li W, Wang P,2018. Reemergence of human plague in Yunnan China in 2016[J]. PLoS ONE,13(6):e0198067. DOI:10.1371/journal.pone.0198067 . | |
| Shkoporov A N, Clooney A G, Sutton T D S, Ryan F J, Daly K M, Nolan J A, McDonnell S A, Khokhlova E V, Draper L A, Forde A, Guerin E, Velayudhan V, Ross R P, Hill C,2019. The human gut virome is highly diverse,stable,and individual specific[J]. Cell Host & Microbe,26(4):527‑541. DOI:10.1016/j.chom.2019.09.009 . | |
| Su Hang, Xu Chunyu, Jia Xiuqi, Yu Chengwen, Wang Shuyan, Jin Zhimin,2021. Advances and prospect for rodents carry viruses researches[J]. Biological Chemical Engineering,7(1):123‑126.(in Chinese with English abstract) | |
| Sun Yue, Tu Changchun, Zhu Guoqiang, He Biao,2021. Viral Metagenomic analyses of Himalayan Marmots in Kashgar,Xinjiang,China[J]. Chinese Journal of Virology,37(2):404‑414. DOI:10.13242/j.cnki.bingduxuebao.003894. (in Chinese with English abstract) | |
| Taylor L H, Latham S M, Woolhouse M E,2001. Risk factors for human disease emergence[J]. Philos Trans R SocLond B Biol Sci,356(1411):983‑989. DOI:10.2307/3066689 . | |
| Tu C, Crameri G, Kong X, Chen J, Sun Y, Yu M, Xiang H, Xia X, Liu S, Ren T, Yu Y, Eaton B T, Xuan H, Wang L F,2004. Antibodies to SARS coronavirus in civets[J]. Emerging Infectious Diseases,10(12):2244‑2248. DOI:10.3201/eid1012.040520 . | |
| Wang W L, Xu S Y, Ren Z G, Tao L, Jiang J W, Zheng S S,2015. Application of metagenomics in the human gut microbiome[J]. World J Gastroenterol,21(3):803‑814. | |
| Zhao L, Rosario K, Breitbart M, Duffy S,2019. Eukaryotic circular repencoding single‑stranded DNA(CRESS DNA)virses:Ubiquitous viruses with small genomes and a diverse host rang[J]. Adv Virus Res,103:71‑133. DOI:10.1016/bs.aivir. 2018.10.001 . | |
| Zhang Cuiyuan,2016. Application of high‑throughput sequencing for virus detection and analysis in wild animalfeces[D]. Changsha:Hunan Normal University.(in Chinese with English abstract) | |
| 孙悦,涂长春,朱国强,何彪,2021. 新疆喀什喜马拉雅旱獭病毒组研究[J]. 病毒学报,37(2):404‑414. DOI:10.13242/j.cnki.bingduxuebao.003894 . | |
| 苏航,徐春雨,贾修岐,于成文,王述艳,金志民,2021. 啮齿目动物携带病毒的研究进展与展望[J]. 生物化工,7(1):123‑126. | |
| 张翠媛,2016. 高通量测序技术应用于野生动物粪便中病毒的发现与分析[D]. 长沙:湖南师范大学. |
| [1] | HUANG Xiaolong, WANG Zhijun, WU Jiasong, LIU Lixian . THE BREEDING BIOLOGICAL CHARACTERISTICS OF Marmota himalayana IN RESHUITAN AND WULANNAOTAN, HAIYAN COUNTY, QINGHAI PROVINCE [J]. , 1986, 6(4): 307-311. |
| Viewed | ||||||
|
Full text |
|
|||||
|
Abstract |
|
|||||
青公网安备 63010402000199号
青ICP备05000010号-2