The estrogen receptor locus is associated with a major gene influencing litter size in pigs

Genetic polymorphisms in ESR and FSHβ genes and their association with litter traits in Large White pigs

The estrogen receptor (ESR) gene and follicle-stimulating hormone β (FSHβ) gene are responsible for litter traits. The present study aimed to verify the polymorphisms of ESR and FSHβ and assess their effects on the litter traits in 201 Large White pigs. Four SNPs (g.C669T, g.A1296G, g.C1665T and g.A1755G) were found in ESR. The TT genotype at g.C1665T locus and AA genotype at g.A1755G locus could significantly increase the total litter size of the first litter of American Large White pigs (p < 0.05). Eight SNPs were found in exon 3 of FSHβ. The AA genotype at g.A511G locus, AA and AG genotypes at g.A617G locus, CC and CT genotypes at g.C630T locus, CT and TT genotypes at g.C652T locus, CT and TT genotypes at g.C735T locus, AA and AG genotypes at g.A746G, AA and AG genotypes at g.A921G and CT genotype at g.C678T could significantly increase the litter size of different strains of Large White pigs (p < 0.05). Our study revealed that the genetic variations of ESR and FSHβ were closely related to the litter trait of Large White pigs. Therefore, ESR and FSHβ genes could be used as molecular markers for the genetic selection of Large White pigs.

Effect of NLR family pyrin domain containing 9 gene polymorphism on litter size in large white pigs

NLR family pyrin domain containing 9 (NLRP9) is a mammalian reproduction-related gene. In this study, we researched the associations between polymorphisms located in the coding sequence (CDS) of the NLRP9 gene, and both the total number of piglets born per litter (TNB) and the number of piglets born alive per litter (NBA) in Canada Large White pigs (CLW). We detected a single nucleotide polymorphism (SNP) within exon 3 (g.10910C > T). The allele frequencies at the NLRP9 locus were 0.474 for the C allele and 0.526 for the T allele. Three genotypes, CC, CT, and TT, occurred with frequencies of 0.216, 0.515, and 0.269, respectively. Sows with the CC genotype had the largest TNB and NBA, sows with TT had the smallest, and those with CT were in-between. This difference was statistically significant (p < 0.05). Furthermore, CC females grew faster than CT or TT females, and there was a significant relationship between NLRP9 polymorphism and the average daily gain (p < 0.05). Here, we provide the first evidence for a novel SNP in NLRP9 associated with litter size in CLW sows, which could be used as a genetic marker to improve litter size in pig breeding and production.

Identification and characterization of circRNAs in peri-implantation endometrium between Yorkshire and Erhualian pigs

BACKGROUND

One of the most critical periods for the loss of pig embryos is the 12th day of gestation when implantation begins. Recent studies have shown that non-coding RNAs (ncRNAs) play important regulatory roles during pregnancy. Circular RNAs (circRNAs) are a kind of ubiquitously expressed ncRNAs that can directly regulate the binding proteins or regulate the expression of target genes by adsorbing micro RNAs (miRNA).

RESULTS

We used the Illumina Novaseq6,000 technology to analyze the circRNA expression profile in the endometrium of three Erhualian (EH12) and three Yorkshire (YK12) pigs on day 12 of gestation. Overall, a total of 22,108 circRNAs were identified. Of these, 4051 circRNAs were specific to EH12 and 5889 circRNAs were specific to YK12, indicating a high level of breed specificity. Further analysis showed that there were 641 significant differentially expressed circRNAs (SDEcircRNAs) in EH12 compared with YK12 (FDR < 0.05). Functional enrichment of differential circRNA host genes revealed many pathways and genes associated with reproduction and regulation of embryo development. Network analysis of circRNA-miRNA interactions further supported the idea that circRNAs act as sponges for miRNAs to regulate gene expression. The prediction of differential circRNA binding proteins further explored the potential regulatory pathways of circRNAs. Analysis of SDEcircRNAs suggested a possible reason for the difference in embryo survival between the two breeds at the peri-implantation stage.

CONCLUSIONS

Together, these data suggest that circRNAs are abundantly expressed in the endometrium during the peri-implantation period in pigs and are important regulators of related genes. The results of this study will help to further understand the differences in molecular pathways between the two breeds during the critical implantation period of pregnancy, and will help to provide insight into the molecular mechanisms that contribute to the establishment of pregnancy and embryo loss in pigs.

The transcription factor RUNX1 affects the maturation of porcine oocytes via the BMP15/TGF-β signaling pathway

Bone morphogenetic protein 15 (BMP15) is specifically expressed in oocytes in pigs at all stages from early stages to ovulation and has an important role in oocyte maturation. However, there are few reports on the molecular mechanisms by which BMP15 affects oocyte maturation. In this study, we identified the core promoter region of BMP15 using a dual luciferase activity assay and successfully predicted the DNA binding motif of the transcription factor RUNX1. The effect of BMP15 and RUNX1 on oocyte maturation was examined using the first polar body extrusion rate, a reactive oxygen species (ROS) assay and total glutathione (GSH) content at three time points of 12, 24 and 48 h of in vitro culture of porcine isolated oocytes. Subsequently, the effect of the transcription factor RUNX1 on the TGF-β signaling pathway (BMPR1B and ALK5) was further verified using RT-qPCR and Western blotting. We found that the overexpression of BMP15 significantly increased the first polar body extrusion rate (P < 0.01) and total glutathione content of oocytes cultured in vitro for 24 h and decreased reactive oxygen levels (P < 0.01), whereas interference with BMP15 decreased the first polar body extrusion rate (P < 0.01), increased reactive oxygen levels in oocytes cultured in vitro for 24 h (P < 0.01), and decreased glutathione content (P < 0.01). The dual luciferase activity assay and online software prediction showed that RUNX1 is a potential transcription factor binding to the core promoter region (-1203/-1423 bp) of BMP15. Overexpression of RUNX1 significantly increased the expression of BMP15 and oocyte maturation rate, while inhibition of RUNX1 decreased the expression of BMP15 and the oocyte maturation rate. Moreover, the expression of BMPR1B and ALK5 in the TGF-β signaling pathway increased significantly after overexpression of RUNX1, whereas their expression decreased after inhibition of RUNX1. Overall, our results suggest that the transcription factor RUNX1 positively regulates the expression of BMP15 and influences oocyte maturation through the TGF-β signaling pathway. This study provides a theoretical basis for further complementing the BMP15/TGF-β signaling pathway to regulate mammalian oocyte maturation.

Identification of the genetic basis of sow pelvic organ prolapse

Introduction: Pelvic organ prolapse (POP) is one contributor to recent increases in sow mortality that have been observed in some populations and environments, leading to financial losses and welfare concerns. Methods: With inconsistent previous reports, the objective here was to investigate the role of genetics on susceptibility to POP, using data on 30,429 purebred sows, of which 14,186 were genotyped (25K), collected from 2012 to 2022 in two US multiplier farms with a high POP incidence of 7.1% among culled and dead sows and ranging from 2% to 4% of all sows present by parity. Given the low incidence of POP for parities 1 and >6, only data from parities 2 to 6 were retained for analyses. Genetic analyses were conducted both across parities, using cull data (culled for POP versus another reason), and by parity, using farrowing data. (culled for POP versus culled for another reason or not culled). Results and Discussion: Estimates of heritability from univariate logit models on the underlying scale were 0.35 ± 0.02 for the across-parity analysis and ranged from 0.41 ± 0.03 in parity 2 to 0.15 ± 0.07 in parity 6 for the by-parity analyses. Estimates of genetic correlations of POP between parities based on bivariate linear models indicated a similar genetic basis of POP across parities but less similar with increasing distance between parities. Genome wide association analyses revealed six 1 Mb windows that explained more than 1% of the genetic variance in the across-parity data. Most regions were confirmed in several by-parity analyses. Functional analyses of the identified genomic regions showed a potential role of several genes on chromosomes 1, 3, 7, 10, 12, and 14 in susceptibility to POP, including the Estrogen Receptor gene. Gene set enrichment analyses showed that genomic regions that explained more variation for POP were enriched for several terms from custom transcriptome and gene ontology libraries. Conclusion: The influence of genetics on susceptibility to POP in this population and environment was confirmed and several candidate genes and biological processes were identified that can be targeted to better understand and mitigate the incidence of POP.

