Genome-wide study on intramuscular fat in Italian Large White pig breed using the PorcineSNP60 BeadChip

Smartphone-based digital phenotyping for genome-wide association study of intramuscular fat traits in longissimus dorsi muscle of pigs

Intramuscular fat (IMF) content and distribution significantly contribute to the eating quality of pork. However, the current methods used for measuring these traits are complex, time-consuming and costly. To simplify the measurement process, this study developed a smartphone application (App) called Pork IMF. This App serves as a rapid and portable phenotyping tool for acquiring pork images and extracting the image-based IMF traits through embedded deep-learning algorithms. Utilizing this App, we collected the IMF traits of the longissimus dorsi muscle in a crossbred population of Large White × Tongcheng pigs. Genome-wide association studies detected 13 and 16 SNPs that were significantly associated with IMF content and distribution, respectively, highlighting NR2F2, MCTP2, MTLN, ST3GAL5, NDUFAB1 and PID1 as candidate genes. Our research introduces a user-friendly digital phenotyping technology for quantifying IMF traits and suggests candidate genes and SNPs for genetic improvement of IMF traits in pigs.

Genome-Wide Association Study and Identification of Candidate Genes for Intramuscular Fat Fatty Acid Composition in Ningxiang Pigs

Ningxiang pigs exhibit a diverse array of fatty acids, making them an intriguing model for exploring the genetic underpinnings of fatty acid metabolism. We conducted a genome-wide association study using a dataset comprising 50,697 single-nucleotide polymorphisms (SNPs) and samples from over 600 Ningxiang pigs. Our investigation yielded novel candidate genes linked to five saturated fatty acids (SFAs), four monounsaturated fatty acids (MUFAs), and five polyunsaturated fatty acids (PUFAs). Significant associations with SFAs, MUFAs, and PUFAs were found for 37, 21, and 16 SNPs, respectively. Notably, some SNPs have significant PVE, such as ALGA0047587, which can explain 89.85% variation in Arachidic acid (C20:0); H3GA0046208 and DRGA0016063 can explain a total of 76.76% variation in Elaidic Acid (C18:1n-9(t)), and the significant SNP ALGA0031262 of Arachidonic acid (C20:4n-6) can explain 31.76% of the variation. Several significant SNPs were positioned proximally to previously reported genes. In total, we identified 11 candidate genes (hnRNPU, CEPT1, ATP1B1, DPT, DKK1, PRKG1, EXT2, MEF2C, IL17RA, ITGA1 and ALOX5), six candidate genes (ALOX5AP, MEDAG, ISL1, RXRB, CRY1, and CDKAL1), and five candidate genes (NDUFA4L2, SLC16A7, OTUB1, EIF4E and ROBO2) associated with SFAs, MUFAs, and PUFAs, respectively. These findings hold great promise for advancing breeding strategies aimed at optimizing meat quality and enhancing lipid metabolism within the intramuscular fat (IMF) of Ningxiang pigs.

Genome-wide mapping of signatures of selection using a high-density array identified candidate genes for growth traits and local adaptation in chickens

BACKGROUND

Availability of single nucleotide polymorphism (SNP) genotyping arrays and progress in statistical analyses have allowed the identification of genomic regions and genes under selection in chicken. In this study, SNP data from the 600 K Affymetrix chicken array were used to detect signatures of selection in 23 local Italian chicken populations. The populations were categorized into four groups for comparative analysis based on live weight (heavy vs light) and geographical area (Northern vs Southern Italy). Putative signatures of selection were investigated by combining three extended haplotype homozygosity (EHH) statistical approaches to quantify excess of haplotype homozygosity within (iHS) and between (Rsb and XP-EHH) groups. Presence of runs of homozygosity (ROH) islands was also analysed for each group.

RESULTS

After editing, 541 animals and 313,508 SNPs were available for statistical analyses. In total, 15 candidate genomic regions that are potentially under selection were detected among the four groups: eight within a group by iHS and seven by combining the results of Rsb and XP-EHH, which revealed divergent selection between the groups. The largest overlap between genomic regions identified to be under selection by the three approaches was on chicken chromosome 8. Twenty-one genomic regions were identified with the ROH approach but none of these overlapped with regions identified with the three EHH-derived statistics. Some of the identified regions under selection contained candidate genes with biological functions related to environmental stress, immune responses, and disease resistance, which indicate local adaptation of these chicken populations.

CONCLUSIONS

Compared to commercial lines, local populations are predominantly reared as backyard chickens, and thus, may have developed stronger resistance to environmental challenges. Our results indicate that selection can play an important role in shaping signatures of selection in local chicken populations and can be a starting point to identify gene mutations that could have a useful role with respect to climate change.

RNA sequencing analysis of the longissimus dorsi to identify candidate genes underlying the intramuscular fat content in Anqing Six-end-white pigs

