Genetic Diversity of Indigenous Pigs from South China Area Revealed by SNP Array

Exploring the Interplay between the Hologenome and Complex Traits in Bovine and Porcine Animals Using Genome-Wide Association Analysis

Genome-wide association studies (GWAS) significantly enhance our ability to identify trait-associated genomic variants by considering the host genome. Moreover, the hologenome refers to the host organism's collective genetic material and its associated microbiome. In this study, we utilized the hologenome framework, called Hologenome-wide association studies (HWAS), to dissect the architecture of complex traits, including milk yield, methane emissions, rumen physiology in cattle, and gut microbial composition in pigs. We employed four statistical models: (1) GWAS, (2) Microbial GWAS (M-GWAS), (3) HWAS-CG (hologenome interaction estimated using COvariance between Random Effects Genome-based restricted maximum likelihood (CORE-GREML)), and (4) HWAS-H (hologenome interaction estimated using the Hadamard product method). We applied Bonferroni correction to interpret the significant associations in the complex traits. The GWAS and M-GWAS detected one and sixteen significant SNPs for milk yield traits, respectively, whereas the HWAS-CG and HWAS-H each identified eight SNPs. Moreover, HWAS-CG revealed four, and the remaining models identified three SNPs each for methane emissions traits. The GWAS and HWAS-CG detected one and three SNPs for rumen physiology traits, respectively. For the pigs' gut microbial composition traits, the GWAS, M-GWAS, HWAS-CG, and HWAS-H identified 14, 16, 13, and 12 SNPs, respectively. We further explored these associations through SNP annotation and by analyzing biological processes and functional pathways. Additionally, we integrated our GWA results with expression quantitative trait locus (eQTL) data using transcriptome-wide association studies (TWAS) and summary-based Mendelian randomization (SMR) methods for a more comprehensive understanding of SNP-trait associations. Our study revealed hologenomic variability in agriculturally important traits, enhancing our understanding of host-microbiome interactions.

Population genetics reveals new introgression in the nucleus herd of min pigs

OBJECTIVE

Min pigs are a unique genetic resource among local pig breeds in China. They have more excellent characteristics in cold and stress resistance, good meat quality, and a high reproductive rate. However, the genetic structure and driving factors remain unclear in the nucleus herd. In this study, the genetic diversity of Min pigs was studied to reveal the formation mechanism of its unique genetic structure. We hope to protect and develop the genetic resources of Min pigs.

METHODS

We analyzed different types of genes to identify the genetic structure and gene introgression pattern of Min pigs. The nuclear DNA dataset includes information on 21 microsatellite loci and 6 Y-chromosome genes, and the mitochondrial D-loop gene is selected to represent maternal lineages. The above genes are all from the nucleus herd of Min pigs.

RESULTS

The results of genetic structure identification and analysis of potential exogenous gene introgression patterns indicate that the nucleus herd of Min pigs maintains a high level of genetic diversity (polymorphism information content = 0.713, expected heterozygosity = 0.662, observed heterozygosity = 0.612). Compared with other Asian pig breeds, the formation of Min pig breeds is more special. Gene introgression from European pig breeds to Min pigs has occurred, which is characterized by complete introgression of paternal genes and incomplete introgression of maternal genes.

CONCLUSION

Gene introgression caused by cross-breeding is not the main factor leading to the formation of the current genetic structure of Min pigs, but this process has increased the level of genetic diversity in the nucleus herd. Compared with the influence of gene introgression, our research suggest that artificial selection and environmental adaptive evolution make Min pigs form unique genetic characteristics.

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

Evaluation of the Genetic Diversity, Population Structure and Selection Signatures of Three Native Chinese Pig Populations

