Revisiting the Evolutionary History of Pigs via De Novo Mutation Rate Estimation in A Three-generation Pedigree

Construction of the porcine genome mobile element variations and investigation of its role in population diversity and gene expression

BACKGROUND

Mobile element variants (MEVs) have a significant and complex impact on genomic diversity and phenotypic traits. However, the quantity, distribution, and relationship with gene expression and complex traits of MEVs in the pig genome remain poorly understood.

RESULTS

We constructed the most comprehensive porcine MEV library based on high-depth whole genome sequencing (WGS) data from 747 pigs across 59 breeds worldwide. This database identified a total of 147,993 polymorphic MEVs, including 121,099 short interspersed nuclear elements (SINEs), 26,053 long interspersed nuclear elements (LINEs), 802 long terminal repeats (LTRs), and 39 other transposons, among which 54% are newly discovered. We found that MEVs are unevenly distributed across the genome and are strongly influenced by negative selection effects. Importantly, we identified 514, 530, and 584 candidate MEVs associated with population differentiation, domestication, and breed formation, respectively. For example, a significantly differentiated MEV is located in the ATRX intron between Asian and European pigs, whereas ATRX is also differentially expressed between Asian and European pigs in muscle tissue. In addition, we identified 4,169 expressed MEVs (eMEVs) significantly associated with gene expression and 6,914 splicing MEVs (sMEVs) associated with gene splicing based on RNA-seq data from 266 porcine liver tissues. These eMEVs and sMEVs explain 6.24% and 9.47%, respectively, of the observed cis-heritability and highlight the important role of MEVs in the regulation of gene expression. Finally, we provide a high-quality SNP-MEV reference haplotype panel to impute MEV genotypes from genome-wide SNPs. Notably, we identified a candidate MEV significantly associated with total teat number, demonstrating the functionality of this reference panel.

CONCLUSIONS

The present investigation demonstrated the importance of MEVs in pigs in terms of population diversity, gene expression and phenotypic traits, which may provide useful resources and theoretical support for pig genetics and breeding.

Characterization and distribution of de novo mutations in the zebra finch

Germline de novo mutations (DNMs) provide the raw material for evolution. The DNM rate varies considerably between species, sexes and chromosomes. Here, we identify DNMs in the zebra finch (Taeniopygia guttata) across 16 parent-offspring trios using two genome assemblies of different quality. Using an independent genotyping assay, we validate 82% of the 150 candidate DNMs. DNM rates are consistent between both assemblies, with estimates of 6.14 × 10-9 and 6.36 × 10-9 per site per generation. We observe a strong paternal bias in DNM rates (male-to-female ratio ɑ ≈ 4), but this bias is in transition mutations only, leading to a transition-to-transversion ratio of 3.18 and 3.57. Finally, we find that DNMs tend to be randomly distributed across chromosomes, not associated with recombination hotspots or genic regions. However, the sex chromosome chrZ shows a roughly fourfold increased DNM rate compared to autosomes, which is more than the expected increase due to chrZ spending two-thirds of its time in males. Overall, our results further enhance our understanding of DNMs in passerine songbirds.

A comprehensive atlas of nuclear sequences of mitochondrial origin (NUMT) inserted into the pig genome

BACKGROUND

The integration of nuclear mitochondrial DNA (mtDNA) into the mammalian genomes is an ongoing, yet rare evolutionary process that produces nuclear sequences of mitochondrial origin (NUMT). In this study, we identified and analysed NUMT inserted into the pig (Sus scrofa) genome and in the genomes of a few other Suinae species. First, we constructed a comparative distribution map of NUMT in the Sscrofa11.1 reference genome and in 22 other assembled S. scrofa genomes (from Asian and European pig breeds and populations), as well as the assembled genomes of the Visayan warty pig (Sus cebifrons) and warthog (Phacochoerus africanus). We then analysed a total of 485 whole genome sequencing datasets, from different breeds, populations, or Sus species, to discover polymorphic NUMT (inserted/deleted in the pig genome). The insertion age was inferred based on the presence or absence of orthologous NUMT in the genomes of different species, taking into account their evolutionary divergence. Additionally, the age of the NUMT was calculated based on sequence degradation compared to the authentic mtDNA sequence. We also validated a selected set of representative NUMT via PCR amplification.

RESULTS

