Haplotypes of the porcine peroxisome proliferator-activated receptor delta gene are associated with backfat thickness

Meta-analysis of genome-wide association studies uncovers shared candidate genes across breeds for pig fatness trait

Background

Average backfat thickness (BFT) is a critical complex trait in pig and an important indicator for fat deposition and lean rate. Usually, genome-wide association study (GWAS) was used to discover quantitative trait loci (QTLs) of BFT in a single population. However, the power of GWAS is limited by sample size in a single population. Alternatively, meta-analysis of GWAS (metaGWAS) is an attractive method to increase the statistical power by integrating data from multiple breeds and populations. The aim of this study is to identify shared genetic characterization of BFT across breeds in pigs via metaGWAS. 

Results

In this study, we performed metaGWAS on BFT using 15,353 pigs (5,143 Duroc, 7,275 Yorkshire, and 2,935 Landrace) from 19 populations. We detected 40 genome-wide significant SNPs (Bonferroni corrected P < 0.05) and defined five breed-shared QTLs in across-breed metaGWAS. Markers within the five QTL regions explained 7 ~ 9% additive genetic variance and showed strong heritability enrichment. Furthermore, by integrating information from multiple bioinformatics databases, we annotated 46 candidate genes located in the five QTLs. Among them, three important (MC4R, PPARD, and SLC27A1) and seven suggestive candidate genes (PHLPP1, NUDT3, ILRUN, RELCH, KCNQ5, ITPR3, and U3) were identified.

Conclusion

QTLs and candidate genes underlying BFT across breeds were identified via metaGWAS from multiple populations. Our findings contribute to the understanding of the genetic architecture of BFT and the regulating mechanism underlying fat deposition in pigs.

Expression and Genetic Effects of GLI Pathogenesis-Related 1 Gene on Backfat Thickness in Pigs

Backfat thickness (BFT) is an important carcass composition trait and regarded as a breeding focus. Our initial transcriptome analysis of pig BFT identified GLI pathogenesis-related 1 (GLIPR1) as one of the promising candidate genes. This study was conducted to identify the expression profiles, polymorphisms, and genetic effects of the GLIPR1 gene on BFT in pigs. The expression of the GLIPR1 gene existed in every detected tissue, and there was a significantly higher expression in spleen and adipose tissue than others (p < 0.05). At the different ages of pig, the expression of the GLIPR1 gene was low at an early age, increased with growth, and reached the highest level at 180 days. Genetic polymorphism analysis was detected in 553 individuals of the Large White × Minzhu F2 population. Four SNPs in the promoter significantly associated with 6−7 rib BFT (p < 0.05) were predicted to alter the transcription factor binding sites (TFBS), and the mutations of g.38758089 T>G and g.38758114 G>C were predicted to change the TFs associated with the regulation of adipogenesis. Haplotypes were formed by the detected SNPs, and one block showed a strong association with BFT (p < 0.05). In summary, our results indicate that the expression profiles and genetic variants of GLIPR1 affected the BFT of pigs. To our knowledge, this study is the first to demonstrate the biological function and genetic effects of the GLIPR1 gene on the BFT of pig and provide genetic markers to optimize breeding for BFT in pigs.

Screening of Genes Related to Growth, Development and Meat Quality of Sahan Crossbred F1 Sheep Based on RNA-Seq Technology