Indigenous pig populations in Hainan Province live in tropical climate conditions and a relatively closed geographical environment, which has contributed to the formation of some excellent characteristics, such as heat tolerance, strong disease resistance and excellent meat quality. Over the past few decades, the number of these pig populations has decreased sharply, largely due to a decrease in growth rate and poor lean meat percentage. For effective conservation of these genetic resources (such as heat tolerance, meat quality and disease resistance), the whole-genome sequencing data of 78 individuals from 3 native Chinese pig populations, including Wuzhishan (WZS), Tunchang (TC) and Dingan (DA), were obtained using a 150 bp paired-end platform, and 25 individuals from two foreign breeds, including Landrace (LR) and Large White (LW), were downloaded from a public database. A total of 28,384,282 SNPs were identified, of which 27,134,233 SNPs were identified in native Chinese pig populations. Both genetic diversity statistics and linkage disequilibrium (LD) analysis indicated that indigenous pig populations displayed high genetic diversity. The result of population structure implied the uniqueness of each native Chinese pig population. The selection signatures were detected between indigenous pig populations and foreign breeds by using the population differentiation index (F_ST) method. A total of 359 candidate genes were identified, and some genes may affect characteristics such as immunity (IL-2, IL-21 and ZFYVE16), adaptability (APBA1), reproduction (FGF2, RNF17, ADAD1 and HIPK4), meat quality (ABCA1, ADIG, TLE4 and IRX5), and heat tolerance (VPS13A, HSPA4). Overall, the findings of this study will provide some valuable insights for the future breeding, conservation and utilization of these three Chinese indigenous pig populations.

A web tool for the global identification of pig breeds

BACKGROUND

Natural and artificial selection for more than 9000 years have led to a variety of domestic pig breeds. Accurate identification of pig breeds is important for breed conservation, sustainable breeding, pork traceability, and local resource registration.

RESULTS