We have constructed an atlas of 418 NUMT regions, 70 of which were not present in any assembled genomes. We identified ancient NUMT regions (older than 55 million years ago, Mya) and NUMT that appeared at different time points along the Suinae evolutionary lineage. We identified very recent polymorphic NUMT (private to S. scrofa, with < 1 Mya), and more ancient polymorphic NUMT (3.5-10 Mya) present in various Sus species. These latest polymorphic NUMT regions, which segregate in European and Asian pig breeds and populations, are likely the results of interspecies admixture within the Sus genus.

CONCLUSIONS

This study provided a first comprehensive analysis of NUMT present in the Sus scrofa genome, comparing them to NUMT found in other species within the order Cetartiodactyla. The NUMT-based evolutionary window that we reconstructed from NUMT integration ages could be useful to better understand the micro-evolutionary events that shaped the modern pig genome and enriched the genetic diversity of this species.

Physiological basis for xenotransplantation from genetically modified pigs to humans

The collective efforts of scientists over multiple decades have led to advancements in molecular and cellular biology-based technologies including genetic engineering and animal cloning that are now being harnessed to enhance the suitability of pig organs for xenotransplantation into humans. Using organs sourced from pigs with multiple gene deletions and human transgene insertions, investigators have overcome formidable immunological and physiological barriers in pig-to-nonhuman primate (NHP) xenotransplantation and achieved prolonged pig xenograft survival. These studies informed the design of Revivicor's (Revivicor Inc, Blacksburg, VA) genetically engineered pigs with 10 genetic modifications (10 GE) (including the inactivation of 4 endogenous porcine genes and insertion of 6 human transgenes), whose hearts and kidneys have now been studied in preclinical human xenotransplantation models with brain-dead recipients. Additionally, the first two clinical cases of pig-to-human heart xenotransplantation were recently performed with hearts from this 10 GE pig at the University of Maryland. Although this review focuses on xenotransplantation of hearts and kidneys, multiple organs, tissues, and cell types from genetically engineered pigs will provide much-needed therapeutic interventions in the future.

Naturally occurring genetic diseases caused by de novo variants in domestic animals

With the advent of next-generation sequencing, an increasing number of cases of de novo variants in domestic animals have been reported in scientific literature primarily associated with clinically severe phenotypes. The emergence of new variants at each generation is a crucial aspect in understanding the pathology of early-onset diseases in animals and can provide valuable insights into similar diseases in humans. With the aim of collecting deleterious de novo variants in domestic animals, we searched the scientific literature and compiled reports on 42 de novo variants in 31 genes in domestic animals. No clear disease-associated phenotype has been established in humans for three of these genes (NUMB, ANKRD28 and KCNG1). For the remaining 28 genes, a strong similarity between animal and human phenotypes was recognized from available information in OMIM and OMIA, revealing the importance of comparative studies and supporting the use of domestic animals as natural models for human diseases, in line with the One Health approach.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Exploration of Genome-Wide Recombination Rate Variation Patterns at Different Scales in Pigs

Meiotic recombination is a prevalent process in eukaryotic sexual reproduction organisms that plays key roles in genetic diversity, breed selection, and species evolution. However, the recombination events differ across breeds and even within breeds. In this study, we initially computed large-scale population recombination rates of both sexes using approximately 52 K SNP genotypes in a total of 3279 pigs from four different Chinese and Western breeds. We then constructed a high-resolution historical recombination map using approximately 16 million SNPs from a sample of unrelated individuals. Comparative analysis of porcine recombination events from different breeds and at different resolutions revealed the following observations: Firstly, the 1Mb-scale pig recombination maps of the same sex are moderately conserved among different breeds, with the similarity of recombination events between Western pigs and Chinese indigenous pigs being lower than within their respective groups. Secondly, we identified 3861 recombination hotspots in the genome and observed medium- to high-level correlation between historical recombination rates (0.542~0.683) and estimates of meiotic recombination rates. Third, we observed that recombination hotspots are significantly far from the transcription start sites of pig genes, and the silico-predicted PRDM9 zinc finger domain DNA recognition motif is significantly enriched in the regions of recombination hotspots compared to recombination coldspots, highlighting the potential role of PRDM9 in regulating recombination hotspots in pigs. Our study analyzed the variation patterns of the pig recombination map at broad and fine scales, providing a valuable reference for genomic selection breeding and laying a crucial foundation for further understanding the molecular mechanisms of pig genome recombination.

The Genetic Selection of HSPD1 and HSPE1 Reduce Inflammation of Liver and Spleen While Restraining the Growth and Development of Skeletal Muscle in Wuzhishan Pigs

