A genome-wide association analysis for carcass traits in a commercial Duroc pig population

Gestational Caloric Restriction Alters Adipose Tissue Methylome and Offspring's Metabolic Profile in a Swine Model

Limited nutrient supply to the fetus results in physiologic and metabolic adaptations that have unfavorable consequences in the offspring. In a swine animal model, we aimed to study the effects of gestational caloric restriction and early postnatal metformin administration on offspring's adipose tissue epigenetics and their association with morphometric and metabolic variables. Sows were either underfed (30% restriction of total food) or kept under standard diet during gestation, and piglets were randomly assigned at birth to receive metformin (n = 16 per group) or vehicle treatment (n = 16 per group) throughout lactation. DNA methylation and gene expression were assessed in the retroperitoneal adipose tissue of piglets at weaning. Results showed that gestational caloric restriction had a negative effect on the metabolic profile of the piglets, increased the expression of inflammatory markers in the adipose tissue, and changed the methylation of several genes related to metabolism. Metformin treatment resulted in positive changes in the adipocyte morphology and regulated the methylation of several genes related to atherosclerosis, insulin, and fatty acids signaling pathways. The methylation and gene expression of the differentially methylated FASN, SLC5A10, COL5A1, and PRKCZ genes in adipose tissue associated with the metabolic profile in the piglets born to underfed sows. In conclusion, our swine model showed that caloric restriction during pregnancy was associated with impaired inflammatory and DNA methylation markers in the offspring's adipose tissue that could predispose the offspring to later metabolic abnormalities. Early metformin administration could modulate the size of adipocytes and the DNA methylation changes.

Genomic architecture of carcass and pork traits and their association with immune capacity

Carcass and pork traits have traditionally been considered of prime importance in pig breeding programmes. However, the changing conditions in modern farming, coupled with antimicrobial resistance issues, are raising the importance of health and robustness-related traits. Here, we explore the genetic architecture of carcass and pork traits and their relationship with immunity phenotypes in a commercial Duroc pig population. A total of nine traits related to fatness, lean content and meat pH were measured at slaughter (∼190 d of age) in 378 pigs previously phenotyped (∼70 d of age) for 36 immunity-related traits, including plasma concentrations of immunoglobulins, acute-phase proteins, leukocytes subpopulations and phagocytosis. Our study showed medium to high heritabilities and strong genetic correlations between fatness, lean content and meat pH at 24 h postmortem. Genetic correlations were found between carcass and pork traits and white blood cells. pH showed strong positive genetic correlations with leukocytes and eosinophils, and strong negative genetic correlations with haemoglobin, haematocrit and cytotoxic T cell proportion. In addition, genome-wide association studies (GWASs) pointed out four significantly associated genomic regions for lean meat percentages in different muscles, ham fat, backfat thickness, and semimembranosus pH at 24 h. The functional annotation of genes located in these regions reported a total of 14 candidate genes, with BGN, DPP10, LEPR, LEPROT, PDE4B and SLC6A8 being the strongest candidates. After performing an expression GWAS for the expression of these genes in muscle, two signals were detected in cis for the BGN and SLC6A8 genes. Our results indicate a genetic relationship between carcass fatness, lean content and meat pH with a variety of immunity-related traits that should be considered to improve immunocompetence without impairing production traits.

Assessing the relationship between the in silico predicted consequences of 97 missense mutations mapping to 68 genes related to lipid metabolism and their association with porcine fatness traits

In general, the relationship between the predicted functional consequences of missense mutations mapping to genes known to be involved in human diseases and the severity of disease manifestations is weak. In this study, we tested in pigs whether missense single nucleotide polymorphisms (SNPs), predicted to have consequences on the function of genes related to lipid metabolism are associated with lipid phenotypes. Association analysis demonstrated that nine out of 72 nominally associated SNPs were classified as "highly" or "very highly consistent" in silico-predicted functional mutations and did not show association with lipid traits expected to be affected by inactivation of the corresponding gene. Although the lack of endophenotypes and the limited sample size of certain genotypic classes might have limited to some extent the reach of the current study, our data indicate that present-day bioinformatic tools have a modest ability to predict the impact of missense mutations on complex phenotypes.

Integrated Omics analysis of pig muscle metabolism under the effects of dietary Chlorella vulgaris and exogenous enzymes

Monogastric feeding is dependent on costly conventional feedstuffs. Microalgae such as Chlorella vulgaris are a sustainable alternative; however, its recalcitrant cell wall hinders monogastric digestion. Carbohydrate Active Enzyme (CAZyme) supplementation is a possible solution. The objective of this work was to evaluate the effect of 5% dietary C. vulgaris (CV) and enzymatic supplementation (CV + R-Rovabio® Excel AP; CV + M-four CAZyme mix) on muscle transcriptome and proteome of finishing pigs, in an integrated approach. Control pigs increased the abundance of contractile apparatus (MYH1, MYH2, MYH4) and energy metabolism (CKMT1, NDUFS3) proteins, demonstrating increased nutrient availability. They had increased expression of SCD, characteristic of increased glucose availability, via the activation of SREBP-1c and ChREBP. CV and CV + R pigs upregulated proteolytic and apoptotic genes (BAX, DDA1), whilst increasing the abundance of glucose (UQCRFS1) and fatty acid catabolism (ACADS) proteins. CV + R pigs upregulated ACOT8 and SIRT3 genes as a response to reduced nutrient availability, maintaining energy homeostasis. The cell wall specific CAZyme mix, CV + M, was able to comparatively reduce Omics alterations in the muscle, thereby reducing endogenous nutrient catabolism compared to the CV + R and CV.

