Integration of ATAC-seq and RNA-seq analysis identifies key genes affecting intramuscular fat content in pigs

Genetic basis of phenotypic convergence in pig terminal sires using pathway-based selection signature detection methods

The primary purpose of genetic improvement in lean pig breeds is to enhance production performance. Owing to their similar breeding directions, Duroc and Pietrain pigs are ideal models for investigating the phenotypic convergence underlying artificial selection. However, most important economic traits are controlled by a polygenic basis, so traditional strategies for detecting selection signatures may not fully reveal the genetic basis of complex traits. The pathway-based gene network analysis method utilizes each pathway as a unit, overcoming the limitations of traditional strategies for detecting selection signatures by revealing the selection of complex biological processes. Here, we utilized 13 122 398 high-quality SNPs from whole-genome sequencing data of 48 Pietrain pigs, 156 Duroc pigs and 36 European wild boars to detect selective signatures. After calculating FST and iHS scores, we integrated the pathway information and utilized the r/bioconductor graphite and signet packages to construct gene networks, identify subnets and uncover candidate genes underlying selection. Using the traditional strategy, a total of 47 genomic regions exhibiting parallel selection were identified. The enriched genes, including INO80, FZR1, LEPR and FAF1, may be associated with reproduction, fat deposition and skeletal development. Using the pathway-based selection signatures detection method, we identified two significant biological pathways and eight potential candidate genes underlying parallel selection, such as VTN, FN1 and ITGAV. This study presents a novel strategy for investigating the genetic basis of complex traits and elucidating the phenotypic convergence underlying artificial selection, by integrating traditional selection signature methods with pathway-based gene network analysis.

Identification of Heilongjiang crossbred beef cattle pedigrees and reveals functional genes related to economic traits based on whole-genome SNP data

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies.

Integrated ATAC-seq and RNA-seq Analysis of In Vitro Cultured Skeletal Muscle Satellite Cells to Understand Changes in Cell Proliferation

Skeletal muscle satellite cells, the resident stem cells in pig skeletal muscle, undergo proliferation and differentiation to enable muscle tissue repair. The proliferative and differentiative abilities of these cells gradually decrease during in vitro cultivation as the cell passage number increases. Despite extensive research, the precise molecular mechanisms that regulate this process are not fully understood. To bridge this knowledge gap, we conducted transcriptomic analysis of skeletal muscle satellite cells during in vitro cultivation to quantify passage number-dependent changes in the expression of genes associated with proliferation. Additionally, we explored the relationships between gene transcriptional activity and chromatin accessibility using transposase-accessible chromatin sequencing. This revealed the closure of numerous open chromatin regions, which were primarily located in intergenic regions, as the cell passage number increased. Integrated analysis of the transcriptomic and epigenomic data demonstrated a weak correlation between gene transcriptional activity and chromatin openness in expressed genic regions; although some genes (e.g., GNB4 and FGD5) showed consistent relationships between gene expression and chromatin openness, a substantial number of differentially expressed genes had no clear association with chromatin openness in expressed genic regions. The p53-p21-RB signaling pathway may play a critical regulatory role in cell proliferation processes. The combined transcriptomic and epigenomic approach taken here provided key insights into changes in gene expression and chromatin openness during in vitro cultivation of skeletal muscle satellite cells. These findings enhance our understanding of the intricate mechanisms underlying the decline in cellular proliferation capacity in cultured cells.

Translatomics reveals the role of dietary calcium addition in regulating muscle fat deposition in pigs

Intramuscular fat (IMF) in pork holds significant importance for economic performance within the pig industry and dietary calcium supplementation enhances the accumulation of intramuscular fat. Additionally, calcium ions inhibit translation and reduce protein synthesis. However, the mechanism by which calcium regulates IMF deposition in muscle through translation remains largely unknown. In this study, we compared the ribosome profiles of the longissimus dorsi muscles of Duroc × Landrace × Large white pigs from the normal calcium (NC) group or calcium supplement (HC) group by Ribo-seq, and RNA-seq. By integrating multiple-omics analysis, we further discovered 437 genes that were transcriptionally unchanged but translationally altered and these genes were significantly enriched in the oxidative phosphorylation signaling pathway. Furthermore, experimental data showed that inhibiting the expression of COX10 and mtND4L increased triglyceride accumulation in C2C12 cells, providing new targets for intramuscular fat deposition. Finally, this work links dietary calcium, translation regulation and IMF deposition, providing a new strategy for both meat quality and economic performance within the pig industry.

