Porcine SPP1 gene polymorphism association with phenotypic traits in the Landrace × Jeju (Korea) black pig F2 population

Genome-Wide Association Study for Weight Loss at the End of Dry-Curing of Hams Produced from Purebred Heavy Pigs

Dissecting the genetics of production traits in livestock is of outmost importance, both to understand biological mechanisms underlying those traits and to facilitate the design of selection programs incorporating that information. For the pig industry, traits related to curing are key for protected designation of origin productions. In particular, appropriate ham weight loss after dry-curing ensures high quality of the final product and avoids economic losses. In this study, we analyzed data (N = 410) of ham weight loss after approximately 20 months of dry-curing. The animals used for ham production were purebred pigs belonging to a commercial line. A genome-wide association study (GWAS) of 29,844 SNP markers revealed the polygenic nature of the trait: 221 loci explaining a small percentage of the variance (0.3-1.65%) were identified on almost all Sus scrofa chromosomes. Post-GWAS analyses revealed 32 windows located within regulatory regions and 94 windows located in intronic regions of specific genes. In total, 30 candidate genes encoding receptors and enzymes associated with ham weight loss (MTHFD1L, DUSP8), proteolysis (SPARCL1, MYH8), drip loss (TNNI2), growth (CDCA3, LSP1, CSMD1, AP2A2, TSPAN4), and fat metabolism (AGPAT4, IGF2R, PTDSS2, HRAS, TALDO1, BRSK2, TNNI2, SYT8, GTF2I, GTF2IRD1, LPCAT3, ATN1, GNB3, CMIP, SORCS2, CCSER1, SPP1) were detected.

Characterization, mRNA expression profile, subcellular distribution and association analysis with piglet diarrhea of porcine matrix metallopeptidase 7 (pMMP7)

Matrix metalloproteinase 7 (MMP7) is involved in the degradation of extracellular matrix in disease processes and therefore plays an important role in host disease resistance/susceptibility. To better understanding the effects of porcine MMP7 (pMMP7) on piglets diarrhea trait, we characterized pMMP7 gene, identified genetic variations in pMMP7 and explored the relationship between pMMP7 polymorphisms and piglets diarrhea in Min pig and Landrace populations. The complete coding sequence of pMMP7 is 804 bp encoding a protein of 267 amino acids. Sequence alignment showed that the identity between pMMP7 and human MMP7 was approximately 80%. The expression of pMMP7 in the gut of healthy piglets were weak and the distribution of the pMMP7-EGFP fusion protein was observed mainly in the cytoplasm. After the identification of 21 genetic variations in 5' flanking region and exons, Hae III and Eco72 Ⅰ PCR-RFLP were established to genotype SNP rs327380117 and rs329429922, respectively. Statistical analysis indicated that Landrace piglets with a TT genotype at rs327380117 had a lower diarrhea score and day-14 wt than TC piglets (p < 0.05); the diarrhea score of AA Landrace animals with rs329429922 was lower than that of GG individuals (p < 0.05). The findings presented here contribute to the understanding of the biological function of pMMP7 and may provide new molecular markers for pig breeding.

Associations between maternal milk protein genotypes with preweaning calf growth traits in beef cattle

The objectives of this study were to investigate milk casein polymorphisms in dams and to determine the impacts of maternal casein genotypes on growth traits of their sucking calves. Milk samples from 433 dams of the breeds German Angus (GA) and German Simmental (GS) were typed at the milk protein loci α s1-casein (αs1-CN), β-casein (β-CN), α s2-casein (αs2-CN), and κ-casein (κ-CN) via isoelectric focusing. Associations between casein genotypes in maternal milk with growth traits of their 1,872 calves were analyzed until the age of weaning using linear mixed models, considering either genotypes of individual casein loci (model 1) or composite α s1-β-α s2-κ-CN genotypes within the casein cluster (model 2). Besides environmental effects such as sex, age of the dam, and calving year-season, genetic effects (breed group and maternal and paternal effects) were considered in statistical models. The composite casein genotype BBǀA2A2ǀAAǀAB (order of genes on bovine chromosome 6: α s1-ǀβ-ǀα s2-ǀκ-CN) was associated with greater average daily weight gains (ADG) and heavier age-adjusted weaning weights (WW) of calves (P < 0.05). The effects of composite genotypes on birth weight of calves were similar (P > 0.05; model 2). With regard to individual casein loci, greater ADG and WW were observed for calves from dams with the genotypes κ-CN BB and α s1-CN BB, respectively (P < 0.05; model 1). Age-adjusted WW was largest for calves from dams carrying the κ-CN genotype BB (215 kg) compared with calves representing the maternal AB and AA genotypes (both 204 kg). Results from the present study suggested selectable casein genotypes due to their nutritional value of milk (value in terms of offspring performances), offering new perspectives for breeding strategies in beef cattle to improve preweaning calf performance.