A unique 15-bp InDel in the first intron of BMPR1B regulates its expression in Taihu pigs

BACKGROUND

BMPR1B (Bone morphogenetic protein receptor type-1B) is a receptor in the bone morphogenetic protein (BMP) family and has been identified as a candidate gene for reproductive traits in pigs. Our previous study in Taihu pigs found a specific estrogen response element (ERE) in the first intron of the BMPR1B gene that is associated with the number born alive trait. However, little is known about the mechanism by which the ERE regulates the expression of BMPR1B in the endometrium.

RESULTS

Here, a 15-bp InDel (insertion/deletion) (AGCCAGAAAGGAGGA) was identified as a unique variation in Taihu pigs, and was shown to be responsible for the binding of the type I receptor of estrogen (ESR1) to the ERE using dual-luciferase assays. Four BMPR1B transcripts (T1, T2, T3, and T4) were identified by 5' RACE in endometrial tissue. Expression of T3 and T4 in the endometrium of Meishan pigs was significantly higher than in Duroc pigs during pregnancy. Luciferase assays showed that three distinct BMPR1B promoters may drive expression of T1, T3, and T4. Interestingly, ERE-mediated enhancement of T4 promoter activity significantly increased expression of Transcript T4 in the endometrium of Taihu pigs (P < 0.05). In contrast, the ERE inhibited activity of the T3 promoter and decreased expression of the T3 transcript in the Duroc background (P < 0.05). In summary, we identified a 15-bp InDel in the Taihu ERE that can be used as a molecular marker for the number born alive trait, characterized the 5' untranslated regions (UTRs) of BMPR1B transcripts in the endometrium, and determined how the transcripts are processed by alternative splicing events.

CONCLUSIONS

Our results provide a foundation for understanding the transcriptional regulation of BMPR1B and its contributions to the unique breeding prolificacy characteristics of Taihu pigs.

Exploring a Possible Link between the Fecal Microbiota and the Production Performance of Pigs

The backfat thickness of pigs not only affects the physical properties and taste of meat, but it also closely relates to the reproduction performance of sows. Accumulating evidence indicates that, apart from genetic factors, gut microbiota can also modulate the fat deposition and muscle growth. However, the differential microbiota in pigs with different backfat thickness, and whether microbiota affects backfat thickness, remains elusive. Firstly, 16S ribosomal RNA (16S rRNA) gene sequencing was performed on 62 fecal samples from pigs with different backfat thicknesses, and the compositions of microbiota among different groups with different backfat thicknesses were different. The abundance of Lactobacillus. reuteri (L. reuteri) and Prevotella sp RS2 was significantly higher in pigs with low-backfat thickness than that in pigs with middle and high-backfat thickness; meanwhile, the abundance of Desulfovibrio piger was significantly lower (p < 0.05) in pigs with low-backfat thickness. Furthermore, the functional profiling of microbial communities suggested that the abundance of isoquinoline alkaloid biosynthesis and styrene degradation were significantly lower (p < 0.05) in the low-backfat thickness group than that in middle and high-backfat thickness groups. Finally, L. reuteri fed to Meishan piglets was capable of improving the production performance and had the potential to reduce backfat thickness. This study provides new evidence that microbiota can regulate the phenotype of the host, and dietary supplementation with L. reuteri can improve the production performance of piglets.

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.

Challenging Sustainable and Innovative Technologies in Cheese Production: A Review

It is well known that cheese yield and quality are affected by animal genetics, milk quality (chemical, physical, and microbiological), production technology, and the type of rennet and dairy cultures used in production. Major differences in the same type of cheese (i.e., hard cheese) are caused by the rennet and dairy cultures, which affect the ripening process. This review aims to explore current technological advancements in animal genetics, methods for the isolation and production of rennet and dairy cultures, along with possible applications of microencapsulation in rennet and dairy culture production, as well as the challenge posed to current dairy technologies by the preservation of biodiversity. Based on the reviewed scientific literature, it can be concluded that innovative approaches and the described techniques can significantly improve cheese production.

Whole‐genome analysis reveals the hybrid formation of Chinese indigenous DHB pigs following human migration

Hybridization is widespread in nature and is a valuable tool in domestic breeding. The DHB (DaHuaBai) pig in South China is the product of such a breeding strategy, resulting in increased body weight compared with other pigs in the surrounding area. We analyzed genomic data from 20 Chinese pig breeds and investigated the genomic architecture after breed formation of DHB. The breed showed inconsistency in genotype and body weight phenotype, in line with selection after hybridization. By quantifying introgression with a haplotype‐based approach, we proposed a two‐step introgression from large‐sized pigs into small‐sized pigs to produce DHB, consistent with the human migration events in Chinese history. Combining with gene prioritization and allele frequency analysis, we identify candidate genes that showed selection after introgression and that may affect body weight, such as IGF1R, SRC, and PCM1. Our research provides an example of a hybrid formation of domestic breeds along with human migration patterns.

Allelic variability of estrogen receptor (ESR) gene and its effect on litter traits of Doom pigs

The present study was conducted to identify the polymorphism of estrogen receptor (ESR) gene and its biological association with litter traits (litter size at birth, litter size at weaning, litter weight at birth and litter weight at weaning) of Doom pigs native to Assam. A total of 50 adult pigs (12 males and 38 females) chosen randomly from three different herds (Herd I, Herd II, and Herd III) were utilized in the present study. Detection of polymorphism of ESR gene was done by means of PCR-RFLP method. The amplified PCR product was digested with Pvu II restriction enzyme. PCR-RFLP analysis of ESR gene revealed polymorphic banding pattern. Two genotypes viz. AA and AB were identified. AA genotype yielded one fragment (120 bp) and AB genotype yielded three fragments (120, 65, and 55 bp). BB genotype was not found in the population under study. The frequencies of A and B alleles were found to be 0.650 and 0.350, respectively, and the genotypic frequencies of ESR gene were found to be 0.300 and 0.700 for AA and AB genotypes, respectively. There was no significant (P > 0.05) effect of ESR genotype on litter traits and the population under study was not in Hardy-Weinberg Equilibrium for ESR gene. Clustal W Multiple alignment of partial sequence of ESR gene revealed single nucleotide changes at 33, 65, 70, 83, and 92^nd nucleotide positions. The presence of Pvu II polymorphism and identification of single nucleotide variation of ESR gene opens interesting prospects for improvement of litter traits in Doom pig through selective breeding program, especially based on marker-assisted selection.

Transcriptomic Data Analyses Reveal That Sow Fertility-Related lincRNA NORFA Is Essential for the Normal States and Functions of Granulosa Cells

