African Suid Genomes Provide Insights into the Local Adaptation to Diverse African Environments.

African wild suids consist of several endemic species that represent ancient members of the family Suidae and have colonized diverse habitats on the African continent. However, limited genomic resources for African wild suids hinder our understanding of their evolution and genetic diversity. In this study, we assembled high-quality genomes of a common warthog (Phacochoerus africanus), a red river hog (Potamochoerus porcus), as well as an East Asian Diannan small-ear pig (Sus scrofa). Phylogenetic analysis showed that common warthog and red river hog diverged from their common ancestor around the Miocene/Pliocene boundary, putatively predating their entry into Africa. We detected species-specific selective signals associated with sensory perception and interferon signaling pathways in common warthog and red river hog, respectively, which contributed to their local adaptation to savannah and tropical rainforest environments, respectively. The structural variation and evolving signals in genes involved in T-cell immunity, viral infection, and lymphoid development were identified in their ancestral lineage. Our results provide new insights into the evolutionary histories and divergent genetic adaptations of African suids.

Genetic Architecture Underlying Nascent Speciation-The Evolution of Eurasian Pigs under Domestication.

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between "good species." Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky-Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.

Whole genome resequencing reveals an association of ABCC4 variants with preaxial polydactyly in pigs.

BACKGROUND: Polydactyly is one of the most common congenital limb dysplasia in many animal species. Although preaxial polydactyly (PPD) has been comprehensively studied in humans as a common abnormality, the genetic variations in other animal species have not been fully understood. Herein, we focused on the pig, as an even-toed ungulate mammal model with its unique advantages in medical and genetic researches, two PPD families consisting of four affected and 20 normal individuals were sequenced. RESULTS: Our results showed that the PPD in the sampled pigs were not related to previously reported variants. A strong association was identified at ABCC4 and it encodes a transmembrane protein involved in ciliogenesis. We found that the affected and normal individuals were highly differentiated at ABCC4, and all the PPD individuals shared long haplotype stretches as compared with the unaffected individuals. A highly differentiated missense mutation (I85T) in ABCC4 was observed at a residue from a transmembrane domain highly conserved among a variety of organisms. CONCLUSIONS: This study reports ABCC4 as a new candidate gene and identifies a missense mutation for PPD in pigs. Our results illustrate a putative role of ciliogenesis process in PPD, coinciding with an earlier observation of ciliogenesis abnormality resulting in pseudo-thumb development in pandas. These results expand our knowledge on the genetic variations underlying PPD in animals.

Analysis of the genetic variation in mitochondrial DNA, Y-chromosome sequences, and MC1R sheds light on the ancestry of Nigerian indigenous pigs.

BACKGROUND: The history of pig populations in Africa remains controversial due to insufficient evidence from archaeological and genetic data. Previously, a Western ancestry for West African pigs was reported based on loci that are involved in the determination of coat color. We investigated the genetic diversity of Nigerian indigenous pigs (NIP) by simultaneously analyzing variation in mitochondrial DNA (mtDNA), Y-chromosome sequence and the melanocortin receptor 1 (MC1R) gene. RESULTS: Median-joining network analysis of mtDNA D-loop sequences from 201 NIP and previously characterized loci clustered NIP with populations from the West (Europe/North Africa) and East/Southeast Asia. Analysis of partial sequences of the Y-chromosome in 57 Nigerian boars clustered NIP into lineage HY1. Finally, analysis of MC1R in 90 NIP resulted in seven haplotypes, among which the European wild boar haplotype was carried by one individual and the European dominant black by most of the other individuals (93%). The five remaining unique haplotypes differed by a single synonymous substitution from European wild type, European dominant black and Asian dominant black haplotypes. CONCLUSIONS: Our results demonstrate a European and East/Southeast Asian ancestry for NIP. Analyses of MC1R provide further evidence. Additional genetic analyses and archaeological studies may provide further insights into the history of African pig breeds. Our findings provide a valuable resource for future studies on whole-genome analyses of African pigs.

A human-specific insertion promotes cell proliferation and migration by enhancing TBC1D8B expression.