Intramuscular fat (IMF) is a significant marker for pork quality. The Anqing Six-end-white pig has the characteristics of high meat quality and IMF content. Owing to the influence of European commercial pigs and a late start in resource conservation, the IMF content within local populations varies between individuals. This study analyzed the longissimus dorsi transcriptome of purebred Anqing Six-end-white pigs with varying IMF content to recognize differentially expressed genes. We identified 1528 differentially expressed genes between the pigs with high (H) and low (L) IMF content. Based on these data, 1775 Gene Ontology terms were significantly enriched, including lipid metabolism, modification and storage, and regulation of lipid biosynthesis. Pathway analysis revealed 79 significantly enriched pathways, including the Peroxisome proliferator-activated receptor and mitogen-activated protein kinase signaling pathways. Moreover, gene set enrichment analysis indicated that the L group had increased the expression of genes related to ribosome function. Additionally, the protein-protein interaction network analyses revealed that VEGFA, KDR, LEP, IRS1, IGF1R, FLT1 and FLT4 were promising candidate genes associated with the IMF content. Our study identified the candidate genes and pathways involved in IMF deposition and lipid metabolism and provides data for developing local pig germplasm resources.

Genome-Wide Association Studies for Flesh Color and Intramuscular Fat in (Duroc × Landrace × Large White) Crossbred Commercial Pigs

The intuitive impression of pork is extremely important in terms of whether consumers are enthusiastic about purchasing it. Flesh color and intramuscular fat (IMF) are indispensable indicators in meat quality assessment. In this study, we determined the flesh color and intramuscular fat at 45 min and 12 h after slaughter (45 mFC, 45 mIMF, 12 hFC, and 12 hIMF) of 1518 commercial Duroc × Landrace × Large White (DLY) pigs. We performed a single nucleotide polymorphism (SNP) genome-wide association study (GWAS) analysis with 28,066 SNPs. This experiment found that the correlation between 45 mFC and 12 hFC was 0.343. The correlation between 45 mIMF and 12 hIMF was 0.238. The heritability of the traits 45 mFC, 12 hFC, 45 mIMF, and 12 hIMF was 0.112, 0.217, 0.139, and 0.178, respectively, and we identified seven SNPs for flesh color and three SNPs for IMF. Finally, several candidate genes regulating these four traits were identified. Three candidate genes related to flesh color were provided: SNCAIP and PRR16 on SSC2, ST3GAL4 on SSC5, and GALR1 on SSC1. A total of three candidate genes related to intramuscular fat were found, including ABLIM3 on SSC2, DPH5 on SSC4, and DOCK10 on SSC15. Furthermore, GO and KEGG analysis revealed that these genes are involved in the regulation of apoptosis and are implicated in functions such as pigmentation and skeletal muscle metabolism. This study applied GWAS to analyze the scoring results of flesh color and IMF in different time periods, and it further revealed the genetic structure of flesh color and IMF traits, which may provide important genetic loci for the subsequent improvement of pig meat quality traits.

A composite strategy of genome-wide association study and copy number variation analysis for carcass traits in a Duroc pig population

BACKGROUND

Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs.

RESULTS

Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function.

CONCLUSIONS

We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations.

Genome-wide identification of functional enhancers and their potential roles in pig breeding

Background

The pig is an economically important livestock species and is a widely applied large animal model in medical research. Enhancers are critical regulatory elements that have fundamental functions in evolution, development and disease. Genome-wide quantification of functional enhancers in the pig is needed.

Results

We performed self-transcribing active regulatory region sequencing (STARR-seq) in the porcine kidney epithelial PK15 and testicular ST cell lines, and reliably identified 2576 functional enhancers. Most of these enhancers were located in repetitive sequences and were enriched within silent and lowly expressed genes. Enhancers poorly overlapped with chromatin accessibility regions and were highly enriched in chromatin with the repressive histone modification H3K9me3, which is different from predicted pig enhancers detected using ChIP-seq for H3K27ac or/and H3K4me1 modified histones. This suggests that most pig enhancers identified with STARR-seq are endogenously repressed at the chromatin level and may function during cell type-specific development or at specific developmental stages. Additionally, the PPP3CA gene is associated with the loin muscle area trait and the QKI gene is associated with alkaline phosphatase activity that may be regulated by distal functional enhancers.

Conclusions

In summary, we generated the first functional enhancer map in PK15 and ST cells for the pig genome and highlight its potential roles in pig breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00726-y.

Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds

Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.

Genome-wide associations for immune traits in two maternal pig lines

Background

