Muscle glycogen level and occurrence of acid meat in commercial hybrid pigs are regulated by two low-frequency causal variants with large effects and multiple common variants with small effects

Profiling of m6A methylation in porcine intramuscular adipocytes and unravelling PHKG1 represses porcine intramuscular lipid deposition in an m6A-dependent manner

Intramuscular fat (IMF) content is mainly determined by intramuscular preadipocyte adipogenesis. Epigenetic modifications are known to have a regulatory effect on IMF. As N6-methyladenosine (m6A) is the most abundant epigenetic modification in eukaryotic RNAs. In the present study, we used m6A methylation and RNA sequencing (seq) to identify the m6A-modified RNAs associated with the adipogenic differentiation of intramuscular preadipocytes. Among them, the expression and m6A level of phosphorylase kinase subunit G1 (PHKG1) were found to be significantly changed during adipogenesis. Further studies revealed that knockdown of the methylase METTL3 decreased the m6A methylation of PHKG1 and led to a reduction in PHKG1. Moreover, knockdown of PHKG1 promoted adipogenic differentiation by upregulating the expression of adipogenic genes. In addition, we found that the IMF content in the longissimus thoracis (LT) of Bamei (BM) pigs was greater than that in Large White (LW) pigs, whereas the m6A and PHKG1 expression levels were lower in BM pigs. These findings indicate that the m6A level and expression of PHKG1 were significantly correlated with IMF content and meat quality. In conclusion, this study sheds light on the mechanism by which m6A modification regulates IMF deposition.

Comparison of Growth Performance and Plasma Metabolomics between Two Sire-Breeds of Pigs in China

The Yorkshire pigs, renowned for their remarkable growth rate, low feed conversion ratio (FCR), and high meat production, emerge as a novel preference for paternal breeding. In this study, we found that purebred paternal Yorkshire pigs (PY) surpass the purebred Duroc breed in terms of growth rate. Specifically, purebred PY attain a weight of 100 kg at an earlier age compared to purebred Duroc (Male, 145.07 vs. 162.91; Female, 145.91 vs. 167.57; p-value < 0.01). Furthermore, different hybrid combinations suggest that offspring involving purebred PY exhibit superior growth performance. Compared with purebred Duroc, the offspring of purebred PY have an earlier age in days (173.23 vs. 183.54; p-value < 0.05) at the same slaughter weight. The changes of plasma metabolites of 60-day-old purebred boars in the two sire-breeds showed that 1335 metabolites in plasma were detected. Compared with Duroc, 28 metabolites were down-regulated and 49 metabolites were up-regulated in PY. Principal component analysis (PCA) discerned notable dissimilarities in plasma metabolites between the two sire-breeds of pigs. The levels of glycerol 3-phosphate choline, cytidine, guanine, and arachidonic acid increased significantly (p-value < 0.05), exerting an impact on their growth and development. According to our results, PY could be a new paternal option as a terminal sire in three-way cross system.

Carcass and meat quality traits and their relationships in Duroc × Landrace × Yorkshire barrows slaughtered at various seasons

To understand characteristics of carcass traits and meat quality in pig population, 22 indicators of carcass characteristics and meat quality traits were measured on 278 Duroc × Landrace × Yorkshire barrows that were slaughtered in different seasons (spring, summer, autumn and winter). The effects of body weight and season on carcass characteristics and meat quality were analyzed by GLM procedure, followed the Bonferroni multiple test. The phenotypic correlations among those traits were calculated by employing the CORR procedure. In addition, the linear regression equations were constructed by stepwise regression model in REG procedure. The results showed that pigs slaughtered in spring had the heaviest body weight among the four seasons (P < 0.05), pigs slaughtered in summer had the lowest backfat depth and shear force (P < 0.05), and pigs slaughtered in winter had the lowest drip loss (P < 0.05). The results showed more variation in backfat depth, drip loss, intramuscular fat content, and shear force, compared with other indicators across pigs. Body weight had a significant association with loin eye area, average backfat depth and L^⁎_{24 h} (P < 0.05). Furthermore, regression equations for drip loss, cooking loss, shear force, and intramuscular fat content were constructed using more accessible indicators. Collectively, this study provided an overall view of carcass and meat quality traits in a commercial pig population in China, and illustrated that season significantly affected carcass characteristics and meat quality traits independently of body weight.