NORFA, the first lincRNA associated with sow fertility, has been shown to control granulosa cell (GC) functions and follicular atresia. However, the underlying mechanism is not fully understood. In this study, RNA-seq was performed and we noticed that inhibition of NORFA led to dramatic transcriptomic alterations in porcine GCs. A total of 1,272 differentially expressed transcripts were identified, including 1167 DEmRNAs and 105 DEmiRNAs. Furthermore, protein–protein interaction, gene-pathway function, and TF–miRNA–mRNA regulatory networks were established and yielded four regulatory modules with multiple hub genes, such as AR, ATG5, BAK1, CENPE, NR5A1, NFIX, WNT5B, ssc-miR-27b, and ssc-miR-126. Functional assessment showed that these hub DEGs were mainly enriched in TGF-β, PI3K-Akt, FoxO, Wnt, MAPK, and ubiquitin pathways that are essential for GC states (apoptosis and proliferation) and functions (hormone secretion). In vitro, we also found that knockdown of NORFA in porcine GCs significantly induced cell apoptosis, impaired cell viability, and suppressed 17β-estradiol (E2) synthesis. Notably, four candidate genes for sow reproductive traits (INHBA, NCOA1, TGFβ-1, and TGFBR2) were also identified as potential targets of NORFA. These findings present a panoramic view of the transcriptome in NORFA-reduced GCs, highlighting that NORFA, a candidate lincRNA for sow fertility, is crucial for the normal states and functions of GCs.

Evaluation of candidate genes related to litter traits in Indian pig breeds

Improvement in litter traits is the key to profitable pig farming that directly enhances the economic standing of the farmers in developing countries. The present study aimed to explore oestrogen receptor (ESR), epidermal growth factor (EGF), follicle-stimulating hormone beta subunit (FSHβ), prolactin receptor (PRLR) and retinol-binding protein 4 (RBP4) genes as possible candidate genetic markers for litter traits in indigenous pigs of India. The breeds included in the study were Ghungroo, Mali, Niang Megha and Tenyi Vo, and the reproductive traits considered were litter size at birth (LSB), number born alive (NBA), litter weight at birth (LWB), litter size at weaning (LSW) and litter weight at weaning (LWW) at their first parity. PCR-RFLP and primer-based mutation detection methods were used to identify polymorphism, and associations between the genotypes and the traits were analysed using a general linear model. The Ghungroo pigs recorded the best litter performances among the breeds (p < .05, LWB p < .01). Different alleles and genotypes of the genes under study were detected. Short interspersed nuclear element (SINE) -/- genotype of FSHβ revealed significantly higher litter traits (p < .05, LSB p < .01). The LWW was also found to be significantly influenced by ESR BB and AB, EGF AB and BB, and PRLR CC genotypes (p < .05). Although we did not find statistically significant and consistently superior litter traits with respect to different genotypes of other studied genes than genotype SINE -/- of the FSHβ, PRLR CC genotype demonstrated superior performances for all the litter traits. Our study revealed the FSHβ as a potential candidate genetic marker for litter traits in indigenous pig breeds of India.

SMARCA2 is regulated by NORFA-miR-29c, a novel pathway that controls granulosa cell apoptosis and is related to female fertility

SMARCA2, an evolutionarily conserved catalytic ATPase subunit of SWI/SNF complexes, has been implicated in development and diseases; however, its role in mammalian ovarian function and female fertility is unknown. Here, we identified and characterized the 3'-UTR of the porcine SMARCA2 gene and identified a novel adenylate number variation. Notably, this mutation was significantly associated with sow litter size traits and SMARCA2 levels, due to its influence on the stability of SMARCA2 mRNA in ovarian granulosa cells (GCs). Immunohistochemistry and functional analysis showed that SMARCA2 is involved in the regulation of follicular atresia by inhibiting GC apoptosis. In addition, miR-29c, a pro-apoptotic factor, was identified as a functional miRNA that targets SMARCA2 in GCs and mediates regulation of SMARCA2 expression via the NORFA-SMAD4 axis. Although a potential miR-29c-responsive element was identified within NORFA, negative regulation of miR-29c expression by NORFA was not due to activity as a competing endogenous RNA. In conclusion, our findings demonstrate that SMARCA2 is a candidate gene for sow litter size traits, because it regulates follicular atresia and GC apoptosis. Additionally, we have defined a novel candidate pathway for sow fertility, the NORFA-TGFBR2-SMAD4-miR-29c-SMARCA2 pathway.This article has an associated First Person interview with the first author of the paper.

Sicilian Black Pig: An Overview

Simple Summary

The conservation of the genetic variability of animals used for food production and non-food raw materials and services is a problem of primary importance at a global level. In recent years, conservation of biodiversity in livestock species has been favoring the need to preserve genetic variability of the autochthonous breeds, exploiting them in the context of production systems. In this context, a precious genetic reserve is represented by autochthonous breeds used for the production of typical products used in Italian gastronomic traditions, of which some organoleptic properties of their meats that could disappear due to severe selection programs are being recovered. Currently, the survival of autochthonous breeds is linked to various reasons such as their rusticity, i.e., the adaptability to difficult environmental conditions, and to the higher market value of their productions obtained according to traditional methods compared to the industrial production types. As information on autochthonous Italian pigs is limited, further research aims at making better use of these breeds and at increasing the knowledge of their genetic variability.

Abstract

The Sicilian black pig (SB) (Nero Siciliano), also known as the Nero dei Nebrodi, Nero delle Madonie, or Nero dell’Etna pig ecotype, is an autochthonous Italian breed. The origins of this breed date back to Greek and Carthaginian dominations. In ancient times, its breeding was fairly common throughout Sicily, registering only a temporary reduction during the Arab domination. This breed is known primarily for its distinctive black coat, although some individuals display wattles and a partially or wholly white face. The SB pig has a birth rate with an average per sow of 7.6 piglets, each of 1.4 kg live body weight, showing an average daily gain (ADG) of 346 g/day during the fattening period. Slaughter generally takes place at an average age of 390 days, with an average live weight of 95 kg. This breed also appears to withstand adverse climatic conditions and resist disease. The purpose of this manuscript is to offer a general overview regarding the Sicilian Black pig and to consider the recent findings related to genome investigation. The recent application of Next Generation Sequencing (NGS) technologies in the study of the genome of autochthonous breeds showed that polymorphisms of some candidate genes for production performance and phenotypic traits represent important information for selection processes. The protection of autochthonous breeds, intended as sources of genomic diversity for the further improvements of pigs for commercial use, constitutes a valuable opportunity to create new sustainable pig chains.

Effects of follicle-stimulating hormone beta subunit and nuclear receptor coactivator 1 gene polymorphisms and expressions on pink-eyed white mink reproductive traits

The present study was designed to investigate comparative expressions of follicle-stimulating hormone beta subunit (FSHβ) and nuclear receptor coactivator 1 (NCOA1) genes by real-time polymerase chain reaction, using polymerase chain reaction single-strand conformation polymorphism methods to investigate the effects of gene polymorphisms on reproductive traits, including total number of kits born (TNB) and number of born alive (NBA) in pink-eyed white mink. Four single-nucleotide polymorphisms were identified in the FSHβ and NCOA1 genes. The g.1228G>A polymorphism of FSHβ was associated with NBA and TNB (P < 0.01). The g.151536T>C polymorphism of NCOA1 was associated with NBA and TNB (P < 0.01). The NCOA1 mRNA levels in hypothalamus, ovary, and uterus during the first half of gestation were higher than during the middle term and last half of gestation (P < 0.01). The FSHβ mRNA levels in the hypothalamus and uterus were higher during the first half of gestation than during the middle term and last half of gestation (P < 0.05). In conclusion, the g.1866T>C polymorphism of FSHβ and the g.151536T>C polymorphism of NCOA1 could be molecular markers for reproductive traits, and expressions of FSHβ and NCOA1 might be involved in the regulation of embryo attachment mechanisms in pink-eyed white mink breeding.

Identification of Ovarian Circular RNAs and Differential Expression Analysis between MeiShan and Large White Pigs