HOXA9 gene inhibits proliferation and differentiation and promotes apoptosis of bovine preadipocytes

BACKGROUND

Hox gene family is an important transcription factor that regulates cell process, and plays a role in the process of adipocytes differentiation and fat deposition. Previous transcriptome sequencing studies have indicated that the Homeobox A9 gene (HOXA9) is a candidate gene for regulating the process of bovine lipid metabolism, but the function and specific mechanism of action remain unclear. Therefore, this study aims to explore the role of HOXA9 in the proliferation, differentiation and apoptosis of bovine preadipocytes through gain-of-function and lose-of-function.

RESULT

It found HOXA9 highly expressed in bovine adipose tissue, and its expression level changed significantly during adipocytes differentiation process. It gave a hint that HOXA9 may be involved in the process of bovine lipid metabolism. The results of HOXA9 gain-of-function experiments indicated that HOXA9 appeared to act as a negative regulator not only in the differentiation but also in the proliferation of bovine preadipocytes, which is mainly reflected that overexpression of HOXA9 down-regulate the mRNA and protein expression level of PPARγ, CEBPα and FABP4 (P < 0.05). The mRNA expression level of CDK1, CDK2, PCNA, CCNA2, CCNB1, CCND1 and CCNE2, as well as the protein expression of CDK2 also significantly decreased. The decrease of lipid droplets content was the main characteristic of the phenotype (P < 0.01), which further supported the evidence that HOXA9 was a negative regulator of preadipocytes differentiation. The decrease of cell proliferation rate and EdU positive rate, as well as the limitation of transition of preadipocytes from G0/G1 phase to S phase also provided evidence for the inhibition of proliferation. Apart from this above, we noted an interesting phenomenon that overexpression of HOXA9 showed in a significant upregulation of both mRNA and protein level of apoptosis markers, accompanied by a significant increase in cell apoptosis rate. These data led us not to refute the fact that HOXA9 played an active regulatory role in apoptosis. HOXA9 loss-of-function experiments, however, yielded the opposite results. Considering that HOXA9 acts as a transcription factor, we predicted its target genes. Dual luciferase reporter assay system indicated that overexpression of HOXA9 inhibits activity of PCNA promoter.

CONCLUSION

Taken together, we demonstrated for the first time that HOXA9 played a role as a negative regulatory factor in the differentiation and proliferation of preadipocytes, but played a positive regulatory role in apoptosis, and it may play a regulatory role by targeting PCNA. This study provides basic data for further exploring the regulatory network of intramuscular fat deposition in bovine.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BACKGROUND

Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.

RESULTS

Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.

CONCLUSIONS

Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

Verification of Key Target Molecules for Intramuscular Fat Deposition and Screening of SNP Sites in Sheep from Small-Tail Han Sheep Breed and Its Cross with Suffolk

Intramuscular fat (IMF) is vital for meat tenderness and juiciness. This study aims to explore the IMF deposition mechanism and the related molecular markers in sheep. Two populations, Small-tail Han Sheep (STH) and STH × Suffolk (SFK) F1 (SFK × STH), were used as the research object. Histological staining techniques compared the differences in the longissimus dorsi muscle among populations. A combination of transcriptome sequencing and biological information analysis screened and identified IMF-related target genes. Further, sequencing technology was employed to detect SNP loci of target genes to evaluate their potential as genetic markers. Histological staining revealed that the muscle fiber gap in the SFK × STH F1 was larger and the IMF content was higher. Transcriptome analysis revealed that PIK3R1 and PPARA were candidate genes. Histological experiments revealed that the expressions of PIK3R1 mRNA and PPARA mRNA were lower in SFK × STH F1 compared with the STH. Meanwhile, PIK3R1 and PPARA proteins were located in intramuscular adipocytes and co-located with the lipid metabolism marker molecule (FASN). SNP locus analysis revealed a mutation site in exon 7 of the PIK3R1 gene, which served as a potential genetic marker for IMF deposition. This study's findings will provide a new direction for meat quality breeding in sheep.

Genome-wide detection of selection signatures in Jianli pigs reveals novel cis-regulatory haplotype in EDNRB associated with two-end black coat color

BACKGROUND

Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown.

RESULTS

In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function.

CONCLUSION

Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs.

Characteristics of Transcriptome and Metabolome Concerning Intramuscular Fat Content in Beijing Black Pigs