Human-specific insertions play important roles in human phenotypes and diseases. Here we reported a 446-bp insertion (Insert-446) in intron 11 of the TBC1D8B gene, located on chromosome X, and traced its origin to a portion of intron 6 of the EBF1 gene on chromosome 5. Interestingly, Insert-446 was present in the human Neanderthal and Denisovans genomes, and was fixed in humans after human-chimpanzee divergence. We have demonstrated that Insert-446 acts as an enhancer through binding transcript factors that promotes a higher expression of human TBC1D8B gene as compared with orthologs in macaques. In addition, over-expression TBC1D8B promoted cell proliferation and migration through "a dual finger" catalytic mechanism (Arg538 and Gln573) in the TBC domain in vitro and knockdown of TBC1D8B attenuated tumorigenesis in vivo. Knockout of Insert-446 prevented cell proliferation and migration in cancer and normal cells. Our results reveal that the human-specific Insert-446 promotes cell proliferation and migration by upregulating the expression of TBC1D8B gene. These findings provide a significant insight into the effects of human-specific insertions on evolution.

Target capture sequencing for the first Nigerian genotype I ASFV genome.

African swine fever (ASF) is a contagious viral disease that affects domestic pigs and wild boars, causing significant economic losses globally. After the first Nigerian outbreak in 1997, there have been frequent reports of ASF in pig-producing regions in the country. To facilitate control, it is important to understand the genotype and phylogenetic relationship of ASF viruses (ASFVs). Recent genetic analysis of Nigerian ASFV isolates has revealed the presence of both genotypes I and II; this is based on analysis of a few selected genes. Phylogenetic analysis of ASFV whole genomes highlights virus origins and evolution in greater depth. However, there is currently no information on the ASFV genome from Nigerian isolates. Two ASFV-positive samples were detected during a random survey of 150 Nigerian indigenous pig samples collected in 2016. We assembled near-complete genomes of the two ASFV-positive samples using in-solution hybrid capture sequencing. The genome-wide phylogenetic tree assigned these two genomes into p72 genotype I, particularly close to the virulent Benin 97/1 strain. The two ASFVs share 99.94 and 99.92 % genomic sequence identity to Benin97/1. This provides insight into the origin and relationship of ASFV strains from Nigeria and Italy. The study reports for the first time the determination of near-complete genomes of ASFV using in-solution hybrid capture sequencing, which represents an important advance in understanding the global evolutionary landscape of ASFVs.

CaptureProbe: a java tool for designing probes for capture Hi-C applications.

Many functional elements associated with traits and diseases are located in non-coding regions and act on distant target genes via chromatin looping and folding, making it difficult for scientists to reveal the genetic regulatory mechanisms. Capture Hi-C is a newly developed chromosome conformation capture technology based on hybridization capture between probes and target genomic regions. It can identify interactions among target loci and all other loci in a genome with low cost and high resolution. Here, we developed CaptureProbe, a user-friendly, graphical Java tool for the design of capture probes across a range of target sites or regions. Numerous parameters helped to achieve and optimize the designed probes. Design testing of CaptureProbe showed high efficiency in the design success ratio of target loci and probe specificity. Hence, this program will help scientists conduct genome spatial interaction research. CaptureProbe and source code are available at https://sourceforge.net/projects/captureprobe/.

OrthReg: a tool to predict cis-regulatory elements based on cross-species orthologous sequence conservation.

Cis-regulatory elements play an important role in the development of traits and disease in organisms (Ma et al., 2020; Woolfe et al., 2005) and their annotation could facilitate genetic studies. The Encyclopedia of DNA Elements (ENCODE) (Davis et al., 2018) and Functional Annotation of Animal Genomes (FAANG) (FAANG Consortium et al., 2015) offer pioneering data on regulatory elements in several species. Currently, however, regulatory element annotation data remain limited for most organisms. In this study, we developed a tool (OrthReg) for annotating conserved orthologous cis-regulatory elements in targeted genomes using an annotated reference genome. Cross-species validation of this annotation tool using human and mouse ENCODE data confirmed the robustness of this strategy. To explore the efficiency of the tool, we annotated the pig genome and identified more than 28 million regulatory annotation records using the reference human ENCODE data. With this regulatory annotation, some putative regulatory non-coding variants were identified within domestication sweeps in European and East Asian pigs. Thus, this tool can utilize data produced by ENCODE, FAANG, and similar projects, and can be easily extended to customized experimental data. The extensive application of this tool will help to identify informative single nucleotide polymorphisms (SNPs) in post-genome-wide association studies and resequencing analysis of organisms with limited regulatory annotation data.