Comparison of meat quality and glycolysis potential of two hybrid pigs in three-way hybrid model

With the improvement of consumers' requirements for pork quality, the method of crossbreeding with excellent local pig breeds to improve meat quality is popular. Saba pig has high reproduction rate, good meat quality and high utilization rate of roughage, but its excellent characteristics have not been fully developed and utilized. To promote the development and utilization of Saba pigs and production of high-quality pork, the meat quality traits and glycolysis potential of Duroc × (Landrace × Yorkshire) (DLY), Berkshire × (Duroc × Saba) (BDS), and Duroc × (Berkshire × Saba) (DBS) three-way crossbred pigs were compared. The results showed that DLY had the highest live weight, carcass weight, lean meat percentage, drip loss, glycolysis potential, muscle diameter, and relative mRNA expression levels of type IIb muscle fibers as well as the lowest ultimate pH (p < 0.05). The lightness value of DBS was the highest (p < 0.05). Among the three crossbred pigs, myristic, arachidic, palmitoleic, and eicosenoic acids were the highest in BDS. These results indicated that the carcass traits of local crossbred pigs were worse than those of DLY pigs, but meat quality was markedly higher, with BDS showing the best meat quality.

Genome-Wide Association Analysis Reveals Novel Loci Related with Visual Score Traits in Nellore Cattle Raised in Pasture-Based Systems

Body conformation traits assessed based on visual scores are widely used in Zebu cattle breeding programs. The aim of this study was to identify genomic regions and biological pathways associated with body conformation (CONF), finishing precocity (PREC), and muscling (MUSC) in Nellore cattle. The measurements based on visual scores were collected in 20,807 animals raised in pasture-based systems in Brazil. In addition, 2775 animals were genotyped using a 35 K SNP chip, which contained 31,737 single nucleotide polymorphisms after quality control. Single-step GWAS was performed using the BLUPF90 software while candidate genes were identified based on the Ensembl Genes 69. PANTHER and REVIGO platforms were used to identify key biological pathways and STRING to create gene networks. Novel candidate genes were revealed associated with CONF, including ALDH9A1, RXRG, RAB2A, and CYP7A1, involved in lipid metabolism. The genes associated with PREC were ELOVL5, PID1, DNER, TRIP12, and PLCB4, which are related to the synthesis of long-chain fatty acids, lipid metabolism, and muscle differentiation. For MUSC, the most important genes associated with muscle development were SEMA6A, TIAM2, UNC5A, and UIMC1. The polymorphisms identified in this study can be incorporated in commercial genotyping panels to improve the accuracy of genomic evaluations for visual scores in beef cattle.

The Genetic Architecture of Meat Quality Traits in a Crossbred Commercial Pig Population

Meat quality is of importance in consumer acceptance and purchasing tendency of pork. However, the genetic architecture of pork meat quality traits remains elusive. Herein, we conducted genome-wide association studies to detect single nucleotide polymorphisms (SNPs) and genes affecting meat pH and meat color (L*, lightness; a*, redness; b*, yellowness) in 1518 three-way crossbred pigs. All individuals were genotyped using the GeneSeek Porcine 50K BeadChip. In sum, 30 SNPs and 20 genes are found to be associated with eight meat quality traits. Notably, we detect one significant quantitative trait locus (QTL) on SSC15 with a 143 kb interval for meat pH (pH_12h), together with the most promising candidate TNS1. Interestingly, two newly identified SNPs located in the TTLL4 gene demonstrate the highest phenotypic variance of pH_12h in this QTL, at 2.67%. The identified SNPs are useful for the genetic improvement of meat quality traits in pigs by assigning higher weights to associated SNPs in genomic selection.

Enhancer-promoter interaction maps provide insights into skeletal muscle-related traits in pig genome

Background

