Cpt1a gene expression in peripheral blood mononuclear cells as an early biomarker of diet-related metabolic alterations

DNA Methylation Near CPT1A and Changes in Triglyceride-rich Lipoproteins in Response to Weight-loss Diet Interventions

CONTEXT

Carnitine palmitoyltransferase-1A, encoded by the CPT1A gene, plays a key role in the oxidation of long-chain fatty acids in the mitochondria and may be important in triglyceride metabolism. Previous work has shown that high fat intake was negatively associated with CPT1A methylation and positively associated with CPT1A expression.

OBJECTIVE

We aim to investigate the association of DNA methylation (DNAm) at the CPT1A gene with reductions in triglycerides and triglyceride-rich lipoproteins (TRLs) in response to weight-loss diet interventions.

METHODS

The current study included 538 White participants, who were randomly assigned to 1 of 4 diets varying in macronutrient components. We defined the regional DNAm at CPT1A as the average methylation level over CpGs within 500 bp of the 3 triglyceride-related DNAm sites.

RESULTS

Dietary fat intake significantly modified the association between baseline DNAm at CPT1A and 2-year changes in total plasma triglycerides, independent of concurrent weight loss. Among participants assigned to a low-fat diet, a higher regional DNAm level at CPT1A was associated with a greater reduction in total plasma triglycerides at 2 years (P = .01), compared with those assigned to a high-fat diet (P = .64) (P interaction = .018). Further investigation on lipids and apolipoproteins in very low-density lipoprotein (VLDL) revealed similar interaction patterns for 2-year changes in VLDL-triglycerides, VLDL-cholesterol, and VLDL-apolipoprotein B (P interaction = .009, .002, and .016, respectively), but not for VLDL-apoC-III (P interaction = .36).

CONCLUSION

Participants with a higher regional DNAm level at CPT1A benefit more in long-term improvement in triglycerides, particularly in the TRLs and related apolipoproteins when consuming a low-fat weight-loss diet.

Mitochondrial reprogramming in peripheral blood mononuclear cells of patients with glycogen storage disease type Ia

BACKGROUND

Glycogen storage disease type Ia (GSDIa) is an inborn metabolic disorder caused by the deficiency of glucose-6-phospatase-α (G6Pase-α) leading to mitochondrial dysfunction. It remains unclear whether mitochondrial dysfunction is present in patients' peripheral blood mononuclear cells (PBMC) and whether dietary treatment can play a role. The aim of this study was to investigate mitochondrial function in PBMC of GSDIa patients.

METHODS

Ten GSDIa patients and 10 age-, sex- and fasting-time matched controls were enrolled. Expression of genes involved in mitochondrial function and activity of key fatty acid oxidation (FAO) and Krebs cycle proteins were assessed in PBMC. Targeted metabolomics and assessment of metabolic control markers were also performed.

RESULTS

Adult GSDIa patients showed increased CPT1A, SDHB, TFAM, mTOR expression (p < 0.05) and increased VLCAD, CPT2 and citrate synthase activity in PBMC (p < 0.05). VLCAD activity directly correlated with WC (p < 0.01), BMI (p < 0.05), serum malonycarnitine levels (p < 0.05). CPT2 activity directly correlated with BMI (p < 0.05).

CONCLUSION

Mitochondrial reprogramming is detectable in PBMC of GSDIa patients. This feature may develop as an adaptation to the liver enzyme defect and may be triggered by dietary (over)treatment in the frame of G6Pase-α deficiency. PBMC can represent an adequate mean to assess (diet-induced) metabolic disturbances in GSDIa.

Gene expression associations with body mass index in the Multi-Ethnic Study of Atherosclerosis

BACKGROUND/OBJECTIVES

Obesity, defined as excessive fat accumulation that represents a health risk, is increasing in adults and children, reaching global epidemic proportions. Body mass index (BMI) correlates with body fat and future health risk, yet differs in prediction by fat distribution, across populations and by age. Nonetheless, few genetic studies of BMI have been conducted in ancestrally diverse populations. Gene expression association with BMI was assessed in the Multi-Ethnic Study of Atherosclerosis (MESA) in four self-identified race and ethnicity (SIRE) groups to identify genes associated with obesity.

SUBJECTS/METHODS

