Comparison of backfat thickness and serum estradiol levels between OVX and Sham groups
We observed that the backfat thickness of female pigs in the OVX group was remarkably higher than that in the Sham group (
P<0.05), while the estradiol level was significantly lower than that in the sham-operated pigs (472.49±27.45 ng/L VS. 238.71±11.30 ng/L) (
P<0.05, Fig 1), which indicates that removal of ovary lead to subcutaneous fat accumulation in pigs. Given that the liver was essential to regulate glucose and lipid metabolism, alterations in hepatic physiological functions may be responsible for obesity of pigs.
Comparison of hepatic lipid metabolism indexes and histopathology between OVX and sham groups
Compared with sham-operated controls, OVX pigs exhibited markedly increased levels of TC and TG
and reduced levels of HDL-C (
P<0.05, Fig 2). These results were consistent with previous studies in mice and rats (
Zhu et al., 2023). High TG caused by any reason can lead to hepatic steatosis (
Miller 2000). It has been reported that the plasma TC was increased in the aging humans and rodents since their physiological function in eliminating cholesterol decreased (
Lei et al., 2021). Several studies have demonstrated that ectopic accumulation of lipid within liver could specifically cause hepatic insulin resistance in humans and rodents (
Fabbrini et al., 2009). A previous study found that ovariectomized mice had increased monounsaturated fatty acid levels of hepatic triglyceride in the liver (
Jackson et al., 2011).
Histomorphologic analysis showed that OVX females had mor1e lipid drops accumulated in the livers compared to Sham pigs (Fig 3). OVX mice commonly demonstrate increased body weight owing to excessive fat accumulation, as well as increased liver lipid accumulation, leading to fatty liver (
Nanashima et al., 2020). Hepatic cholesterol accumulation is driven by a deeply deranged cellular cholesterol homeostasis, characterized by elevated cholesterol synthesis and uptake from circulating lipoproteins and by a reduced cholesterol excretion (
Musso et al., 2013).
Transcriptomic analysis of liver tissues of OVX and Sham groups
In the present study, average number of raw reads, clean reads and mapped reads in the liver of Sham and OVX groups were shown in Supplementary Table 1. By alignment analysis, 21,832 genes commonly expressed in the liver tissues of the OVX and Sham female pigs. A total of 414 genes were identified as differentially expressed genes (DEGs) between two groups under the thresholds of
P<0.05 and |log2Foldchange| >1, of which 130 and 314 DEGs were up-regulated and down-regulated in the OVX group. Hierarchical clustering showed a very clear separation of DEGs between the OVX and Sham groups (Fig 4).
Gene ontology (GO) enrichment and KEGG pathway analysis of all DEGs were carried out between the two groups (Fig 5). GO enrichment analysis revealed that main enriched biological processes were cell-cell signaling, G protein-coupled adenosine receptor signaling pathway, cholesterol biosynthetic process. For GO cellular components, DEGs involved in the extracellular matrix, perinuclear region of cytoplasm, integral component of plasma membrane and integrin complex were most strongly represented. In the category of GO molecular function, DEGs were mainly associated with metallocarboxy peptidase activity, calcium ion binding, creatine kinase activity
and protein kinase activator activity. KEGG pathway analysis showed that PPAR signaling pathway, insulin signaling pathway, hippo signaling pathway, bile secretion seemed to be the obvious core metabolic pathways. Furthermore, we screened out twelve DEGs involved in glucose and lipid metabolism, including
PCK1,
FDPS,
HMGCR,
HMGCS1,
HMGCS2,
PPP1R3B,
PPP1R3C,
ACAT2,
SIK1,
OGDHL,
SOCS2 and
IGFBP1 (Table 2). The qRT-PCR analysis revealed that mRNA expression levels of these genes were consistent with the results of RNA-Seq (
P<0.05, Fig 6).
PPAR are ligand-activated transcriptional modulators with well-documented functions in hepatic, whole-body energy homeostasis, lipid and glucose metabolism
and inflammatory responses (
Berthier et al., 2021). A study by
Hu et al., (2021) highlighted that hepatic lipid metabolism in the hyperlipidemia model rats could be regulated by atorvastatin ester through the PPAR signaling pathway and
HMGCR expression. In the present study, we observed the OVX females had reduced expression levels of
PCK1 and
HMGCS2, whereas increased levels of
HMGCR and
HMGCS1 compared with the Sham females. Accumulating evidence indicates that
PCK1, the rate-limiting enzyme in gluconeogenesis, regulates not only glucose homeostasis but also lipogenesis by activating sterol-regulatory element-binding proteins (
Ye et al., 2023). A recent study showed that
PCK1 deficiency stimulated lipogenic gene expression and lipid synthesis and activated the RhoA/PI3K/AKT pathway by increasing intracellular GTP levels, increasing secretion of platelet-derived growth factor-AA
and promoting hepatic stellate cell activation (
Ye et al., 2023).
The steady state of the cholesterol metabolism depends on a complex network involving cholesterol uptake, synthesis, transport and excretion (
Musso et al., 2013). Many clinical and basic research have confirmed that lipid metabolism is bound up with cholesterol synthesis (
Luo et al., 2023). Both
HMGCS and
HMGCR are the two key rate-limiting enzymes in cholesterol synthesis (
Hu et al., 2021). A study conducted on OVX mice model suggested that the expression level of
HMGCR related with cholesterol metabolism also significantly increased in the liver of ovariectomized mice compared with sham mice (
Lei et al., 2021). Besides, several noncoding RNA has been found to regulating hepatic cholesterol and lipid metabolism by targeting
HMGCR (
Sun et al., 2015). Additionally, in cholesterol metabolism, ACAT2 promotes the secretion of cholesteryl ester-enriched very low-density lipoproteins by the liver and lacking hepatic ACAT2 can prevent dietary cholesterol-driven hepatic steatosis in mice (
Alger et al., 2010). Hence, OVX females had elevated
ACAT2 expression may be related to increase of hapetic TC levels.
In the liver, two major subunits, protein phosphatase 1 regulatory subunit 3B(
PPP1R3B) and regulatory subunit 3C (
PPP1R3C) are expressed at approximately equivalent levels and together facilitate the storage of hepatic glycogen. It has been confirmed that overexpression of
PPP1R3B lead to increase of liver glycogen in human cell lines and in mice (
Agius, 2015). Furthermore,
López-Soldado et al., (2015) found that
PPP1R3C transgenic mice increased hepatic glycogen and improved glucose and insulin sensitivity. A survey proved that hepatic
PPP1R3C mRNA expression in diabetic mice was increased compared with that in control mice, suggesting that the abnormally high expression of
PPP1R3C in livers may be one of the main causes of hyperglycemia in diabetes (
Ji et al., 2019). Therefore, the expression levels of
PPP1R3B and
PPP1R3C in the livers of OVX females were upregulated in this study, which may imply ovariectomy affects liver glucose metabolism.
Salt-inducible kinase 1 (
SIK1) belongs to the serine-threonine kinase family
and also is an AMPK-related protein kinase (
Darling and Cohen, 2021). A recent study reported that the overexpression of
SIK1 in the liver of the rats could downregulate the expression of gluconeogenic genes (
Song et al., 2019). In primary mouse hepatocytes, forced
SIK1 expression not only promoted gluconeogenesis but also suppressed lipogenesis (
Wang et al., 2020). Consistent with this, data from
Zhang et al., (2017) demonstrate that
SIK1 expression is inversely correlated with the expression of lipogenic molecules and lipid biosynthesis. In the present study, downregulation of
SIK1 in the OVX livers may facilitate gluconeogenic process.
Oxoglutarate dehydrogenase-like (
OGDHL) participates in regulating the degradation of glucose and glutamate (
Bunik and Degtyarev 2008). A study conducted on hepatocellular carcinoma suggested that silencing of
OGDHL result in decrease of lipogenesis (
Dai et al., 2020), so upregulation of
OGDHL may be associated with elevated triglyceride levels.
Insulin-like growth factor binding protein 1 (
IGFBP1), a binding protein with a high affinity to insulin-like growth factors, mainly produced by hepatocytes and is proved to be associated with glucose regulation and insulin resistance.
IGFBP1 treatment significantly ameliorated hepatic steatosis by interacting with
ITGB1 and suppressed inflammation by inhibiting NF-κB and ERK signaling pathways (
Pan et al., 2021). Previous research has indicated that reduced levels of serum IGFBP1 are commonly observed in obesity, hyperinsulinemia
and nonalcoholic fatty liver disease (NAFLD) (
Graffigna et al., 2009). Furthermore, as one of the suppressors of cytokine signaling (
SOCS) family, it has been clarified that
SOCS2 participated in hepatic steatosis and NAFLD (
Yuan et al., 2016). Recent research published demonstrated that
SOCS2 plays a role in inhibiting inflammation and apoptosis in macrophages during nonalcoholic steatohepatitis through the NF-κB and inflammasome signaling pathways (
Li et al., 2021). Therefore, decreased expression levels of
IGFBP1 and
SOCS2 in the OVX livers may increase the risk of hepatic steatosis.