Genome-wide association analysis reveals genomic regions on Chromosome 13 affecting litter size and candidate genes for uterine horn length in Erhualian pigs

Genetic diversity, growth and heart function of Auckland Island pigs, a potential source for organ xenotransplantation

One of the prerequisites for successful organ xenotransplantation is a reasonable size match between the porcine organ and the recipient's organ to be replaced. Therefore, the selection of a suitable genetic background of source pigs is important. In this study, we investigated body and organ growth, cardiac function, and genetic diversity of a colony of Auckland Island pigs established at the Center for Innovative Medical Models (CiMM), LMU Munich. Male and female Auckland Island pig kidney cells (selected to be free of porcine endogenous retrovirus C) were imported from New Zealand, and founder animals were established by somatic cell nuclear transfer (SCNT). Morphologically, Auckland Island pigs have smaller body stature compared to many domestic pig breeds, rendering their organ dimensions well-suited for human transplantation. Furthermore, echocardiography assessments of Auckland Island pig hearts indicated normal structure and functioning across various age groups throughout the study. Single nucleotide polymorphism (SNP) analysis revealed higher runs of homozygosity (ROH) in Auckland Island pigs compared to other domestic pig breeds and demonstrated that the entire locus coding the swine leukocyte antigens (SLAs) was homozygous. Based on these findings, Auckland Island pigs represent a promising genetic background for organ xenotransplantation.

Genomic Scanning of Inbreeding Depression for Litter Size in Two Varieties of Iberian Pigs

Inbreeding depression is expected to be more pronounced in fitness-related traits, such as pig litter size. Recent studies have suggested that the genetic determinism of inbreeding depression may be heterogeneous across the genome. Therefore, the objective of this study was to conduct a genomic scan of the whole pig autosomal genome to detect the genomic regions that control inbreeding depression for litter size in two varieties of Iberian pigs (Entrepelado and Retinto). The datasets consisted of 2069 (338 sows) and 2028 (327 sows) records of litter size (Total Number Born and Number Born Alive) for the Entrepelado and Retinto varieties. All sows were genotyped using the Geneseek GGP PorcineHD 70 K chip. We employed the Unfavorable Haplotype Finder software to extract runs of homozygosity (ROHs) and conducted a mixed-model analysis to identify highly significant differences between homozygous and heterozygous sows for each specific ROH. A total of eight genomic regions located on SSC2, SSC5, SSC7, SSC8, and SSC13 were significantly associated with inbreeding depression, housing some relevant genes such as FSHR, LHCGR, CORIN, AQP6, and CEP120.

Effects of Varying Levels of Wheat Bran Dietary Fiber on Growth Performance, Fiber Digestibility and Gut Microbiota in Erhualian and Large White Pigs