In recent years, animal welfare and health has become more and more important in pig breeding. So far, numerous parameters have been considered as important biomarkers, especially in the immune reaction and inflammation. Previous studies have shown moderate to high heritabilities in most of these traits. However, the genetic background of health and robustness of pigs needs to be extensively clarified. The objective of this study was to identify genomic regions with a biological relevance for the immunocompetence of piglets. Genome-wide Association Studies (GWAS) in 535 Landrace (LR) and 461 Large White (LW) piglets were performed, investigating 20 immune relevant traits. Besides the health indicators of the complete and differential blood count, eight different cytokines and haptoglobin were recorded in all piglets and their biological dams to capture mediating processes and acute phase reactions. Additionally, all animals were genotyped using the Illumina PorcineSNP60v2 BeadChip.

Results

In summary, GWAS detected 25 genome-wide and 452 chromosome-wide significant SNPs associated with 17 immune relevant traits in the two maternal pig lines LR and LW. Only small differences were observed considering the maternal immune records as covariate within the statistical model. Furthermore, the study identified across- and within-breed differences as well as relevant candidate genes. In LR more significant associations and related candidate genes were detected, compared with LW. The results detected in LR and LW are partly in accordance with previously identified quantitative trait loci (QTL) regions. In addition, promising novel genomic regions were identified which might be of interest for further detailed analysis. Especially putative pleiotropic regions on SSC5, SSC12, SSC15, SSC16 and SSC17 are of major interest with regard to the interacting structure of the immune system. The comparison with already identified QTL gives indications on interactions with traits affecting piglet survival and also production traits.

Conclusion

In conclusion, results suggest a polygenic and breed-specific background of immune relevant traits. The current study provides knowledge about regions with biological relevance for health and immune traits. Identified markers and putative pleiotropic regions provide first indications in the context of balancing a breeding-based modification of the porcine immune system.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07997-1).

Genomic Analyses Revealed the Genetic Difference and Potential Selection Genes of Growth Traits in Two Duroc Lines

Duroc pigs are famous for their high growth rate, feed conversion efficiency, and lean meat percentage. Given that they have been subjected to artificial selection and breeding in multiple countries, various lines with obvious differences in production performance have formed. In this study, we genotyped 3,770 American Duroc (AD) pigs and 2,098 Canadian Duroc (CD) pigs using the GeneSeek Porcine SNP50 Beadchip to dissect the genetic differences and potential selection genes of growth traits in these two Duroc pig lines. Population structure detection showed that there were significant genetic differences between the two Duroc pig lines. Hence, we performed F _ST and cross-population extended haplotype homozygosity (XP-EHH) analyses between the two lines. As a result, we identified 38 annotated genes that were significantly enriched in the gland development pathway in the AD line, and 61 annotated genes that were significantly enriched in the immune-related pathway in the CD line. For three growth traits including backfat thickness (BFT), loin muscle depth (LMD), and loin muscle area (LMA), we then performed selection signature detection at 5 and 10% levels within the line and identified different selected regions and a series of candidate genes that are involved in lipid metabolism and skeletal muscle development or repair, such as IRX3, EBF2, WNT10B, TLR2, PITX3, and SGCD. The differences in selected regions and genes between the two lines may be the cause of the differences in growth traits. Our study suggests significant genetic differences between the AD and CD lines, which provide a theoretical basis for selecting different Duroc lines as sires for different needs.

Genome-wide association studies of preweaning growth and in vivo carcass composition traits in Esme sheep

Sheep are considered as a major contributor of global food security. Moreover, sheep preweaning growth traits as well as in vivo carcass composition traits such as ultrasonic measurements of Longissimus dorsi muscle depth (UMD) and back-fat thickness (UFD) are crucially important indicators of meat yield and hot carcass composition. Despite their relative importance for productivity and profitability of a sheep production system, detected QTL for these traits are quite scarce. Therefore, we implemented GWAS for these traits using animal mixed model-based association approach provided by GenABEL in Esme sheep. Three genome-wide and 14 individual chromosome-wide associated SNPs were discovered. As a result, ESRP1, LOC105613082, ZNF641, DUSP5, TEAD1, SMOX, PTPRT, RALYL, POM121C, PHIP, LOC101106051, ZIM3, PEG3, TRPC7, FBXL4, LOC105610397, LOC105616489 and DNAAF2 were suggested as candidates. Some of the discovered genes and involved pathways were already annotated to contribute growth and development in various species including human, mice and cattle. All in all, the results of this study are expected to strongly contribute to shed a light on the underlying molecular mechanisms behind growth and carcass composition traits, with potential implications on studies aiming faster genetic improvement, targeted low-resolution SNP panel designs and genome-editing studies.

Evidence Against the Causal Relationship Between a Putative Cis-Regulatory Variant of MYH3 and Intramuscular Fat Content in Pigs