Convergent genomic signatures of high-altitude adaptation among domestic mammals.

Abundant and diverse domestic mammals living on the Tibetan Plateau provide useful materials for investigating adaptive evolution and genetic convergence. Here, we used 327 genomes from horses, sheep, goats, cattle, pigs and dogs living at both high and low altitudes, including 73 genomes generated for this study, to disentangle the genetic mechanisms underlying local adaptation of domestic mammals. Although molecular convergence is comparatively rare at the DNA sequence level, we found convergent signature of positive selection at the gene level, particularly the EPAS1 gene in these Tibetan domestic mammals. We also reported a potential function in response to hypoxia for the gene C10orf67, which underwent positive selection in three of the domestic mammals. Our data provide an insight into adaptive evolution of high-altitude domestic mammals, and should facilitate the search for additional novel genes involved in the hypoxia response pathway.

Maternal genomic variability of the wild boar (Sus scrofa) reveals the uniqueness of East-Caucasian and Central Italian populations.

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East-Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.

Population Genomics Analysis Revealed Origin and High-altitude Adaptation of Tibetan Pigs.

Tibetan pig is native to the Qinghai-Tibet Plateau and has adapted to the high-altitude environmental condition such as hypoxia. However, its origin and genetic mechanisms underlying high-altitude adaptation still remain controversial and enigmatic. Herein, we analyze 229 genomes of wild and domestic pigs from Eurasia, including 63 Tibetan pigs, and detect 49.6 million high-quality variants. Phylogenomic and structure analyses show that Tibetan pigs have a close relationship with low-land domestic pigs in China, implying a common domestication origin. Positively selected genes in Tibetan pigs involved in high-altitude physiology, such as hypoxia, cardiovascular systems, UV damage, DNA repair. Three of loci with strong signals of selection are associated with EPAS1, CYP4F2, and THSD7A genes, related to hypoxia and circulation. We validated four non-coding mutations nearby EPAS1 and CYP4F2 showing reduced transcriptional activity in Tibetan pigs. A high-frequency missense mutation is found in THSD7A (Lys561Arg) in Tibetan pigs. The selective sweeps in Tibetan pigs was found in association with selection against non-coding variants, indicating an important role of regulatory mutations in Tibetan pig evolution. This study is important in understanding the evolution of Tibetan pigs and advancing our knowledge on animal adaptation to high-altitude environments.

Evolution of conserved secondary structures and their function in transcriptional regulation networks.

BACKGROUND: Many conserved secondary structures have been identified within conserved elements in the human genome, but only a small fraction of them are known to be functional RNAs. The evolutionary variations of these conserved secondary structures in human populations and their biological functions have not been fully studied. RESULTS: We searched for polymorphisms within conserved secondary structures and identified a number of SNPs within these elements even though they are highly conserved among species. The density of SNPs in conserved secondary structures is about 65% of that of their flanking, non-conserved, sequences. Classification of sites as stems or as loops/bulges revealed that the density of SNPs in stems is about 62% of that found in loops/bulges. Analysis of derived allele frequency data indicates that sites in stems are under stronger evolutionary constraint than sites in loops/bulges. Intergenic conserved secondary structures tend to associate with transcription factor-encoding genes with genetic distance being the measure of regulator-gene associations. A substantial fraction of intergenic conserved secondary structures overlap characterized binding sites for multiple transcription factors. CONCLUSION: Strong purifying selection implies that secondary structures are probably important carriers of biological functions for conserved sequences. The overlap between intergenic conserved secondary structures and transcription factor binding sites further suggests that intergenic conserved secondary structures have essential roles in directing gene expression in transcriptional regulation networks.

Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs.