This study aimed to identify genes related to sheep growth, development and meat quality. Small-tailed Han sheep (STH), and small-tailed Han sheep and Suffolk crossbred F1 (STH×SFK), were selected to determine the growth performance, slaughter performance, and meat quality. The longissimus dorsi muscle was selected for transcriptome sequencing, and the target gene was screened based on bioinformatics analysis; real-time fluorescent quantitative PCR (RT-PCR) and western blotting (WB) were conducted to verify the target gene. Locations of genes in tissues were confirmed via immunofluorescence. The results showed that the pre-slaughter live weight, bust circumference, slaughter performance, and marbling score of the STH×SFK population were significantly higher than those of the STH population (P < 0.01). Sequencing results showed that 560 differentially expressed genes (DEGs) were identified in the STH×SFK population, of which 377 exhibited up-regulated and 183 exhibited down-regulated expression levels. GO annotation revealed that DEGs could be classified into 13 cell components, 10 molecular functions, and 22 biological processes. The KEGG enrichment analysis showed that DEGs were mainly enriched in the Rap1 signaling pathway, Ras signaling pathway, and other pathways related to growth and meat quality. Based on the GO and KEGG analyses, four candidate genes related to sheep growth and meat quality, namely myostain (MSTN), interferon-related developmental regulator 1 (IFRD1), peroxisome proliferator activator receptor delta (PPARD), and myosin light chain 2 (MLC2 or MYL2), were screened. The expression levels of genes and proteins were verified via RT-PCR and WB, and the results were consistent with the trend of transcriptome sequencing. Immunofluorescence results showed that IFRD1 was expressed in the cytoplasm and nucleus, and MYL2 was expressed in the cytoplasm. This study revealed the mechanism of gene regulation of sheep growth and development at the molecular level and provided a theoretical basis for studying sheep genetics and breeding.

Identification of candidate genes that specifically regulate subcutaneous and intramuscular fat deposition using transcriptomic and proteomic profiles in Dingyuan pigs

Subcutaneous fat and intramuscular fat (IMF) deposition are closely related to meat production and pork quality. Dingyuan pig, as a native pig breed in China, low selection leads to obvious genetic and phenotypic differences in the population. Individuals with extreme fat content in the population are ideal models for studying the mechanism of fat deposition. In this study, we used RNA-Seq and tandem mass tags-based (TMT) proteomics to analyze the key pathways and genes that specifically regulate subcutaneous fat and IMF deposition in Dingyuan pigs. We identified 191 differentially expressed genes (DEGs) and 61 differentially abundant proteins (DAPs) in the high and low back fat thickness (HBF, LBF) groups, 85 DEGs and 12 DAPs were obtained in the high and low intramuscular fat (HIMF, LIMF) groups. The functional analysis showed that the DEGs and DAPs in the backfat groups were mainly involved in carbohydrates, amino acids, and fatty acids metabolism, whereas the IMF groups were involved in the insulin pathway, longevity, and some disease-related pathways. We found 40 candidate genes that might tissue-specifically lipids deposition for subcutaneous and intramuscular fat. Our research provides theoretical reference materials for the improvement of fat deposition traits of local pig breeds in my country.

Genome-Wide Assessment of Runs of Homozygosity and Estimates of Genomic Inbreeding in a Chinese Composite Pig Breed

The primary purpose of the current study was to assess the genetic diversity, runs of homozygosity (ROH) and ROH islands in a Chinese composite pig and explore hotspot regions for traces of selection. First, we estimated the length, number, and frequency of ROH in 262 Xidu black pigs using the Porcine SNP50 BeadChip and compared the estimates of inbreeding coefficients, which were calculated based on ROHs (F_ROH) and homozygosity (F_HOM). Our result shows that a total of 7,248 ROH exceeding 1Mb were detected in 262 pigs. In addition, Sus scrofa chromosome (SSC) 8 and SSC10, respectively, has the highest and lowest chromosome coverage by ROH. These results suggest that inbreeding estimation based on total ROH may be a useful method, especially for crossbreed or composite populations. We also calculated an inbreeding coefficient of 0.077 from the total ROH. Eight ROH islands were found in this study. These ROH islands harbored genes associated with fat deposition, muscular development, reproduction, ear shape, and adaptation, such as TRAF7, IGFBP7, XPO1, SLC26A8, PPARD, and OR1F1. These findings may help to understand the effects of environmental and artificial selection on the genome structure of composite pigs. Our results provide a basis for subsequent genomic selection (GS), and provides a reference for the hybrid utilization of other pig breeds.

Discovery of selection-driven genetic differences of Duroc, Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.

Genes involved in angiogenesis and circulatory system development are differentially expressed in porcine epithelial oviductal cells during long-term primary in vitro culture – a transcriptomic study