Improving meat quality has become the main goal of modern pig breeding. Intramuscular fat content (IMF) is an important trait influencing meat quality of livestock, but the molecular mechanism behind this trait is still unclear. Recently, Cho et al. reported the discovery of the first causal mutation affecting IMF and red flesh color (a*) in pigs, namely XM_013981330.2:g.−1805_−1810del, a 6-bp deletion variant in the porcine MYH3 promoter region. The objective of this study was to reassess the causality of this mutation for its potential commercial application. By Sanger sequencing, we firstly identified several new variants (including a 4-bp deletion) at or near the 6-bp deletion site, which formed four haplotypes in multiple breeds. Unexpectedly, the 6-bp deletion allele, previously determined as the MYH3 Q allele because of its significantly positive effect on IMF and a*, was found not only in Chinese indigenous breeds, but also in four western commercial breeds with relatively lower IMF levels, including Duroc, Large White, Landrace and Pietrain. More surprisingly, we found that the MYH3 Q allele and the haplotypes harboring it had no significant effects on IMF, marbling and color score in three large-scale divergent pig populations: the heterogeneous F6 and F7 pigs and commercial crossbred Duroc × (Landrace × Yorkshire) pigs. Transient transfection analysis in porcine satellite cells showed that the 6-bp deletion variants had a negligible effect on transcription of reporter gene, but could attenuate the MRF (myogenesis regulatory factors)-induced increase in luciferase activity of the MYH3 promoter vector. The MYH3 protein level in muscle did not differ significantly among the haplotype groups. Therefore, our results cannot support the causal relationship between the 6-bp deletion in MYH3 and IMF trait, suggesting that the causal mutation for the IMF QTL on SSC12 needs to be further identified.

Gut Microbiota Influence Lipid Metabolism of Skeletal Muscle in Pigs

Gut microbiota is recognized as a strong determinant of host physiology including fat metabolism and can transfer obesity-associated phenotypes from donors to recipients. However, the relationship between gut microbiota and intramuscular fat (IMF) is still largely unknown. Obese Jinhua pigs (JP) have better meat quality that is associated with higher IMF content than lean Landrace pigs (LP). The present study was conducted to test the contribution of gut microbiota to IMF properties by transplanting fecal microbiota of adult JP and LP to antibiotics-treated mice. Similar to JP donors, the mice receiving JP's microbiota (JM) had elevated lipid and triglyceride levels and the lipoprotein lipase activity, as well as reduced mRNA level of angiopoietin-like 4 (ANGPTL4) in the gastrocnemius muscles, compared to those in mice receiving LP's microbiota (LM). High-throughput 16S rRNA sequencing confirmed that transplantation of JP and LP feces differently reconstructed the gut microbiota in both jejunum and colon of mouse recipients. In colonic samples, we observed an elevated ratio of Firmicutes to Bacteroidetes and increased abundance of genus Romboutsia in JM, which were positively correlated with obesity. Furthermore, the abundance of Akkermansia decreased in JM, which is positively correlated with lean. Colonic concentrations of acetate (P = 0.047) and butyrate (P = 0.014) were significantly lower in JM than in LM, and consistently, the terminal genes for butyrate synthesis, butyryl CoA: acetate CoA transferase were less abundant in colonic microbiota of JM. Taken together, these gut microbiota of obese JP intrinsically promotes IMF accumulation and can transfer the properties to mouse recipients. Manipulation of intestinal microbiota will, therefore, have the potential to improve the meat quality and flavor of pigs and even to ameliorate the metabolic syndrome in human.

Dissecting the Gene Expression Networks Associated with Variations in the Major Components of the Fatty Acid Semimembranosus Muscle Profile in Large White Heavy Pigs

Simple Summary

The amount and fatty acid composition of intramuscular fat are important features for the qualitative characteristics of processed and fresh meat products, but the knowledge of the key molecular drivers controlling these traits is still scant. To this aim, the present study investigated the co-expression networks of genes related to variations in the major fatty acids deposited in pig Semimembranosus muscle. Palmitic and palmitoleic acid contents were associated with a downregulation of genes involved in autophagy, mitochondrial fusion, and mitochondrial activity, suggesting that the deposition of these fatty acids may be enhanced in muscles with a reduced mitochondrial function. A higher proportion of oleic acid and a reduction in the percentages of n-6 and n-3 polyunsaturated fatty acids were related to changes in the mRNA levels of genes involved in Mitogen-Activated Protein Kinase (MAPK) signaling. The obtained results indicated gene expression networks and new candidate genes associated with the studied traits. Further studies are needed to confirm our results and identify in the discussed genes molecular markers for future selection schemes aimed at improving pork nutritional and technological quality. Furthermore, as pigs are considered reliable animal models for several human conditions, the obtained results may also be of interest for improving the knowledge of the molecular pathways associated with obesity and diabetes.

Abstract