Gene expression programs are intimately linked to the interplay of active cis regulatory elements mediated by chromatin contacts and associated RNAs. Genome-wide association studies (GWAS) have identified many variants in these regulatory elements that can contribute to phenotypic diversity. However, the functional interpretation of these variants remains nontrivial due to the lack of chromatin contact information or limited contact resolution. Furthermore, the distribution and role of chromatin-associated RNAs in gene expression and chromatin conformation remain poorly understood. To address this, we first present a comprehensive interaction map of nuclear dynamics of 3D chromatin-chromatin interactions (H3K27ac BL-HiChIP) and RNA-chromatin interactions (GRID-seq) to reveal genomic variants that contribute to complex skeletal muscle traits.

Results

In a genome-wide scan, we provide systematic fine mapping and gene prioritization from GWAS leading signals that underlie phenotypic variability of growth rate, meat quality, and carcass performance. A set of candidate functional variants and 54 target genes previously not detected were identified, with 71% of these candidate functional variants choosing to skip over their nearest gene to regulate the target gene in a long-range manner. The effects of three functional variants regulating KLF6 (related to days to 100 kg), MXRA8 (related to lean meat percentage), and TAF11 (related to loin muscle depth) were observed in two pig populations. Moreover, we find that this multi-omics interaction map consists of functional communities that are enriched in specific biological functions, and GWAS target genes can serve as core genes for exploring peripheral trait-relevant genes.

Conclusions

Our results provide a valuable resource of candidate functional variants for complex skeletal muscle-related traits and establish an integrated approach to complement existing 3D genomics by exploiting RNA-chromatin and chromatin-chromatin interactions for future association studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01322-2.

Identification of two novel pathogenic variants of the NR1H4 gene in intrahepatic cholestasis of pregnancy patients

Background

Intrahepatic cholestasis of pregnancy (ICP) can cause adverse pregnancy outcomes, such as spontaneous preterm delivery and stillbirth. It is a complex disease influenced by multiple factors, including genetics and the environment. Previous studies have reported that functioning nuclear receptor subfamily 1 group H member 4 (NR1H4) plays an essential role in bile acid (BA) homeostasis. However, some novel variants and their pathogenesis have not been fully elucidated. Therefore, this research aimed to investigate the genetic characteristics of the NR1H4 gene in ICP.

Methods

In this study, we sequenced the entire coding region of NR1H4 in 197 pregnant women with ICP disease. SIFT and PolyPhen2 were used to predict protein changes. Protein structure modelling and comparisons between NR1H4 reference and modified protein structures were performed by SWISS-MODEL and Chimera 1.14rc, respectively. T-tests were used to analyse the potential significant differences between NR1H4 mutations and wild types for 29 clinical features. Fisher’s test was conducted to test the significance of differences in mutation frequencies between ICP and the three databases.

Results

We identified four mutations: two novel missense mutations, p.S145F and p.M185L; rs180957965 (A230S); and rs147030757 (N275N). The two novel missense mutations were absent in 1029 controls and three databases, including the 1000 Genomes Project (1000G_ALL), Exome Aggregation Consortium (ExAC) and ChinaMAP. Two web-available tools, SIFT and PolyPhen2, predicted that these mutations are harmful to the function of the protein. Moreover, compared to the wild-type protein structure, the NR1H4 p.S145F and p.M185L protein structure showed a slight change in the chemical bond in two zinc finger structures. Combined clinical data indicate that the mutation group had higher levels of total bile acid (TBA) than the wild-type group. Therefore, we hypothesized that these two mutations altered the protein structure of NR1H4, which impaired the function of NR1H4 itself and its target gene and caused an increase in TBA.

Conclusions

To our knowledge, this is the first study to identify the novel p.S145F and p.M185L mutations in 197 ICP patients. Our present study provides new insights into the genetic architecture of ICP involving the two novel NR1H4 mutations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01240-w.

Genome-wide association study reveals a quantitative trait locus and two candidate genes on Sus scrofa chromosome 5 affecting intramuscular fat content in Suhuai pigs