We evaluated the performance of four selectors and six classifiers for breed identification using a wide range of pig breeds (N = 91). The internal cross-validation and external independent testing showed that partial least squares regression (PLSR) was the most effective selector and partial least squares-discriminant analysis (PLS-DA) was the most powerful classifier for breed identification among many breeds. Five-fold cross-validation indicated that using PLSR as the selector and PLS-DA as the classifier to discriminate 91 pig breeds yielded 98.4% accuracy with only 3K single nucleotide polymorphisms (SNPs). We also constructed a reference dataset with 124 pig breeds and used it to develop the web tool iDIGs ( http://alphaindex.zju.edu.cn/iDIGs_en/ ) as a comprehensive application for global pig breed identification. iDIGs allows users to (1) identify pig breeds without a reference population and (2) design small panels to discriminate several specific pig breeds.

CONCLUSIONS

In this study, we proved that breed identification among a wide range of pig breeds is feasible and we developed a web tool for such pig breed identification.

Lost pigs of Angola: Whole genome sequencing reveals unique regions of selection with emphasis on metabolism and feed efficiency

Angola, in the western coast of Africa, has been through dramatic social events that have led to the near-disappearance of native swine populations, and the recent introduction of European exotic breeds has also contributed to the erosion of this native swine repertoire. In an effort to investigate the genetic basis of native pigs in Angola (ANG) we have generated whole genomes from animals of a remote local pig population in Huambo province, which we have compared with 78 genomes of European and Asian pig breeds as well as European and Asian wild boars that are currently in public domain. Analyses of population structure showed that ANG pigs grouped within the European cluster and were clearly separated from Asian pig breeds. Pairwise F_ST ranged from 0.14 to 0.26, ANG pigs display lower levels of genetic differentiation towards European breeds. Finally, we have identified candidate regions for selection using a complementary approach based on various methods. All results suggest that selection towards feed efficiency and metabolism has occurred. Moreover, all analysis identified CDKAL1 gene, which is related with insulin and cholesterol metabolism, as a candidate gene overlapping signatures of selection unique to ANG pigs. This study presents the first assessment of the genetic relationship between ANG pigs and other world breeds and uncovers selection signatures that may indicate adaptation features unique to this important genetic resource.

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.

Genetic diversity and selection signatures of four indigenous pig breeds from eastern China

Chinese indigenous pig breeds have been undergoing selection for thousands of years, and have become invaluable genetic sources over the world. To investigate the population structure and genetic diversity of Jinhua (JH), Longyou Black (LYW), Shengxian Spotted (SXH), and Lanxi Spotted (LXH) breeds, a total of 200 pigs belonging to 10 diverse population were genotyped using SNP chips. The results showed that LYW pigs exhibited higher level of heterozygosity than the other indigenous pigs. In addition, gene introgression from intensively reared commercial pig breeds to LYW pigs was detected. Moreover, selection signature analysis revealed the possibility of differences between Chinese indigenous and intensively reared commercial pig breeds were mainly present for meat and carcass traits. Furthermore, we found that ANXA13, DISP1, and SRSF6 were the nearest genes located around the common selection signatures detected between each indigenous pig breed and Chinese wild boars. Our findings provide new insights into the selection signatures of Chinese indigenous pigs, and may contribute to future pig breeding.

Genome-Wide Assessment of a Korean Composite Pig Breed, Woori-Heukdon

A Korean synthetic pig breed, Woori-Heukdon (WRH; F3), was developed by crossing parental breeds (Korean native pig [KNP] and Korean Duroc [DUC]) with their crossbred populations (F1 and F2). This study in genome-wide assessed a total of 2,074 pigs which include the crossbred and the parental populations using the Illumina PorcineSNP60 BeadChip. After quality control of the initial datasets, we performed population structure, genetic diversity, and runs of homozygosity (ROH) analyses. Population structure analyses showed that crossbred populations were genetically influenced by the parental breeds according to their generation stage in the crossbreeding scheme. Moreover, principal component analysis showed the dispersed cluster of WRH, which might reflect introducing a new breeding group into the previous one. Expected heterozygosity values, which were used to assess genetic diversity, were .365, .349, .336, .330, and .211 for WRH, F2, F1, DUC, and KNP, respectively. The inbreeding coefficient based on ROH was the highest in KNP (.409), followed by WRH (.186), DUC (.178), F2 (.107), and F1 (.035). Moreover, the frequency of short ROH decreased according to the crossing stage (from F1 to WRH). Alternatively, the frequency of medium and long ROH increased, which indicated recent inbreeding in F2 and WRH. Furthermore, gene annotation of the ROH islands in WRH that might be inherited from their parental breeds revealed several interesting candidate genes that may be associated with adaptation, meat quality, production, and reproduction traits in pigs.

ITGB6 inhibits the proliferation of porcine skeletal muscle satellite cells

The formation of embryonic muscle fibers determines the amount of postnatal muscles and is regulated by a variety of signaling pathways and transcription factors. Previously, by using chromatin immunoprecipitation-sequencing and RNA-Seq techniques, we identified a large number of genes that are regulated by H3K27me3 in porcine embryonic skeletal muscles. Among these genes, we found that ITGB6 is regulated by H3K27me3. However, its function in muscle development is unknown. In this study, we first verified that ITGB6 was differentially regulated by H3K27me3 and that its expression levels were upregulated in porcine skeletal muscles at embryonic Days 33, 65, and 90. Then, we performed gain- or loss-of-function studies on porcine skeletal muscle satellite cells to study the role of ITGB6 in porcine skeletal muscle development. The proliferation of porcine skeletal muscle satellite cells was studied through real-time polymerase chain reaction, Cell Counting Kit-8, 5-ethynyl-2'-deoxyuridine staining, Western blot, and flow cytometry analyses. We found that the ITGB6 gene was regulated by H3K27me3 during muscle development and had an inhibitory effect on the proliferation of porcine skeletal muscle satellite cells.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Copy Number Variants in Four Italian Turkey Breeds

Simple Summary

Hybrid Turkey selection is focusing on meat production traits characterized by high genetic heritability; the strong directional selection is well known to produce a constant loss in genetic diversity. Genetic characterization is one of the essential activities in the management of populations at risk of extinction. In addition, the genetic structure at the population level and the relationships between individuals are nowadays analysable at the genomic level. In this paper, the genome of 4 different Italian turkey breeds included in the Autochthonous Italian Poultry Breeds Register are analysed in order to obtain a genome-wide Copy Number Variant scan to ameliorate the existing knowledge of the genomic structure of Italian local turkey breeds. Differences have been described at genomic level for physiological, reproductive, and behavioral traits. The analyzed breeds are clearly distinguishable at the genomic level, and their relationships are clearly linked to their geographical origin and to the history of the rural structure of their developing regions. Genome information based on Copy Number Variant (CNV) detection has generated important information in this study concerning the uniqueness of the Italian local turkey breeds.

Abstract

Heritage breeds can be considered a genetic reservoir of genetic variability to be conserved and valorized considering their historical, cultural, and adaptive characteristics and possibly for their high potential in commercial hybrid genetic improvement by gene introgression. The aim of the present research is to investigate via Copy Number Variant (CNVs) the genomic makeup of 4 Italian autochthonous turkey breeds (Bronzato Comune—BrCI, 24; Ermellinato di Rovigo—ErRo, 24; Parma e Piacenza—PrPc, 25; Romagnolo—RoMa, 29). CNVs detection was performed using two different software and an interbreed CNVs comparison was carried out. A total of 1077 CNVs were identified in 102 turkeys, summarized into 519 CNV regions (CNVRs), which resulted after merging in 101 and 18 breed and shared regions. Biodiversity was analyzed using the effective information supplied by CNVs analysis, and BrCI and ErRo were characterized by a low mapped CNV number. Differences were described at a genomic level related to physiological, reproductive, and behavioral traits. The comparison with other three Italian turkey breeds (Brianzolo, Colle Euganei, and Nero Italiano) using a CNV data set available in the literature showed high clustering properties at the genomic level, and their relationships are strictly linked to the geographical origin and to the history of the rural structure of their native regions.

Assessment of the Growth and Reproductive Performance of Cloned Pietrain Boars

How to maximize the use of the genetic merits of the high-ranking boars (also called superior ones) is a considerable question in the pig breeding industry, considering the money and time spent on selection. Somatic cell nuclear transfer (SCNT) is one of the potential ways to answer the question, which can be applied to produce clones with genetic resources of superior boar for the production of commercial pigs. For practical application, it is essential to investigate whether the clones and their progeny keep behaving better than the "normal boars", considering that in vitro culture and transfer manipulation would cause a series of harmful effects to the development of clones. In this study, 59,061 cloned embryos were transferred into 250 recipient sows to produce the clones of superior Pietrain boars. The growth performance of 12 clones and 36 non-clones and the semen quality of 19 clones and 28 non-clones were compared. The reproductive performance of 21 clones and 25 non-clones were also tested. Furthermore, we made a comparison in the growth performance between 466 progeny of the clones and 822 progeny of the non-clones. Our results showed that no significant difference in semen quality and reproductive performance was observed between the clones and the non-clones, although the clones grew slower and exhibited smaller body size than the non-clones. The F1 progeny of the clones showed a greater growth rate than the non-clones. Our results demonstrated through the large animal population showed that SCNT manipulation resulted in a low growth rate and small body size, but the clones could normally produce F1 progeny with excellent growth traits to bring more economic benefits. Therefore, SCNT could be effective in enlarging the merit genetics of the superior boars and increasing the economic benefits in pig reproduction and breeding.

Identification and Expression Pattern of EZH2 in Pig Developing Fetuses

The proper methylation status of histones is essential for appropriate cell lineage and organogenesis. EZH2, a methyltransferase catalyzing H3K27me3, has been abundantly studied in human and mouse embryonic development. The pig is an increasing important animal model for molecular study and pharmaceutical research. However, the transcript variant and temporal expression pattern of EZH2 in the middle and late porcine fetus are still unknown. Here, we identified the coding sequence of the EZH2 gene and characterized its expression pattern in fetal tissues of Duroc pigs at 65- and 90-day postcoitus (dpc). Our results showed that the coding sequence of EZH2 was 2241 bp, encoding 746 amino acids. There were 9 amino acid insertions and an amino acid substitution in this transcript compared with the validated reference sequence in NCBI. EZH2 was ubiquitously expressed in the fetal tissues of two time points with different expression levels. These results validated a different transcript in pigs and characterized its expression profile in fetal tissues of different gestation stages, which indicated that EZH2 played important roles during porcine embryonic development.

Whole genome sequence analysis reveals genetic structure and X-chromosome haplotype structure in indigenous Chinese pigs

Chinese indigenous pigs exhibit considerable phenotypic diversity, but their population structure and the genetic basis of agriculturally important traits need further exploration. Here, we sequenced the whole genomes of 24 individual pigs representing 22 breeds distributed throughout China. For comparison with European and commercial breeds (one pig per breed), we included seven published pig genomes with our new genomes for analyses. Our results showed that breeds grouped together based on morphological classifications are not necessarily more genetically similar to each other than to breeds from other groups. We found that genetic material from European pigs likely introgressed into five Chinese breeds. We have identified two new subpopulations of domestic pigs that encompass morphology-based criteria in China. The Southern Chinese subpopulation comprises the classical South Chinese Type and part of the Central China Type. In contrast, the Northern Chinese subpopulation comprises the North China Type, the Lower Yangtze River Basin Type, the Southwest Type, the Plateau Type, and the remainder of the Central China Type. Eight haplotypes and two recombination sites were identified within a conserved 40.09 Mb linkage-disequilibrium (LD) block on the X chromosome. Potential candidate genes (LEPR, FANCC, COL1A1, and PCCA) influencing body size were identified. Our findings provide insights into the phylogeny of Chinese indigenous pig breeds and benefit gene mining efforts to improve major economic traits.

The fine-scale genetic structure and selection signals of Chinese indigenous pigs

Genome-wide SNP profiling has yielded insights into the genetic structure of China indigenous pigs, but has focused on a limited number of populations. Here, we present an analysis of population structure and signals of positive selection in 42 Chinese pig populations that represent the most extensive pig phenotypic diversity in China, using genotype data of 1.1 million SNPs on customized Beadchips. This unravels the fine-scale genetic diversity, phylogenic relationships, and population structure of these populations, which shows remarkably concordance between genetic clusters and geography with few exceptions. We also reveal the genetic contribution to North Chinese pigs from European modern pigs. Furthermore, we identify possible targets of selection in the Tibetan pig, including the well-characterized hypoxia gene (EPAS1) and several previously unrecognized candidates. Intriguingly, the selected haplotype in the EPAS1 gene is associated with higher hemoglobin contents in Tibetan pigs, which is different from the protective role of EPAS1 in the high-altitude adaptation in Tibetan dogs and their owners. Additionally, we present evidence for the causality between EDNRB variants and the two-end-black (TEB) coat color phenotype in all Chinese pig populations except the Jinhua pig. We hypothesize that distinct targets have been independently selected for the formation of the TEB phenotype in Chinese pigs of different geographic origins. This highlights the importance of characterizing population-specific genetic determinants for heritable phenotype in diverse pig populations.

Study on Hematological and Biochemical Characters of Cloned Duroc Pigs and Their Progeny

Simple Summary

Cloning is the most promising technique for passing the excellent phenotypes of the best individuals in the population. Here we studied the effects of cloning on Duroc pig, which is the most popular sire used in pig production due to its good growth and meat quality. Understanding the changes of cloned Duroc pigs and their progenies is of great importance for animal breeding and public acceptance. The results of this study suggested that there were no difference in blood parameters between the cloned Duroc and the conventionally bred Duroc and their progenies.

Abstract

To increase public understanding in cloned animals produced by somatic cell nuclear transfer technology, our previous study investigated the carcass trait and meat quality of the clones (paper accepted), and this study we further evaluate differences by investigating the blood parameters in cloned pigs and their progeny. We collected blood samples from the clones and conventionally bred non-clones and their progeny, and investigated their hematological and blood biochemical characters. Our results supported the hypothesis that there was no significant difference between clones and non-clones, or their progeny. Taken together, the data demonstrated that the clones or their progeny were similar with their controls in terms of blood parameters, although there were still other kinds of disorders, such as abnormal DNA methylation or histone modifications that needs further investigation. The data in this study agreed that cloning technique could be used to preserve and enlarge the genetics of the superior boars in pig breeding industry, especially in facing of the deadly threat of African Swine fever happened in China.