To date, high-throughput technology such as RNA-sequencing has been successfully applied in livestock sciences to investigate molecular networks involved in complex traits, such as meat quality. Pork quality depends on several organoleptic, technological, and nutritional characteristics, and it is also influenced by the fatty acid (FA) composition of intramuscular fat (IMF). To explore the molecular networks associated with different IMF FA compositions, the Semimembranosus muscle (SM) from two groups of Italian Large White (ILW) heavy pigs divergent for SM IMF content was investigated using transcriptome analysis. After alignment and normalization, the obtained gene counts were used to perform the Weighted Correlation Network Analysis (WGCNA package in R environment). Palmitic and palmitoleic contents showed association with the same gene modules, comprising genes significantly enriched in autophagy, mitochondrial fusion, and mitochondrial activity. Among the key genes related to these FAs, we found TEAD4, a gene regulating mitochondrial activity that seems to be a promising candidate for further studies. On the other hand, the genes comprised in the modules associated with the IMF contents of oleic, n-6, and n-3 polyunsaturated FAs (PUFAs) were significantly enriched in Mitogen-Activated Protein Kinase (MAPK) signaling, in agreement with previous studies suggesting that several MAPK players may have a primary role in regulating lipid deposition. These results give an insight into the molecular cascade associated with different IMF FA composition in ILW heavy pigs. Further studies are needed to validate the results and confirm whether some of the identified key genes may be effective candidates for pork quality.

Investigating the Features of PDO Green Hams during Salting: Insights for New Markers and Genomic Regions in Commercial Hybrid Pigs

Protected Designation of Origin (PDO) dry-cured hams production is greatly dependent on raw meat quality. This study was performed to identify genetic markers associated with the quality of dry-cured ham. Carcass traits of 229 heavy pigs belonging to three commercial genetic lines were registered (weight, EUROP classification). Phenotypic traits (Semimembranosus muscle ultimate pH, ham weight and lean meat content, adsorbed salt) of the corresponding thighs, undergone PDO ham process in three different plants, were measured, using a fast and non-invasive technology. Green ham weight and lean meat percentage influenced the estimated salt content and the weight loss during salting, even if the processing plant greatly affected the variability of the measured ham traits. The genomic data were obtained with the GeneSeek Genomic Profiler (GGP) 70k HD Porcine Array, using the slaughter day and the sex of the animals in the statistical analyses. The phenotypic traits were associated with the genotypes through GenAbel software. The results showed that 18 SNPs located on nine porcine chromosomes were found to be associated with nine phenotypic traits, mainly related to ham weight loss during salting. New associations were found between markers in the genes Neural Precursor Cell Expressed Developmentally Down-Regulated 9 (NEDD9, SSC7), T-Cell Lymphoma Invasion and Metastasis 2 (TIAM2, SSC1), and the ham quality traits. After validation, these SNPs may be useful to improve the quality of thighs for the production of PDO dry-cured hams.

Genomic Tools for the Conservation and Genetic Improvement of a Highly Fragmented Breed-The Ramo Grande Cattle from the Azores

Simple Summary

Inbreeding control is a key concern in managing local endangered breeds, which often have developed unique adaptation features. Ramo Grande is a local cattle breed raised in the Azores archipelago under very harsh conditions, with a census of about 1300 cows dispersed by various islands. This fragmentation is a challenge when the goal is to keep inbreeding under control. Currently, panels of genetic markers are available which enable the assessment of inbreeding and the occurrence of previous bottlenecks in a population. These panels also allow the identification of genes associated with specific production traits, if reliable phenotypic information is available. We used a panel of genetic markers and estimated that the degree of inbreeding was approaching a level of concern, while some exotic gene inflow may have occurred in the past. We were able to identify genetic markers significantly associated with longevity, which reflects the ability of these cattle to remain productive under severe environmental conditions. Genetic markers were also identified as significantly associated with age at first calving and calf growth rate. The results indicate that genomic information can be used to control inbreeding and to implement genomic selection in Ramo Grande cattle to enhance adaptation and production traits.

Abstract

Ramo Grande is a local cattle breed raised in the archipelago of Azores, with a small and dispersed census, where inbreeding control is of utmost importance. A single nucleotide polymorphism (SNP) Beadchip array was used to assess inbreeding, by analysis of genomic regions harboring contiguous homozygous genotypes named runs of homozygosity (ROH), and to estimate past effective population size by analysis of linkage disequilibrium (LD). Genetic markers associated with production traits were also investigated, exploiting the unique genetic and adaptation features of this breed. A total of 639 ROH with length >4 Mb were identified, with mean length of 14.96 Mb. The mean genomic inbreeding was 0.09, and long segments of ROH were common, indicating recent inbred matings. The LD pattern indicates a large effective population size, suggesting the inflow of exotic germplasm in the past. The genome-wide association study identified novel markers significantly affecting longevity, age at first calving and direct genetic effects on calf weight. These results provide the first evidence of the association of longevity with genes related with DNA recognition and repair, and the association of age at first calving with aquaporin proteins, which are known to have a crucial role in reproduction.