RNA-sequencing was performed on 1096 MESA participants (37.8% white, 24.3% Hispanic, 28.4% African American, and 9.5% Chinese American) and linear models were used to assess the association of expression from each gene for its effect on BMI, adjusting for age, sex, sequencing center, study site, five expression and four genetic principal components in each self-identified race group. Sample-size-weighted meta-analysis was performed to identify genes with BMI-associated expression across ancestry groups.

RESULTS

Within individual SIRE groups, there were zero to three genes whose expression is significantly (p < 1.97 × 10^-6) associated with BMI. Across all groups, 45 genes were identified by meta-analysis whose expression was significantly associated with BMI, explaining 29.7% of BMI variation. The 45 genes are expressed in a variety of tissues and cell types and are enriched for obesity-related processes including erythrocyte function, oxygen binding and transport, and JAK-STAT signaling.

CONCLUSIONS

We have identified genes whose expression is significantly associated with obesity in a multi-ethnic cohort. We have identified novel genes associated with BMI as well as confirmed previously identified genes from earlier genetic analyses. These novel genes and their biological pathways represent new targets for understanding the biology of obesity as well as new therapeutic intervention to reduce obesity and improve global public health.

The potential of DNA methylation as a biomarker for obesity and smoking

DNA methylation is an epigenetic modification that has consistently been shown to be linked with a variety of human traits and diseases. Because DNA methylation is dynamic and potentially reversible in nature and can reflect environmental exposures and predict the onset of diseases, it has piqued interest as a potential disease biomarker. DNA methylation patterns are more stable than transcriptomic or proteomic patterns, and they are relatively easy to measure to track exposure to different environments and risk factors. Importantly, technologies for DNA methylation quantification have become increasingly cost effective-accelerating new research in the field-and have enabled the development of novel DNA methylation biomarkers. Quite a few DNA methylation-based predictors for a number of traits and diseases already exist. Such predictors show potential for being more accurate than self-reported or measured phenotypes (such as smoking behavior and body mass index) and may even hold potential for applications in clinics. In this review, we will first discuss the advantages and challenges of DNA methylation biomarkers in general. We will then review the current state and future potential of DNA methylation biomarkers in two human traits that show rather consistent alterations in methylome-obesity and smoking. Lastly, we will briefly speculate about the future prospects of DNA methylation biomarkers, and possible ways to achieve them.

Role of cellular retinol-binding protein, type 1 and retinoid homeostasis in the adult mouse heart: A multi-omic approach

The main active metabolite of Vitamin A, all-trans retinoic acid (RA), is required for proper cellular function and tissue organization. Heart development has a well-defined requirement for RA, but there is limited research on the role of RA in the adult heart. Homeostasis of RA includes regulation of membrane receptors, chaperones, enzymes, and nuclear receptors. Cellular retinol-binding protein, type 1 (CRBP1), encoded by retinol-binding protein, type 1 (Rbp1), regulates RA homeostasis by delivering vitamin A to enzymes for RA synthesis and protecting it from non-specific oxidation. In this work, a multi-omics approach was used to characterize the effect of CRBP1 loss using the Rbp1^-/- mouse. Retinoid homeostasis was disrupted in Rbp1^-/- mouse heart tissue, as seen by a 33% and 24% decrease in RA levels in the left and right ventricles, respectively, compared to wild-type mice (WT). To further inform on the effect of disrupted RA homeostasis, we conducted high-throughput targeted metabolomics. A total of 222 metabolite and metabolite combinations were analyzed, with 33 having differential abundance between Rbp1^-/- and WT hearts. Additionally, we performed global proteome profiling to further characterize the impact of CRBP1 loss in adult mouse hearts. More than 2606 unique proteins were identified, with 340 proteins having differential expression between Rbp1^-/- and WT hearts. Pathway analysis performed on metabolomic and proteomic data revealed pathways related to cellular metabolism and cardiac metabolism were the most disrupted in Rbp1^-/- mice. Together, these studies characterize the effect of CRBP1 loss and reduced RA in the adult heart.

Genome-Wide Expression Profiling of mRNAs, lncRNAs and circRNAs in Skeletal Muscle of Two Different Pig Breeds

Simple Summary