An oviduct is an essential organ for gamete transport, oocyte maturation, fertilization, spermatozoon capacitation and early embryo development. The epithelium plays an important role in oviduct functioning. The products of secretory cells provide an optimal environment and influence gamete activities and embryonic development. The oviduct physiology changes during the female cycle, thus, the ratio of the secreted molecules in the oviduct fluid differs between phases. In this study, a differential gene expression in porcine oviduct epithelial cells was examined during the long-term primary in vitro culture. The microarray expression analysis revealed 2552 genes, 1537 of which were upregulated and 995 were downregulated after 7 days of culture, with subsequent changes in expression during 30 day-long culture. The obtained genes were classified into 8 GO BP terms, connected with angiogenesis and circulatory system development, extracted by DAVID software. Among all genes, 10 most up-regulated and 10 most down-regulated genes were selected for further investigation. Interactions between genes were indicated by STRING software and REACTOME FIViz application to the Cytoscape 3.6.0 software. Most of the genes belonged to more than one ontology group. Although studied genes are mostly responsible for angiogenesis and circulatory system development, they can also be found to be expressed in processes connected with fertilization and early embryo development. The latter function is focused on more, considering the fact that these genes were expressed in epithelial cells of the fallopian tube which is largely responsible for reproductive processes.

A genome-wide detection of selection signatures in conserved and commercial pig breeds maintained in Poland

BACKGROUND

Identification of selection signatures can provide a direct insight into the mechanism of artificial selection and allow further disclosure of the candidate genes related to the animals' phenotypic variation. Domestication and subsequent long-time selection have resulted in extensive phenotypic changes in domestic pigs, involving a number of traits, like behavior, body composition, disease resistance, reproduction and coat color. In this study, based on genotypes obtained from PorcineSNP60 Illumina assay we attempt to detect both diversifying and within-breed selection signatures in 530 pigs belonging to four breeds: Polish Landrace, Puławska, Złotnicka White and Złotnicka Spotted, of which the last three are a subject of conservative breeding and substantially represent the native populations.

RESULTS

A two largely complementary statistical methods were used for signatures detection, including: pairwise F_ST and relative extended haplotype homozygosity (REHH) test. Breed-specific diversifying selection signals included several genes involved in processes connected with fertility, growth and metabolism which are potentially responsible for different phenotypes of the studied breeds. The diversifying selection signals also comprised PPARD gene that was previously found to have a large effect on the shape of the external ear in pigs or two genes encoding neuropeptide Y receptors (Y2 and Y5) involved in fat deposition and stress response which are important features differentiating the studied breeds. REHH statistics allowed detecting several within-breed selection signatures overlapping with genes connected with a range of functions including, among others: metabolic pathways, immune system response or implantation and development of the embryo.

CONCLUSIONS

The study provides many potential candidate genes with implication for traits selected in the individual breeds and gives strong basis for further studies aiming at identification of sources of variation among the studied pig breeds.

The Functional SNPs in the 5' Regulatory Region of the Porcine PPARD Gene Have Significant Association with Fat Deposition Traits

Peroxisome proliferator-activated receptor delta (PPARD) is a key regulator of lipid metabolism, insulin sensitivity, cell proliferation and differentiation. In this study, we identified two Single Nucleotide Polymorphisms (SNPs, g.1015 A>G and g.1018 T>C) constituting four haplotypes (GT, GC, AC and AT) in the 5’ regulatory region of porcine PPARD gene. Functional analysis of the four haplotypes showed that the transcriptional activity of the PPARD promoter fragment carrying haplotype AC was significantly lower than that of the other haplotypes in 3T3-L1, C2C12 and PK-15 cells, and haplotype AC had the lowest binding capacities to the nuclear extracts. Transcription factor 7-like 2 (TCF7L2) enhanced the transcription activities of promoter fragments of PPARD gene carrying haplotypes GT, GC and AT in C2C12 and 3T3-L1 cells, and increased the protein expression of PPARD gene in C2C12 myoblasts. TCF7L2 differentially bound to the four haplotypes, and the binding capacity of TCF7L2 to haplotype AC was the lowest. There were significant associations between -655A/G and fat deposition traits in three pig populations including the Large White × Meishan F₂ pigs, France and American Large White pigs. Pigs with genotype GG had significantly higher expression of PPARD at both mRNA and protein level than those with genotype AG. These results strongly suggested that the SNPs in 5’ regulatory region of PPARD genes had significant impact on pig fat deposition traits.