Transcriptome Characterization of Repressed Embryonic Myogenesis Due to Maternal Calorie Restriction

Fetal malnutrition decreases skeletal myofiber number and muscle mass in neonatal mammals, which increases the risk of developing obesity and diabetes in adult life. However, the associated molecular mechanisms remain unclear. Here, we investigated how the nutrient (calorie) availability affects embryonic myogenesis using a porcine model. Sows were given a normal or calorie restricted diet, following which skeletal muscle was harvested from the fetuses at 35, 55, and 90 days of gestation (dg) and used for histochemical analysis and high-throughput sequencing. We observed abrupt repression of primary myofiber formation following maternal calorie restriction (MCR). Transcriptome profiling of prenatal muscles revealed that critical genes and muscle-specific miRNAs associated with increased proliferation and myoblast differentiation were downregulated during MCR-induced repression of myogenesis. Moreover, we identified several novel miRNA-mRNA interactions through an integrative analysis of their expression profiles, devising a putative molecular network involved in the regulation of myogenesis. Interestingly, NC_010454.3_1179 was identified as a novel myogenic miRNA that can base-pair with sequences in the 3'-UTR of myogenic differentiation protein 1 (MyoD1). And we found that this UTR inhibited the expression of a linked reporter gene encoding a key myogenic regulatory factor, resulting in suppression of myogenesis. Our results greatly increase our understanding of the mechanisms underlying the nutrient-modulated myogenesis, and may also serve as a valuable resource for further investigation of fundamental developmental processes or assist in rational target selection ameliorating repressed myogenesis under fetal malnutrition.

Effect of single-nucleotide polymorphisms on specific reproduction parameters in Hungarian Large White sows

The aim of this study was to reveal the effect of single-nucleotide polymorphisms (SNPs) on the total number of piglets born (TNB), the litter weight born alive (LWA), the number of piglets born dead (NBD), the average litter weight on the 21st day (M21D) and the interval between litters (IBL). Genotypes were determined on a high-density Illumina Porcine SNP 60K BeadChip. Data screening and data identification were performed by a multi-locus mixed-model. Statistical analyses were carried out to find associations between individual genotypes of 290 Hungarian Large White sows and the investigated reproduction parameters. According to the analysis outcome, three SNPs were identified to be associated with TNB. These loci are located on chromosomes 1, 6 and 13 (-log₁₀P = 6.0, 7.86 and 6.22, the frequencies of their minor alleles, MAF, were 0.298, 0.299 and 0.364, respectively). Two loci showed considerable association (-log₁₀P = 10.35 and 10.46) with LWA on chromosomes 5 and X, the MAF were 0.425 and 0.446, respectively. Seven loci were found to be associated with NBD. These loci are located on chromosomes 5, 6, 13, 14, 15, 16 and 18 (-log₁₀P = 10.95, 5.43, 8.29, 6.72, 6.81, 5.90, and 5.15, respectively). One locus showed association (-log₁₀P = 5.62) with M21D on chromosome 1 (the MAF was 0.461). Another locus was found to be associated with IBL on chromosome 8 (-log₁₀P = 7.56; the MAF was 0.438). The above-mentioned loci provide a straightforward possibility to assist selection by molecular tools and, consequently, to improve the competitiveness of the Hungarian Large White (HLW) breed.

A note on cloning and analysis of the osteopontin gene with egg quality traits in two Chinese domestic ducks