Variation exists in muscle-related traits, such as muscle growth and meat quality, between obese and lean pigs. In this study, the transcriptome profiles of skeletal muscle between Beijing Blackand Yorkshire pigs were characterized to explore the molecular mechanism underlying skeletal muscle-relatedtraits. Gene Ontology (GO) and KEGG pathway enrichment analyses showed that differentially expressed mRNAs, lncRNAs, and circRNAs involved in skeletal muscle development and fatty acid metabolism played a key role in the determination of muscle-related traits between different pig breeds. These results provide candidate genes responsible for muscle phenotypic variation and are valuable for pig breeding.

Abstract

RNA-Seq technology is widely used to analyze global changes in the transcriptome and investigate the influence on relevant phenotypic traits. Beijing Black pigs show differences in growth rate and meat quality compared to western pig breeds. However, the molecular mechanisms responsible for such phenotypic differences remain unknown. In this study, longissimus dorsi muscles from Beijing Black and Yorkshire pigs were used to construct RNA libraries and perform RNA-seq. Significantly different expressions were observed in 1051 mRNAs, 322 lncRNAs, and 82 circRNAs. GO and KEGG pathway annotation showed that differentially expressed mRNAs participated in skeletal muscle development and fatty acid metabolism, which determined the muscle-related traits. To explore the regulatory role of lncRNAs, the cis and trans-target genes were predicted and these lncRNAswere involved in the biological processes related to skeletal muscle development and fatty acid metabolismvia their target genes. CircRNAs play a ceRNA role by binding to miRNAs. Therefore, the potential miRNAs of differentially expressed circRNAs were predicted and interaction networks among circRNAs, miRNAs, and key regulatory mRNAs were constructed to illustrate the function of circRNAs underlying skeletal muscle development and fatty acid metabolism. This study provides new clues for elucidating muscle phenotypic variation in pigs.

Use of human PBMC to analyse the impact of obesity on lipid metabolism and metabolic status: a proof-of-concept pilot study

Peripheral blood mononuclear cells (PBMC) are widely used as a biomarker source in nutrition/obesity studies because they reflect gene expression profiles of internal tissues. In this pilot proof-of-concept study we analysed in humans if, as we previously suggested in rodents, PBMC could be a surrogate tissue to study overweight/obesity impact on lipid metabolism. Pre-selected key lipid metabolism genes based in our previous preclinical studies were analysed in PBMC of normoglycemic normal-weight (NW), and overweight-obese (OW-OB) subjects before and after a 6-month weight-loss plan. PBMC mRNA levels of CPT1A, FASN and SREBP-1c increased in the OW-OB group, according with what described in liver and adipose tissue of humans with obesity. This altered expression pattern was related to increased adiposity and early signs of metabolic impairment. Greater weight loss and/or metabolic improvement as result of the intervention was related to lower CPT1A, FASN and SREBP-1c gene expression in an adjusted linear mixed-effects regression analysis, although no gene expression recovery was observed when considering mean comparisons. Thus, human PBMC reflect lipid metabolism expression profile of energy homeostatic tissues, and early obesity-related alterations in metabolic at-risk subjects. Further studies are needed to understand PBMC usefulness for analysis of metabolic recovery in weigh management programs.

Lipid Phenotypes and DNA Methylation: a Review of the Literature

PURPOSE OF REVIEW

Epigenetic modifications via DNA methylation have previously been linked to blood lipid levels, dyslipidemias, and atherosclerosis. The purpose of this review is to discuss current literature on the role of DNA methylation on lipid traits and their associated pathologies.

RECENT FINDINGS

Candidate gene and epigenome-wide approaches have identified differential methylation of genes associated with lipid traits (particularly CPT1A, ABCG1, SREBF1), and novel approaches are being implemented to further characterize these relationships. Moreover, studies on environmental factors have shown that methylation variations at lipid-related genes are associated with diet and pollution exposure. Further investigation is needed to elucidate the directionality of the associations between the environment, lipid traits, and epigenome. Future studies should also seek to increase the diversity of cohorts, as European and Asian ancestry populations are the predominant study populations in the current literature.

CUN-BAE Index as a Screening Tool to Identify Increased Metabolic Risk in Apparently Healthy Normal-Weight Adults and Those with Obesity

BACKGROUND

Imbalanced dietary intake is related to increased adiposity, which is linked to increased metabolic risk even in the absence of obesity. BMI is traditionally used to classify body fatness and weight range, but it only considers body weight and height. The Clínica Universidad de Navarra-Body Adiposity Estimator (CUN-BAE) equation has appeared as an additional tool to estimate adiposity considering also other relevant parameters, i.e., sex and age.