To evaluate the tolerance of a high-fiber diet in Erhualian pigs (Er-HL), the present investigation systematically investigated the ramifications of varying wheat bran fiber levels, specified as total dietary fiber (TDF) values of 14.07%, 16.32%, 17.99%, and 18.85%, on growth performance, fiber digestibility and gut microbiota in Er-HL, large Large White pigs (L-LW, the same physiological stage as the Er-HL) and small Large White pigs (S-LW, the same body weight as the Er-HL). Our results revealed that fiber levels exerted no discernable impact on growth performance (average daily feed intake (ADFI), and average daily gain (ADG)) of Er-HL (p > 0.05). Conversely, L-LW exhibited a decrease in ADFI and ADG with increasing fiber levels (p < 0.05). Notably, the apparent total tract digestibility (ATTD) of various fiber components, including neutral detergent fiber (NDF), acid detergent fiber (ADF), hemicellulose, TDF and insoluble dietary fiber (IDF), in Er-HL were significantly higher than those in S-LW and L-LW irrespective of diets (p < 0.05). The ATTD of cellulose and hemicellulose in Er-HL significantly decreased with increasing fiber levels (p < 0.05), yet remained statistically indifferent when comparing the 7%-wheat-bran-replaced diet (7% WRB, TDF 16.32%) to the basal diet (TDF 14.07%) (p > 0.05). The cecal microbiota of Er-HL had higher richness estimators (Chao1 and ACE) than those of S-LW and L-LW irrespective of diets (p < 0.01). Breed serves as a pivotal determinant in shaping swine gut microbiota. Thirteen genera were selected as the key bacteria related to high fiber digestibility of Er-HL. Further functional examination of these key genera elucidated an enrichment of pathways pertinent to carbohydrate metabolism in Er-HL samples compared with S-LW and L-LW samples. In summary, Er-HL exhibited high-fiber tolerance both in terms of growth performance and fiber digestibility compared with Large White pigs. Specifically, the ATTD of NDF, ADF, hemicellulose, IDF and TDF were significantly higher in Er-HL compared with L-LW and S-LW, irrespective of diets. Fiber level exerted no discernable impact on growth performance (ADFI, ADG) and the ATTD of fiber (NDF, ADF, IDF and TDF) in Er-HL. The optimum fiber level of the Er-HL was identified as 7% WRB (TDF 16.32%). Thirteen genera were ascertained to significantly contribute to high fiber digestibility of Er-HL, correlating with an enhancement of carbohydrate metabolism pathways.

Genome-Wide Association Analysis Identified Variants Associated with Body Measurement and Reproduction Traits in Shaziling Pigs

With the increasing popularity of genomic sequencing, breeders pay more attention to identifying the crucial molecular markers and quantitative trait loci for improving the body size and reproduction traits that could affect the production efficiency of pig-breeding enterprises. Nevertheless, for the Shaziling pig, a well-known indigenous breed in China, the relationship between phenotypes and their corresponding genetic architecture remains largely unknown. Herein, in the Shaziling population, a total of 190 samples were genotyped using the Geneseek Porcine 50K SNP Chip, obtaining 41857 SNPs for further analysis. For phenotypes, two body measurement traits and four reproduction traits in the first parity from the 190 Shaziling sows were measured and recorded, respectively. Subsequently, a genome-wide association study (GWAS) between the SNPs and the six phenotypes was performed. The correlation between body size and reproduction phenotypes was not statistically significant. A total of 31 SNPs were found to be associated with body length (BL), chest circumference (CC), number of healthy births (NHB), and number of stillborns (NSB). Gene annotation for those candidate SNPs identified 18 functional genes, such as GLP1R, NFYA, NANOG, COX7A2, BMPR1B, FOXP1, SLC29A1, CNTNAP4, and KIT, which exert important roles in skeletal morphogenesis, chondrogenesis, obesity, and embryonic and fetal development. These findings are helpful to better understand the genetic mechanism for body size and reproduction phenotypes, while the phenotype-associated SNPs could be used as the molecular markers for the pig breeding programs.

Allelic to genome wide perspectives of swine genetic variation to litter size and its component traits

Litter size is a complex and sex limited trait that depends on various biological, managemental and environmental factors. Owing to its low heritability it is inefficaciously selected by traditional methods. However, due to higher heritability of ovulation rate and embryo survival, selection based on component traits of litter size is advocated. QTL analysis and candidate gene approach are among the various supplementary/alternate strategies for selection of litter size. QTL analysis is aimed at identifying genomic regions affecting trait of interest significantly. Candidate gene approach necessitates identification of genes potentially affecting the trait. There are various genes that significantly affect litter size and its component traits viz. ESR, LEP, BF, IGFBP, RBP4, PRLR, CTNNAL1, WNT10B, TCF12, DAZ, and RNF4. These genes affect litter size in a complex interacting manner. Lately, genome wide association study (GWAS) have been utilized to unveil the genetic and biological background of litter traits, and elucidate the genes governing litter size. Favorable SNPs in these genes have been identified and offers a scope for inclusion in selection programs thereby increasing breeding efficiency and profit in pigs. The review provides a comprehensive coverage of investigations carried out globally to unravel the genetic variation in litter size and its component traits in pigs, both at allelic and genome wide level. It offers a current perspective on different strategies including the profiling of candidate genes, QTLs, and genome wide association studies as an aid to efficient selection for litter size and its component traits.

Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation

