Genome-wide association study for growth and fatness traits in Chinese Sujiang pigs

Potential selection for lipid kinase activity and spermatogenesis in Henan native pig breeds and growth shaping by introgression of European genes

BACKGROUND

China has one third of the worldwide indigenous pig breeds. The Henan province is one of the earliest pig domestication centers of China (about 8000 years ago). However, the precise genetic characteristics of the Henan local pig breeds are still obscure. To understand the origin and the effects of selection on these breeds, we performed various analyses on lineage composition, genetic structure, and detection of selection sweeps and introgression in three of these breeds (Queshan, Nanyang and Huainan) using genotyping data on 125 Queshan, 75 Nanyang, 16 Huainan pigs and 878 individuals from 43 Eurasian pig breeds.

RESULTS

We found no clear evidence of ancestral domestic pig DNA lineage in the Henan local breeds, which have an extremely complicated genetic background. Not only do they share genes with some northern Chinese pig breeds, such as Erhualian, Hetaodaer, and Laiwu, but they also have a high admixture of genes from foreign pig breeds (33-40%). Two striking selection sweeps in small regions of chromosomes 2 and 14 common to the Queshan and Nanyang breeds were identified. The most significant enrichment was for lipid kinase activity (GO:0043550) with the genes FII, AMBRA1, and PIK3IP1. Another interesting 636.35-kb region on chromosome 14 contained a cluster of spermatogenesis genes (OSBP2, GAL3ST1, PLA2G3, LIMK2, and PATZ1), a bisexual sterility gene MORC2, and a fat deposition gene SELENOM. Reproduction and growth genes LRP4, FII, and ARHGAP1 were present in a 238.05-kb region on SSC2 under selection. We also identified five loci associated with body length (P = 0.004) on chromosomes 1 and 12 that were introgressed from foreign pig breeds into the Henan breeds. In addition, the Chinese indigenous pig breeds fell into four main types instead of the previously reported six, among which the Eastern type could be divided into two subgroups.

CONCLUSIONS

Admixture of North China, East China and foreign pigs contributed to high genetic diversity of Henan local pigs. Ontology terms associated with lipid kinase activity and spermatogenesis and growth shaping by introgression of European genes in Henan pigs were identified through selective sweep analyses.

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.

Genome-Wide Association Analysis Identifies Genomic Regions and Candidate Genes for Growth and Fatness Traits in Diannan Small-Ear (DSE) Pigs

In the livestock industry, the growth and fatness traits are directly related to production efficiency and economic profits. As for Diannan small-ear (DSE) pigs, a unique indigenous breed, the genetic architecture of growth and fatness traits is still elusive. The aim of this study was to search the genetic loci and candidate genes associated with phenotypic traits in DSE pigs using GWAS based on the Geneseek Porcine 50K SNP Chip data. A total of 22,146 single nucleotide polymorphisms (SNPs) were detected in 265 DSE pigs and used for Genome-wide association studies (GWAS) analysis. Seven SNPs were found to be associated with back height, chest circumference, cannon bone circumference, and backfat thickness at the suggestive significance level. Based on gene annotation results, these seven SNPs were, respectively, mapped to the following candidate genes, VIPR2, SLC10A2, NUCKS1, MCT1, CHCHD3, SMOX, and GPR1, which are mainly involved with adipocyte differentiation, lipid metabolism, skeletal muscle development, and average daily weight gain. Our work offers novel insights into the genetic architecture of economically important traits in swine and may play an important role in breeding using molecular markers in the DSE breed.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle

Locating the genetic variation of important livestock and poultry economic traits is essential for genetic improvement in breeding programs. Identifying the candidate genes for the productive ability of Huaxi cattle was one crucial element for practical breeding. Based on the genotype and phenotype data of 1,478 individuals and the RNA-seq data of 120 individuals contained in 1,478 individuals, we implemented genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and Fisher's combined test (FCT) to identify the candidate genes for the carcass trait, the weight of longissimus dorsi muscle (LDM). The results indicated that GWAS, TWAS, and FCT identified seven candidate genes for LDM altogether: PENK was located by GWAS and FCT, PPAT was located by TWAS and FCT, and XKR4, MTMR3, FGFRL1, DHRS4, and LAP3 were only located by one of the methods. After functional analysis of these candidate genes and referring to the reported studies, we found that they were mainly functional in the progress of the development of the body and the growth of muscle cells. Combining advanced breeding techniques such as gene editing with our study will significantly accelerate the genetic improvement for the future breeding of Huaxi cattle.