To study the characteristics of genes and metabolites related to intramuscular fat (IMF) content with less influence by breed background and individual differences, the skeletal muscle samples from 40 Beijing black pigs with either high or low IMF content were used to perform transcriptome and metabolome analyses. About 99 genes (twofold-change) were differentially expressed. Up-regulated genes in the high IMF pigs were mainly related to fat metabolism. The key genes in charge of IMF deposition are ADIPOQ, CIDEC, CYP4B1, DGAT2, LEP, OPRL1, PLIN1, SCD, and THRSP. KLHL40, TRAFD1, and HSPA6 were novel candidate genes for the IMF trait due to their high abundances. In the low IMF pigs, the differentially expressed genes involved in virus resistance were up-regulated. About 16 and 18 differential metabolites (1.5 fold-change) were obtained in the positive and negative modes, respectively. Pigs with low IMF had weaker fatty acid oxidation due to the down-regulation of various carnitines. Differentially expressed genes were more important in determining IMF deposition than differential metabolites because relatively few differential metabolites were obtained, and they were merely the products under the physiological status of diverged IMF content. This study provided valuable information for further studies on IMF deposition.

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.

An Integrative ATAC-Seq and RNA-Seq Analysis of the Endometrial Tissues of Meishan and Duroc Pigs

Meishan pigs are a well-known indigenous pig breed in China characterized by a high fertility. Notably, the number of endometrial grands is significantly higher in Meishan pigs than Duroc pigs. The characteristics of the endometrial tissue are related to litter size. Therefore, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-sequencing (RNA-seq) to analyze the mechanisms underlying the differences in fecundity between the breeds. We detected the key transcription factors, including Double homeobox (Dux), Ladybird-like homeobox gene 2 (LBX2), and LIM homeobox 8 (Lhx8), with potentially pivotal roles in the regulation of the genes related to endometrial development. We identified the differentially expressed genes between the breeds, including SOX17, ANXA4, DLX3, DMRT1, FLNB, IRF6, CBFA2T2, TFCP2L1, EFNA5, SLIT2, and CYFIP2, with roles in epithelial cell differentiation, fertility, and ovulation. Interestingly, ANXA4, CBFA2T2, and TFCP2L1, which were upregulated in the Meishan pigs in the RNA-seq analysis, were identified again by the integration of the ATAC-seq and RNA-seq data. Moreover, we identified genes in the cancer or immune pathways, FoxO signaling, Wnt signaling, and phospholipase D signaling pathways. These ATAC-seq and RNA-seq analyses revealed the accessible chromatin and potential mechanisms underlying the differences in the endometrial tissues between the two types of pigs.

Single-nucleus and bulk RNA sequencing reveal cellular and transcriptional mechanisms underlying lipid dynamics in high marbled pork

Pork is the most consumed meat in the world, and its quality is associated with human health. Intramuscular fat (IMF) deposition (also called marbling) is a key factor positively correlated with various quality traits and lipo-nutritional values of meat. However, the cell dynamics and transcriptional programs underlying lipid deposition in highly marbled meat are still unclear. Here, we used Laiwu pigs with high (HLW) or low (LLW) IMF contents to explore the cellular and transcriptional mechanisms underlying lipid deposition in highly-marbled pork by single-nucleus RNA sequencing (snRNA-seq) and bulk RNA sequencing. The HLW group had higher IMF contents but less drip loss than the LLW group. Lipidomics results revelled the changes of overall lipid classes composition (e.g., glycerolipids including triglycerides, diglycerides, and monoglycerides; sphingolipids including ceramides and monohexose ceramide significantly increased) between HLW and LLW groups. SnRNA-seq revealed nine distinct cell clusters, and the HLW group had a higher percentage of adipocytes (1.40% vs. 0.17%) than the LLW group. We identified 3 subpopulations of adipocytes, including PDE4D⁺/PDE7B⁺ (in HLW and LLW), DGAT2⁺/SCD⁺ (mostly in HLW) and FABP5⁺/SIAH1⁺ cells (mostly in HLW). Moreover, we showed that fibro/adipogenic progenitors could differentiate into IMF cells and contribute to 43.35% of adipocytes in mice. In addition, RNA-seq revealed different genes involved in lipid metabolism and fatty acid elongation. Our study provides new insights into the cellular and molecular signatures of marbling formation; such knowledge may facilitate the development of new strategies to increase IMF deposition and the lipo-nutritional quality of high marbled pork.