It has been established that through binding to bone morphogenetic proteins (BMPs), bone morphogenetic protein receptor I B (BMPR1B) can mediate transforming growth factor β (TGF-β) signal transduction, and is involved in the regulation of several biological processes, such as bone and muscle formation and homeostasis, as well as folliculogenesis. Also known as FecB, BMPR1B has been reported as the major gene for sheep prolificacy. A number of previous studies have analyzed the relationship between single nucleotide polymorphisms (SNPs) in this gene and its related performance. In recent years, with the illustration of the effect of retrotransposon insertion on the expression of the proximal genes or phenotypic variation, retrotransposon insertion polymorphisms (RIPs) have been used as a novel type of molecular marker in the evaluation of evolution, population structure and breeding of plant and domestic animals. In this study, the RIPs in porcine BMPR1B gene were excavated, and thereafter verified using a comparative genome and polymerase chain reaction (PCR). The potential effects of phenotype, gene expression and functions related to RIPs were also explored. The results showed that 13 distinct RIPs were identified in introns of porcine BMPR1B. Among these, only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 displayed a close relationship with the growth traits of Large White pigs. Moreover, the total number of BMPR1B-SINE+/+-RIP9 individuals born was found to be significantly higher than that of SINE−/− (p < 0.05). These two RIPs showed an obvious distribution pattern among Chinese indigenous breeds and Western commercial breeds. The expression of BMPR1B in ovaries of adult BMPR1B-SINE+/+-RIP9 Sushan pigs was found to be significantly higher in comparison to those of BMPR1B-SINE−/−-RIP9 (p < 0.05). SINE insertion of BMPR1B-SINE-RIP9 and LINE insertion of BMPR1B-LINE-RIP13 were observed to significantly increase the activity of Octamer binding transcription factor 4 (OCT4) minipromoter in CHO and C2C12 cells (p < 0.01). Therefore, these two RIPs could serve as useful molecular markers for modulating the growth or reproductive traits in assisted selection of pig breeding, while the mechanisms of the insertion function should be studied further.

Genome‐wide analysis of hybridization in wild boar populations reveals adaptive introgression from domestic pig

The admixture of domestic pig into French wild boar populations has been monitored since the 1980s thanks to the existence of a cytogenetic difference between the two sub‐species. The number of chromosomes is 2n = 36 in wild boar and 2n = 38 in pig, respectively. This difference makes it possible to assign the “hybrid” status to wild boar individuals controlled with 37 or 38 chromosomes. However, it does not make it possible to determine the timing of the hybridization(s), nor to guarantee the absence of domestic admixture in an animal with 2n = 36 chromosomes. In order to analyze hybridization in greater detail and to avoid the inherent limitations of the cytogenetic approach, 362 wild boars (WB) recently collected in different French geographical areas and in different environments (farms, free ranging in protected or unprotected areas, animals with 2n = 36, 37 or 38 chromosomes) were genotyped on a 70K SNP chip. Principal component analyses allowed the identification of 13 “outliers” (3.6%), for which the proportion of the genome of “domestic” origin was greater than 40% (Admixture analyses). These animals were probably recent hybrids, having Asian domestic pig ancestry for most of them. For the remaining 349 animals studied, the proportion of the genome of “wild” origin varied between 83% and 100% (median: 94%). This proportion varied significantly depending on how the wild boar populations were managed. Local ancestry analyses revealed adaptive introgression from domestic pig, suggesting a critical role of genetic admixture in improving the fitness and population growth of WB. Overall, our results show that the methods used to monitor the domestic genetic contributions to wild boar populations should evolve in order to limit the level of admixture between the two gene pools.

Identification of Genomic Regions and Candidate Genes for Litter Traits in French Large White Pigs Using Genome-Wide Association Studies

