Genome-wide identification and comparative analysis of endogenous retroviruses across five breed-representative pig assemblies.
5 affiliations
- 1College of Animal Science and Technology, Yunnan Agricultural University, Kunming, Yunnan, China.
- 2Honghe Prefecture Feed Inspection Station, Mengzi, Yunnan, China.
- 3CP Kunming Co., Ltd, Kunming, Yunnan, China.
- 4College of Animal Science and Technology, Yunnan Agricultural University, Kunming, Yunnan, China. [email protected].
- 5College of Veterinary Medicine, Yunnan Agricultural University, Kunming, Yunnan, China. [email protected].
Abstract
Background
Porcine xenotransplantation offers a promising solution to human organ shortages, yet cross-species transmission of endogenous retroviruses (ERVs) remains a major safety barrier. Current risk assessments predominantly rely on a single reference genome, creating a distinct blind spot regarding ERV diversity and structural variation across diverse genetic backgrounds. We profiled ERV loci across five chromosome-level pig genome assemblies (Bama miniature pig, Duroc, Meishan, Ningxiang, and USMARC), characterizing structural completeness, phylogenetic relationships, relative insertion ages, and host-gene/QTL genomic context. ERVs were identified with the LTRharvest/LTRdigest pipeline, annotated by Pfam domain screening, and classified as complete, near-complete, or partial. Phylogenetic analysis used maximum-likelihood inference, insertion age was estimated from paired long terminal repeat (LTR) divergence, and ERV-proximal genes within 10 kb were analyzed for Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment and tissue expression.Results
Across the five chromosome-level assemblies, the paired-LTR-anchored catalog comprised 23,767 ERV candidates in the main autosome/chrX analysis set, including 75 structurally complete, 322 near-complete, and 23,370 partial loci. An additional 86 chrY records were detected in assemblies with chrY and were all classified as partial. Complete-ERV insertion-age profiles showed assembly-level variation, although Kruskal-Wallis tests did not detect a significant overall difference among assemblies. The five-assembly shared ERV-proximal host-gene set contained 278 genes at the 10 kb window; this set decreased to 229 genes at 5 kb and increased to 358 and 618 genes at 20 and 50 kb, respectively. Chromosome-stratified randomization showed that the shared 10 kb gene set overlapped pig QTL intervals more often than expected by chance (56.47% observed versus 43.23% random mean; empirical P < 0.0001).Conclusions
Across five breed-representative pig assemblies, we identified 23,767 paired-LTR ERV candidates, of which only 75 (0.32%) met the structural-completeness criteria, confirming that the detected porcine ERV catalog is dominated by degraded elements. Complete loci showed assembly-level variation in retroviral affinities and relative insertion-age profiles, while the five-assembly shared 10 kb ERV-proximal gene set showed non-random QTL colocalization (56.47% observed versus 43.23% random mean; empirical P < 0.0001). Benchmarking against published porcine endogenous retrovirus (PERV) and retroviral insertion polymorphism (RIP) resources linked 9 complete loci to reconstructed Sus scrofa ERV6A/B full-length ERVs (SscERV6A/B Fl-ERVs) and left 66 complete loci outside that focused set. This catalog provides concrete ERV loci and host-gene contexts for donor-line evaluation and future population-scale, T2T-informed studies of porcine ERV diversity.
Related
- Structural Insights into Endogenous Retroviruses: Evolution and Co-Opted Functions.
- HIF1A recruits primate-specific endogenous retroviruses into the human hypoxic and immune responses
- Purification-free cryo-EM reveals endogenous assemblies from minute tardigrade samples
- Identification of Synchronously Active Cells (ISAC) to study the activity of cell assemblies
- A comprehensive comparative genomic analysis of the endogenous Florendovirus-1 in Prunus
BibTeX
@article{zhao2026genomewide,
title = {{Genome-wide identification and comparative analysis of endogenous retroviruses across five breed-representative pig assemblies.}},
author = {Xiaotong Zhao and Yongping Hu and Zhidong Zhang and Xun Xiang and Hua Chang},
journal = {BMC genomics},
year = {2026},
month = jul,
doi = {10.1186/s12864-026-13211-x},
url = {https://europepmc.org/article/MED/42897448},
note = {CC BY-NC-ND},
}
RIS
TY - JOUR
TI - Genome-wide identification and comparative analysis of endogenous retroviruses across five breed-representative pig assemblies.
AU - Xiaotong Zhao
AU - Yongping Hu
AU - Zhidong Zhang
AU - Xun Xiang
AU - Hua Chang
PY - 2026
DA - 2026/07/01
JO - BMC genomics
DO - 10.1186/s12864-026-13211-x
UR - https://europepmc.org/article/MED/42897448
AB - Background Porcine xenotransplantation offers a promising solution to human organ shortages, yet cross-species transmission of endogenous retroviruses (ERVs) remains a major safety barrier. Current risk assessments predominantly rely on a single reference genome, creating a distinct blind spot regarding ERV diversity and structural variation across diverse genetic backgrounds. We profiled ERV loci across five chromosome-level pig genome assemblies (Bama miniature pig, Duroc, Meishan, Ningxiang, and USMARC), characterizing structural completeness, phylogenetic relationships, relative insertion ages, and host-gene/QTL genomic context. ERVs were identified with the LTRharvest/LTRdigest pipeline, annotated by Pfam domain screening, and classified as complete, near-complete, or partial. Phylogenetic analysis used maximum-likelihood inference, insertion age was estimated from paired long terminal repeat (LTR) divergence, and ERV-proximal genes within 10 kb were analyzed for Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment and tissue expression.Results Across the five chromosome-level assemblies, the paired-LTR-anchored catalog comprised 23,767 ERV candidates in the main autosome/chrX analysis set, including 75 structurally complete, 322 near-complete, and 23,370 partial loci. An additional 86 chrY records were detected in assemblies with chrY and were all classified as partial. Complete-ERV insertion-age profiles showed assembly-level variation, although Kruskal-Wallis tests did not detect a significant overall difference among assemblies. The five-assembly shared ERV-proximal host-gene set contained 278 genes at the 10 kb window; this set decreased to 229 genes at 5 kb and increased to 358 and 618 genes at 20 and 50 kb, respectively. Chromosome-stratified randomization showed that the shared 10 kb gene set overlapped pig QTL intervals more often than expected by chance (56.47% observed versus 43.23% random mean; empirical P < 0.0001).Conclusions Across five breed-representative pig assemblies, we identified 23,767 paired-LTR ERV candidates, of which only 75 (0.32%) met the structural-completeness criteria, confirming that the detected porcine ERV catalog is dominated by degraded elements. Complete loci showed assembly-level variation in retroviral affinities and relative insertion-age profiles, while the five-assembly shared 10 kb ERV-proximal gene set showed non-random QTL colocalization (56.47% observed versus 43.23% random mean; empirical P < 0.0001). Benchmarking against published porcine endogenous retrovirus (PERV) and retroviral insertion polymorphism (RIP) resources linked 9 complete loci to reconstructed Sus scrofa ERV6A/B full-length ERVs (SscERV6A/B Fl-ERVs) and left 66 complete loci outside that focused set. This catalog provides concrete ERV loci and host-gene contexts for donor-line evaluation and future population-scale, T2T-informed studies of porcine ERV diversity.
SN - 1471-2164
ER -