Genome-wide association study identifying genetic variants associated with carcass backfat thickness, lean percentage and fat percentage in a four-way crossbred pig population using SLAF-seq technology

Background

Carcass backfat thickness (BFT), carcass lean percentage (CLP) and carcass fat percentage (CFP) are important to the commercial pig industry. Nevertheless, the genetic architecture of BFT, CLP and CFP is still elusive. Here, we performed a genome-wide association study (GWAS) based on specific-locus amplified fragment sequencing (SLAF-seq) to analyze seven fatness-related traits, including five BFTs, CLP, and CFP on 223 four-way crossbred pigs.

Results

A total of 227, 921 highly consistent single nucleotide polymorphisms (SNPs) evenly distributed throughout the genome were used to perform GWAS. Using the mixed linear model (MLM), a total of 20 SNP loci significantly related to these traits were identified on ten Sus scrofa chromosomes (SSC), of which 10 SNPs were located in previously reported quantitative trait loci (QTL) regions. On SSC7, two SNPs (SSC7:29,503,670 and rs1112937671) for average backfat thickness (ABFT) exceeded 1% and 10% Bonferroni genome-wide significance levels, respectively. These two SNP loci were located within an intron region of the COL21A1 gene, which was a protein-coding gene that played an important role in the porcine backfat deposition by affecting extracellular matrix (ECM) remodeling. In addition, based on the other three significant SNPs on SSC7, five candidate genes, ZNF184, ZNF391, HMGA1, GRM4 and NUDT3 were proposed to influence BFT. On SSC9, two SNPs for backfat thickness at 6–7 ribs (67RBFT) and one SNP for CLP were in the same locus region (19 kb interval). These three SNPs were located in the PGM2L1 gene, which encoded a protein that played an indispensable role in glycogen metabolism, glycolysis and gluconeogenesis as a key enzyme. Finally, one significant SNP on SSC14 for CLP was located within the PLBD2 gene, which participated in the lipid catabolic process.

Conclusions

A total of two regions on SSC7 and SSC9 and eight potential candidate genes were found for fatness-related traits in pigs. The results of this GWAS based on SLAF-seq will greatly advance our understanding of the genetic architecture of BFT, CLP, and CFP traits. These identified SNP loci and candidate genes might serve as a biological basis for improving the important fatness-related traits of pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08827-8.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in Muscle Development of Wu'an Goat Longissimus Dorsi Tissues

The growth and development of skeletal muscle is a physiological process regulated by a variety of genes and signaling pathways. As a posttranscriptional regulatory factor, circRNA plays a certain regulatory role in the development of animal skeletal muscle in the form of a miRNA sponge. However, the role of circRNAs in muscle development and growth in goats is still unclear. In our study, apparent differences in muscle fibers in Wu'an goats of different ages was firstly detected by hematoxylin-eosin (HE) staining, the circRNA expression profiles of longissimus dorsi muscles from 1-month-old (mon1) and 9-month-old (mon9) goats were screened by RNA-seq and verified by RT-qPCR. The host genes of differentially expressed (DE) circRNAs were predicted, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses (KEGG) of host genes with DE circRNAs were performed to explore the functions of circRNAs. The circRNA-miRNA-mRNA networks were then constructed using Cytoscape software. Ten significantly differentially expressed circRNAs were also verified in the mon1 and mon9 groups by RT-qPCR. Luciferase Reporter Assay was used to verify the binding site between circRNA and its targeted miRNA. The results showed that a total of 686 DE circRNAs were identified between the mon9 and mon1 groups, of which 357 were upregulated and 329 were downregulated. Subsequently, the 467 host genes of DE circRNAs were predicted using Find_circ and CIRI software. The circRNA-miRNA-mRNA network contained 201 circRNAs, 85 miRNAs, and 581 mRNAs; the host mRNAs were associated with "muscle fiber development" and "AMPK signaling pathway" and were enriched in the FoxO signaling pathway. Competing endogenous RNA (ceRNA) network analysis showed that novel_circ_0005314, novel_circ_0005319, novel_circ_0009256, novel_circ_0009845, novel_circ_0005934 and novel_circ_0000134 may play important roles in skeletal muscle growth and development between the mon9 and mon1 groups. Luciferase Reporter Assay confirmed the combination between novel_circ_0005319 and chi-miR-199a-5p, novel_circ_0005934 and chi-miR-450-3p and novel_circ_0000134 and chi-miR-655. Our results provide specific information related to goat muscle development and a reference for the goat circRNA profile.