Indel mutations of sheep PLAG1 gene and their associations with growth traits

Pleiomorphic adenoma gene 1 (PLAG1) is mainly expressed in embryonic development, and it is reported to take an effect on the growth performance of mice, cattle, pigs, and sheep. To explore how conservative the PLAG1 is in different sheep breeds, the effects of the two indel variants on the growth traits of the Chinese Luxi blackhead (LXBH) sheep were firstly detected. The P2-del 30 bp and P4-del 45 bp indel loci of the sheep PLAG1 gene were significantly related to 15 growth traits (P < 0.05). Genotype ID and genotype II were dominant for the P2-del 30 bp and P4-del 45 bp loci, respectively. The above findings indicated that the two indel mutations in the ovine PLAG1 gene were suggested to become the molecular markers for the selection of economic traits in sheep.

PHARP: a pig haplotype reference panel for genotype imputation

Pigs not only function as a major meat source worldwide but also are commonly used as an animal model for studying human complex traits. A large haplotype reference panel has been used to facilitate efficient phasing and imputation of relatively sparse genome-wide microarray chips and low-coverage sequencing data. Using the imputed genotypes in the downstream analysis, such as GWASs, TWASs, eQTL mapping and genomic prediction (GS), is beneficial for obtaining novel findings. However, currently, there is still a lack of publicly available and high-quality pig reference panels with large sample sizes and high diversity, which greatly limits the application of genotype imputation in pigs. In response, we built the pig Haplotype Reference Panel (PHARP) database. PHARP provides a reference panel of 2012 pig haplotypes at 34 million SNPs constructed using whole-genome sequence data from more than 49 studies of 71 pig breeds. It also provides Web-based analytical tools that allow researchers to carry out phasing and imputation consistently and efficiently. PHARP is freely accessible at http://alphaindex.zju.edu.cn/PHARP/index.php . We demonstrate its applicability for pig commercial 50 K SNP arrays, by accurately imputing 2.6 billion genotypes at a concordance rate value of 0.971 in 81 Large White pigs (~ 17 × sequencing coverage). We also applied our reference panel to impute the low-density SNP chip into the high-density data for three GWASs and found novel significantly associated SNPs that might be casual variants.

An imputation-based genome-wide association study for growth and fatness traits in Sujiang pigs

Sujiang pigs are a synthetic breed derived from Jiangquhai, Fengjing, and Duroc pigs. In this study, we sequenced the genome of 62 pigs with a coverage depth of 10× to 20×, including 27 Sujiang and 35 founder breed pigs, and we collected 360 global pigs' genome sequence data from public databases including 39 Duroc pigs. We obtained a high-quality variant dataset of 365 Sujiang pigs by imputing the porcine 80 K single nucleotide polymorphism (SNP) Beadchip to the whole-genome scale with a total of 422 pigs as a reference panel. A dataset of 365 imputated Sujiang pigs was used to perform single-trait genome-wide association study (GWAS) and meta-analyses for growth and fatness traits. Single-trait GWAS identified 1 907, 18, and 14 SNPs surpassing the suggestively significant threshold for backfat thickness, chest circumference, and chest width, respectively. Meta-analyses identified 2 400 genome-wide significant SNPs and 520 suggestively significant SNPs for backfat thickness and chest circumference, and 719 genome-wide significant SNPs and 1 225 suggestively significant SNPs for all seven traits. According to the meta-analysis of backfat thickness and chest circumference, a remarkable region of 2.69 Mb on Sus scrofa chromosome 4 containing FAM110B, IMPAD1, LYN, MOS, PENK, PLAG1, SDR16C5 and XKR4 was identified as a candidate region. The haplotype heat map of the 2.69 Mb region verified that Sujiang pigs were derived from Duroc and Chinese indigenous pigs, especially Jiangquhai pigs. The Kruskal-Wallis test showed that haplotypes of the 2.69 Mb region significantly affected backfat thickness and chest circumference traits. We then focused on PLAG1, an important growth-related gene, and identified two synonymous SNPs with obvious differences among different breeds in the PLAG1 gene. We then performed genotyping of 365 Sujiang, 150 Duroc, 95 Jiangquhai, and 100 Fengjing pigs to confirm the above result and verified that the two variants significantly affected phenotypes of growth and fatness traits. Our findings not only provide insights into the genetic architecture of porcine growth and fatness traits but also provide potential markers for selective breeding of these traits in Sujiang pigs.