OBJECTIVES

We aimed to determine whether the CUN-BAE index could estimate adiposity-related metabolic risk in apparently healthy, normoglycemic adults.

METHODS

In this case-control study, men and women (18-45 y old) were classified as normal-weight (NW) [n = 20; BMI (in kg/m2) <25] or overweight-obese (OW-OB) (n = 34; BMI ≥25). The primary outcome was body fat content and clinical circulating parameters to assess by correlation analysis CUN-BAE's usefulness as a predictor of metabolic risk. In addition, transcriptomic biomarkers of lipid metabolism were analyzed in peripheral blood mononuclear cells (PBMCs) as secondary outcome indicators of metabolic impairment. Data were analyzed by correlation analysis and comparison of means.

RESULTS

CUN-BAE values correlated directly with body fatness obtained by DXA (r = 0.89, P < 0.01), with classical molecular biomarkers of metabolic risk, and with PBMC gene expression of carnitine palmitoyltransferase 1A (CPT1A), sterol regulatory element binding transcription factor 1c (SREBP-1c), and fatty acid synthase (FASN), early markers of metabolic impairment (P < 0.05). Moreover, CUN-BAE allowed identification of NW individuals with excessive body fatness, who were not yet presenting obesity-related molecular alterations. In these subjects, visceral fat correlated directly with circulating glucose, triglycerides, and total and LDL cholesterol, and with triglyceride-glucose and fatty liver indexes (P < 0.05). This is indicative of a metabolically obese NW phenotype.

CONCLUSIONS

Data obtained in our cohort of young normoglycemic volunteers support the use of the CUN-BAE index as a tool to estimate accurately body fat mass, but also as a first easy/effective screening tool to identify lean people with increased fat mass and increased metabolic risk.This trial was registered at clinicaltrials.gov as NCT04402697.

Deregulation of lncRNA HIST1H2AG-6 and AIM1-3 in peripheral blood mononuclear cells is associated with newly diagnosed type 2 diabetes

BACKGROUND

Type 2 diabetes mellitus (T2DM) is mainly affected by genetic and environmental factors; however, the correlation of long noncoding RNAs (lncRNAs) with T2DM remains largely unknown.

METHODS

Microarray analysis was performed to identify the differentially expressed lncRNAs and messenger RNAs (mRNAs) in patients with T2DM and healthy controls, and the expression of two candidate lncRNAs (lnc-HIST1H2AG-6 and lnc-AIM1-3) were further validated using quantitative real-time polymerase chain reaction (qRT-PCR). Spearman's rank correlation coefficient was used to measure the degree of association between the two candidate lncRNAs and differentially expressed mRNAs. Furthermore, the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway and GO (Gene Ontology) enrichment analysis were used to reveal the biological functions of the two candidate lncRNAs. Additionally, multivariate logistic regression analysis and receiver operating characteristic (ROC) curve analysis were performed.

RESULTS

The microarray analysis revealed that there were 55 lncRNAs and 36 mRNAs differentially expressed in patients with T2DM compared with healthy controls. Notably, lnc-HIST1H2AG-6 was significantly upregulated and lnc-AIM1-3 was significantly downregulated in patients with T2DM, which was validated in a large-scale qRT-PCR examination (90 controls and 100 patients with T2DM). Spearman's rank correlation coefficient revealed that both lncRNAs were correlated with 36 differentially expressed mRNAs. Furthermore, functional enrichment (KEGG and GO) analysis demonstrated that the two lncRNA-related mRNAs might be involved in multiple biological functions, including cell programmed death, negative regulation of insulin receptor signal, and starch and sucrose metabolism. Multivariate logistic regression analysis revealed that lnc-HIST1H2AG-6 and lnc-AIM1-3 were significantly correlated with T2DM (OR = 5.791 and 0.071, respectively, both P = 0.000). Furthermore, the ROC curve showed that the expression of lnc-HIST1H2AG-6 and lnc-AIM1-3 might be used to differentiate patients with T2DM from healthy controls (area under the ROC curve = 0.664 and 0.769, respectively).

CONCLUSION