Potential functional variants in AHR signaling pathways are associated with age at puberty in swine

Puberty in female pigs is defined as age at first estrus and gilts that have an earlier age at puberty are more likely to have greater lifetime productivity. Because age at puberty is predictive for sow longevity and lifetime productivity, but not routinely measured in commercial herds, it would be beneficial to use genomic or marker-assisted selection to improve these traits. A GWAS at the US Meat Animal Research Center (USMARC) identified several loci associated with age at puberty in pigs. Candidate genes in these regions were scanned for potential functional variants using sequence information from the USMARC swine population founder animals and public databases. In total, 135 variants (SNP and insertion/deletions) in 39 genes were genotyped in 1284 phenotyped animals from a validation population sired by Landrace and Yorkshire industry semen using the Agena MassArray system. Twelve variants in eight genes were associated with age at puberty (P < 0.005) with estimated additive SNP effects ranging from 1.6 to 5.3 days. Nine of these variants were non-synonymous coding changes in AHR, CYP1A2, OR2M4, SDCCAG8, TBC1D1 and ZNF608, two variants were deletions of one and four codons in aryl hydrocarbon receptor, AHR, and the most significant SNP was near an acceptor splice site in the acetyl-CoA carboxylase alpha, ACACA. Several of the loci identified have a physiological and a genetic role in sexual maturation in humans and other animals and are involved in AHR-mediated pathways. Further functional validation of these variants could identify causative mutations that influence age at puberty in gilts and possibly sow lifetime productivity.

Genome-Wide Association Study and Pathway Analysis for Heterophil/Lymphocyte (H/L) Ratio in Chicken

Disease control and prevention have been critical factors in the dramatic growth of the poultry industry. Disease resistance in chickens can be improved through genetic selection for immunocompetence. The heterophil/lymphocyte ratio (H/L) in the blood reflects the immune system status of chickens. Our objective was to conduct a genome-wide association study (GWAS) and pathway analysis to identify possible biological mechanisms involved in H/L traits. In this study, GWAS for H/L was performed in 1317 Cobb broilers to identify significant single-nucleotide polymorphisms (SNPs) associated with H/L. Eight SNPs (p < 1/8068) reached a significant level of association. The significant SNP on GGA 19 (chicken chromosome 19) was in the gene for complement C1q binding protein (C1QBP). The wild-type and mutant individuals showed significant differences in H/L at five identified SNPs (p < 0.05). According to the results of pathway analysis, nine associated pathways (p < 0.05) were identified. By combining GWAS with pathway analysis, we found that all SNPs after QC explained 12.4% of the phenotypic variation in H/L, and 52 SNPs associated with H/L explained as much as 9.7% of the phenotypic variation in H/L. Our findings contribute to understanding of the genetic regulation of H/L and provide theoretical support.

Co-expression network analysis predicts a key role of microRNAs in the adaptation of the porcine skeletal muscle to nutrient supply

BACKGROUND

The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake.

RESULTS

We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation.

CONCLUSIONS

As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

Identifying Selection Signatures for Backfat Thickness in Yorkshire Pigs Highlights New Regions Affecting Fat Metabolism

Identifying the genetic basis of improvement in pigs contributes to our understanding of the role of artificial selection in shaping the genome. Here we employed the Cross Population Extended Haplotype Homozogysity (XPEHH) and the Wright's fixation index (F_ST) methods to detect trait-specific selection signatures by making phenotypic gradient differential population pairs, and then attempted to map functional genes of six backfat thickness traits in Yorkshire pigs. The results indicate that a total of 283 and 466 single nucleotide polymorphisms (SNPs) were identified as trait-specific selection signatures using F_ST and XPEHH, respectively. Functional annotation suggested that the genes overlapping with the trait-specific selection signatures such as OSBPL8, ASAH2, SMCO2, GBE1, and ABL1 are responsible for the phenotypes including fat metabolism, lean body mass and fat deposition, and transport in mouse. Overall, the study developed the methods of gene mapping on the basis of identification of selection signatures. The candidate genes putatively associated with backfat thickness traits can provide important references and fundamental information for future pig-breeding programs.

Genomic Analysis To Identify Signatures of Artificial Selection and Loci Associated with Important Economic Traits in Duroc Pigs

Identifying genetic basis of domestication and improvement in livestock contributes to our understanding of the role of artificial selection in shaping the genome. Here we used whole-genome sequencing and the genotyping by sequencing approach to detect artificial selection signatures and identify the associated SNPs of two economic traits in Duroc pigs. A total of 38 candidate selection regions were detected by combining the fixation index and the Composite Likelihood Ratio methods. Further genome-wide association study revealed seven associated SNPs that were related with intramuscular fat content and feed conversion ratio traits, respectively. Enrichment analysis suggested that the artificial selection regions harbored genes, such as MSTN, SOD2, MC5R and CD83, which are responsible for economic traits including lean muscle mass, fertility and immunization. Overall, this study found a series of candidate genes putatively associated with the breeding improvement of Duroc pigs and the polygenic basis of adaptive evolution, which can provide important references and fundamental information for future breeding programs.