Intramuscular fat content (IFC) is an essential quantitative trait of meat, affecting multiple meat quality indicators. A certain amount of IFC could not only improve the sensory score of pork but also increase the flavour, tenderness, juiciness and shelf-life. To dissect the genetic determinants of IFC, two methods, including genome-wide efficient mixed-model analysis (GEMMA) and linkage disequilibrium adjusted kinships (LDAKs), were used to carry out genome-wide association studies for IFC in Suhuai pig population. A total of 14 and 18 significant single nucleotide polymorphisms (SNPs) were identified by GEMMA and LDAK, respectively. The results of these two methods were highly consistent and all 14 significant SNPs in GEMMA were detected by LDAK. Seven of the 18 SNPs reached the genome-wide significance level (P < 9.85E-07) while 11 cases reached the suggestive significance level (P < 1.77E-05). These significant SNPs were mainly distributed on Sus scrofa chromosome (SSC) 5, 3, and 7. Moreover, one locus resides in a 2.27 Mb (71.37-73.64 Mb) region on SSC5 harbouring 13 significant SNPs associated with IFC, and the lead SNP (rs81302978) also locates in this region. Linkage disequilibrium (LD) analysis showed that there were four pairs of complete LD (r² = 1) among these 13 SNPs, and the remaining 9 SNPs with incomplete LD (r² ≠ 1) were selected for subsequent analyses of IFC. Association analyses showed that 7 out of 9 SNPs were significantly associated with IFC (P < 0.05) in 330 Suhuai pigs, and the other 2 SNPs tended to reach a significant association level with IFC (P < 0.1). The phenotypic variance explained (PVE) range of these 9 SNPs was 0.92-3.55%. Meanwhile, the lead SNP was also significantly associated (rs81302978) with IFC (P < 0.05) in 378 commercial hybrid pigs (Pietrain × Duroc) × (Landrace × Yorkshire) (PDLY), and the PVE was 1.38%. Besides, two lipid metabolism-relevant candidate genes, the leucine rich repeat kinase 2 (LRRK2) and PDZ domain containing ring finger 4 (PDZRN4) were identified in the 2.27 Mb region on SSC5. In conclusion, our results may provide a set of markers useful for genetic improvement of IFC in pigs and will advance the genome selection process of IFC on pig breeding programmes.

Whole-exome sequencing reveals ANO8 as a genetic risk factor for intrahepatic cholestasis of pregnancy

Background

Intrahepatic cholestasis of pregnancy (ICP) is characterized by pruritus and cholestasis in late pregnancy and results in adverse pregnancy outcomes, including preterm delivery and birth weight, which are affected by the genetic and environmental background. However, until now, the genetic architecture of ICP has remained largely unclear.

Methods

Twenty-six clinical data points were recorded for 151 Chinese ICP patients. The data generated from whole-exome sequencing (WES) using the BGISEQ-500 platform were further analyzed by Burrows-Wheeler Aligner (BWA) software, Genome Analysis Toolkit (GATK), ANNOVAR tool, etc. R packages were used to conduct t-test, Fisher’s test and receiver operating characteristic (ROC) curve analyses.

Results

We identified eighteen possible pathogenic loci associated with ICP disease in known genes, covering ABCB4, ABCB11, ATP8B1 and TJP2. The loci Lys386Gln, Gly527Gln and Trp708Ter in ABCB4, Leu589Met, Gln605Pro and Gln1194Ter in ABCB11, and Arg189Ser in TJP2 were novel discoveries. In addition, WES analysis indicated that the gene ANO8 involved in the transport of bile salts is newly identified as associated with ICP. The functional network of the ANO8 gene confirmed this finding. ANO8 contained 8 rare missense mutations that were found in eight patients among the 151 cases and were absent from 1029 controls. Out of the eight SNPs, 3 were known, and the remaining five are newly identified. These variants have a low frequency, ranging from 0.000008 to 0.00001 in the ExAC, gnomAD – Genomes and TOPMED databases. Bioinformatics analysis showed that the sites and their corresponding amino acids were both highly conserved among vertebrates. Moreover, the influences of all the mutations on protein function were predicted to be damaging by the SIFT tool. Combining clinical data, it was found that the mutation group (93.36 µmol/L) had significantly (P = 0.038) higher total bile acid (TBA) levels than the wild-type group (40.81 µmol/L).

Conclusions

To the best of our knowledge, this is the first study to employ WES technology to detect genetic loci for ICP. Our results provide new insights into the genetic basis of ICP and will benefit the final identification of the underlying mutations.