The profiles of lncRNA and mRNA were significantly changed in patients with T2DM. The expression levels of lnc-HIST1H2AG-6 and lnc-AIM1-3 genes were significantly correlated with some features of T2DM, which may be used to distinguish patients with T2DM from healthy controls and may serve as potential novel biomarkers for diagnosis in the future.

Carbohydrate and fat intake associated with risk of metabolic diseases through epigenetics of CPT1A

BACKGROUND

Epigenome-wide association studies identified the cg00574958 DNA methylation site at the carnitine palmitoyltransferase-1A (CPT1A) gene to be associated with reduced risk of metabolic diseases (hypertriglyceridemia, obesity, type 2 diabetes, hypertension, metabolic syndrome), but the mechanism underlying these associations is unknown.

OBJECTIVES

We aimed to elucidate whether carbohydrate and fat intakes modulate cg00574958 methylation and the risk of metabolic diseases.

METHODS

We examined associations between carbohydrate (CHO) and fat (FAT) intake, as percentages of total diet energy, and the CHO/FAT ratio with CPT1A-cg00574958, and the risk of metabolic diseases in 3 populations (Genetics of Lipid Lowering Drugs and Diet Network, n = 978; Framingham Heart Study, n = 2331; and REgistre GIroní del COR study, n = 645) while adjusting for confounding factors. To understand possible causal effects of dietary intake on the risk of metabolic diseases, we performed meta-analysis, CPT1A transcription analysis, and mediation analysis with CHO and FAT intakes as exposures and cg00574958 methylation as the mediator.

RESULTS

We confirmed strong associations of cg00574958 methylation with metabolic phenotypes (BMI, triglyceride, glucose) and diseases in all 3 populations. Our results showed that CHO intake and CHO/FAT ratio were positively associated with cg00574958 methylation, whereas FAT intake was negatively correlated with cg00574958 methylation. Meta-analysis further confirmed this strong correlation, with β = 58.4 ± 7.27, P = 8.98 x 10-16 for CHO intake; β = -36.4 ± 5.95, P = 9.96 x 10-10 for FAT intake; and β = 3.30 ± 0.49, P = 1.48 x 10-11 for the CHO/FAT ratio. Furthermore, CPT1A mRNA expression was negatively associated with CHO intake, and positively associated with FAT intake, and metabolic phenotypes. Mediation analysis supports the hypothesis that CHO intake induces CPT1A methylation, hence reducing the risk of metabolic diseases, whereas FAT intake inhibits CPT1A methylation, thereby increasing the risk of metabolic diseases.

CONCLUSIONS

Our results suggest that the proportion of total energy supplied by CHO and FAT can have a causal effect on the risk of metabolic diseases via the epigenetic status of CPT1A.Study registration at https://www.clinicaltrials.gov/: the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN)-NCT01023750; and the Framingham Heart Study (FHS)-NCT00005121.

Impaired CPT1A Gene Expression Response to Retinoic Acid Treatment in Human PBMC as Predictor of Metabolic Risk

Ex vivo human peripheral blood mononuclear cell (PBMC) systems offer the possibility to test transcriptomic effects of food bioactive compounds with potential health effects. We investigated all-trans retinoic acid (ATRA) effect on mRNA expression of key lipid metabolism and inflammatory genes in PBMCs from normal-weight (NW) and overweight-obese (OW-OB) men with different metabolic syndrome-related features. PBMCs were incubated with 10 µM ATRA and mRNA levels of selected genes were analyzed using real-time RT-qPCR. Human ex vivo PBMCs responded to ATRA treatment, but the response for some genes was dependent on body mass index (BMI), with a lower response in PBMC from OW-OB than from NW donors. Moreover, gene expression response was affected by circulating high-density lipoprotein (HDL)-cholesterol levels. Particularly, the response to ATRA of CPT1A, previously reported as a sensitive metabolic risk predictive biomarker, was dependent on HDL levels and not on BMI, being impaired in those individuals with lower HDL levels, specifically in OW-OB. Thus, PBMCs' insensitivity to ATRA, which can be considered as indicative of impaired metabolism, was observed in individuals with higher metabolic risk (OW-OB with low HDL levels). In conclusion, an ex vivo human PBMC system indicates that ATRA response could be influenced by metabolic syndrome features. Moreover, our study reinforces the role of CPT1A as a marker of metabolic risk and points to plasmatic HDL-cholesterol levels as a parameter to take into consideration when the effects of nutritional factors and/or dietary interventions on humans are under study. Further studies including women are required to detect potential gender differences in the observed effects.