1. Osteopontin (OPN) is a highly phosphorylated acidic glycoprotein that plays a crucial role in eggshell formation. In this study, an 893-bp cDNA sequence of the OPN gene, which encodes 180 amino acids, was obtained. 2. Polymorphisms of the OPN gene were analysed with DNA sequencing and polymerase chain reaction-single-strand conformation polymorphism methods in two Chinese domestic laying ducks (Jingding n = 100, Youxian n = 478, respectively). 3. One polymorphism was identified in exon 7 (NM_ 004676534.1:c.267T>C) of the OPN gene, with three genotypes: TT (both T allels weren't mutated (wild type)), TC (one T allel was mutated to C (heterozygote genotype)) and CC. (both T allels were mutated to C (homozygote 20 genotype)) Association analysis with egg quality traits in the two Chinese domestic laying ducks showed that the ducks with the CC genotype had significantly greater eggshell strength and eggshell thickness (p < 0.05). Hence, the exon 7 267T>C polymorphism of the OPN gene is a potentially valuable genetic marker for laying duck breeding programmes.

Linkage disequilibrium and functional analysis of PRE1 insertion together with SNPs in the promoter region of IGFBP7 gene in different pig breeds

Polymorphisms in regions upstream of transcription initiation site may modify the transcriptional activity of target genes by changing promoter activity. This study aims to determine whether or not polymorphisms at porcine IGFBP7 promoter region affect gene expression. In this study, eight SNPs and one PRE1 insertion in this region were first confirmed. The PRE1 insertion was widespread in 20 Chinese indigenous breeds, but was not observed in three commercial breeds. A perfect linkage disequilibrium, consisting of six of those SNPs and a PRE1, was observed with two haplotypes (h1 and h2) in five pig breeds. The h1 haplotype had an overwhelming superiority distribution in Large White, Landrace, and Bama mini-pig; in turn, the h2 only existed in the PRE1 presence breeds. As the haplotypes and PRE1 were located at gene promoter regions, we further investigated the transfection of plasmids with three different fragments of IGFBP-7 promoter region (H1, H2, RF). The CMV promoter of the pEGFP-N1 was substituted by these three different fragments, respectively. Different transcriptional and translational activities of EGFP in PK-15 cells were observed in these three recombinant plasmids by quantitative real-time PCR and flow cytometric analysis. The results indicated that H1 had the higher transcriptional and translational activities of EGFP as compared to the H2 (P < 0.05, P < 0.05). As compared to the RF group, EGFP mRNA expression level was significantly higher in H1 groups (P < 0.05). The IGFBP-7 promoter polymorphisms detected in this study may be important functional variants and potential genetic markers for pig population genetic study.

Transcriptional profiling of subcutaneous adipose tissue in Italian Large White pigs divergent for backfat thickness

Fat deposition is a widely studied trait in pigs because of its implications with animal growth efficiency, technological and nutritional characteristics of meat products, but the global framework of the biological and molecular processes regulating fat deposition in pigs is still incomplete. This study describes the backfat tissue transcription profile in Italian Large White pigs and reports genes differentially expressed between fat and lean animals according to RNA-seq data. The backfat transcription profile was characterised by the expression of 23 483 genes, of which 54.1% were represented by known genes. Of 63 418 expressed transcripts, about 80% were non-previously annotated isoforms. By comparing the expression level of fat vs. lean pigs, we detected 86 robust differentially expressed transcripts, 72 more highly expressed (e.g. ACP5, BCL2A1, CCR1, CD163, CD1A, EGR2, ENPP1, GPNMB, INHBB, LYZ, MSR1, OLR1, PIK3AP1, PLIN2, SPP1, SLC11A1, STC1) and 14 lower expressed (e.g. ADSSL1, CDO1, DNAJB1, HSPA1A, HSPA1B, HSPA2, HSPB8, IGFBP5, OLFML3) in fat pigs. The main functional categories enriched in differentially expressed genes were immune system process, response to stimulus, cell activation and skeletal system development, for the overexpressed genes, and unfolded protein binding and stress response, for the underexpressed genes, which included five heat shock proteins. Adipose tissue alterations and impaired stress response are linked to inflammation and, in turn, to adipose tissue secretory activity, similar to what is observed in human obesity. Our results provide the opportunity to identify biomarkers of carcass fat traits to improve the pig production chain and to identify genetic factors that regulate the observed differential expression.

Analysis of Stage-Specific Gene Expression Profiles in the Uterine Endometrium during Pregnancy in Pigs

The uterine endometrium plays a critical role in regulating the estrous cycle and the establishment and maintenance of pregnancy in mammalian species. Many studies have investigated the expression and function of genes in the uterine endometrium, but the global expression pattern of genes and relationships among genes differentially expressed in the uterine endometrium during gestation in pigs remain unclear. Thus, this study investigated global gene expression profiles using microarray in pigs. Diverse transcriptome analyses including clustering, network, and differentially expressed gene (DEG) analyses were performed to detect endometrial gene expression changes during the different gestation stages. In total, 6,991 genes were found to be differentially expressed by comparing genes expressed on day (D) 12 of pregnancy with those on D15, D30, D60, D90 and D114 of pregnancy, and clustering analysis of detected DEGs distinguished 8 clusters. Furthermore, several pregnancy-related hub genes such as ALPPL2, RANBP17, NF1B, SPP1, and CST6 were discovered through network analysis. Finally, detected hub genes were technically validated by quantitative RT-PCR. These results suggest the complex network characteristics involved in uterine endometrial gene expression during pregnancy and indicate that diverse patterns of stage-specific gene expression and network connections may play a critical role in endometrial remodeling and in placental and fetal development to establish and maintenance of pregnancy in pigs.

Comprehensive analysis of the whole transcriptomes from two different pig breeds using RNA-Seq method

Next-generation sequencing RNA-Seq technology is a powerful tool that creates new possibilities for whole-transcriptome analysis. In our study, the RNA-Seq method was applied to analyze global changes in transcriptome from muscle tissue (m. semimembranosus) in two pig breeds (Pietrain and Polish Landrace, PL). The breeds differ in terms of muscularity, growth rate and reproduction traits. Using three different approaches (deseq, cufflinks and edger) and taking into account the most restrictive criteria, 35 genes differentially expressed between Pietrain and PL pigs were identified. In both breeds, the most abundant were transcripts encoding ribosomal and cytoskeletal proteins (TPM3, TCAP, TMOD4, TPM2, TNNC1) and calcium-binding proteins involved in muscle contraction, calcium-mediated signaling or cation transport (CASQ1, MLC2V, SLC25A4, MYL3). In PL pigs, we identified up-regulation of several genes that play crucial roles in reproduction: female gamete generation (BDP1, PTPN21, USP9X), fertilization (EGFR) and embryonic development (CPEB4). In the Pietrain breed, only seven genes were over-expressed (CISH, SPP1, TUBA8, ATP6V1C2, IGKC, predicted LOC100510960 and LOC100626400), and they play important roles in, for example, negative regulation of apoptosis, immune response, cell-cell signaling, cell growth and migration as well as the metabolic process. The functions of the majority of selected genes were consistent with phenotypic variation in investigated breeds; thus, we proposed a new panel of candidate genes that can be associated with economically important pig traits.

Genome wide screening of candidate genes for improving piglet birth weight using high and low estimated breeding value populations

Birth weight is an economically important trait in pig production because it directly impacts piglet growth and survival rate. In the present study, we performed a genome wide survey of candidate genes and pathways associated with individual birth weight (IBW) using the Illumina PorcineSNP60 BeadChip on 24 high (HEBV) and 24 low estimated breeding value (LEBV) animals. These animals were selected from a reference population of 522 individuals produced by three sires and six dam lines, which were crossbreds with multiple breeds. After quality-control, 43,257 SNPs (single nucleotide polymorphisms), including 42,243 autosomal SNPs and 1,014 SNPs on chromosome X, were used in the data analysis. A total of 27 differentially selected regions (DSRs), including 1 on Sus scrofa chromosome 1 (SSC1), 1 on SSC4, 2 on SSC5, 4 on SSC6, 2 on SSC7, 5 on SSC8, 3 on SSC9, 1 on SSC14, 3 on SSC18, and 5 on SSCX, were identified to show the genome wide separations between the HEBV and LEBV groups for IBW in piglets. A DSR with the most number of significant SNPs (including 7 top 0.1% and 31 top 5% SNPs) was located on SSC6, while another DSR with the largest genetic differences in F_ST was found on SSC18. These regions harbor known functionally important genes involved in growth and development, such as TNFRSF9 (tumor necrosis factor receptor superfamily member 9), CA6 (carbonic anhydrase VI) and MDFIC (MyoD family inhibitor domain containing). A DSR rich in imprinting genes appeared on SSC9, which included PEG10 (paternally expressed 10), SGCE (sarcoglycan, epsilon), PPP1R9A (protein phosphatase 1, regulatory subunit 9A) and ASB4 (ankyrin repeat and SOCS box containing 4). More importantly, our present study provided evidence to support six quantitative trait loci (QTL) regions for pig birth weight, six QTL regions for average birth weight (ABW) and three QTL regions for litter birth weight (LBW) reported previously by other groups. Furthermore, gene ontology analysis with 183 genes harbored in these 27 DSRs suggested that protein, metal, ion and ATP binding, viral process and innate immune response present important pathways for deciphering their roles in fetal growth or development. Overall, our study provides useful information on candidate genes and pathways for regulating birth weight in piglets, thus improving our understanding of the genetic mechanisms involved in porcine embryonic or fetal development.