miR-34a regulates adipogenesis in porcine intramuscular adipocytes by targeting ACSL4

BACKGROUND

Intramuscular fat (IMF) content is an important factor in porcine meat quality. Previously, we showed that miR-34a was less abundant in liver tissue from pigs with higher backfat thickness, compared to pigs with lower backfat thickness. The purpose of this present study was to explore the role of miR-34a in adipogenesis.

RESULT

Bioinformatics analysis identified Acyl-CoA synthetase long chain family member 4 (ACSL4) as a putative target of miR-34a. Using a luciferase reporter assay, we verified that miR-34a binds the ACSL4 mRNA at the 3'UTR. To examine the role of the miR-34a-ACSL4 interaction in IMF deposition in the pig, mRNA and protein expression of the ACSL4 gene was measured in primary intramuscular preadipocytes transfected with miR-34a mimic and inhibitor. Our results showed that ACSL4 is expressed throughout the entire differentiation process in pig preadipocytes, similar to the lipogenesis-associated genes PPARγ and aP2. Transfection with miR-34a mimic reduced lipid droplet formation during adipogenesis, while miR-34a inhibitor increased lipid droplet accumulation. Transfection with miR-34a mimic also reduced the mRNA and protein expression of ACSL4 and lipogenesis genes, including PPARγ, aP2, and SREBP-1C, but increased the expression of steatolysis genes such as ATGL and Sirt1. In contrast, the miR-34a inhibitor had the opposite effect on gene expression. Further, knockdown of ACSL4 decreased lipid droplet accumulation.

CONCLUSIONS

Our results support the hypothesis that miR-34a regulates intramuscular fat deposition in porcine adipocytes by targeting ACSL4.

Integrating genome-wide co-association and gene expression to identify putative regulators and predictors of feed efficiency in pigs

BACKGROUND

Feed efficiency (FE) has a major impact on the economic sustainability of pig production. We used a systems-based approach that integrates single nucleotide polymorphism (SNP) co-association and gene-expression data to identify candidate genes, biological pathways, and potential predictors of FE in a Duroc pig population.

RESULTS

We applied an association weight matrix (AWM) approach to analyse the results from genome-wide association studies (GWAS) for nine FE associated and production traits using 31K SNPs by defining residual feed intake (RFI) as the target phenotype. The resulting co-association network was formed by 829 SNPs. Additive effects of this SNP panel explained 61% of the phenotypic variance of RFI, and the resulting phenotype prediction accuracy estimated by cross-validation was 0.65 (vs. 0.20 using pedigree-based best linear unbiased prediction and 0.12 using the 31K SNPs). Sixty-eight transcription factor (TF) genes were identified in the co-association network; based on the lossless approach, the putative main regulators were COPS5, GTF2H5, RUNX1, HDAC4, ESR1, USP16, SMARCA2 and GTF2F2. Furthermore, gene expression data of the gluteus medius muscle was explored through differential expression and multivariate analyses. A list of candidate genes showing functional and/or structural associations with FE was elaborated based on results from both AWM and gene expression analyses, and included the aforementioned TF genes and other ones that have key roles in metabolism, e.g. ESRRG, RXRG, PPARGC1A, TCF7L2, LHX4, MAML2, NFATC3, NFKBIZ, TCEA1, CDCA7L, LZTFL1 or CBFB. The most enriched biological pathways in this list were associated with behaviour, immunity, nervous system, and neurotransmitters, including melatonin, glutamate receptor, and gustation pathways. Finally, an expression GWAS allowed identifying 269 SNPs associated with the candidate genes' expression (eSNPs). Addition of these eSNPs to the AWM panel of 829 SNPs did not improve the accuracy of genomic predictions.

CONCLUSIONS

Candidate genes that have a direct or indirect effect on FE-related traits belong to various biological processes that are mainly related to immunity, behaviour, energy metabolism, and the nervous system. The pituitary gland, hypothalamus and thyroid axis, and estrogen signalling play fundamental roles in the regulation of FE in pigs. The 829 selected SNPs explained 61% of the phenotypic variance of RFI, which constitutes a promising perspective for applying genetic selection on FE relying on molecular-based prediction.

Single-Locus and Multi-Locus Genome-Wide Association Studies for Intramuscular Fat in Duroc Pigs