The reproductive traits of sows are one of the important economic traits in pig production, and their performance directly affects the economic benefits of the entire pig industry. In this study, a total of 895 French Large White pigs were genotyped by GeneSeek Porcine 50K SNP Beadchip and four phenotypic traits of 1407 pigs were recorded, including total number born (TNB), number born alive (NBA), number healthy piglets (NHP) and litter weight born alive (LWB). To identify genomic regions and genes for these traits, we used two approaches: a single-locus genome-wide association study (GWAS) and a single-step GWAS (ssGWAS). Overall, a total of five SNPs and 36 genomic regions were identified by single-locus GWAS and ssGWAS, respectively. Notably, fourof all five significant SNPs were located in 10.72–11.06 Mb on chromosome 7, were also identified by ssGWAS. These regions explained the highest or second highest genetic variance in the TNB, NBA and NHP traits and harbor the protein coding gene ENSSSCG00000042180. In addition, several candidate genes associated with litter traits were identified, including JARID2, PDIA6, FLRT2 and DICER1. Overall, these novel results reflect the polygenic genetic architecture of the litter traits and provide a theoretical reference for the following implementation of molecular breeding.

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.

Discovery of Differentially Expressed MicroRNAs in Porcine Ovaries With Smaller and Larger Litter Size

The number of live births in a litter is an important reproductive trait, and is one of the main indicators which reflect the production level and economic benefit of a pig farm. The ovary is an important reproductive organ of the sow, and it undergoes a series of biological processes during each estrous cycle. A complex transcriptional network containing coding and non-coding RNAs in the ovary closely regulates the reproductive capability of sows. However, the molecular regulation mechanisms affecting sow litter size are still unclear. We investigated the expression profiles of microRNAs (miRNAs) in porcine ovaries from sows with smaller than average litter sizes (SLS) and those with larger litter sizes (LLS). In total, 411 miRNAs were identified, and of these 17 were significantly down-regulated and 16 miRNAs were up-regulated when comparing sows with LLS and SLS, respectively. We further characterized the role of miR-183 which was one of the most up-regulated miRNAs. CCK-8, EdU incorporation and western blotting assays demonstrated that miR-183 promoted the proliferation of granulosa cells (GCs) in pig ovaries. Moreover, miR-183 inhibited the synthesis of estradiol in GCs and promoted the synthesis of progesterone. These results will help in gaining understanding of the role of miRNAs in regulating porcine litter size.

Differences in gene expression and variable splicing events of ovaries between large and small litter size in Chinese Xiang pigs

BACKGROUND

Although lots of quantitative trait loci (QTLs) and genes present roles in litter size of some breeds, the information might not make it clear for the huge diversity of reproductive capability in pig breeds. To elucidate the inherent mechanisms of heterogeneity of reproductive capability in litter size of Xiang pig, we performed transcriptome analysis for the expression profile in ovaries using RNA-seq method.

RESULTS

We identified 1,419 up-regulated and 1,376 down-regulated genes in Xiang pigs with large litter size. Among them, 1,010 differentially expressed genes (DEGs) were differently spliced between two groups with large or small litter sizes. Based on GO and KEGG analysis, numerous members of genes were gathered in ovarian steroidogenesis, steroid biosynthesis, oocyte maturation and reproduction processes.

CONCLUSIONS

Combined with gene biological function, twelve genes were found out that might be related with the reproductive capability of Xiang pig, of which, eleven genes were recognized as hub genes. These genes may play a role in promoting litter size by elevating steroid and peptide hormones supply through the ovary and facilitating the processes of ovulation and in vivo fertilization.

FoxA2 and p53 regulate the transcription of HSD17B1 in ovarian granulosa cells of pigs