Population structure of 3907 worldwide pigs and the introgression of Chinese indigenous pigs by European pigs

With the improvement in sequencing technology and the decrease in sequencing cost, increasing amounts of genomic data for pigs have been uploaded to public databases. However, no researchers have to date integrated all currently available data to uncover the global genetic status of pigs. Meanwhile, little is known about the introgression from European to Chinese pigs and its underlying influences. Therefore, we integrated the effective genotype data of 3907 pigs from 193 populations worldwide using population genetic analysis, gene flow analysis and a sharing-IBD study. These findings illustrate not only the population structure of 59 Chinese native breeds and others but also the amounts of gene flow and introgression that have occurred between Western and Chinese pigs. In addition, we demonstrate the presence of introgressed European haplotypes in Chinese indigenous breeds and identify relevant introgressed regions that contain genes associated with growth and feed efficiency. Moreover, we compare the introgression patterns of Western and Chinese pigs and further discuss possible explanations for why the level of introgression differs between Chinese pig breeds and Western modern breeds. Collectively, this study provides a fine global population structure analysis of pigs and presents evidence of European pigs being interbred with local breeds in China.

Effects of age on subcutaneous adipose tissue proteins in Chinese indigenous Ningxiang pig by TMT-labeled quantitative proteomics

Adipose tissue not only affects meat quality and animal productivity, but also participates in inflammation and immunity. Ningxiang pig is famous for their excellent meat quality, disease resistance and tolerance of roughage. It is not yet well known how proteins in adipose tissue is dynamically regulated during the growth of Ningxiang pig. This report studies the proteomic changes in subcutaneous adipose tissue in Ningxiang pigs to gain a better understanding of the molecular mechanism of fat development during the growth period. By TMT-labeled quantitative proteomic analysis of subcutaneous adipose tissue of 9 purebred Ningxiang pigs of different ages, we identified 2533 unique proteins and 716 differentially abundant proteins (DAPs), of which more than half of the DAPs were concentrated in the 90d-210d period. Retrograde endocannabinoid signaling was only significantly enriched in DAPs of N90d vs N30d, Alcoholism and Graft-versus-host disease were only significantly enriched in DAPs of N210d vs N90d. Proteins related to dilated cardiomyopathy was found to be an important pathway in fat development and lipid metabolism. A variety of novel DAPs involved in maintaining mitochondrial function and cell viability, such as NDUFS6, SDHB, COX5A, ATP5D and TNNT1, which play a role in controlling the prediction networks, may indirectly regulate the development and functional maintenance of adipocytes. SIGNIFICANCE: These age-dependent DAPs discovered in this study may help expand the understanding of the molecular mechanisms of the development, function maintenance and transformation of adipose tissue in Ningxiang pig for developing new strategies for improving meat quality and pig breeding in the future.

Genetic Evaluation and Population Structure of Jiangsu Native Pigs in China Revealed by SINE Insertion Polymorphisms

Simple Summary