MeiShan and Large White pigs differ in their female fecundity. However, the mechanisms behind the gene expression and regulation that cause these differences remain unclear. In this study, we profiled circRNAs and identified 5,879 circRNAs from the ovaries of MeiShan and Large White pigs. Eighty-five circRNAs were differentially expressed between the two pig breeds. Of these, 37 were up-regulated and 48 were down-regulated in MeiShan pigs. Gene ontology enrichment analysis suggested that the differentially expressed circRNA were involved in the hormone-mediated signaling pathway. We verified that circSCIN and its parent gene, scinderin (SCIN), were differentially expressed by reverse transcription and quantitative PCR (RT-qPCR). Luciferase assays demonstrated that circSCIN can target and sponge miR-133 and miR-148a/b. The identification of differentially expressed circRNAs (DECs) and their regulatory functions increased our understanding of the differences in reproductive efficiency between MeiShan and Large White pigs.

Lifetime cancer prevalence and life history traits in mammals

Background

Cancer is a common diagnosis in many mammalian species, yet they vary in their vulnerability to cancer. The factors driving this variation are unknown, but life history theory offers potential explanations to why cancer defense mechanisms are not equal across species.

Methodology

Here we report the prevalence of neoplasia and malignancy in 37 mammalian species, representing 11 mammalian orders, using 42 years of well curated necropsy data from the San Diego Zoo and San Diego Zoo Safari Park. We collected data on life history components of these species and tested for associations between life history traits and both neoplasia and malignancy, while controlling for phylogenetic history.

Results

These results support Peto's paradox, in that we find no association between lifespan and/or body mass and the prevalence of neoplasia or malignancy. However, a positive relationship exists between litter size and prevalence of malignancy (P = 0.005, Adj. R2 = 0.212), suggesting that a species' life history strategy may influence cancer vulnerabilities. Lastly, we tested for the relationship between placental invasiveness and malignancy. We find no evidence for an association between placental depth and malignancy prevalence (P = 0.618, Adj. R2 = 0.068).

Conclusions

Life history theory offers a powerful framework to understand variation in cancer defenses across the tree of life. These findings provide insight into the relationship between life history traits and cancer vulnerabilities, which suggest a trade-off between reproduction and cancer defenses.

Lay summary

Why are some mammals more vulnerable to cancer than others? We test whether life history trade-offs may explain this variation in cancer risk. Bigger, longer-lived animals do not develop more cancer compared to smaller, shorter-lived animals. However, we find a positive association between litter size and cancer prevalence in mammals.

The Interaction of lncRNA XLOC-2222497, AKR1C1, and Progesterone in Porcine Endometrium and Pregnancy

The endometrium is an important tissue for pregnancy and plays an important role in reproduction. In this study, high-throughput transcriptome sequencing was performed in endometrium samples of Meishan and Yorkshire pigs on days 18 and 32 of pregnancy. Aldo-keto reductase family 1 member C1 (AKR1C1) was found to be a differentially expressed gene, and was identified by quantitative real-time PCR (qRT-PCR) and Western blot. Immunohistochemistry results revealed the cellular localization of the AKR1C1 protein in the endometrium. Luciferase activity assay demonstrated that the AKR1C1 core promoter region was located in the region from -706 to -564, containing two nuclear factor erythroid 2-related factor 2 (NRF2) binding sites (antioxidant response elements, AREs). XLOC-2222497 was identified as a nuclear long non-coding RNA (lncRNA) highly expressed in the endometrium. XLOC-2222497 overexpression and knockdown have an effect on the expression of AKR1C1. Endocrinologic measurement showed the difference in progesterone levels between Meishan and Yorkshire pigs. Progesterone treatment upregulated AKR1C1 and XLOC-2222497 expression in porcine endometrial epithelial cells. In conclusion, transcriptome analysis revealed differentially expressed transcripts during the early pregnancy process. Further experiments demonstrated the interaction of XLOC-2222497/AKR1C1/progesterone in the endometrium and provided new potential targets for pregnancy maintenance and its control.

Genome-Wide Analysis Reveals Human-Mediated Introgression from Western Pigs to Indigenous Chinese Breeds

Genetic variations introduced via introgression from Western to Chinese pigs have contributed to the performance of Chinese breeds in traits such as growth rate and feed conversion efficiency. However, little is known about the underlying genomic changes that occurred during introgression and the types of traits affected by introgression. To address these questions, 525 animals were characterized using an SNP array to detect genomic regions that had been introgressed from European to indigenous Chinese breeds. The functions of genes located in introgressed regions were also investigated. Our data show that five out of six indigenous Chinese breeds show evidence of introgression from Western pigs, and eight introgressed genome regions are shared by five of the Chinese breeds. A region located on chr13: 12.8-13.1 M was affected by both introgression and artificial selection, and this region contains the glucose absorption related gene, OXSM, and the sensory related gene, NGLY. The results provide a foundation for understanding introgression from Western to indigenous Chinese pigs.

Integrated Analysis of miRNA-mRNA Network Reveals Different Regulatory Patterns in the Endometrium of Meishan and Duroc Sows during Mid-Late Gestation

Simple Summary

Meishan pigs have a lower fetal loss rate during mid-late gestation compared to Duroc pigs. Differentially expressed mRNAs and miRNAs detected in endometrial tissue from Meishan and Duroc sows at mid-late gestation are involved in regulating hormone and oxygen levels, blood vessel development, and developmental processes affecting reproduction. In addition, ssc-miR-503 and ssc-miR-671-5p were shown to target the EGF and ESR1 genes, respectively. These findings provided an important resource for studying embryonic mortality during mid-late gestation in pigs.

Abstract

Embryo loss is a major factor affecting profitability in the pig industry. Embryonic mortality occurs during peri-implantation and mid-late gestation in pigs. Previous investigations have shown that the embryo loss rate in Meishan pigs is significantly lower than in commercial breeds. Most studies have focused on embryonic mortality during early gestation, but little is known about losses during mid-late gestation. In this study, we performed a transcriptome analysis of endometrial tissue in mid-late gestation sows (gestation days 49 and 72) sampled from two breeds (Meishan (MS) and Duroc (DU)) that have different embryo loss rates. We identified 411, 1113, 697, and 327 differentially expressed genes, and 14, 36, 57, and 43 differentially expressed miRNAs in four comparisons (DU49 vs. DU72, DU49 vs. MS49, DU72 vs. MS72, and MS49 vs. MS72), respectively. Subsequently; seven differentially expressed mRNAs and miRNAs were validated using qPCR. Functional analysis suggested the differentially expressed genes and miRNAs target genes mainly involved in regulation of hormone levels, blood vessel development, developmental process involved in reproduction, embryonic placenta development, and the immune system. A network analysis of potential miRNA-gene interactions revealed that differentially expressed miRNAs in Meishan pigs are involved in the response to estradiol and oxygen levels, and affect angiogenesis and blood vessel development. The binding site on ssc-miR-503 for epidermal growth factor (EGF) and the binding site on ssc-miR-671-5p for estrogen receptor α (ESR1) were identified using a dual luciferase assay. The results of this study will enable further exploration of miRNA-mRNA interactions important in pig pregnancy and will help to uncover molecular mechanisms affecting embryonic mortality in pigs during mid-late gestation.

Relationship between mitochondrial DNA haplogroup and litter size in the pig

Mitochondrial DNA (mtDNA) has been widely associated with complex traits in farm animals. The present study evaluated the effects of mtDNA on litter size in pigs. Mitogenome sequencing of 1017 sows distinguished 232 variations, including 229 single nucleotide polymorphisms and three indels, which constituted 11 haplotypes and further clustered into two haplogroups that differed significantly (P&lt;0.05) in litter size. In order to explain the associations between the effect of haplogroup on litter size and different maternal origins, extant mitogenome sequences were used for phylogenetic or principal component analyses. The results of these analyses led to the identification of two groups, representing Chinese and European origins. The haplotypes corresponding to high litter size were all in the Chinese cluster, whereas haplotypes corresponding to low litter size were all in the European cluster. The results of this study suggest that the effect of haplogroup on litter size in the pig could be caused by diverse maternal origins, and that mtDNA haplogroup may be a marker for genetic selection for pig litter size.