Wuzhishan (WZS) pigs, which are minipigs native to Hainan Province in China, are characterized by strong resistance to extreme hot temperatures and humidity. The relationship between their immune response and growth still needs to be clarified. In this study, we used whole genome sequencing (WGS) to detect variations within 37 WZS pigs, 32 Large White (LW) pigs, and 22 Xiangxi black (XXB) pigs, and ~2.49 GB of SNPs were obtained. These data were combined with those of two other pig breeds, and it was found that most of the genes detected (354) were located within the distinct genetic regions between WZS pigs and LW pigs. The network that was constructed using these genes represented a center including 12 hub genes, five of which had structural variations (SVs) within their regulatory regions. Furthermore, RNA-seq and RT-qPCR data for 12 genes were primarily consistent in liver, spleen, and LDM tissues. Notably, the expression of HSPs (HSPD1 and HSPE1) was higher while that of most genes involved in the JAK3-STAT pathway were lower in liver tissue of WZS pigs, compared with LW pigs. This likely not only reduced inflammation-related immune response but also impaired their growth. Our findings demonstrated the role of HSPs in the connection between inflammation and growth rate, while also providing the fundamental genetic selection of the adaptability of WZS pigs.

Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences

High-quality whole-genome resequencing in large-scale pig populations with pedigree structure and multiple breeds would enable accurate construction of haplotype and robust selection-signature detection. Here, we sequence 740 pigs, combine with 149 of our previously published resequencing data, retrieve 207 resequencing datasets, and form a panel of worldwide distributed wild boars, aboriginal and highly selected pigs with pedigree structures, amounting to 1096 genomes from 43 breeds. Combining with their haplotype-informative reads and pedigree structure, we accurately construct a panel of 1874 haploid genomes with 41,964,356 genetic variants. We further demonstrate its valuable applications in GWAS by identifying five novel loci for intramuscular fat content, and in genomic selection by increasing the accuracy of estimated breeding value by 36.7%. In evolutionary selection, we detect MUC13 gene under a long-term balancing selection, as well as NPR3 gene under positive selection for pig stature. Our study provides abundant genomic variations for robust selection-signature detection and accurate haplotypes for deciphering complex traits in pigs.

X-linked hydrocephalus genes: Their proximity to telomeres and high A + T content compared to Parkinson's disease

Proximity to telomeres (i) and high adenine and thymine (A + T) content (ii) are two factors associated with high mutation rates in human chromosomes. We have previously shown that >100 human genes when mutated to cause congenital hydrocephalus (CH) meet either factor (i) or (ii) at 91% matching, while two factors are poorly satisfied in human genes associated with familial Parkinson's disease (fPD) at 59%. Using the sets of mouse, rat, and human chromosomes, we found that 7 genes associated with CH were located on the X chromosome of mice, rats, and humans. However, genes associated with fPD were in different autosomes depending on species. While the contribution of proximity to telomeres in the autosome was comparable in CH and fPD, high A + T content played a pivotal contribution in X-linked CH (43% in all three species) than in fPD (6% in rodents or 13% in humans). Low A + T content found in fPD cases suggests that PARK family genes harbor roughly 3 times higher chances of methylations in CpG sites or epigenetic changes than X-linked genes.

Warthog Genomes Resolve an Evolutionary Conundrum and Reveal Introgression of Disease Resistance Genes

African wild pigs have a contentious evolutionary and biogeographic history. Until recently, desert warthog (Phacochoerus aethiopicus) and common warthog (P. africanus) were considered a single species. Molecular evidence surprisingly suggested they diverged at least 4.4 million years ago, and possibly outside of Africa. We sequenced the first whole-genomes of four desert warthogs and 35 common warthogs from throughout their range. We show that these two species diverged much later than previously estimated, 400,000–1,700,000 years ago depending on assumptions of gene flow. This brings it into agreement with the paleontological record. We found that the common warthog originated in western Africa and subsequently colonized eastern and southern Africa. During this range expansion, the common warthog interbred with the desert warthog, presumably in eastern Africa, underlining this region’s importance in African biogeography. We found that immune system–related genes may have adaptively introgressed into common warthogs, indicating that resistance to novel diseases was one of the most potent drivers of evolution as common warthogs expanded their range. Hence, we solve some of the key controversies surrounding warthog evolution and reveal a complex evolutionary history involving range expansion, introgression, and adaptation to new diseases.