The oestrogens have been highly implicated in the fertility of female animals. It is widely known that the oestrogens are primarily synthetized by the ovarian granulosa cells (GCs), and the final and essential step of this process is to catalyse the oestrone to the more active oestradiol by the protein coded by hydroxysteroid 17-beta dehydrogenase 1 (HSD17B1) gene. However, the molecular mechanism regarding the transcription of HSD17B1 remains to be fully elucidated in ovarian GCs. In this study, the 5'-deletion, luciferase assay and chromatin immunoprecipitation (ChIP) were utilized to explore the molecular regulation of transcription of HSD17B1 with the porcine ovarian GCs as the cellular model. After the deletions with -2105 to -1754 bp, -1753 to -1429 bp, -1430 to -1081 bp and -1082 to -730 bp, the relative luciferase activity of HSD17B1 promoter did not change significantly, but the deletion of -731 to -332 bp significantly increased the relative luciferase activity of HSD17B1 promoter, and an insertion (GTTT) that might raise the transcription of HSD17B1 was identified at -401 bp of HSD17B1. These findings suggested the region from -731 to +38 bp was the core promoter of HSD17B1, and the region between -731 to -332 bp might be a silence element for HSD17B1. Furthermore, the forkhead box A2 (FoxA2) directly bound at -412 to -401 bp to negatively but p53 bound at -383 to -374 bp to positively regulate the transcription and translation of HSD17B1 in ovarian GCs. These findings will improve our understanding on HSD17B1-mediated oestrogens and provide useful information for further investigations into fertility of females.

Survey of SNPs Associated with Total Number Born and Total Number Born Alive in Pig

Reproductive productivity depend on a complex set of characteristics. The number of piglets at birth (Total number born, Litter size, TNB) and the number of alive piglets at birth (Total number born alive, NBA) are the main indicators of the reproductive productivity of sows in pig breeding. Great hopes are pinned on GWAS (Genome-Wide Association Studies) to solve the problems associated with studying the genetic architecture of reproductive traits of pigs. This paper provides an overview of international studies on SNP (Single nucleotide polymorphism) associated with TNB and NBA in pigs presented in PigQTLdb as "Genome map association". Currently on the base of Genome map association results 306 SNPs associated with TNB (218 SNPs) and NBA (88 SNPs) have been identified and presented in the Pig QTLdb database. The results are based on research of pigs such as Large White, Yorkshire, Landrace, Berkshire, Duroc and Erhualian. The presented review shows that most SNPs found in chromosome areas where candidate genes or QTLs (Quantitative trait locus) have been identified. Further research in the given direction will allow to obtain new data that will become an impulse for creating breakthrough breeding technologies and increase the production efficiency in pig farming.

Porcine LIF gene polymorphisms and their association with litter size traits in four pig breeds

Leukemia inhibitory factor (LIF) is an important productivity-related gene in pigs. We found two polymorphisms — g.6646C>T and g.6988C>T — in exon 3 of LIF in pigs by using DNA sequencing and polymerase chain reaction-restriction fragment length polymorphism. Three genotypes were obtained and associated with litter size traits in Anqing Six-end-white (AQ), Wei (W), Wannan Black (WNB), and Large White (LW) pigs. At locus g.6646C>T, the g.6646C allele frequency variation was 0.6869 (AQ), 0.7473 (W), 1 (WNB), and 0.6852 (LW). In AQ pigs, sows with the TT genotype had higher total number of piglets born (TNB) and number of piglets born alive (NBA) in the first parity and multiparities (P < 0.01). In W and LW pigs, sows with the CC genotype had higher TNB and NBA in multiparities (P < 0.01). At locus g.6988C>T, the g.6988C allele frequency variation was 1 (AQ), 0.6154 (W), 1 (WNB), and 0.6667 (LW). The CC genotype significantly differed from CT or TT genotypes (P < 0.01) for TNB and NBA in W and LW pigs. Thus, LIF was shown to have a significant influence on litter size. Therefore, g.6646C>T and g.6988C>T loci of LIF could be potential marker-assisted selection tools for improving litter size in pig production.

Association of Rs339939442 in the AHR Gene with Litter Size are Inconsistent among Chinese Indigenous Pigs and Western Commercial Pigs