In a previous study, 30 SINE-RIPs were applied for population genetic analysis in 7 Chinese miniature pig populations and approved effectively in the genetic distances and breed-relationships between these populations. There are abundant indigenous pigs famous across the world for their prolificacy in the Jiangsu Province of eastern China, such as Meishan, Erhualian. Since pork production relies on limited commercial breeds such as Landrace, Large White, and Duroc pigs, characterized by maximizing productivity in intensive production systems, these indigenous pigs are nowadays decrease sharply. The genetic characterizations of these resources are essential requirements for the development of conservation, selection, and sustainable utilizations. Therefore, SINE-RIPs were selected to evaluate the genetic variation and population structure of Jiangsu pig populations and the results may assist with the conservation and utilization of these native pig populations.

Abstract

Short interspersed nuclear elements (SINEs), one type of retrotransposon, are considered to be ideal molecular markers due to their wide distribution in the genome, high copy number, and high polymorphism. Preliminary studies have identified more than 35,000 SINE-retrotransposon insertion polymorphisms (RIPs) in the pig genome. In this study, 18 SINE-RIPs were used to evaluate the genetic variation and population structure of seven native pig populations and two crossbreeds in the Jiangsu Province of China. Two commercial pig breeds (Duroc and Large White) and one Italian native breed (Sicilian Black pig) were selected as the control. The results showed that all 18 SINE-RIPs were polymorphic among these pigs. The Jiangsu native pig populations (Erhualian, Fengjing, Middle Meishan, Mi, Shawutou, Small Meishan, and Huai) were shown to be more polymorphic than the crossbreeds (Sushan and Sujiang) and external breeds (Sicilian Black pig, Large White, and Duroc) based on the expected heterozygosity and polymorphic information content values. Some native pigs, including Small Meishan, Mi, Middle Meishan, and Erhualian, had a higher degree of inbreeding according to the F_IS values. Based on the neighbor-joining tree, all of the Jiangsu native pig populations formed one branch, while the three external pig breeds formed the other branches, with the two crossbreeds containing more than 50% external pig ancestry. The Huai pigs were independent of the other Jiangsu native pigs but shared a common ancestor with Sujiang and Mi. The results provide a new perspective on the population structure of these native pig breeds and will assist with the conservation and utilization of Chinese native pigs.

Genome-Wide Selection Sweep between Wild and Local Pigs from Europe for the Investigation of the Hereditary Characteristics of Domestication in Sus Scrofa

The phenotypic characteristics of existing domestic pigs (DPs) greatly differ from those of wild boar (WB) populations thousands of years ago. After thousands of years of human domestication, WB and DP have diverged greatly in terms of genetics. Theoretically, worldwide local pigs have independent contributions from their local WBs at the beginning of Sus scrofa domestication. The investigation of the vicissitude of the heredity material between domestic populations and their wild ancestors will help in further understanding the domestication history of domestic animals. In the present study, we performed a genome-wide association scan (GWSA) and phylogeny estimation with a total of 1098 public European Illumina 60K single nucleotide polymorphism data, which included 650 local DPs and 448 WBs. The results revealed that the phylogenetic relationship of WBs corresponds to their geographical distribution and carries large divergence with DPs, and all WB breeds (e.g., HRWB, SBWB, and TIWB) presents a closely linkage with the middle WB (e.g., HRWB, and PLWB). In addition, 64 selected candidate genes (e.g., IDH2, PIP5K1B, SMARCA2, KIF5C, and TJP2) were identified from GWSA. A total of 63 known multiple biological functional pathways were annotated by 22 genes, and ubiquinone and other terpenoid-quinone biosynthesis pathways that belong to the metabolism of cofactors and vitamins were significantly enriched (p < 0.05). The most frequent (28.57%) pathways were classified under metabolism. We confirmed that the middle European WB has made an important genetic contribution to the entire European WB populations. A series of selected genes discovered from this study provides the scientific community with a deeper understanding of the heredity performance of metabolism and emotion and the real purpose behind domestication.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

Identify known and novel candidate genes associated with backfat thickness in Duroc pigs by large-scale genome-wide association analysis