A genome scan for selection signatures in pigs

Identifying signatures of selection can provide a straightforward insight into the mechanism of artificial selection and further uncover the causal genes related to the phenotypic variation. Based on Illumina Porcine60KSNP chip data, four complementary methods, Long-Range Haplotype (LRH), Tajima’s D, Cross Population Extend Haplotype Homozygosity Test (XPEHH) and F_ST, were implemented in this study to detect the selection signatures in the whole genome of one typical Chinese indigenous breed, Rongchang, one Chinese cultivated breed, Songliao, and two western breeds, Landrace and Yorkshire. False Discovery Rate (FDR) was implemented to control the false positive rates. In our study, a total of 159, 127, 179 and 159 candidate selection regions with average length of 0.80 Mb, 0.73 Mb, 0.78 Mb and 0.73 Mb were identified in Landrace, Rongchang, Songliao and Yorkshire, respectively, that span approximately 128.00 Mb, 92.38 Mb, 130.30 Mb and 115.40 Mb and account for approximately 3.74–5.33% of genome across all autosomes. The selection regions of 11.52 Mb shared by Landrace and Yorkshire were the longest when chosen pairs from the pool of the four breeds were examined. The overlaps between Yorkshire and Songliao, approximately 9.20 Mb, were greater than those of Yorkshire and Rongchang. Meanwhile, the overlaps between Landrace and Songliao were greater than those of Landrace and Rongchang but less than those of Songliao and Ronchang. Bioinformatics analysis showed that the genes/QTLs relevant to fertility, coat color, and ear morphology were found in candidate selection regions. Some genes, such as LEMD3, MC1R, KIT, TRHR etc. that were reported under selection, were confirmed in our study, and this analysis also demonstrated the diversity of breeds.

Polymorphism in 3' untranslated region of the pig PPARA gene influences its transcript level and is associated with adipose tissue accumulation

The PPARA (peroxisome proliferator-activated receptor-α) gene encodes a nuclear receptor that plays an important role in fatty acid catabolism by transcriptional regulation of genes involved in fatty acid oxidation and can be considered as a candidate gene for fatness traits in the pig. The aim of the study was to search for a functional polymorphism in 3' untranslated region (UTR), their association with production traits, and postnatal PPARA transcript level in 2 skeletal muscles (longissimus and semimembranosus) of 5 commercial pig breeds (Polish Landrace [PL], Polish Large White [PLW], Duroc, Pietrain, and Pulawska). Altogether, 9 novel polymorphisms (8 SNP and 1 indel) were found in the 3' UTR. The in silico analysis revealed 6 putative microRNA target sequences in the analyzed region. The c.*636A>G substitution was widely distributed across breeds and located near the putative target sequence for miR-224. The relative PPARA transcript level was higher (P < 0.05) in LM of AA than in those of GG homozygous animals for SNP c.*636A>G. The luciferase assay revealed that miR-224 probably acts as a negative regulator of the PPARA expression in pig adipocytes (P = 2.9 × 10(-7)), but we did not observe the effect of the A or G alleles on the interaction between miR-224 and its putative target sequence. We hypothesize that the 2 predominant haplotypes, differing at 4 sites (including c.*636A>G), present different architecture of its 3' UTR and it could affect the level of the transcript. The c.*636A>G SNP, analyzed in PL and PLW, was significantly associated with backfat thickness at 3 points (P < 0.05) and intramuscular fat content (P < 0.01) in PL. Suggestive associations were found between 4 SNP (c.*321A>C, c.*324G>C, c.*626T>C, and c.*636A>G) and fatty acid contents in LM and subcutaneous and visceral fat tissue of PL, PLW, Duroc and Pietrain pigs. The PPARA mRNA level was higher in semimembranosus muscle than in LM (P = 8.38 × 10(-12)) in a general comparison and the same trend was found in most breeds (except for PL) and at all tested days of age (60, 90, 120, 150, 180, and 210 d). The effect of breed was highly significant in a general comparison (P = 0.48 × 10(-8)), but there was no common expression pattern in both muscles among different age groups. We conclude that the c.*636A>G SNP in the PPARA gene can be considered in PL breed as a useful genetic marker for adipose tissue accumulation.