Whole genome SNPs discovery in Nero Siciliano pig

Autochthonous pig breeds represent an important genetic reserve to be utilized mainly for the production of typical products. To explore its genetic variability, here we present for the first time whole genome sequencing data and SNPs discovered in a male domestic Nero Siciliano pig compared to the last pig reference genome Sus scrofa11.1.A total of 346.8 million paired reads were generated by sequencing. After quality control, 99.03% of the reads were mapped to the reference genome, and over 11 million variants were detected.Additionally, we evaluated sequence diversity in 21 fitness-related loci selected based on their biological function and/or their proximity to relevant QTLs. We focused on genes that have been related to environmental adaptation and reproductive traits in previous studies regarding local breeds. A total of 6,747 variants were identified resulting in a rate of 1 variant every ~276 bases. Among these variants 1,132 were novel to the dbSNP151 database. This study represents a first step in the genetic characterization of Nero Siciliano pig and also provides a platform for future comparative studies between this and other swine breeds.

A genome-wide association study of reproduction traits in four pig populations with different genetic backgrounds

Objective

Genome-wide association study and two meta-analysis based on GWAS performed to explore the genetic mechanism underlying variation in pig number born alive (NBA) and total number born (TNB).

Methods

Single trait GWAS and two meta-analysis (single-trait meta analysis and multi-trait meta analysis) were used in our study for NBA and TNB on 3,121 Yorkshires from 4 populations, including three different American Yorkshire populations (n = 2,247) and one British Yorkshire populations (n = 874).

Results

The result of single trait GWAS showed that no significant associated single nucleotide polymorphisms (SNPs) were identified. Using single-trait meta analysis and multi-trait meta analysis within populations, 11 significant loci were identified associated with target traits. Spindlin 1, vascular endothelial growth factor A, forkhead box Q1, msh homeobox 1, and LHFPL tetraspan submily member 3 are five functionally plausible candidate genes for NBA and TNB. Compared to the single population GWAS, single-trait Meta analysis can improve the detection power to identify SNPs by integrating information of multiple populations. The multiple-trait analysis reduced the power to detect trait-specific loci but enhanced the power to identify the common loci across traits.

Conclusion

In total, our findings identified novel genes to be validated as candidates for NBA and TNB in pigs. Also, it enabled us to enlarge population size by including multiple populations with different genetic backgrounds and increase the power of GWAS by using meta analysis.

Methods for reproductive tract scoring as a tool for improving sow productivity

Improving sow lifetime productivity is essential for maximizing farm profitability. Study objectives were to determine the accuracy for different vulva scoring methods in a commercial production system and to assess whether gilt reproductive tract scoring [evaluated by vulva width (VW)] prior to puberty could serve as useful gilt selection criteria. To accomplish this objective, 958 prepubertal replacement gilts in a commercial system were evaluated at approximately 15 wk of age. Gilt body weight (BW) was recorded in addition to 4 different methods to evaluate VW. Methods for VW assessment included digital caliper measurement (mm), visual evaluation and scoring by trained farm personnel [Farm Score (FS)], and 2 methods using scoring tools [Vulva Score Method A and B (VSA and VSB, respectively)] specifically calibrated from the VW distribution measured on gilts from previous studies. The VSA and FS methods assigned gilts to one of 3 categories (S, M, L, and 1, 2, 3, respectively) whereas VSB classified gilts vulvas using a 5-point scoring system (1 to 5). At 15 wk of age, a low proportion of variability in vulva size (27.8 ± 0.1 mm) could be explained by BW (62.2 ± 0.2 kg; R² = 0.05). All 3 scoring methods were effective in categorizing gilts based upon VW, as the measured VW size within methods differed by score (P < 0.01). The proportion of gilts achieving their first parity increased with score for VSA (64.7%, 73.2%, and 84.4%; P = 0.02), VSB (66.0%, 71.7%, 79.2%, 76.4%, and 84.2%; P = 0.02), and FS (67.2%, 75.0%, and 88.8%; P = 0.03), but VSA, VSB, and FS did not influence percentage of gilts achieving their second parity (P = 0.32, 0.29, and 0.30, respectively). Litter performance of gilts scored as M or L using VSA improved with an increased total born over 2 parities compared to those scored as S (23.96 vs. 26.38 pigs; P < 0.01) as well as born alive (21.13 vs. 23.05 pigs; P < 0.05). Results were similar for VSB, where scores 2 to 5 had greater total born (23.97 vs. 26.33 pigs; P < 0.01) and born alive (21.11 vs. 23.02 pigs; P < 0.05) through 2 parities compared to gilts scored 1. Using the FS method, total born pigs tended to be increased (P = 0.06) through 2 parities for gilts having a 2 or 3 vulva score compared to those scored as a 1. Collectively, assessing VW at approximately 15 wk of age may identify sows with improved productivity through 2 parities as breeding herd females.

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.

Genetic effect and combined genotype effect of ESR, FSHβ, CTNNAL1 and miR-27a loci on litter size in a Large White population

To select new Large White line with high number of piglets born, genotypes of estrogen receptor (ESR), the follicle stimulating hormone β subunit (FSHβ), catenin alpha like 1 (CTNNAL1) and miR-27a were tested in 472 Large White sows. The associations of different genotypes with litter size traits were also studied. The results showed ESR^BB and FSHβ^BB sows produced 0.41-1.49 more pigs per litter (p < .05) for total number born (TNB) and number born alive (NBA) than did other corresponding genotypes. TNB of CTNNAL1^CG sows is 0.50 more pigs per litter (p < .05) than that of CTNNAL1^GG sows with the dominance effect of 0.25 pigs per litter (p < .05). miR-27a^BB sows had a less estimated breeding value (EBV) to TNB and had a more number of mummified pigs (NM) than did miR-27a^AA or miR-27a^AB sows (p < .05). Therefore, ESR^B, FSHβ^B, CTNNAL1^G, miR-27a^A allele was favorable for litter size traits. Furthermore, combined genetic effect analysis showed ESR^AAFSHβ^BB, ESR^AACTNNAL1^CG, ESR^AAmiR-27a^AA, FSHβ^BBCTNNAL1^CC, FSHβ^BBmiR-27a^AA and CTNNAL1^CG miR-27a^AB was the favorable combined genotype for litter size traits. These results identified favorable alleles and genotypes for litter size traits and suggested a potential selection scheme for litter size in Large White pigs.

Whole-genome resequencing reveals candidate mutations for pig prolificacy

Changes in pig fertility have occurred as a result of domestication, but are not understood at the level of genetic variation. To identify variations potentially responsible for prolificacy, we sequenced the genomes of the highly prolific Taihu pig breed and four control breeds. Genes involved in embryogenesis and morphogenesis were targeted in the Taihu pig, consistent with the morphological differences observed between the Taihu pig and others during pregnancy. Additionally, excessive functional non-coding mutations have been specifically fixed or nearly fixed in the Taihu pig. We focused attention on an oestrogen response element (ERE) within the first intron of the bone morphogenetic protein receptor type-1B gene (BMPR1B) that overlaps with a known quantitative trait locus (QTL) for pig fecundity. Using 242 pigs from 30 different breeds, we confirmed that the genotype of the ERE was nearly fixed in the Taihu pig. ERE function was assessed by luciferase assays, examination of histological sections, chromatin immunoprecipitation, quantitative polymerase chain reactions, and western blots. The results suggest that the ERE may control pig prolificacy via the cis-regulation of BMPR1B expression. This study provides new insight into changes in reproductive performance and highlights the role of non-coding mutations in generating phenotypic diversity between breeds.