The Bama Xiang pig (BMX) is a famous early-maturing Chinese indigenous breed with a two-end black coat. To uncover the genetic basis of the BMX phenotype, we conducted comparative genomic analyses between BMX and East Asian wild boars and Laiwu pigs, respectively. Genes under positive selection were enriched in pathways associated with gonadal hormone and melanin synthesis, consistent with the phenotypic changes observed during development in BMX pigs. We also performed differentially expressed gene analysis based on RNA-seq data from pituitary tissues of BMX and Large White pigs. The CTTNBP2NL, FRS2, KANK4, and KATNAL1 genes were under selection and exhibited expressional changes in the pituitary tissue, which may affect BMX pig puberty. Our study demonstrated the positive selection of early maturity in the development of BMX pigs and advances our knowledge on the role of regulatory elements in puberty evolution in pigs.

PigVar: a database of pig variations and positive selection signatures.

Database URL: http://www.ibiomedical.net/pigvar/.

Genetic variations associated with six-white-point coat pigmentation in Diannan small-ear pigs.

A common phenotypic difference among domestic animals is variation in coat color. Six-white-point is a pigmentation pattern observed in varying pig breeds, which seems to have evolved through several different mechanistic pathways. Herein, we re-sequenced whole genomes of 31 Diannan small-ear pigs from China and found that the six-white-point coat color in Diannan small-ear pigs is likely regulated by polygenic loci, rather than by the MC1R locus. Strong associations were observed at three loci (EDNRB, CNTLN, and PINK1), which explain about 20 percent of the total coat color variance in the Diannan small-ear pigs. We found a mutation that is highly differentiated between six-white-point and black Diannan small-ear pigs, which is located in a conserved noncoding sequence upstream of the EDNRB gene and is a putative binding site of the CEBPB protein. This study advances our understanding of coat color evolution in Diannan small-ear pigs and expands our traditional knowledge of coat color being a monogenic trait.

DNA methylation signatures of long intergenic noncoding RNAs in porcine adipose and muscle tissues.

Long intergenic noncoding RNAs (lincRNAs) are one of the major unexplored components of genomes. Here we re-analyzed a published methylated DNA immunoprecipitation sequencing (MeDIP-seq) dataset to characterize the DNA methylation pattern of pig lincRNA genes in adipose and muscle tissues. Our study showed that the methylation level of lincRNA genes was higher than that of mRNA genes, with similar trends observed in comparisons of the promoter, exon or intron regions. Different methylation pattern were observed across the transcription start sites (TSS) of lincRNA and protein-coding genes. Furthermore, an overlap was observed between many lincRNA genes and differentially methylated regions (DMRs) identified among different breeds of pigs, which show different fat contents, sexes and anatomic locations of tissues. We identify a lincRNA gene, linc-sscg3623, that displayed differential methylation levels in backfat between Min and Large White pigs at 60 and 120 days of age. We found that a demethylation process occurred between days 150 and 180 in the Min and Large White pigs, which was followed by remethylation between days 180 and 210. These results contribute to our understanding of the domestication of domestic animals and identify lincRNA genes involved in adipogenesis and muscle development.

Domestication genomics: evidence from animals.

Animal domestication has far-reaching significance for human society. The sequenced genomes of domesticated animals provide critical resources for understanding the genetic basis of domestication. Various genomic analyses have shed a new light on the mechanism of artificial selection and have allowed the mapping of genes involved in important domestication traits. Here, we summarize the published genomes of domesticated animals that have been generated over the past decade, as well as their origins, from a phylogenomic point of view. This review provides a general description of the genomic features encountered under a two-stage domestication process. We also introduce recent findings for domestication traits based on results from genome-wide association studies and selective-sweep scans for artificially selected genomic regions. Particular attention is paid to issues relating to the costs of domestication and the convergent evolution of genes between domesticated animals and humans.