Backfat thickness (BFT) is complex and economically important traits in the pig industry, since it reflects fat deposition and can be used to measure the carcass lean meat percentage in pigs. In this study, all 6,550 pigs were genotyped using the Geneseek Porcine 50K SNP Chip to identify SNPs related to BFT and to search for candidate genes through genome-wide association analysis in two Duroc populations. In total, 80 SNPs, including 39 significant and 41 suggestive SNPs, and 6 QTLs were identified significantly associated with the BFT. In addition, 9 candidate genes, including a proven major gene MC4R, 3 important candidate genes (RYR1, HMGA1, and NUDT3) which were previously described as related to BFT, and 5 novel candidate genes (SIRT2, NKAIN2, AMH, SORCS1, and SORCS3) were found based on their potential functional roles in BFT. The functions of candidate genes and gene set enrichment analysis indicate that most important pathways are related to energy homeostasis and adipogenesis. Finally, our data suggest that most of the candidate genes can be directly used for genetic improvement through molecular markers, except that the MC4R gene has an antagonistic effect on growth rate and carcass lean meat percentage in breeding. Our results will advance our understanding of the complex genetic architecture of BFT traits and laid the foundation for additional genetic studies to increase carcass lean meat percentage of pig through marker-assisted selection and/or genomic selection.

A novel 31bp deletion within the CDKL5 gene is significantly associated with growth traits in Dezhou donkey

The discovery of molecular markers which associate with livestock economic traits is of great significance for livestock breeding. Selective analysis has found a potential correlation between CDKL5 and growth traits, but there is still a lack of experimental proof. In this study, a 31-bp deletion (g.176595_176626delATGTCACATGTGGTACTGCCATGTGGAATTT) of CDKL5 gene was found by sequencing. The 31-bp indel was then genotyped in 380 individuals of Dezhou donkeys by polyacrylamide gel electrophoresis and there were three genotypes in this population. After the association analysis between growth traits and genotypes, it was found that this 31-bp indel polymorphism was significantly associated with the chest circumference of Dezhou donkeys (p < 0.05), and body length, chest depth and rump width (p < 0.01). In addition, all individuals with DD genotype were better than those with other genotypes in growth traits. This study revealed that a newly identified polymorphic locus in the CDKL5 gene is related to growth traits, which provides a molecular marker for genetic improvement of Dezhou donkey and may lay a solid foundation for the breeding of Dezhou donkey.

Whole-genome SNP markers reveal conservation status, signatures of selection, and introgression in Chinese Laiwu pigs

Laiwu pigs are a Chinese indigenous breed that is renowned for its exceptionally high intramuscular fat content (average greater than 6%), providing an excellent genetic resource for the genetic improvement of meat quality of modern commercial pigs. To uncover genetic diversity, population structure, signature of selection, and potential exotic introgression in this breed, we sampled 238 Laiwu pigs from a state-supported conservation population and genotyped these individuals using GeneSeek 80K SNP BeadChip. We then conducted in-depth population genetics analyses for the Laiwu pig in a context of 1,116 pigs from 42 Eurasian diverse breeds. First, we show that the current Laiwu population has more abundant genetic diversity than the population of 18 years ago likely due to gene flow from European commercial breeds. Both neighbor-joining (NJ) and principal component analyses indicate the introgression of European haplotypes into Laiwu pigs. The admixture analysis reveals that an average 26.66% of Laiwu genetic components are of European origin. Then, we assigned the tested individuals to different families according to their clustering patterns in the NJ tree and proposed a family-based conservation strategy to reduce the risk of inbreeding depression in Laiwu pigs. Next, we explored three statistics (ROH and iHS and EigenGWAS) to identify a list of candidate genes for fat deposition, reproduction, and growth in Laiwu pigs. Last, we detected a strong signature of introgression from European pigs into Laiwu pigs at the GPC6 locus that regulates the growth of developing long bones. Further association analyses indicate that the introgressed GPC6 haplotype likely contributed to the improvement of growth performance in Laiwu pigs. Altogether, this study not only benefits the better conservation of the Laiwu pig, but also advances our knowledge of the poorly understood effect of human-mediated introgression on phenotypic traits in Chinese indigenous pigs.