SLA-11 mutations are associated with litter size traits in Large White and Chinese DIV pigs

Litter size is an important economic traits in pigs. SLA-11 gene is a member of SLA (swine leukocyte antigen) complex. In our previous study, the SLA-11 gene was differentially expressed in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows. Here, we identified two mutations (c.754-132 T > C and c.1421 + 38 T > C) in SLA-11 gene and analyzed the associations of two SNPs with litter size traits in Large White (n = 263) and DIV (n = 117) sows. The results showed that in Large White pigs, SLA-11 c.754-132 CC sows produced 0.74 and 0.87 more pigs per litter for TNB and NBA of all parities than did TT sows (p < .05); In DIV pigs, SLA-11 c.754-132 CC sows produced 1.17 more pigs per litter for TNB of all parities than did TC sows (p < .05). In Large White pigs, SLA-11 c.1421 + 38 CC sows produced 0.9 more pigs per litter for TNB of all parities than did TT sows (p < .05), while in DIV pigs SLA-11 c.1421 + 38 CC sows produced 0.84 and 0.7 less pigs per litter for TNB and NBA of all parities than did TT sows (p < .05). Our research indicated that SLA-11 mutations were potential molecular markers for improving the litter size traits in pigs.

Genome-wide detection of selection signatures in Chinese indigenous Laiwu pigs revealed candidate genes regulating fat deposition in muscle

Background

Currently, genome-wide scans for positive selection signatures in commercial breed have been investigated. However, few studies have focused on selection footprints of indigenous breeds. Laiwu pig is an invaluable Chinese indigenous pig breed with extremely high proportion of intramuscular fat (IMF), and an excellent model to detect footprint as the result of natural and artificial selection for fat deposition in muscle.

Result

In this study, based on GeneSeek Genomic profiler Porcine HD data, three complementary methods, F_ST, iHS (integrated haplotype homozygosity score) and CLR (composite likelihood ratio), were implemented to detect selection signatures in the whole genome of Laiwu pigs. Totally, 175 candidate selected regions were obtained by at least two of the three methods, which covered 43.75 Mb genomic regions and corresponded to 1.79% of the genome sequence. Gene annotation of the selected regions revealed a list of functionally important genes for feed intake and fat deposition, reproduction, and immune response. Especially, in accordance to the phenotypic features of Laiwu pigs, among the candidate genes, we identified several genes, NPY1R, NPY5R, PIK3R1 and JAKMIP1, involved in the actions of two sets of neurons, which are central regulators in maintaining the balance between food intake and energy expenditure.

Conclusions

Our results identified a number of regions showing signatures of selection, as well as a list of functionally candidate genes with potential effect on phenotypic traits, especially fat deposition in muscle. Our findings provide insights into the mechanisms of artificial selection of fat deposition and further facilitate follow-up functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0622-y) contains supplementary material, which is available to authorized users.

Effects of pituitary-specific overexpression of FSHα/β on reproductive traits in transgenic boars

BACKGROUND

Follicle-stimulating hormone (FSH) is a gonadotropin synthesized and secreted by the pituitary gland. FSH stimulates follicle development and maturation in females. It also plays an important role in spermatogenesis in males, including humans and mice. However, the effects of FSH on male pigs are largely unknown. In this study, we generated transgenic pigs to investigate the effects of FSHα/β overexpression on reproductive traits in boars.

RESULTS

After five transgenic F0 founders were crossed with wide-type pigs, 193 F1 animals were obtained. Of these, 96 were confirmed as transgenic. FSHα and FSHβ mRNAs were detected only in pituitary tissue. Transgenic boars exhibited significantly higher levels of FSHα and FSHβ mRNA, serum FSH, and serum testosterone, compared to full-sib non-transgenic boars. Significant increases in testis weight, vas deferens diameter, seminiferous tubule diameter, and the number of Leydig cells were observed, suggesting that the exogenous FSHα/β affects reproductive traits. Finally, transgenic and non-transgenic boars had similar growth performance and biochemical profiles.

CONCLUSIONS

Pituitary-specific overexpression of FSHα/β genes is likely to impact reproductive traits positively, as indicated by enhancements in serum testosterone level, testis weight, the development of vas deferens, seminiferous tubules, and Leydig cells in transgenic boars. A high level of serum FSH induces secretion of serum testosterone, possibly by boosting the number of Leydig cells, which presumably increases the libido and the frequency of sexual activity in transgenic boars. Our study provides a preliminary foundation for the genetic improvement of reproductive traits in male pigs.

Whole-genome association and genome partitioning revealed variants and explained heritability for total number of teats in a Yorkshire pig population

Objective

The study was designed to perform a genome-wide association (GWA) and partitioning of genome using Illumina’s PorcineSNP60 Beadchip in order to identify variants and determine the explained heritability for the total number of teats in Yorkshire pig.

Methods

After screening with the following criteria: minor allele frequency, MAF≤0.01; Hardy-Weinberg equilibrium, HWE≤0.000001, a pair-wise genomic relationship matrix was produced using 42,953 single nucleotide polymorphisms (SNPs). A genome-wide mixed linear model-based association analysis (MLMA) was conducted. And for estimating the explained heritability with genome- or chromosome-wide SNPs the genetic relatedness estimation through maximum likelihood approach was used in our study.

Results

The MLMA analysis and false discovery rate p-values identified three significant SNPs on two different chromosomes (rs81476910 and rs81405825 on SSC8; rs81332615 on SSC13) for total number of teats. Besides, we estimated that 30% of variance could be explained by all of the common SNPs on the autosomal chromosomes for the trait. The maximum amount of heritability obtained by partitioning the genome were 0.22±0.05, 0.16±0.05, 0.10±0.03 and 0.08±0.03 on SSC7, SSC13, SSC1, and SSC8, respectively. Of them, SSC7 explained the amount of estimated heritability along with a SNP (rs80805264) identified by genome-wide association studies at the empirical p value significance level of 2.35E-05 in our study. Interestingly, rs80805264 was found in a nearby quantitative trait loci (QTL) on SSC7 for the teat number trait as identified in a recent study. Moreover, all other significant SNPs were found within and/or close to some QTLs related to ovary weight, total number of born alive and age at puberty in pigs.

Conclusion

The SNPs we identified unquestionably represent some of the important QTL regions as well as genes of interest in the genome for various physiological functions responsible for reproduction in pigs.

Genome-wide SNP data unveils the globalization of domesticated pigs

BACKGROUND

Pigs were domesticated independently in Eastern and Western Eurasia early during the agricultural revolution, and have since been transported and traded across the globe. Here, we present a worldwide survey on 60K genome-wide single nucleotide polymorphism (SNP) data for 2093 pigs, including 1839 domestic pigs representing 122 local and commercial breeds, 215 wild boars, and 39 out-group suids, from Asia, Europe, America, Oceania and Africa. The aim of this study was to infer global patterns in pig domestication and diversity related to demography, migration, and selection.

RESULTS

A deep phylogeographic division reflects the dichotomy between early domestication centers. In the core Eastern and Western domestication regions, Chinese pigs show differentiation between breeds due to geographic isolation, whereas this is less pronounced in European pigs. The inferred European origin of pigs in the Americas, Africa, and Australia reflects European expansion during the sixteenth to nineteenth centuries. Human-mediated introgression, which is due, in particular, to importing Chinese pigs into the UK during the eighteenth and nineteenth centuries, played an important role in the formation of modern pig breeds. Inbreeding levels vary markedly between populations, from almost no runs of homozygosity (ROH) in a number of Asian wild boar populations, to up to 20% of the genome covered by ROH in a number of Southern European breeds. Commercial populations show moderate ROH statistics. For domesticated pigs and wild boars in Asia and Europe, we identified highly differentiated loci that include candidate genes related to muscle and body development, central nervous system, reproduction, and energy balance, which are putatively under artificial selection.