Intramuscular fat (IMF) is an important quantitative trait of meat, which affects the associated sensory properties and nutritional value of pork. To gain a better understanding of the genetic determinants of IMF, we used a composite strategy, including single-locus and multi-locus association analyses to perform genome-wide association studies (GWAS) for IMF in 1,490 Duroc boars. We estimated the genomic heritability of IMF to be 0.23 ± 0.04. A total of 30 single nucleotide polymorphisms (SNPs) were found to be significantly associated with IMF. The single-locus mixed linear model (MLM) and multiple-locus methods multi-locus random-SNP-effect mixed linear model (mrMLM), fast multi-locus random-SNP-effect efficient mixed model association (FASTmrEMMA), and integrative sure independence screening expectation maximization Bayesian least absolute shrinkage and selection operator model (ISIS EM-BLASSO) analyses identified 5, 9, 8, and 21 significant SNPs, respectively. Interestingly, a novel quantitative trait locus (QTL) on SSC 7 was found to affect IMF. In addition, 10 candidate genes (BDKRB2, GTF2IRD1, UTRN, TMEM138, DPYD, CASQ2, ZNF518B, S1PR1, GPC6, and GLI1) were found to be associated with IMF based on their potential functional roles in IMF. GO analysis showed that most of the genes were involved in muscle and organ development. A significantly enriched KEGG pathway, the sphingolipid signaling pathway, was reported to be associated with fat deposition and obesity. Identification of novel variants and functional genes will advance our understanding of the genetic mechanisms of IMF and provide specific opportunities for marker-assisted or genomic selection in pigs. In general, such a composite single-locus and multi-locus strategy for GWAS may be useful for understanding the genetic architecture of economic traits in livestock.

Transcriptome Analysis Reveals the Effect of Long Intergenic Noncoding RNAs on Pig Muscle Growth and Fat Deposition

Muscle growth and fat deposition are the two important biological processes in the development of pigs which are closely related to the pig production performance. Long intergenic noncoding RNAs (lincRNAs), with lack of coding potential and the length of at least 200nt, have been extensively studied to play important roles in many biological processes. However, the importance and molecular regulation mechanism of lincRNAs in the process of muscle growth and fat deposition in pigs are still to be further studied comprehensively. In our study, we used the data, including liver, abdominal fat, and longissimus dorsi muscle of 240 days' age of two F2 full-sib female individuals from the white Duroc and Erhualian crossbreed, to identify 581 putative lincRNAs associated with pig muscle growth and fat deposition. The 581 putative lincRNAs shared many common features with other mammalian lincRNAs, such as fewer exons, lower expression levels, and shorter transcript lengths. Cross-tissue comparisons showed that many transcripts were tissue-specific and were involved in the important biological processes in their corresponding tissues. Gene ontology and pathway analysis revealed that many potential target genes (PTGs) of putative lincRNAs were involved in pig muscle growth and fat deposition-related processes, including muscle cell proliferation, lipid metabolism, and fatty acid degradation. In Quantitative Trait Locus (QTLs) analysis, some PTGs were screened from putative lincRNAs, MRPL12 is associated with muscle growth, GCGR and SLC25A10 were associated with fat deposition, and PPP3CA, DPYD, and FGGY were related not only to muscle growth but also to fat deposition. Therefore, it implied that these lincRNAs might participate in the biological processes related to muscle growth or fat deposition through homeostatic regulation of PTGs, but the detailed molecular regulatory mechanisms still needed to be further explored. This study lays the molecular foundation for the in-depth study of the role of lincRNAs in the pig muscle growth and fat deposition and further provides the new molecular markers for understanding the complex biological mechanisms of pig muscle growth and fat deposition.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

BACKGROUND

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

RESULTS

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

CONCLUSIONS

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

A whole-genome sequence based association study on pork eating quality traits and cooking loss in a specially designed heterogeneous F6 pig population

To determine the genetic basis of pork eating quality traits and cooking loss, we herein performed a genome-wide association study (GWAS) for tenderness, juiciness, oiliness, umami, overall liking and cooking loss by using whole genome sequences of heterogeneous stock F6 pigs which were generated by crossing 4 typical western pig breeds (Duroc, Landrace, Large White and Pietrain) and 4 typical Asian pig breeds (Erhualian, Laiwu, Bamaxiang and Tibetan). We identified 50 associated loci (QTLs) and most of them are novel. Seven loci also showed pleiotropic associations with different traits. In addition, we identified multiple promising candidate genes for these traits, including PAK1 and AQP11 for cooking loss, EP300 for tenderness, SDK1 for juiciness, FITM2 and 5-linked MYH genes for oiliness, and TNNI2 and TNNT3 for overall liking. Our results provide not only a better understanding of the genetic basis for meat quality, but also a potential application in future breeding for these complex traits.

Extent of linkage disequilibrium and effective population size of the Landrace population in Korea

Objective

The genetic diversity of the Landrace population, a representative maternal pig breed in Korea, is important for genetic improvement. Previously, the effective population size (Ne) has been used to infer the genetic diversity of a population of interest. In this study, we aimed to use single nucleotide polymorphism (SNP) data to characterize linkage disequilibrium (LD) and the Ne of the Korean Landrace population.

Methods