Effects of Breastfeeding and Formula Feeding on the Expression Level of FTO, CPT1A and PPAR-α Genes in Healthy Infants

PURPOSE

The study aimed to investigate the effect of breastfeeding, formula feeding and mix feeding (breastfed plus formula-fed) on the expression level of obesity-predisposing genes including fat mass and obesity-associated (FTO), carnitine palmitoyltransferase 1A (CPT1A), and peroxisome proliferator-activated receptor-α (PPAR-α) in 5- to 6-month-old infants.

PATIENTS AND METHODS

A total of 150 infants participated in this case-control study. All subjects were healthy infants aged 5-6 months that divided into 3 groups: breastfed, formula-fed, and mix-fed. The expression level of FTO, CPT1A, and PPAR-α genes in peripheral blood mononuclear cells (PBMC) was evaluated in each group using reverse transcription-polymerase chain reaction (RT-PCR) method.

RESULTS

Our findings showed that the current weight, height, and head circumference of infants in the formula feeding and mix feeding groups were significantly higher than those in the exclusive breastfeeding group. The expression level of FTO and CPT1A genes in formula-fed and mix-fed infants was significantly higher (p<0.001) than that in breastfed infants, while the expression level of PPAR-α gene was significantly lower (p<0.05).

CONCLUSION

Breastfeeding showed modulatory effects on the expression level of obesity-predisposing genes and can protect against obesity and subsequent non-communicable diseases. However, more investigations are required to explain the epigenetic effects of breast milk.

Novel Markers of the Metabolic Impact of Exogenous Retinoic Acid with A Focus on Acylcarnitines and Amino Acids

Treatment with all-trans retinoic acid (ATRA), the carboxylic form of vitamin A, lowers body weight in rodents by promoting oxidative metabolism in multiple tissues including white and brown adipose tissues. We aimed to identify novel markers of the metabolic impact of ATRA through targeted blood metabolomics analyses, with a focus on acylcarnitines and amino acids. Blood was obtained from mice treated with a high ATRA dose (50 mg/kg body weight/day, subcutaneous injection) or placebo (controls) during the 4 days preceding collection. LC-MS/MS analyses with a focus on acylcarnitines and amino acids were conducted on plasma and PBMC. Main results showed that, relative to controls, ATRA-treated mice had in plasma: increased levels of carnitine, acetylcarnitine, and longer acylcarnitine species; decreased levels of citrulline, and increased global arginine bioavailability ratio for nitric oxide synthesis; increased levels of creatine, taurine and docosahexaenoic acid; and a decreased n-6/n-3 polyunsaturated fatty acids ratio. While some of these features likely reflect the stimulation of lipid mobilization and oxidation promoted by ATRA treatment systemically, other may also play a causal role underlying ATRA actions. The results connect ATRA to specific nutrition-modulated biochemical pathways, and suggest novel mechanisms of action of vitamin A-derived retinoic acid on metabolic health.

Transcriptome Analysis of Three Sheep Intestinal Regions reveals Key Pathways and Hub Regulatory Genes of Large Intestinal Lipid Metabolism

The large intestine, also known as the hindgut, is an important part of the animal digestive system. Recent studies on digestive system development in ruminants have focused on the rumen and the small intestine, but the molecular mechanisms underlying sheep large intestine metabolism remain poorly understood. To identify genes related to intestinal metabolism and to reveal molecular regulation mechanisms, we sequenced and compared the transcriptomes of mucosal epithelial tissues among the cecum, proximal colon and duodenum. A total of 4,221 transcripts from 3,254 genes were identified as differentially expressed transcripts. Between the large intestine and duodenum, differentially expressed transcripts were found to be significantly enriched in 6 metabolism-related pathways, among which PPAR signaling was identified as a key pathway. Three genes, CPT1A, LPL and PCK1, were identified as higher expression hub genes in the large intestine. Between the cecum and colon, differentially expressed transcripts were significantly enriched in 5 lipid metabolism related pathways, and CEPT1 and MBOAT1 were identified as hub genes. This study provides important information regarding the molecular mechanisms of intestinal metabolism in sheep and may provide a basis for further study.