CONCLUSIONS

Key events related to domestication, dispersal, and mixing of pigs from different regions are reflected in the 60K SNP data, including the globalization that has recently become full circle since Chinese pig breeders in the past decades started selecting Western breeds to improve local Chinese pigs. Furthermore, signatures of ongoing and past selection, acting at different times and on different genetic backgrounds, enhance our insight in the mechanism of domestication and selection. The global diversity statistics presented here highlight concerns for maintaining agrodiversity, but also provide a necessary framework for directing genetic conservation.

Genome-wide genetic structure and differentially selected regions among Landrace, Erhualian, and Meishan pigs using specific-locus amplified fragment sequencing

As typical Chinese indigenous pig breeds, Erhualian and Meishan have been widely used to produce new strain or breed in the world. However, the genetic basis of characteristics of these pig breeds is still limited. Moreover, considering cost and output of sequencing, it is necessary to further develop cost-effective method for pig genome screening. To contribute on this issue, we developed a SLAF-seq (specific-locus amplified fragment sequencing) method for pigs and applied it to analyze the genetic difference among Landrace, Erhualian, and Meishan pigs. A total of 453.75 million reads were produced by SLAF-seq. After quality-control, 165,670 SNPs (single nucleotide polymorphisms) were used in further analysis. The results showed that Landrace had distinct genetic relationship compared to Erhualian (F_ST = 0.5480) and Meishan (F_ST = 0.5800), respectively, while Erhualian and Meishan held the relatively close genetic relationship (F_ST = 0.2335). Furthermore, a genome-wide scanning revealed 268 differentially selected regions (DSRs) with 855 genes and 256 DSRs with 347 genes between Landrace and the two Chinese indigenous pig breeds and between Erhualian and Meishan, respectively. This study provides a new cost-effective method for pig genome study and might contribute to a better understanding on the formation mechanism of genetic difference among pigs with different geographical origins.

A novel synonymous SNP (A47A) of the <i>TMEM95</i> gene is significantly associated with the reproductive traits related to testis in male piglets

Abstract. Transmembrane protein 95 (TMEM95) is located on the acrosomal membrane of the sperm head involved in the acrosome reaction; thus, it is regarded as affecting spermatogenesis and reproduction traits. The aim of this study was to explore the novel single nucleotide polymorphisms (SNPs) within the pig TMEM95 gene as well as to evaluate their associations with the testicular sizes in male Landrace (LD) and Large White (LW) breeds. After pool sequencing and bioinformatics analysis, only one novel coding SNP was found in exon 1, namely NC_010454.3: g.341T > C, resulting in a synonymous mutation (A47A). This SNP could be genotyped using the StuI polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) assay. The minor allelic frequencies (MAFs) were 0.259 and 0.480 in the LD and LW breeds. Their polymorphism information content (PIC) values were 0.310 and 0.375. The LW population was at the Hardy–Weinberg equilibrium (HWE) (p > 0.05), whereas the LD population was not (p < 0.05). Association analyses demonstrated that a significant relationship was found between this A47A polymorphism and testis weight at 40 days of age in the LW population (p  =  0.047), and the heterozygote individuals showed lower testis weight than those with other genotypes. Moreover, this SNP was significantly associated with three testis measurement traits at 15 days of age in the LW population (p < 0.05); the individuals with genotypes TT and TC showed consistently superior testis measurement traits than those with genotype CC. These findings demonstrate that the A47A polymorphism had a significant effect on testis measurement traits, suggesting that the TMEM95 gene could be a candidate gene associated with reproductive traits. These results could contribute to breeding and genetics programs in the pig industry via DNA marker-assisted selection (MAS).

Genomics and metabolomics of post-weaning return to estrus

The weaning-to-estrus interval is a multifaceted trait that has the potential to substantially improve production efficiency in today's global swine industry, if variation in this measure can be reduced. Systems-biology approaches should help close the knowledge gap and increase selection tools and management strategies-such as gilt development programs, farrowing, and lactation feeding programs-to decrease the weaning-to-estrus interval. Metabolomics, the study of small compounds within biofluids and tissues, provides links between genotype and phenotype. Given the complexity and influence of the environment on the weaning-to-estrus interval, incorporating metabolomics data will provide valuable insight and guidance for future physiological as well as genetic and genomic strategies to reduce this interval, thereby improving sow productivity.

Individual variation in sow aggressive behavior and its relationship with sow welfare

This study examined the relationships between individual sow aggressive behavior and sow welfare, based on aggression, skin injuries, and stress, in a total of 275 pregnant domestic sows. Over 4 time replicates, sows were randomly mixed into groups of 10 (floor space of 1.8 m/sow) within 7 d of insemination in both their first and second gestations (200 sows per gestation with 126 sows observed in both gestations). Measurements were taken on aggression (both delivered and received) at feeding, skin injuries, and plasma cortisol concentrations at d 2, 9, and 51 after mixing. Live weight gain, nonreproductive removals, litter size (born alive, total born, and stillborn piglets), and farrowing rate were also recorded. In both the first and the second gestations, sows were classified at d 2 after mixing as "submissive" (delivered little or no aggression at feeding relative to aggression received), "subdominant" (received more aggression at feeding than delivered), and "dominant" (delivered more aggression at feeding than received). In both gestations, sows classified as dominant at d 2 subsequently delivered more (gestation 1, < 0.01; gestation 2, < 0.01) and received less (gestation 1, < 0.01; gestation 2, < 0.01) aggression and gained the most weight (gestation 1, < 0.01; gestation 2, < 0.01). Dominant sows had the least skin injuries throughout gestation 1 ( = 0.04), and although submissive sows sustained the most skin injuries at d 9 and 51 of gestation 2, at d 2 the classifications did not differ in skin injuries ( < 0.01). Subdominant sows had the highest cortisol concentrations at d 2 of gestation 2, but there were no differences between classifications at d 9 and 51 in either gestation (gestation 1, > 0.05; gestation 2, = 0.02). There were no significant relationships between aggression classification and reproduction and nonreproductive removals ( > 0.05). In conclusion, sows classified as dominant at feeding at d 2 subsequently received less aggression at feeding, sustained fewer skin injuries, and had higher live weight gain. Submissive and subdominant sows in groups are likely to benefit from the provision of increased resources such as space and access to feed.

On the efficiency of marker-assisted introgression

The efficiency of marker-assisted introgression programmes, expressed as genetic lag relative to a commercial population under continuous selection, was investigated using analytical methods. A genetic model was assumed for which the genetic variance in the introgression population was a function of the within-breed genetic variance and the initial breed difference. It was found that most of the genetic lag occurs in the latter stages of an introgression programme, when males and females which are heterozygous for the alíele to be introgressed are mated to produce homozygous individuals. Reducing genetic lag through selection on genomie proportion by using genetic markers throughout the genome, i.e. by selecting heterozygous individuals which resemble the recipient (commercial) population most, was effective if the initial breed difference was very large (e.g. 20 within-breed phenotypic standard deviations). In that case, selection solely on genetic markers could be practised to speed up genome recovery of the commercial line. If the initial breed difference is small, phenotypic or best linear unbiased prediction (BLUP) selection is superior in reducing genetic lag under the assumed genetic model.

Animal breeding in the (post-) genomic era