Eastern and Southern Chinese pigs have been imported to Western countries to improve economic traits including fertility in Western pig breeds by intensive selecting Chinese advantage genes. It was reported that the selected Asian-derived non-synonymous mutations including rs339939442 (G > T) in the aryl hydrocarbon receptor (AHR) gene could increase litter size in multiple European commercial lines. The objective of this study is to identify whether rs339939442 in the AHR gene is polymorphic and has an influence on the litter size in 10 pig populations including five Chinese indigenous breeds, one cultivated breed, one lean-type breed, two North American lean-type breeds, and one European lean-type breed. We found that rs339939442 had polymorphism in all 10 populations, whereas rs339939442 was associated with litter size only in French Yorkshire (FRA-Y) and Chinese cultivated Suhuai (SH) pigs containing approximately 75% British Yorkshire pigs ancestry. Our results indicated that rs339939442 in the AHR gene was a potential marker to improve litter size in European commercial lines and the pigs containing ancestries of European commercial lines, whereas this locus maybe not a causal mutation affecting the litter size but only in linkage disequilibrium with the causal mutation for litter size.

A haplotype variant of Hu sheep follicle-stimulating hormone receptor promoter region decreases transcriptional activity

In sheep, increased expression of the follicle-stimulating hormone receptor (FSHR) in the ovary is a common feature of ewes that carry the prolific allele. In this study, we demonstrated that polymorphisms in the core promoter of the FSHR gene are associated with the reproductive performance of Hu sheep and are involved in the transcriptional activity of FSHR. An approximately 1.5-kb region of the 5' flanking sequence of the Hu sheep FSHR gene was isolated and characterized, and its core promoter was located in the 5' regulatory region, from nucleotides -580 to -342. Four variants (c.-518T>C, c.-466C>T, c.-414A>G and c.-365C>T) were detected in this region, and six genotypes and three haplotypes were found in the Hu sheep population (n = 245). An association analysis revealed that these polymorphisms are associated with the litter size of Hu ewes. Furthermore, a luciferase assay showed that the T-C-A-C- and C-T-G-T-type core promoters have higher transcriptional activity than does the T-C-G-C type. Notably, the putative binding site for the transcription factor Yin Yang 1 (YY1) was present at the A allele of nucleotide -414, but YY1 can significantly increase the transcriptional activity of the FSHR core promoter, which contains three different haplotypes. Taken together, our results establish that these variants might be involved in regulating the transcriptional activity of FSHR and litter size in Hu ewes and may provide a novel candidate marker for marker-assisted selection in sheep breeding.

Transcriptome analysis of the endometrium from Chinese Erhualian sows that differ in calcium ion concentration and litter size

Embryonic survival rate, an important factor in the fecundity of sows, is affected by endometrium-secreting histotroph. A higher concentration of calcium ion has been observed in the uterus of highly prolific Erhualian sows (EH) compared with those of less prolific (EL) sows. This suggests that EH sows have better establishment and maintenance of pregnancies, thus increasing embryonic survival rate during the peri-implantation period. To understand the mechanisms of how the endometrium-secreting histotroph affects embryonic survival rate during the Erhualian peri-implantation period, the expression patterns of endometrial mRNA in the EH and EL sows on day 12 of gestation were analyzed using RNA sequencing technology. A total of 164 differentially expressed genes (DEGs) were identified (P_adj  < 0.05, |log₂ (FC)| ≥ 1), including 46 upregulated and 118 downregulated genes in EH compared to EL. Gene Ontology enrichment indicated that a subset of DEGs was involved in calcium ion binding and cell adhesion. Solute carrier family 8 member A3 and solute carrier family 24 member 4, identified as upregulated genes (P_adj  < 0.05) in EH, were considered key candidate genes expressed in the endometrium affecting embryonic survival rate during the peri-implantation period. The results improve understanding of the genetic mechanism underlying the variation in litter size of Erhualian pigs during the peri-implantation period.