We genotyped 1,128 Landrace individuals from three representative Korean major grand-grand-parent (GGP) farms using the Illumina PorcineSNP60 version2 BeadChip, which covers >61,565 SNPs located across all autosomes and mitochondrial and sex chromosomes. We estimated the expected LD and current Ne, as well as ancestral Ne.

Results

In the Korean Landrace population, the mean LD (r²) of 3.698 million SNP pairs was 0.135±0.204. The mean r² decreased slowly with as the distance between SNPs increased, and remained constant beyond 3 Mb. According to the r² calculations, 8,085 of 3.698 million SNP pairs were in complete LD. The current Ne (±standard deviation) of the Korean Landrace population is approximately 92.27 [79.46; 105.07] individuals. The ancestral Ne exhibited a slow and steady decline from 186.61 to 92.27 over the past 100 generations. Additionally, we observed more a rapid Ne decrease from the past 20 to 10 generations ago, compared with other intervals.

Conclusion

We have presented an overview of LD and the current and ancestral Ne values in the Korean Landrace population. The mean LD and current Ne for the Korean Landrace population confirm the genetic diversity and reflect the history of this pig population in Korea.

RNA-seq based detection of differentially expressed genes in the skeletal muscle of Duroc pigs with distinct lipid profiles

We have used a RNA-seq approach to investigate differential expression in the skeletal muscle of swine (N = 52) with divergent lipid profiles i.e. HIGH (increased intramuscular fat and muscle saturated and monounsaturated fatty acid contents, higher serum lipid concentrations and fatness) and LOW pigs (leaner and with an increased muscle polyunsaturated fatty acid content). The number of mRNAs and non-coding RNAs (ncRNAs) expressed in the porcine gluteus medius muscle were 18,104 and 1,558, respectively. At the nominal level of significance (P-value ≤ 0.05), we detected 1,430 mRNA and 12 non-coding RNA (ncRNA) transcripts as differentially expressed (DE) in the gluteus medius muscle of HIGH vs LOW pigs. This smaller contribution of ncRNAs to differential expression may have biological and technical reasons. We performed a second analysis, that was more stringent (P-value ≤ 0.01 and fold-change ≥ 1.5), and only 96 and 0 mRNA-and ncRNA-encoding genes happened to be DE, respectively. The subset of DE mRNA genes was enriched in pathways related with lipid (lipogenesis and triacylglycerol degradation) and glucose metabolism. Moreover, HIGH pigs showed a more lipogenic profile than their LOW counterparts.

Study on quantitative expression of PPARγ and ADRP in muscle and its association with intramuscular fat deposition of pig

Background

Intramuscular fat (intramuscular fat, IMF) is one of the important traits of pork quality. How to reasonably improve the intramuscular fat content is the most focus researchers. Some possible regulation of intramuscular fat deposition of candidate genes to cause the attention of people. The objective of this study was to elucidate the relationship between peroxisome proliferator-activated receptor γ (PPARγ) and adipose differentiation-related protein (ADRP) mRNA expression and intramuscular fat (IMF) deposition in the muscle tissue of three breeds of pig: Laiwu (LW), Lulai Black (LL), and Large White (LY).

Results

qPCR analysis of the PPARγ and ADRP genes in the three breeds of pig revealed PPARγ and ADRP mRNA expression profiles of LW > LL > LY and LL > LW > LY, respectively. PPARγ mRNA expression was significantly and positively correlated with IMF deposition (p < 0.05). There were significant correlations between PPARγ and ADRP mRNA expression levels (p < 0.01).

Conclusions

These results suggest correlations between PPARγ and ADRP in fat deposition and regulation in pigs, PPARγ gene may be a main effector of IMF content and play an important role during adipocyte differentiation in pigs, thereby providing new information to further elucidate molecular mechanisms associated with intramuscular fat deposition in Laiwu pigs and provides new data for further molecular studies of mechanisms underlying intramuscular fat deposition in human obesity. The continued elucidation of specific genetic mechanisms between PPARγ and ADRP warrants further studies.

A gene and protein expression study on four porcine genes related to intramuscular fat deposition

Intramuscular fat (IMF) content has a prominent role in meat quality, affecting sensory attributes such as flavour and texture. In the present research, we studied in samples of porcine Semimembranosus muscle four genes related to lipid metabolism and whose gene expressions have been associated to IMF deposition: FASN, SCD, LIPE and LPL. We analysed both mRNA and protein expressions in two groups of Italian Large White pigs divergent for Semimembranosus IMF deposition, with the aim of comparing the levels of four genes and enzymes between the two groups and identifying possible coexpression links. The obtained results suggest a prominent role of LIPE enzyme in IMF hydrolysis, as the samples with low IMF deposition show a significantly higher amount of this lipase. Finally, a poorly known correlation was found between LIPE and FASN enzymes only in female individuals. These results provide new information for the understanding of IMF deposition.