One of the benefits of the genomics revolution for animal production will be knowledge of genes that can be used to select more profitable livestock. Although it is possible to use genetic markers linked to genes of economic importance, tests for the genes themselves will be much more successful. Consequently finding genes of economic importance to livestock will be a major research aim for the future. Most traits of economic importance are quantitative traits affected by many genes. Mutations at many genes (e.g. 500) and at many positions within a gene (e.g. 1000 coding and non-coding bases) can affect a typical quantitative trait. The effect of these mutations on phenotype is usually small (e.g. 0·1 standard deviation) but occasionally large. Many mutations are lost from the population through genetic drift and selection, so that polymorphisms exist at only a subset of the relevant genes (e.g. 100 genes). Finding these genes, that have relatively small effects, is more difficult than finding genes for a classical Mendellian trait but, as the genomic tools become more powerful, it is becoming feasible and some successes have already occurred. The standard approach is to map a quantitative trait loci (QTL) to a chromosome region using linkage and linkage disequilibrium. Then test polymorphisms in positional candidate genes for an effect on the trait. Tools such as genomic sequence, EST collections and comparative maps make this approach feasible. Candidate genes can be selected based on functional data such as gene expression obtained from microarrays. At present the gain in rate of genetic improvement from use of DNA-based tests for QTL is small, because selection without them is already quite accurate, not enough QTL have been identified and genotyping is too expensive. However, in the future, with many QTL identified and inexpensive genotyping combined with decreased generation intervals, large gains are possible.

Genome-wide association analyses using a Bayesian approach for litter size and piglet mortality in Danish Landrace and Yorkshire pigs

BACKGROUND

Litter size and piglet mortality are important traits in pig production. The study aimed to identify quantitative trait loci (QTL) for litter size and mortality traits, including total number of piglets born (TNB), litter size at day 5 (LS5) and mortality rate before day 5 (MORT) in Danish Landrace and Yorkshire pigs by genome-wide association studies (GWAS).

METHODS

The phenotypic records and genotypes were available in 5,977 Landrace pigs and 6,000 Yorkshire pigs born from 1998 to 2014. A linear mixed model (LM) with a single SNP regression and a Bayesian mixture model (BM) including effects of all SNPs simultaneously were used for GWAS to detect significant QTL association. The response variable used in the GWAS was corrected phenotypic value which was obtained by adjusting original observations for non-genetic effects. For BM, the QTL region was determined by using a novel post-Gibbs analysis based on the posterior mixture probability.

RESULTS

The detected association patterns from LM and BM models were generally similar. However, BM gave more distinct detection signals than LM. The clearer peaks from BM indicated that the BM model has an advantage in respect of identifying and distinguishing regions of putative QTL. Using BM and QTL region analysis, for the three traits and two breeds a total of 15 QTL regions were identified on SSC1, 2, 3, 6, 7, 9, 13 and 14. Among these QTL regions, 6 regions located on SSC2, 3, 6, 7 and 13 were associated with more than one trait.

CONCLUSION

This study detected QTL regions associated with litter size and piglet mortality traits in Danish pigs using a novel approach of post-Gibbs analysis based on posterior mixture probability. All of the detected QTL regions overlapped with regions previously reported for reproduction traits. The regions commonly detected in different traits and breeds could be resources for multi-trait and across-bred selection. The proposed novel QTL region analysis method would be a good alternative to detect and define QTL regions.

Detection of one large insertion/deletion (indel) and two novel SNPs within the <i>SPEF2</i> gene and their associations with male piglet reproduction traits

Abstract. The sperm flagella 2 (SPEF2) gene is essential for normal sperm tail development and male fertility. To fully characterize the structure of the mutation and to further study the function of the pig SPEF2 gene, we explored the insertion/deletion (indel) and novel single-nucleotide polymorphisms (SNPs) within the pig SPEF2 gene, and tested their associations with the testicular sizes in male Large White (LW) and Landrace (LD) pigs from China. Herein, a large insertion located at the SPEF2 gene in chromosome 16 was found, and two alleles of "I" (insertion) and "D" (deletion) were designated. Allele "D" was dominant in all analyzed pigs. Two novel SNPs (namely (NC_010458) g.19642G > A, resulting in AfaI aCRS PCR–PFLP, and g.19886C > G, resulting in EcoRI aCRS PCR–PFLP) were found in LW and LD pigs. Association testing revealed that g.19886C > G was significantly associated with the testis long circumference (TLC) in LW pigs (P < 0.05), suggesting that this SNP would be the DNA marker for the marker-assisted selection (MAS) in reproduction traits. This preliminary result indicates that the pig SPEF2 gene had significant effects on male reproduction traits. These findings could not only extend the spectrum of genetic variations in the pig SPEF2 gene but also contribute to implementing MAS in genetics and breeding in pigs.

Expression of genes involved in progesterone receptor paracrine signaling and their effect on litter size in pigs

Background

Embryonic mortality during the period of implantation strongly affects litter size in pigs. Progesterone receptor (PGR) paracrine signaling has been recognized to play a significant role in embryonic implantation. IHH, NR2F2, BMP2, FKBP4 and HAND2 were proved to involve in PGR paracrine signaling. The objective of this study was to evaluate the expression of IHH, NR2F2, BMP2, FKBP4 and HAND2 in endometrium of pregnant sows and to further investigate these genes’ effect on litter size in pigs. Real-time PCR, western blot and immunostaining were used to study target genes/proteins expression in endometrium in pigs. RFLP-PCR was used to detect single nucleotide polymorphisms (SNPs) of target genes.

Results

The results showed that the mRNA and protein expression levels of IHH, NR2F2 and BMP2 were up-regulated during implantation period (P < 0.05 or P < 0.01). All target proteins were mainly observed in luminal epithelium and glandular epithelium. Interestingly, the staining of NR2F2 and HAND2 was also strong in stroma. SNPs detection revealed that there was a -204C > A mutation in promoter region of NR2F2 gene. Three genotypes were found in Large White, Landrace and Duroc sows. A total of 1847 litter records from 625 sows genotyped at NR2F2 gene were used to analyze the total number born (TNB) and number born alive (NBA). The study of the effect on litter size suggested that sows with genotype CC tend to have higher litter size.

Conclusions

These results showed the expression patterns of genes/proteins involved in PGR paracrine signaling over implantation time. And the candidate gene for litter size was identified from genes involved in this signaling. This study could be a resource for further studies to identify the roles of these genes for embryonic implantation in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s40104-016-0090-z) contains supplementary material, which is available to authorized users.

Identification of Differentially Expressed Genes Associated with Litter Size in Berkshire Pig Placenta

Improvement in litter size has become of great interest in the pig industry because fecundity is directly related to sow reproductive life. Improved reproduction has thus been achieved by elucidating the molecular functions of genes associated with fecundity. In the present study, we identified differentially expressed genes (DEGs) via transcriptomic analysis using RNA-sequencing (RNA-Seq) in Berkshire pig placentas from larger (LLG, mean litter size >12) and smaller (SLG, mean litter size < 6.5) litter size groups. In total 588 DEGs were identified (p < 0.05, > 1.5-fold change), of which 98 were upregulated, while 490 were downregulated in the LLG compared with the SLG. Gene Ontology (GO) enrichment was also performed. We concluded that 129 of the 588 DEGs were closely related to litter size according to reproduction related genes selected based on previous reports, as 110 genes were downregulated and 19 upregulated in the LLG compared with the SLG. RT-qPCR utilizing specific primers targeting the early growth response 2 (EGR2), pheromaxein c subunit (PHEROC) and endothelial lipase (LIPG) genes showed high accordance with RNA-Seq results. Furthermore, we investigated the upstream regulators of these three genes in the placenta. We found that WNT9B, a Wnt signaling pathway molecule, and IL-6, known inducers of EGR2 and LIPG, respectively, were significantly increased in LLG compared with SLG. We believe that the induction of IL-6 and LIPG may play an important role in increasing nutrition supply through the placenta from the sow to the piglet during gestation. These results provide novel molecular insights into pig reproduction.