From QTL to QTN: candidate gene set approach and a case study in porcine IGF1-FoxO pathway

Unraveling the genetic background of economic traits is a major goal in modern animal genetics and breeding. Both candidate gene analysis and QTL mapping have previously been used for identifying genes and chromosome regions related to studied traits. However, most of these studies may be limited in their ability to fully consider how multiple genetic factors may influence a particular phenotype of interest. If possible, taking advantage of the combined effect of multiple genetic factors is expected to be more powerful than analyzing single sites, as the joint action of multiple loci within a gene or across multiple genes acting in the same gene set will likely have a greater influence on phenotypic variation. Thus, we proposed a pipeline of gene set analysis that utilized information from multiple loci to improve statistical power. We assessed the performance of this approach by both simulated and a real IGF1-FoxO pathway data set. The results showed that our new method can identify the association between genetic variation and phenotypic variation with higher statistical power and unravel the mechanisms of complex traits in a point of gene set. Additionally, the proposed pipeline is flexible to be extended to model complex genetic structures that include the interactions between different gene sets and between gene sets and environments.

Towards candidate genes affecting body fatness at the SSC7 QTL by expression analyses

A quantitative trait locus (QTL) affecting fatness in a way opposite to expectations based on breed means was mapped to swine chromosome 7 (SSC7) using crosses between Large White (LW) and Meishan (MS) founders. Defining the molecular fatness trait more explicitly would allow deducing positional candidate genes, for which expression differences must be analysed in experimental populations. First, mRNA levels of genes representing sequential steps in adipogenesis or involved in lipid metabolism were studied in backfat of pigs having homozygous LW(QTL7)/LW(QTL7) or heterozygous LW(QTL7)/MS(QTL7) alleles and considered at two ages. mRNA level of DLK1 expressed in preadipocytes was greater in MS(QTL7)/LW(QTL7) pigs than in homozygous pigs at 28 days. Transcript abundances of CEBPA involved in differentiation, the prolipogenic FASN gene and the adipocyte-specific marker FABP4 were lower in MS(QTL7)/LW(QTL7) pigs compared with LW(QTL7)/LW(QTL7) pigs at 150 days. Because these results suggest a lag time in terminal differentiation associated with the MS allele, seven genes in the QTL interval were deduced as promising candidates for the QTL effect by bioinformatics analysis. Among them, PPARD and CDKN1A had lower expression levels in MS(QTL7)/LW(QTL7) pigs at both ages. Genotype-related differences were observed in mRNA levels of PPARD target genes involved in cell differentiation (FZD7) or fatty acid oxidation (ACADL and ACOX1) at 150 days. These results re-evaluate the potential of PPARD to explain part of variation in pig adiposity.

Mapping quantitative trait loci for metabolic and cytological fatness traits of connected F2 crosses in pigs

In the present study 3 connected F(2) crosses were used to map QTL for classical fat traits as well as fat-related metabolic and cytological traits in pigs. The founder breeds were Chinese Meishan, European Wild Boar, and Pietrain with to some extent the same founder animals in the different crosses. The different selection history of the breeds for fatness traits as well as the connectedness of the crosses led to a high statistical power. The total number of F(2) animals varied between 694 and 966, depending on the trait. The animals were genotyped for around 250 genetic markers, mostly microsatellites. The statistical model was a multi-allele, multi-QTL model that accounted for imprinting. The model was previously introduced from plant breeding experiments. The traits investigated were backfat depth and fat area as well as relative number of fat cells with different sizes and 2 metabolic traits (i.e., soluble protein content as an indicator for the level of metabolic turnover and NADP-malate dehydrogenase as an indicator for enzyme activity). The results revealed in total 37 significant QTL on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 14, 17, and 18, with often an overlap of confidence intervals of several traits. These confidence intervals were in some cases remarkably small, which is due to the high statistical power of the design. In total, 18 QTL showed significant imprinting effects. The small and overlapping confidence intervals for the classical fatness traits as well as for the cytological and metabolic traits enabled positional and functional candidate gene identification for several mapped QTL.