Next Article in Journal
Hydroxylation of Progesterone and Its Derivatives by the Entomopathogenic Strain Isaria farinosa KCh KW1.1
Next Article in Special Issue
Revisited Metabolic Control and Reprogramming Cancers by Means of the Warburg Effect in Tumor Cells
Previous Article in Journal
KSK-74: Dual Histamine H3 and Sigma-2 Receptor Ligand with Anti-Obesity Potential
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 13
10.3390/ijms23137014
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lack of TRPV1 Channel Modulates Mouse Gene Expression and Liver Proteome with Glucose Metabolism Changes

by
José Thalles Lacerda
1,
Patrícia R. L. Gomes
2,
Giovanna Zanetti
1,
Nathana Mezzalira
1,
Otoniel G. Lima
1,
Leonardo V. M. de Assis
1,3,
Ali Guler
4,
Ana Maria Castrucci
1,4 and
Maria Nathália Moraes
1,2,*,†
1
Laboratory of Comparative Physiology of Pigmentation, Department of Physiology, Institute of Biosciences, University of São Paulo, São Paulo 05508-090, Brazil
2
Laboratory of Neurobiology, Department of Physiology and Biophysics, Institute of Biomedical Sciences, University of São Paulo, São Paulo 05508-000, Brazil
3
Institute of Neurobiology, University of Lübeck, 23562 Lübeck, Germany
4
Department of Biology, University of Virginia, Charlottesville, VA 22904, USA
*
Author to whom correspondence should be addressed.
†
Permanent address: Departamento de Ciências Biológicas, Instituto de Ciências Ambientais, Químicas e Farmacêuticas, Universidade Federal de São Paulo, Diadema 09972-270, Brazil.
Int. J. Mol. Sci. 2022, 23(13), 7014; https://doi.org/10.3390/ijms23137014
Submission received: 17 May 2022 / Revised: 18 June 2022 / Accepted: 20 June 2022 / Published: 24 June 2022
(This article belongs to the Special Issue Molecular Research on Carbohydrate Metabolism)
Browse Figures
Versions Notes

Abstract

:
To investigate the role of the transient receptor potential channel vanilloid type 1 (TRPV1) in hepatic glucose metabolism, we analyzed genes related to the clock system and glucose/lipid metabolism and performed glycogen measurements at ZT8 and ZT20 in the liver of C57Bl/6J (WT) and Trpv1 KO mice. To identify molecular clues associated with metabolic changes, we performed proteomics analysis at ZT8. Liver from Trpv1 KO mice exhibited reduced Per1 expression and increased Pparα, Pparγ, Glut2, G6pc1 (G6pase), Pck1 (Pepck), Akt, and Gsk3b expression at ZT8. Liver from Trpv1 KO mice also showed reduced glycogen storage at ZT8 but not at ZT20 and significant proteomics changes consistent with enhanced glycogenolysis, as well as increased gluconeogenesis and inflammatory features. The network propagation approach evidenced that the TRPV1 channel is an intrinsic component of the glucagon signaling pathway, and its loss seems to be associated with increased gluconeogenesis through PKA signaling. In this sense, the differentially identified kinases and phosphatases in WT and Trpv1 KO liver proteomes show that the PP2A phosphatase complex and PKA may be major players in glycogenolysis in Trpv1 KO mice.

Graphical Abstract

1. Introduction

The transient receptor potential channel vanilloid type 1 (TRPV1) is a member of the TRP channel superfamily, comprised of a large group of nonselective cation channels activated by heat, protons, capsaicin, and others [1]. They are expressed in many tissues and display a wide variety of physiological functions, including metabolic regulation [1,2,3]. TRP channels have been shown to be involved in the control of appetite and body weight [4,5], pancreatic function [6], glucose metabolism [7,8,9], and the development and progression of diabetes mellitus [10]. In this sense, mice lacking TRPV1 channels exhibit impairment in glucose metabolism, especially evidenced by reduced glucose-induced insulin secretion and high blood glucose levels after the glucose tolerance test [8]. On the other hand, TRPV1 activation by capsaicin in C2C12 cells leads to the recruitment of CAMKK2 and AMPK signaling pathways, which results in increased glucose oxidation and ATP production in an insulin-independent way [9] and increased hepatic glycogen storage in streptozotocin-induced diabetic rats [11]. Pharmacological inhibition of TRPV1 results in decreased CAMKK2 and AMPK activation and, thus, reduced glucose oxidation and ATP production [9]. AMPK is an internal energy sensor that responds to decreased nutrient and energy levels by activating calcium signaling to restore energy levels. Interestingly, the absence of TRPV channels blocks glucose-starvation-induced AMPK activation in the liver [12].
In addition to their metabolic function, TRP channels have been reported to be involved in the regulation of clock genes [13,14]. Our group has demonstrated that the clock molecular machinery of brown adipose tissue (BAT) is affected in the absence of TRP channel types V1 [14] and M8 [13]. Interestingly, in ZEM-2S cells of the teleost Danio rerio, we found that heat shock-induced per2 increase is partially dependent on TRPV1 [15].
The liver performs several functions in different phases of the circadian clock, which is mainly modulated by transcriptional programs [16] and post-translational regulation [17]. During physiological fasting (ZT0 to ZT12), liver metabolism is directed to glycogenolysis, gluconeogenesis, amino acid catabolism, and fatty acid oxidation, and during feeding (ZT12 to ZT24), it switches to glycogenesis, glycolysis, amino acid anabolism, and lipogenesis [18]. The mechanism of glucose regulation involves a variety of liver genes that exhibit robust circadian regulation, suggesting that the circadian clock plays a role in hepatic glucose metabolism [19,20,21]. The activities of the rate-limiting enzymes for glycogenesis and glycogenolysis also show circadian rhythms, and the balance between them contributes to the circadian rhythm of hepatic glycogen content. Accordingly, Clock KO [22] or Bmal1 KO mice [18] display reduced mRNA and protein expression of GYS2 (glycogen synthase 2), which is the rate-limiting enzyme of glycogenesis in the liver. Analysis of the hepatic transcriptome and proteome has shown circadian variation in the expression of many genes and proteins related to energy metabolism, including those regulating gluconeogenesis, lipogenesis, mitochondrial biogenesis, and oxidative metabolism, suggesting that liver circadian clocks are involved in these processes [16,20].
Having this in mind, we hypothesized that hepatic glycogen storage could be compromised in the absence of TRPV1 channels and that this metabolic alteration could be mediated by altered signaling pathways involving genes and proteins related to glucose homeostasis. In fact, we demonstrate that the TRPV1 channel contributes to the regulation of clock genes and genes involved in hepatic metabolism. The TRPV1 channel plays a critical role in the modulation of glycogenolysis and gluconeogenesis. Proteomic analyses indicate that the absence of TRPV1 channels may induce inflammatory features related to oxidative stress and increased gluconeogenesis. We also show that the TRPV1 channel is directly associated with glucagon signaling during the fasting period and with changes in kinases and phosphatases involved in glycogen breakdown. Therefore, we demonstrate that the TRPV1 channel is an important regulator of the liver molecular clock and glucose physiological processes.

2. Results and Discussion

2.1. Clock and Metabolic Hepatic Genes Are Altered in the Absence of TRPV1 Channels

The expression of some TRP channels has been demonstrated in the liver [23]; however, their function has yet to be revealed. Hence, we analyzed Trpv1 expression in WT mice at ZT8 and ZT20, and we found a lower expression of Trpv1 at ZT8 compared to ZT20 (p = 0.0159, Figure 1A). Since TRP channels are involved in energy metabolism [1,7] and we have demonstrated that Tprv1 displays a significant difference in its expression at ZT8 and ZT20, one could speculate that TRPV1 channels are involved in the regulation of hepatic functions. To the best of our knowledge, a functional analysis of calcium influx via TRPV1 in a time-dependent manner has yet to be demonstrated.
One of the most evident rhythms observed in almost all animals is the fasting/feeding cycle in association with the sleep/wake cycle [24,25]. This feature enhances energetic efficiency through the temporal segregation of anabolic and catabolic reactions such as gluconeogenesis and glycolysis. An important player that links metabolism regulation and circadian rhythm is a member of the peroxisome proliferator-activated receptor (PPAR) family that functions as a transcription factor [26,27]. In addition to this and having in mind that the TRPV1 channel is involved in glucose homeostasis [8], we started evaluating the daily expression of Pparα and Pparγ in the liver of mice lacking the TRPV1 channel. Statistical analysis performed by two-way ANOVA pointed to a significant difference in the expression of Pparα related to genotype, time, and interaction of these factors (Table 1). We found that only Trpv1 KO mice displayed daily variation in the Pparα transcript (p = 0.0075) when ZT8 was compared to ZT20 (Figure 1B). Increased Pparα expression at ZT8 was observed in the absence of TRPV1 channels in comparison to WT at the same time point (p = 0.0052, Figure 1B). A similar profile was observed in the expression of Pparγ (Figure 1C); F analysis indicated that genotype and interaction of genotype with time are factors that contribute to the difference observed in the groups (Table 1). Expression of Pparγ did not vary between ZT8 and ZT20 in the liver of WT mice, while it did in the Trpv1 KO liver (p = 0.0184, Figure 1C). Even though liver Pparα and Pparγ from WT did not show temporal variation between ZT8 and ZT20, it is well described that these genes display circadian variation, with a peak at the transition between light and dark phases and in the middle of the light phase, respectively [28,29].
PPARs, which act as transcription factors, regulate the expression of a variety of target genes, especially those related to energy metabolism [30]. It has long been known that metabolism and clock genes are tightly intertwined and that the PPAR family contributes to this process [26]. PPARα and PPARγ are direct regulators of the core clock components Bmal1 and Rev-erbα, and, conversely, PPARα is also a direct Bmal1 target gene [26]. In addition, we have previously demonstrated that in the absence of TRPV1 channels, clock genes are downregulated in brown adipose tissue (BAT) [14]. Hence, our next step was to investigate the contribution of TRPV1 channels to the regulation of hepatic clock genes. Analysis of the F test indicated genotype as the only factor that showed statistical significance (Table 1). At ZT8, we found reduced Per1 expression in the absence of TRPV1 channels when compared to WT mice (p = 0.0244, Figure 1D). It is known that clock genes oscillate in the liver; however, the comparison between ZT8 and 20 in both WT and Trpv1 KO mice was not able to provide such information (Figure 1D), which could be due to the time points analyzed here. In the liver of both WT and Trpv1 KO mice, two-way ANOVA indicated significant differences in Bmal1 expression related to genotype and time but not in the interaction of these two factors (Table 1). In this sense, we observed a daily variation of Bmal1 mRNA in both genotypes (Table 1) when ZT8 was compared to ZT20 (WT p = 0.0008 and Trpv1 KO p = 0.0055). Although the F test points to genotype as a significant factor, no difference in Bmal1 expression was seen in the absence of TRPV1 channels when the genotypes were compared at the same time points (Figure 1E). The antiphase relationship between Per1 and Bmal1, a known feature of synchronized molecular machinery, was observed in both genotypes. In summary, the lack of TRPV1 channels led to the reduced expression of clock genes, at least Per1 and Bmal1, but this reduction was not related to the oscillatory expression profile (day–night variation) of these genes. In BAT, the involvement of TRPV1 and TRPM8 channels with temporal control of the local clock machinery is clear since the expression of clock genes was drastically reduced in Trpv1 and Trpm8 KO mice compared to WT mice [13,14].
Liver circadian gene oscillation is evidenced by the study of Zhang and colleagues [31], who demonstrated that among the 19,788 genes analyzed, 16% showed a circadian variation. It is well known that the circadian system regulates, at gene and protein levels, a series of enzymes involved in glucose and lipid metabolism [32]. Temporal oscillation in triglyceride and cholesterol levels in plasma is known to be high during the feeding period, while triglyceride breakdown and free fatty acid mobilization in the adipose tissue are maintained during the fasting period [33,34].
In addition to the involvement of the TRPV1 channel in the regulation of liver clock genes in our experimental setting, its role in metabolic functions has become an attractive research theme. In fact, Trpv1 KO mice are more insulin resistant after consuming a high-fat diet (HFD) in comparison to HFD-fed WT mice [5]. Following this line, the activation of TRPV1 channels by a diet supplemented with capsaicin, a TRPV1 channel agonist, was shown to decrease lipid accumulation and triglyceride levels in the liver. In addition, TRPV1 activation increased the expression of UCP2, which plays a role in hepatic lipid metabolism and mitochondrial bioenergetics [35]. In accordance with these results, Trpv1 KO [5] and Trpv2 KO mice [36] are more susceptible to HFD-induced weight gain than WT mice and show an increased mass of fat deposits.
The maintenance of blood glucose homeostasis represents one of the fundamental challenges in the regulation of body metabolism, and circadian control has an important role in this scenario [37]. In this sense, it is known that glucose levels and tolerance and insulin action are known to vary throughout the day [24]. Considering that clock and PPAR genes are altered in the absence of TRPV1 channels, we decided to investigate how the TRPV1 channel participates in the regulation of hepatic glucose metabolism genes. One of the important functions of the liver is glucose homeostasis, and GLUT2 is the prevalent glucose transporter in this tissue [38]. We did not find temporal differences in the expression of Glut2 in the liver of WT mice, at least at the time points analyzed here. On the other hand, in the absence of TRPV1 channels, an increase in Glut2 expression was observed at ZT8 compared to WT mice. The increased expression of Glut2 in the light phase resulted in a daytime variation in its transcript (Figure 2A). In fact, the F test points out that genotype and the interaction between time and genotype are significant factors for Glut2 transcript regulation (Table 1).
We subsequently evaluated how the lack of TRPV1 would affect liver physiology by assessing the gene expression of G6pase, Akt, and Gsk3β in mice of both genotypes. In agreement with the literature reports on C57BL/6J mice [22,39], G6pase mRNA was found to vary in the liver of WT and Trpv1 KO mice. The highest expression was identified in the dark phase (ZT20) compared to the light phase (ZT8, Table 1 and Figure 2B). No differences were found between genotypes, as demonstrated by the F test (Table 1). Upon insulin interaction with its receptor, the phosphoinositide 3-kinase (PI3K)–AKT pathway is activated, resulting in AKT phosphorylation. Active AKT signals to multiple downstream pathways to control, among other functions, glycogen synthesis through the inhibition of GSK3β [40]. Intriguingly, the expression of both Akt and Gsk3β was increased at ZT8 in the absence of TRPV1 channels (Figure 2C,D). No temporal variation was seen for both genes (Table 1).
Keeping in mind that genes involved in hepatic glucose homeostasis showed increased transcripts in the fasting phase (ZT8), mainly because we found high expression of Glut2 in the liver of Trpv1 KO mice, we questioned whether hepatic glycogen storage would be altered in the absence of the channel. We found a lower glycogen concentration in the liver of Trpv1 KO mice, only at ZT8, compared to WT mice (Figure 2E). Thus, we investigated if low glycogen storage would be related to high glycogenolysis via PEPCK activation by analyzing the Pepck transcript. We found an increase in its expression at ZT8 (Figure 2F), at the same time point when glycogen storage was reduced in the liver of mice lacking the TRPV1 channel. Given that hepatic glycogen storage was reduced and Pepck transcript was increased at ZT8 in the absence of TRPV1 channels, we next evaluated which protein could be altered in the liver of Trpv1 KO mice by analyzing the proteome of both genotypes.
Overall, we demonstrated that the lack of TRPV1 channels leads to the decreased expression of clock genes and the increased expression of metabolic genes at ZT8 compared to WT mice. The time-dependent functions in the liver have been often observed since some components of glucose hepatic production are increased at the light phase [18,20]. For instance, hepatic glucose, glucose-6-phosphate, and oligosaccharides of glycogen breakdown (represented by maltopentaose levels) accumulate during the light phase [20]. These time-specific changes in glucose metabolism are coincident with high glycogenolysis and gluconeogenesis in the liver during fasting. Thus, our data indicate that reduced expression of Per1 and Bmal1 at ZT 8, observed in the absence of TRPV1 channels, could directly influence these two processes. In addition, the reduced glycogen storage at ZT8 in the liver of Trpv1 KO mice reinforces the high hepatic glucose production in the light phase.
Notably, it was reported that Bmal1 and Glut2 have an antiphase relationship, in which at ZT8, Bmal1 expression is down, while Glut2 displays a peak expression at ZT12 [20]. Interestingly, Glut2 has an antiphase relationship with hepatic glucose production since at the end of the light phase (ZT8), when hepatic glucose is high, the expression of Glut2 is increased to enable the export of hepatic glucose to the circulation. Thus, Glut2 reaches its peak expression at ZT12, the moment when the eating phase starts, and then hepatic glucose goes down [20]. Therefore, we highlight that the gene expression changes related to energy metabolism, seen in the absence of TRPV1 channels, may be associated with reduced Bmal1 expression at ZT8 and/or directly triggered by the lack of the channel in the regulation of metabolic processes. Considering that such changes were not observed in the gene expression of Bmal1, Glut2, Pepck, and G6pase at ZT20 in the liver of Trpv1 KO mice compared to WT mice, we hypothesized that the TRPV1 channel plays an important role under physiological fasting conditions.

2.2. Lack of TRPV1 Channel Affects Proteome and Disrupts Biological Functionality in Liver

To explore the mechanisms underlying the fasting metabolic adaptation in the absence of TRPV1 channels, we used a quantitative relative-based mass spectrometry approach (label-free quantification, LFQ) to investigate changes in the liver proteomic signature of Trpv1 KO mice compared to WT animals at ZT8 (Supplementary Information File S1, Tables S1–S3). We found 1920 proteins in both proteomes, of which 234 and 231 proteins were unique to the liver of WT and Trpv1 KO mice, respectively (Figure 3A). We found that 1455 proteins were shared between the proteomes, with 26 proteins downregulated and 17 upregulated in the liver of Trpv1 KO mice in comparison to WT mice (Figure 3B).
Differential proteins encoded in the liver of WT and Trpv1 KO mice could indicate functional differences induced by the absence of TRPV1 channels. In this regard, we highlighted the downregulation of core components of nucleosomes (H2BC9, H2BC4, and H1-0) and lipid metabolism (DHCR24 and DDX3X) in the absence of TRPV1 channels. It has been shown that nucleosomes serve as a regulator of transcription since their interaction with gene promoters prevents the initiation of transcription [41]. Considering that in the absence of TRPV1 channels, reduced nucleosome proteins were found, we speculated whether transcription regulation could have been affected. In fact, we found increased gene expression of the majority of the hepatic metabolic genes analyzed here (Figure 1 and Figure 2).
In addition, DHCR24, a protein involved in an anti-inflammatory and pro-resolving phenotype (accelerating inflammation resolution) [42], was reduced in the absence of TRPV1 channels (Figure 3B), which could be related to a more inflammatory phenotype. RNA helicase DDX3X is involved in a variety of transcript metabolic regulations [43], including lipid metabolism [44]. Downregulation of DDX3X leads to decreased protein and transcript expression of microsomal triglyceride transfer protein (MTP) and, therefore, may enable the blocking of lipoprotein assembly and the inhibition of lipid export, which is a leading cause of lipid accumulation [44]. We found a reduction of this protein in the absence of TRPV1 channels (Figure 3B). In this sense, we found high expression of Pparα and Pparγ (Figure 1B,C) in the liver of Trpv1 KO mice, two important transcription factors involved in lipid metabolism [26].
Some of the most significantly increased proteins (ALDOC, SSBP1, PLIN3, EHD3) (Figure 3B) participate in the same molecular mechanisms as some decreased proteins in processes that converge to the establishment of the liver phenotype in Trpv1 KO mice. For instance, some increased proteins, in addition to decreased DHCR24 and DDX3X, also suggest a phenotype related to fatty acid metabolism changes in the absence of TRPV1 channels. Among them, ALDOC, an aldolase family member that has a central role not just in glycolysis but also in gluconeogenesis pathways, was upregulated 3.4-fold (log2), indicating glucose and glycerolipid metabolism changes in the liver of Trpv1 KO mice [45,46]. PLIN3 is involved in lipid storage, mobilization, and droplet biogenesis, and it has also been associated with hepatic steatosis [47]. This protein also regulates early-endosome-to-Golgi transport [48], similar to EHD3 [49]. In parallel, SSBP1, which was also increased in Trpv1 KO mice (Figure 3B), is involved in the replication and maintenance of mitochondrial DNA, DNA repair, and cell protection from proteotoxic stresses [50].
To assess in depth the hepatic functional dysregulation caused by the lack of TRPV1 channels, we next conducted a gene ontology biological process (GO-BP) enrichment analysis using exclusively identified proteins of each group. We found that 43 differential GO-BPs were evidenced in both groups (Figure 3C). The differential GO-BP of “regulation of alternative mRNA splicing, via spliceosome” in the liver of Trpv1 KO mice indicated enhanced RNA processing, which is in line with the downregulation of H2BC9, H2BC4, and H1-0 (core components of nucleosomes, Figure 3B). Proteins clustered in this enriched GO-BP may regulate a variety of post-transcriptional processes (e.g., alternative splicing, nucleocytoplasmic transport, translation, and degradation) under stress conditions and cellular reprogramming [51,52]. Of these, THRAP3 promotes the binding of CLOCK-BMAL1 to DNA and links it to the basic transcriptional machinery [53]; it is also required for the DNA damage response, control of mRNA splicing, and nuclear mRNA degradation [54]. In addition, we evidenced “cellular response to oxidative stress”, which might be intimately linked to DNA damage and inflammation, as demonstrated by the participation of TXN2 in previous studies [55,56]. In this context, the immune response mediated by immunoglobulin in Trpv1 KO mice seems to indicate inflammatory status [57,58], which is in line with reduced DHCR24 involved in the anti-inflammatory phenotype (Figure 3B).
In contrast to enriched immune responses mediated by immunoglobulin in the liver of Trpv1 KO mice, we observed that “response to interferon-gamma” was enriched only in the WT liver. It has been demonstrated that interferon-gamma deficiency attenuates hepatic inflammation in a steatohepatitis model [59]. Thus, the lack of proteins related to “response to interferon-gamma”, associated with the compensatory hepatoprotective mechanism shown in the liver of Trpv1 KO mice, may be suggestive of a protective effect displayed against inflammation mediated by the absence of TRPV1 channels. In this line, it was demonstrated that decreased HPGD expression in hepatic endothelial cells is related to hepatoprotective effects in response to liver injury [60]. Indeed, we identified the HPGD protein in the liver of WT mice but not in Trpv1 KO mice (Figure 3B), reinforcing our hypothesis that Trpv1 KO mice display an inflammatory status and subsequent hepatoprotective response. Other proteins that participate in the fatty acid metabolism identified in the liver of WT mice (ACNAT1, FADS1, and NADK2) have been described to prevent inflammation and oxidative stress [61,62,63] but were not found in the Trpv1 KO liver (Figure 3C).
Some of the enriched GO-BPs in the liver of WT mice are related to metabolism responses, including cellular sphingolipid homeostasis, the fatty acid metabolic process, and the oxidation-reduction process (Figure 3C), indicating, mainly, metabolic reprogramming of the lipid metabolism in the absence of TRPV1 channels. Among the associated proteins, mitochondrial enzyme malonate-CoA ligase (ACSF3) induces malonyl-CoA production from malonate, a toxic metabolic product from mitochondrial metabolism [64]. Thus, the fact that ACSF3 was not identified in the liver of Trpv1 KO mice could be related to a possible reduction in malonyl-CoA biosynthesis. Importantly, the acetyl-CoA carboxylase 2 (ACACB, also known as ACC2), a mitochondrial enzyme that catalyzes the carboxylation of acetyl-CoA to form malonyl-CoA [65], was reduced 1.85-fold (log2) in Trpv1 KO liver (Figure 3B). In line with these findings, reduced malonyl-CoA availability would lead to increased CPT1 activity, followed by increased fatty acid oxidation (FAO) and reduced lipogenesis [66]. Next, we identified three fatty-acid-binding proteins (FABP1, FABP4, and FABP5) in the WT proteome and four proteoforms (FABP1, FABP2, FABP4, and FABP5) in the Trpv1 KO proteome. These proteoforms bind long-chain FA with high affinity and are involved in lipid uptake and FA transport [67]. Among these, FABP2 is associated with lipid droplet accumulation in HepG2 cells, being a downstream target of PPARγ as well [68]. It is also related to lipid metabolism changes in the liver of type 2 diabetes mellitus in a nonalcoholic fatty mice model [69].
Hepatic lipid accumulation occurs when the uptake of circulating FA and de novo lipogenesis surpasses the FAO and export of lipids in very-low-density lipoproteins (VLDLs) [70]. Based on our findings in the liver proteomics analysis of Trpv1 KO mice compared with WT mice, as demonstrated by DDX3X (less lipid export), ACACB and ACSF3 (less malonyl-CoA and subsequently higher FAO), FABP2 (high FA uptake) and PLIN3 (lipid availability), we suggest that the liver of Trpv1 KO mice could have greater FA availability to beta-oxidation to support a boost in gluconeogenesis, thus preventing lipid accumulation. In fact, hepatic triglyceride levels in WT did not differ from Trpv1 KO mice [35].
Interestingly, proteins involved in glucose metabolism in Trpv1 KO and WT mice function in opposite directions. UDP galactose-4-epimerase (GALE), an enzyme that participates in the synthesis of glycoproteins and glycolipids, has been found in the liver of Trpv1 KO mice. The hepatic overexpression of GALE has been reported to enhance gluconeogenesis from pyruvate [71]. On the other hand, glucosamine-6-phosphate isomerase 1 (GNPDA1) has been found only in the liver of WT mice. This enzyme provides a source of energy in the form of phosphosugars derived from the catabolism of glycoproteins and glycolipids and transfers them into the glycolytic pathway [72].
Taken together, the above data seem to indicate that the absence of TRPV1 channels alters gene expression and RNA processing and induces inflammatory features, oxidative stress, and an increase in FAO to support more substrates for gluconeogenesis, resulting in changes in glucose metabolism.

2.3. Protein–Protein Interaction Network Correlates the TRPV1 Channel and Hepatic Glucose Metabolism

We constructed a protein–protein interaction network (PPIN) based on a physical interaction complex of all proteins identified in each proteome and by aggregating them to the TRPV1 protein, although this key protein was not found in our analysis. Further, we removed nodes (proteins) disconnected from the rest of the network. The PPINs assembled through WT and Trpv1 KO proteomes, with 1389 proteins/16,394 interactions and 1415 proteins/16,633 interactions, respectively, were used to build a TRPV1-driven subnetwork based on the oriented network method (diffusion) to explore connectivity patterns directly linked to it. The resulting oriented subnetworks had 141 proteins (~10%) from PPINs, of which 100 were overlapped among them. Each cluster of protein interaction was submitted to enrichment analysis using the GO biological process and KEGG pathways, where most of the enriched GO terms and signaling pathways share subnetworks (Figure 4A).
Using this approach, we identified that glucagon signaling was enriched in cluster 4 of the PPIN from the liver of WT mice (Figure 4A). This was an important and unexpected finding since glucagon signaling in Trpv1 KO mice was not identified (Figure 4B) in a similar way as described for WT mice. Cluster 4 from the WT network included “glucagon signaling pathway” (PYGL, PYGM, CALM3, PHKB, PLCB1) and “phosphatidylinositol signaling system” (CALM3, PLCB1, SACM1L, INPP1, IMPA1). Calmodulin-3 (CALM3) is a key protein in both clusters and is connected to TRPV1; however, connections with proteins exclusively identified in the WT proteome (PLCB1 and PHKB) enabled glucagon pathway enrichment. In a downstream pathway, PLCB1 interacts with INPP1, SACM1L, and IMPA1 in phosphatidylinositol signaling. Thus, the relationship shown here proposes that TRPV1 absence leads to PLCB1 signaling impairment and, thus, directly influences associated pathways such as glucose metabolism. Indeed, the relationship between TRPV1 and PLCB1 has been reported in several studies [73,74,75,76].
Given that the glucagon signaling pathway (KEGG ID: mmu04922) is directly related to the TRPV1 channel, we focused on further analysis of this pathway to maximize insight into the effect of TRPV1 absence on glycogenolysis/gluconeogenesis signaling pathways. We found that of the 25 genes mapped for this database, 21 proteins from WT and Trpv1 KO were overlapping. The proteins phospholipase C beta 1 (PLCB1), phosphorylase kinase subunit beta (PHKB), and 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 1 (PFKFB1, known as PFK/FBPase-1) were exclusively identified in the WT proteome, while fructose-1,6-bisphosphatase isozyme 2 (FBP2) was uniquely identified in the Trpv1 KO proteome.
Based on the direct visualization of this pathway and the resulting outputs to glucagon stimulus reported by the KEGG database, we reassembled this signaling pathway using shared proteins among groups and exclusively identified proteins in the liver of WT mice for a better understanding of glycogenolysis and gluconeogenesis in the absence of TRPV1 channels. It was observed that cyclic AMP-dependent protein kinase (PKA) displays a key role in hepatic glucose regulation, phosphorylating PFKFB1 and PHKB proteins (Figure 4C). In response to glucagon, PKA phosphorylates Ser32 of PFKFB1 in the liver, which, in turn, leads to PFK-2 inactivation and FBPase-2 activation to reduce D-fructose 2,6-bisphosphate (F2,6-P2), thus increasing gluconeogenesis [77]. As F2,6-P2 is a powerful allosteric activator of 6-phosphofructo-1-kinase (PFK-1) to stimulate glycolysis and an inhibitor of fructose-1,6-bisphosphatase (FBP), a regulatory enzyme of gluconeogenesis, glucagon decreases the concentration of hepatic F2,6-P2, allowing gluconeogenesis activation [78].
On the other hand, we also investigated the role of the FBP-2 protein and upstream kinases that indirectly regulate its activity in the signaling pathway of Trpv1 KO liver. In this signaling pathway, PKA and AMP-activated protein kinase (AMPK) represent the main regulatory proteins of PFKFB1 and, subsequently, FBP-2 (Figure 4D). Unlike the PKA-mediated mechanism in the liver, activation of AMPK leads to the phosphorylation and activation of PFK-2 (heart isoform), increasing F2,6-P2 and stimulating glycolysis [79]. In addition, AMPK activation has antagonistic mechanisms to PKA activity in the liver. An increase in the intracellular AMP level not only increases AMPK activity but also inhibits adenylyl cyclase, reducing the ability of PKA to promote gluconeogenesis [80]. Interestingly, liver AMPK activity is impaired upon glucose starvation of Trpv1 KO mice [12], which suggests decreased F2,6-P2 levels in this genotype. In addition, the relationship between decreased AMPK activity and the downregulation of ACC2 in the liver of the Trpv1 KO mice shown here can be found in the glucagon signaling pathway from the KEGG database.
Since FBP-2 protein, a regulatory enzyme of gluconeogenesis that converts F1,6-P2 to F6P for hepatic glucose production, was positively regulated in the liver of Trpv1 KO mice, it is suggestive that a high gluconeogenesis rate is found in the absence of TRPV1 channels. Increased gluconeogenesis in Trpv1 KO mice might also be supported by the fact that the lack of TRPV1 channels would lead to deregulated Ca2+ signaling. Ca2+ negatively regulates FBPase [81], and it may inhibit adenylyl cyclases 5 and 6 [82]. In this sense, TRPV1 channel absence and subsequently disrupted Ca2+ signaling could contribute to higher FBPase activation through increased PKA-dependent phosphorylation. In fact, constitutive activation of PKA in the liver results in fasting hyperglycemia and evidence of increased glycogenolysis and gluconeogenesis [83].

2.4. Differential Interactions of Kinases and Phosphatases Underline Changes in Hepatic Glycogen Metabolism in Trpv1 KO Mice

We investigated the participation of exclusively identified kinases and phosphatases in each proteome in the modulation of glycogenolysis/glycogenesis. Using gene ontology annotation from the Uniprot database, we mapped MAPK14, SGK2, and TRPM7 kinases and PGP, PPM1B, LHPP, and DUSP3 phosphatases in Trpv1 KO mice, while in the liver of WT mice, TTN, TWF1, CPNE3, CSNK2A1, MAPK1, TEK, ALPK3, PHKB, and ATM kinases and phosphatase PPP1CB were identified. We also examined proteins that improve or impair the kinase/phosphatase activity (Figure 5A). After that, we used the PhosphoSitePlus database (https://www.phosphosite.org/, accessed on 19 June 2022) with a setting to mouse, KEGG pathway annotation, and literature data to map the relationship between kinases/phosphatases and their upstream regulatory proteins and downstream targets and glycogenolysis/glycogenesis. Using these strategies, 11 proteins related to glycogenolysis/glycogenesis pathways were mapped (Table 2).
For a better understanding of the functionality of these mapped proteins in the glycogenolysis pathways, we assembly a PPIN, remaking connections among the mapped proteins that break down glycogen and are expected to have high activity at ZT8. To avoid a gap of connection in the PPIN, we then inserted upstream proteins of hepatic glycogen (i.e., GYS2 and PYGL) that were identified in both proteomes and proteins not found in the proteome data (PKA, AMPK, GSK3B, and FOXO) that are associated with downstream signaling of the mapped proteins (Figure 5B). Interaction of some kinases with upstream phosphatases was observed, which seems to be required to attenuate kinase activity by negative feedback (e.g., DUSP3 andMAPK14 in Trpv1 KO mice and PPP1CB and PHKB in WT mice). Moreover, our results show that the dephosphorylation and phosphorylation of target proteins are a response to the increased activity of some positive regulatory proteins and, therefore, trigger an important role in downstream signaling.
Taken together, the results obtained with differentially identified kinases, phosphatases, and regulatory proteins in Trpv1 KO and WT mice indicate that changes observed in the Trpv1 KO mice orchestrate an increased breakdown of glycogen. In this context, PKA activity is positively regulated by PRRC1 [84] and ARFGEF2 (also known as BIG2), which acts as an A kinase-anchoring protein (AKAP) [85]. Then, PKA activates protein phosphatase 2A (PP2A complex) [86]. The PP2A phosphatase activator (PTPA) also stimulates PP2A activity in the presence of Mg2+ and Zn2+ [87]. Transient receptor potential cation channel subfamily M member 7 (TRPM7) promotes Zn2+ release from intracellular stores to regulate ROS signaling in oxidative stress responses [88], as indicated by functional enrichment analysis in Trpv1 KO mice. Intriguingly, ENSA, a PP2A inhibitor [89], was also found in the liver of Trpv1 KO mice (Figure 5A), and it seems to act as an “on–off” switch of the phosphatase complex. Once PP2A is active, it mediates the dephosphorylation and inhibition of AMPK and the dephosphorylation and activation of the FoxO1 transcription factor [90]. Among FoxO signaling targeted genes, we observed an increase in the expression of G6pc1 and Pck1 (also known as Pepck) at the transcription level in the absence of TRPV1 channels (Figure 2), which nevertheless did not show significant differences in our proteomic data.
Furthermore, PP2A has also been described as a GSK3B activator [91]. However, it is important to mention that multiple players may regulate GSK3B, as shown in Figure 5B. Thus, the role of the TRPV1 channel in GSK3B phosphorylation is yet unclear due to the many upstream kinases and phosphatases involved in this signaling.
Contrary to Trpv1 KO mice, PPP1CB phosphatase and PHKB kinase were identified in the liver of WT mice, in which PPP1CB acts upstream of PHKB to inhibit it and subsequently promotes the inhibition of PYGL (Figure 5B). PPP1CB is required to regulate glycogen synthesis because it not only inhibits PYGL but also activates GYS2 [92]. We speculate that the PPP1CB may participate in the attenuation of breakdown glycogen magnitude and enable glucose homeostasis at ZT8.
Zhong and collaborators [8] demonstrated that the blood glucose of Trpv1 KO mice has no significant differences in basal conditions when compared to WT mice; however, there was a significantly higher blood glucose rate in Trpv1 KO mice after glucose stimulus. We believe that the initial difference in blood glucose level could not be observed because the daily switch of glycogen metabolism is evidenced by elevated glycogen breakdown, with peaks of glycolytic/gluconeogenic metabolites at ZT8 [18]. Although multiple layers of regulation are required for glycogen metabolism, our results strongly suggest that the lack of TRPV1 channels promotes the increase of glycogenolysis, gluconeogenesis, and FOXO signaling through the PP2A–PKA axis and, thereby, enhances hepatic glucose production.

3. Material and Methods

3.1. In Vivo Procedures

In vivo experimental procedures were conducted on a total of 20 male mice ranging from 3 to 4 months old: 10 C57BL/6J wild-type (WT) and 10 Trpv1 KO in C57BL/6J background mice provided by the Department of Physiology, University of Sao Paulo (USP), originally obtained from Jackson Laboratory (B6.129X1-Trpv1tm1Jul/J, Stock No.: 003770). The lack of Trpv1 channels was confirmed by quantitative PCR (data not shown). Experimental procedures were performed according to the protocol approved by the Ethics Committee for Animal Use (CEUA IB/USP, number 255/16). All methods were carried out in accordance with relevant guidelines and regulations.
Mice were singularly housed in standard propylene cages with access to food and water ad libitum, kept at 22 °C ± 1 under a 12:12 light/dark cycle (800–1000 lux white LED light, ranging from 420 to 750 nm, lights on at 7 AM (ZT 0) and off at 7 PM (ZT 12)). After a 7-day acclimation, the animals were kept in the same conditions for additional 14 days. On the 15th day, the mice were euthanized with CO2 at ZTs 8 and 20, and death was assured by cervical dislocation. At ZT 20, all experimental procedures were carried out with the assistance of night goggles in complete darkness. The liver from 5 WT and 5 Trpv1 KO mice at each time point (ZT8 and ZT20) was removed for RNA extraction, quantitative PCR (qPCR), and glycogen quantification. Proteomics analyses were performed using 4 biological replicates of each group at ZT8, as described below. The time points of euthanasia were chosen based on the fact that peaks of glycolytic/gluconeogenic metabolites and maximal BMAL1-DNA binding occur at ZT8 [18]. To compare time-dependent changes between light and dark phases, we also analyzed gene expression and glycogen storage 12 h later, which corresponds to ZT20.

3.2. Total RNA Extraction and Reverse Transcriptase–Polymerase Chain Reaction (RT–PCR)

Small liver fragments were homogenized in TRIzol (Thermo Fisher Scientific, Waltham, MA, USA), and total RNA was extracted and purified according to the manufacturer’s instructions. Total RNA concentration (OD260) was determined using a spectrophotometer (Nanodrop, Wilmington, DE, USA), and 1 μg was subject to reverse transcription with Superscript III Reverse Transcriptase, random hexamer primers, and other reagents according to the manufacturer’s instructions (Thermo Fisher Scientific, Waltham, MA, USA), generating complementary DNA (cDNA).

3.3. Quantitative PCR (qPCR)

cDNA was subject to quantitative PCR reactions, as previously described [13,14], using a species-specific pair of primers and a probe (Table 3) based on sequences obtained from GenBank (http://www.ncbi.nlm.nih.gov/genbank, accessed on 15 November 2019). Primers and probes were designed as spanning introns by the Primer Express program (Life Technologies, Carlsbad, CA, USA), and synthesized by IDT (Coralville, IA, USA); 18S ribosomal RNA was used to normalize the values of the genes of interest. Prior to this selection, we ascertained that 18S rRNA did not vary among time points (data not shown).
For qPCR using the Taqman multiplex protocol for simultaneous analysis of Per1, Bmal1, and 18S rRNA, the solutions contained primers (300 nM for the genes of interest and 50 nM for 18S rRNA) and probes (200 nM for the genes of interest and 50 nM for 18S rRNA) (Table 3), PowerMix 2X (Bio-Rad Laboratories, Hercules, CA, USA), and cDNA in triplicate for each sample. Individual assays for Trpv1, Akt, G6pase, Glut2, Gsk3β, Pparα, Pparγ, Pepck, or 18S rRNA were prepared with primers at 300 nM (Table 3), KAPA SYBR FAST 2x (Kapa Biosystems, Wilmington, MA, USA), and cDNA in duplicates for each sample. Non-template wells were routinely included.
Reactions were carried out in an i5 thermocycler (Bio-Rad Laboratories, Hercules, CA, USA) in the following conditions: for multiplex assays, 3 min at 95 °C, followed by 45 cycles of 15 s at 95 °C and 60 s at 60 °C; for SYBR Green assays, 10 min at 95 °C, followed by 45 cycles of 15 s at 95 °C, 1 min at 60 °C, and 80 cycles of 10 s, starting at 55 °C with a gradual increase of 0.5 °C.

3.4. Gene Expression Data Analysis

Gene expression assays were quantified according to the 2−ΔΔCT method [93]. All values were normalized by 18S ribosomal RNA, and the resulting minimal mean value obtained from WT mice was subtracted from all other values for both WT and Trpv1 KO mice, obtaining ΔΔCT. Then, ΔΔCT was used as a negative exponential of base 2 (2−ΔΔCT), providing the relative expression of the gene between genotypes and between time points within each genotype.
The log values were obtained from 5 animals, and the median, quantiles, and maximum and minimum expression values of gene transcripts normalized by 18S ribosomal RNA and expressed relative to the minimum mean value of WT mice. were plotted in graphs. Student’s t-test was used to compare Trpv1 gene expression between ZT 8 and ZT 20 in WT mice. To determine the differences between time points of the same genotype or among genotypes at the same time point, logarithmic data were compared by two-way ANOVA, followed by the Bonferroni post-test; p < 0.05 was established to reject the null hypothesis. All statistical analyses were carried out in Prism 7.0 (GraphPad 7.0, San Diego, CA, USA).

3.5. Glycogen Quantification

To measure the amount of liver glycogen, 250 mg of liver was processed as described elsewhere [94]. The absorbance was measured at a wavelength of 650 nm in duplicate for each sample in a plate reader (Biotek Sinergy H1, Winooski, VT, USA). Hepatic glycogen quantification was calculated using the absorbance of the samples in relation to tissue mass (milligrams) and normalized by the glucose concentration curve.

3.6. Sample Preparation for Mass Spectrometry (MS) Analysis

Frozen livers of Trpv1 KO and WT mice (n = 4 per group) euthanized at ZT8 were mechanically homogenized in 100 mM Tris-HCl, pH 7.8, 3% SDS, and 150 mM NaCl supplemented with protease and phosphatase inhibitors (Roche, Indianapolis, IN, USA) using a bead blaster refrigerated homogenizer (Benchmark Scientific, Sayreville, NJ, USA) at a maximum speed of 3650 rpm with three cycles of 30 s on-and-off at 4 °C. Samples were then sonicated for three cycles (20 s bursts with 20 s pauses) and centrifuged at 16,000× g for 10 min at 4 °C. The supernatants were collected, and the protein concentration was measured using the Pierce BCA protein assay in technical duplicates for each sample (Thermo Scientific, Waltham, MA, USA). Aliquots of SDS-lysates containing 300 µg of total protein from each sample were processed according to the filter-aided sample preparation (FASP) method using Microcon 10 kDa centrifugal filter units (Merck Millipore, Burlington, MA, USA) at 10,000× g for 50 min at 20 °C [95,96,97]. To remove SDS, filters containing proteins were washed with urea buffer. Next, Trypsin/LysC Mix (Promega, Madison, WI, USA) was added to the filters at an enzyme-to-protein ratio of 1:100 (w/w) and incubated for 16 h at 37 °C in a thermo-mixer at 600 rpm. Following protein digestion, peptides were filtered through the membrane using 3% acetonitrile in 0.1% trifluoroacetic acid (TFA) and purified with reversed-phase chromatography using C18 micro-pipette tips (TopTipTM, PolyLC, Columbia, MD, USA) according to the manufacturer’s instructions. Peptides were dried in a vacuum concentrator and stored at −20 °C.

3.7. LC–ESI–MS/MS Analysis

Dried peptides were recovered in 20 µL of aqueous 0.1% formic acid and analyzed by online nanoflow LC–ESI–MS/MS using an ACQUITY UPLC system (Waters, Milford, MA, USA) coupled online to a maXis 3G quadrupole time-of-flight (Q-TOF) mass spectrometer (Bruker Daltonics, Bremen, Germany), equipped with a Captive Spray nanoelectrospray source. One µL of the sample was injected, and the peptide mixture was loaded onto a trap column (nanoAcquity UPLC® 2G-V/M Trap 5 µm Symmetry® C18, 180 µm × 20 mm, Waters, Milford, MA, USA) for 3 min at a flow rate of 7 µL/min of aqueous 0.1% formic acid. Peptides were separated on an analytical column (nanoAcquity UPLC® 1.8 µm HSS T3, 75 µm × 200 mm, Waters, Milford, MA, USA) using a 240 min gradient from 2% to 85% of 0.1% formic acid in acetonitrile (1 min at 2%, 209 min 2–30%, 10 min 30–85%, 5 min wash at 85%, 5 min 85–2%, 10 min equilibration at 2%). The flow rate was set at 200 nL/min. Eluted peptides were analyzed in a maXis 3G (Q-TOF) mass spectrometer operated in the positive mode under a data-dependent acquisition manner in the m/z range of 150–2200. Precursor ions were fragmented with collision-induced dissociation (CID).

3.8. MS Data Analysis and Statistical Analysis

The raw files (.d) were loaded into PEAKS Studio 8.5 software (Bioinformatics Solution Inc., Waterloo, ON, Canada), and the PEAKS standard workflow (de novo peptide sequencing, PeaksDB, PeaksPTM, and SPIDER) was applied to the identification of proteins [98,99]. MS/MS spectra were searched against the Mus musculus reference proteome database available at UniprotKB (Proteome ID UP000000589—21,986 protein entries, download one protein sequence per gene, release date Dec 2021). The following parameters were used: precursor mass tolerance of 25 ppm; fragment mass tolerance at 0.02 Da; trypsin was set as the specific enzyme, and up to two missed cleavages were required; carbamidomethylation on cysteine (+57.02 Da) as a fixed modification, oxidation on methionine (+15.99 Da), and acetylation on protein N-terminus (+42.01 Da) were set as variable modifications, with maximal 3 modifications per peptide in SPIDER outcomes. Significance score of −10lgP > 20 (p-value <0.01) for proteins and peptides was applied, and at least 1 unique peptide was required for protein identifications. Proteins were considered to be validated in an experimental group when they were identified in at least two out of the four biological replicates of a given group. The PEAKS Q module was applied (label-free quantification method) to the database search outcome, and the normalization factors based on total-ion count (TIC) were obtained. The normalization factor obtained for each sample was employed to calculate the normalized areas of proteins. Normalized protein areas have been log2-transformed, and the values were subject to Student’s t-test, in which p-values < 0.05, and fold changes ≥1.2 (log2 (FC) ≤ −0.26 or ≥0.26) were selected as cut-off points of significance to the differential relative abundance of proteins.

3.9. MS Data Integration, Network Propagation, and Functional Analysis

To identify biological processes represented by the differentially identified proteins in each proteome, we carried out a functional enrichment analysis of gene ontology biological processes (GO-BPs) using the Database for Annotation, Visualization, and Integrated Discovery (DAVID v.6.8) functional annotation tool [100].
All proteins identified in each proteome were considered for the generation of protein–protein interaction networks (PPINs) of WT and Trpv1 KO mice. UniProt IDs from the protein list were used to perform a multiple protein search in STRING DB v.11.5 (https://string-db.org/, accessed on 19 June 2022) [101]. The chosen protein to set this analysis was based on those that were indicated as part of a physical interaction complex and those that presented a minimum required interaction score of 0.4 (medium confidence score). The resulting interaction network was imported into Cytoscape v.3.8.2 (Boston, MA, USA) [102] for visualization. Then, we conducted a supervised network propagation considering protein–protein interactions and hierarchical relationships using Cytoscape plug-in diffusion [103]. Network propagation is an important method to identify and highlight biological signals within a dataset since it can amplify the association between proteins that lie in network proximity and enable a reduction in the number of paths in the network to guide subsequent analysis focused on targets [104,105]. As the goal of this analysis is to find a subnetwork that is related to the TRPV1 channel and may affect glycogen metabolism, we chose the TRPV1 protein for driving the underlying network to it, using a cut-off point of the top 10% of the most relevant proteins.
Clusters in the subnetwork were identified using the GLay Community Clustering algorithm (default options) incorporated in the Cytoscape plug-in clusterMaker v.2 (Boston, MA, USA) to provide an optimized layout and structured visualization for biological analysis [106,107]. Each cluster was submitted to enrichment analysis in DAVID v.6.8 using the UniProt IDs of clustered proteins, and then the outputs of gene ontology biological process (GO-BP) and annotation protein pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were imported to the Cytoscape plug-in EnrichmentMap (FDR < 0.01) [108].

4. Conclusions

The TRPV1 channel has become a potentially relevant target for metabolic intervention, including insulin resistance, obesity, type 2 diabetes, and nonalcoholic fatty liver, due to its presence in a variety of metabolic tissues. Even though recent studies have greatly contributed to the knowledge regarding the participation of TRPV1 in metabolism [5,8,9,11,12], there are few reports in the literature, thus far, describing the involvement of TRP channels in the mechanisms responsible for hepatic glucose regulation. A phenotype of hyperglycemia is described in mice lacking TRPV1 channels; however, the mechanism underlying this event remains an open question. Our study contributes to addressing this gap in the literature as we have established which molecular components are involved with reduced glycogen storage in Trpv1 KO mice. In the present study, we demonstrate that the lack of TRPV1 channels results in reduced glycogen storage, which is related to increased glycogenolysis but also induces high gluconeogenesis. These alterations were associated with molecular signature changes found in multilayer analysis, from transcriptional modulation of the clock system and metabolic genes to altered proteomes in the absence of TRPV1 channels. We suggest that the increased FAO observed in the liver of Trpv1 KO mice is required to support high hepatic gluconeogenesis since lipid export and malonyl-CoA seem to be reduced. In addition, we found that Per1 and Bmal1 genes, two important components of the molecular machinery responsible for circadian regulation of hepatic glucose metabolism, were altered in the liver of Trpv1 KO mice. Taken together, we demonstrate that pharmacological intervention of signaling pathways, in which the TRPV1 channel is a pivotal component, may modulate glucose metabolism and, therefore, improve liver function. Thus, our results open new avenues to investigate the role of the TRPV1 channel as a key component of hepatic glucose metabolism.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23137014/s1.

Author Contributions

Conceptualization, J.T.L., A.M.C. and M.N.M.; methodology and data acquisition, P.R.L.G., N.M., G.Z., O.G.L. and M.N.M.; formal analysis, J.T.L., P.R.L.G. and M.N.M.; proteomics data acquisition, J.T.L. and M.N.M.; proteomics analysis, J.T.L., with the help of M.N.M.; critical data analysis, all authors; original draft preparation and writing, J.T.L. and M.N.M.; writing, critical review and editing, P.R.L.G., L.V.M.d.A., A.G. and A.M.C. All authors have approved the definitive version of the manuscript and agreed to be accountable for all aspects of the study in ensuring that questions related to the accuracy or integrity of any part of the study are appropriately investigated and resolved. All persons designated as authors qualify for authorship, and all those who qualify for authorship are listed. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the Sao Paulo Research Foundation (FAPESP, grants 2018/14728-0, 2017/24615-5, 2012/50214-4 (A.M.C.), and 2017/26651-9 (M.N.M.)), the National Council of Technological and Scientific Development (CNPq grants 301293/2011-2, 303070/2015-3 (A.M.C.), and 428754/2018-0 (M.N.M.)), and CAPES (Coordination for the Improvement of Higher Education Personnel, Finance Code 001). MN Moraes is a Young Investigator of FAPESP (2019/19005-9 and 2014/16412-9); J.T. Lacerda, P.R.L. Gomes, O.G. Lima, and G. Zanetti are fellows of FAPESP (2020/04524-8, 2021/09923-0, 2018/23915-8, and 2021/01659-2). L.V.M. de Assis and N. Mezzalira were fellows of FAPESP (2018/16511-8 and 2013/02131-5).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Institute of Biosciences, University of São Paulo (protocol CEUA IB/USP, number 255/2016, 14 June 2016).

Data Availability Statement

The raw MS data associated with this manuscript have been submitted to a public repository (the Mass Spectrometry Interactive Virtual Environment–MassIVE, http://massive.ucsd.edu, accessed on 19 June 2022) and deposited to the ProteomeXchange Consortium (http://www.proteomexchange.org/, accessed on 19 June 2022). Both accessed on 21 June 2022. These data are associated with the identifier MassIVE ID MSV000089696 and Proteome Exchange ID PXD034801. Further information on research design is linked to this article. The authors declare that all other data supporting the findings of this study are available within the article and its supplementary information files.

Acknowledgments

We are thankful to Renata Alves dos Santos and Leticia Menezes Vasconcelos (University of São Paulo, São Paulo, Brazil) for providing excellent mouse care.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Christie, S.; Wittert, G.A.; Li, H.; Page, A.J. Involvement of TRPV1 channels in energy homeostasis. Front. Endocrinol. 2018, 9, 420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Dhakal, S.; Lee, Y. Transient receptor potential channels and metabolism. Mol. Cells 2019, 42, 569–578. [Google Scholar] [PubMed]
  3. Ma, W.; Chen, X.; Cerne, R.; Syed, S.K.; Ficorilli, J.V.; Cabrera, O.; Obukhov, A.G.; Efanov, A.M. Catechol estrogens stimulate insulin secretion in pancreatic β-cells via activation of the transient receptor potential A1 (TRPA1) channel. J. Biol. Chem. 2019, 294, 2935–5880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Motter, A.L.; Ahern, G.P. TRPV1-null mice are protected from diet-induced obesity. FEBS Lett. 2008, 582, 2257–2262. [Google Scholar] [CrossRef] [Green Version]
  5. Lee, E.; Jung, D.Y.; Kim, J.H.; Patel, P.R.; Hu, X.; Lee, Y.; Azuma, Y.; Wang, H.F.; Tsitsilianos, N.; Shafiq, U.; et al. Transient receptor potential vanilloid type-1 channel regulates diet-induced obesity, insulin resistance, and leptin resistance. FASEB J. 2015, 29, 3182–3192. [Google Scholar] [CrossRef] [Green Version]
  6. McCoy, D.D.; Zhou, L.; Nguyen, A.K.; Watts, A.G.; Donovan, C.M.; McKemy, D.D. Enhanced insulin clearance in mice lacking TRPM8 channels. Am. J. Physiol. Endocrinol. Metab. 2013, 305, E78–E88. [Google Scholar] [CrossRef] [Green Version]
  7. Page, A.J.; Hatzinikolas, G.; Vincent, A.D.; Cavuoto, P.; Wittert, G.A. The TRPV1 channel regulates glucose metabolism. Am. J. Physiol. Endocrinol. Metab. 2019, 317, E667–E676. [Google Scholar] [CrossRef]
  8. Zhong, B.; Ma, S.; Wang, D.H. TRPV1 mediates glucose-induced insulin secretion through releasing neuropeptides. In Vivo 2019, 33, 1431. [Google Scholar] [CrossRef] [Green Version]
  9. Ferdowsi, P.V.; Ahuja, K.D.K.; Beckett, J.M.; Myers, S. TRPV1 activation by capsaicin mediates glucose oxidation and ATP production independent of insulin signalling in mouse skeletal muscle cells. Cells 2021, 10, 1560. [Google Scholar] [CrossRef]
  10. Suri, A.; Szallasi, A. The emerging role of TRPV1 in diabetes and obesity. Trends Pharmacol. Sci. 2008, 29, 29–36. [Google Scholar] [CrossRef]
  11. Zhang, S.; Tang, L.; Xu, F.; Hui, Y.; Lu, H.; Liu, X. TRPV1 receptor-mediated hypoglycemic mechanism of capsaicin in streptozotocin-induced diabetic rats. Front. Nutr. 2021, 8, 750355. [Google Scholar] [CrossRef] [PubMed]
  12. Li, M.; Zhang, C.S.; Zong, Y.; Feng, J.W.; Ma, T.; Hu, M.; Lin, Z.; Li, X.; Xie, C.; Wu, Y.; et al. Transient receptor potential V channels are essential for glucose sensing by aldolase and AMPK. Cell Metab. 2019, 30, 508–524.e12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Moraes, M.N.; de Assis, L.V.M.; Henriques, F.D.S.; Batista, M.L., Jr.; Güler, A.D.; Castrucci, A.M.L. Cold-sensing TRPM8 channel participates in circadian control of the brown adipose tissue. Biochim. Biophys. Acta Mol. Cell Res. 2017, 1864, 2415–2427. [Google Scholar] [CrossRef] [PubMed]
  14. Moraes, M.N.; Mezzalira, N.; de Assis, L.V.; Menaker, M.; Guler, A.; Castrucci, A.M.L. TRPV1 participates in the activation of clock molecular machinery in the brown adipose tissue in response to light-dark cycle. Biochim. Biophys. Acta Mol. Cell Res. 2017, 1864, 324–335. [Google Scholar] [CrossRef]
  15. Jerônimo, R.; Moraes, M.N.; de Assis, L.V.M.; Ramos, B.C.; Rocha, T.; Castrucci, A.M.L. Thermal stress in Danio rerio: A link between temperature, light, thermo-TRP channels, and clock genes. J. Thermal Biol. 2017, 68, 128–138. [Google Scholar] [CrossRef]
  16. Wang, Y.; Song, L.; Liu, M.; Ge, R.; Zhou, Q.; Liu, W.; Li, R.; Qie, J.; Zhen, B.; Wang, Y.; et al. A proteomics landscape of circadian clock in mouse liver. Nat. Commun. 2018, 9, 1553. [Google Scholar] [CrossRef]
  17. Mauvoisin, D. Circadian rhythms and proteomics: It’s all about posttranslational modifications! Wiley Interdiscip. Rev. Syst. Biol. Med. 2019, 11, e1450. [Google Scholar] [CrossRef]
  18. Greco, C.M.; Koronowski, K.B.; Smith, J.G.; Shi, J.; Kunderfranco, P.; Carriero, R.; Chen, S.; Samad, M.; Welz, P.S.; Zinna, V.M.; et al. Integration of feeding behavior by the liver circadian clock reveals network dependency of metabolic rhythms. Sci. Adv. 2021, 7, eabi7828. [Google Scholar] [CrossRef]
  19. Panda, S.; Antoch, M.P.; Miller, B.H.; Su, A.I.; Schook, A.B.; Straume, M.; Schultz, P.G.; Kay, S.A.; Takahashi, J.S.; Hogenesch, J.B. Coordinated transcription of key pathways in the mouse by the circadian clock. Cell 2002, 109, 307–320. [Google Scholar] [CrossRef] [Green Version]
  20. Koronowski, K.B.; Kinouchi, K.; Welz, P.S.; Smith, J.G.; Zinna, V.M.; Shi, J.; Samad, M.; Chen, S.; Magnan, C.N.; Kinchen, J.M.; et al. Defining the independence of the liver circadian clock. Cell 2019, 177, 1448–1462.e14. [Google Scholar] [CrossRef]
  21. Mukherji, A.; Bailey, S.M.; Staels, B.; Baumert, T.F. The circadian clock and liver function in health and disease. J. Hepatol. 2019, 71, 200–211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Doi, R.; Oishi, K.; Ishida, N. CLOCK regulates circadian rhythms of hepatic glycogen synthesis through transcriptional activation of Gys2. J. Biol. Chem. 2010, 285, 22114–22121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Rychkov, G.Y.; Barritt, G.J. Expression and function of TRP channels in liver cells. Adv. Exp. Med. Biol. 2011, 704, 667–686. [Google Scholar]
  24. Marcheva, B.; Ramsey, K.M.; Peek, C.B.; Affinati, A.; Maury, E.; Bass, J. Circadian clocks and metabolism. Handb. Exp. Pharmacol. 2013, 217, 127–155. [Google Scholar]
  25. de Assis, L.V.M.; Oster, H. The circadian clock and metabolic homeostasis: Entangled networks. Cell. Mol. Life Sci. 2021, 78, 4563–4587. [Google Scholar] [CrossRef]
  26. Chen, L.; Yang, G. PPARs integrate the mammalian clock and energy metabolism. PPAR Res. 2014, 2014, 653017. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, Y.; Nakajima, T.; Gonzalez, F.J.; Tanaka, N. PPARs as metabolic regulators in the liver: Lessons from liver-specific PPAR-null mice. Internat. J. Mol. Sci. 2020, 21, 2061. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, X.; Downes, M.; Yu, R.T.; Bookout, A.L.; He, W.; Straume, M.; Mangelsdorf, D.J.; Evans, R. MNuclear receptor expression links the circadian clock to metabolism. Cell 2006, 126, 801–810. [Google Scholar] [CrossRef] [Green Version]
  29. de Assis, L.V.M.; Moraes, M.N.; Magalhães-Marques, K.K.; Kinker, G.S.; da Silveira Cruz-Machado, S.; Castrucci, A.M.L. Non-metastatic cutaneous melanoma induces chronodisruption in central and peripheral circadian clocks. Int. J. Mol. Sci. 2018, 19, 1065. [Google Scholar] [CrossRef] [Green Version]
  30. Bargut, T.C.L.; Aguila, M.B.; Mandarim-de-Lacerda, C.A. Brown adipose tissue: Updates in cellular and molecular biology. Tissue Cell 2016, 48, 452–460. [Google Scholar] [CrossRef]
  31. Zhang, R.; Lahens, N.F.; Balance, H.I.; Hughes, M.E.; Hogenesch, J.B. A circadian gene expression atlas in mammals: Implications for biology and medicine. Proc. Natl. Acad. Sci. USA 2014, 111, 16219–16224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. McGinnis, G.R.; Young, M.E. Circadian regulation of metabolic homeostasis: Causes and consequences. Nat. Sci. Sleep 2016, 8, 163–180. [Google Scholar] [PubMed] [Green Version]
  33. Pan, X.; Hussain, M.M. Diurnal regulation of microsomal triglyceride transfer protein and plasma lipid levels. J. Biol. Chem. 2007, 282, 24707–24719. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Adamovich, Y.; Rousso-Noori, L.; Zwighaft, Z.; Neufeld-Cohen, A.; Golik, M.; Kraut-Cohen, J.; Wang, M.; Han, X.; Asher, G. Circadian clocks and feeding time regulate the oscillations and levels of hepatic triglycerides. Cell Metab. 2014, 19, 319–330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, L.; Chen, J.; Ni, Y.; Feng, X.; Zhao, Z.; Wang, P.; Sun, J.; Yu, H.; Yan, Z.; Liu, D.; et al. TRPV1 activation prevents nonalcoholic fatty liver through UCP2 upregulation in mice. Pflug. Arch. 2012, 463, 727–732. [Google Scholar] [CrossRef]
  36. Sun, W.; Uchida, K.; Suzuki, Y.; Zhou, Y.; Kim, M.; Takayama, Y.; Takahashi, N.; Goto, T.; Wakabayashi, S.; Kawada, T.; et al. Lack of TRPV2 impairs thermogenesis in mouse brown adipose tissue. EMBO Rep. 2016, 17, 383–399. [Google Scholar] [CrossRef]
  37. Gachon, F.; Loizides-Mangold, U.; Petrenko, V.; Dibner, C. Glucose homeostasis: Regulation by peripheral circadian clocks in rodents and humans. Endocrinology 2017, 158, 1074–1084. [Google Scholar] [CrossRef]
  38. Thorens, B. GLUT2, glucose sensing and glucose homeostasis. Diabetologia 2015, 58, 221–232. [Google Scholar] [CrossRef] [Green Version]
  39. Takeda, Y.; Kang, H.S.; Freudenberg, J.; DeGraff, L.M.; Jothi, R.; Jetten, A.M. Retinoic acid-related orphan receptor γ (RORγ): A novel participant in the diurnal regulation of hepatic gluconeogenesis and insulin sensitivity. PLoS Genet. 2014, 10, e1004331. [Google Scholar] [CrossRef]
  40. Titchenell, P.M.; Lazar, M.A.; Birnbaum, M.J. Unraveling the regulation of hepatic metabolism by insulin. Trends Endocrinol. Metab. 2017, 28, 497–505. [Google Scholar] [CrossRef]
  41. Kornberg, R.D.; Lorch, Y. Primary role of the nucleosome. Mol. Cell 2020, 79, 371–375. [Google Scholar] [CrossRef] [PubMed]
  42. Körner, A.; Zhou, E.; Müller, C.; Mohammed, Y.; Herceg, S.; Bracher, F.; Rensen, P.C.N.; Wang, Y.; Mirakaj, V.; Giera, M. Inhibition of Δ24-dehydrocholesterol reductase activates pro-resolving lipid mediator biosynthesis and inflammation resolution. Proc. Natl. Acad. Sci. USA 2019, 116, 20623–20634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Soto-Rifo, R.; Ohlmann, T. The role of the DEAD-box RNA helicase DDX3 in mRNA metabolism. Wiley Interdiscip. Rev. RNA 2013, 4, 369–385. [Google Scholar] [CrossRef] [PubMed]
  44. Tsai, T.Y.; Wang, W.T.; Li, H.K.; Chen, W.J.; Tsai, Y.H.; Chao, C.H.; Wu Lee, Y.H. RNA helicase DDX3 maintains lipid homeostasis through upregulation of the microsomal triglyceride transfer protein by interacting with HNF4 and SHP. Sci. Rep. 2017, 7, 41452. [Google Scholar] [CrossRef] [PubMed]
  45. Droppelmann, C.A.; Sáez, D.E.; Asenjo, J.L.; Yáñez, A.J.; García-Rocha, M.; Concha, I.I.; Grez, M.; Guinovart, J.J.; Slebe, J.C. A new level of regulation in gluconeogenesis: Metabolic state modulates the intracellular localization of aldolase B and its interaction with liver fructose-1,6-bisphosphatase. Biochem. J. 2015, 472, 225–237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Chang, Y.C.; Yang, Y.C.; Tien, C.P.; Yang, C.J.; Hsiao, M. Roles of aldolase family genes in human cancers and diseases. Trends Endocrinol. Metab. 2018, 29, 549–559. [Google Scholar] [CrossRef]
  47. Ferguson, D.; Zhang, J.; Davis, M.A.; Helsley, R.N.; Vedin, L.L.; Lee, R.G.; Crooke, R.M.; Graham, M.J.; Allende, D.S.; Parini, P.; et al. The lipid droplet-associated protein perilipin 3 facilitates hepatitis C virus-driven hepatic steatosis. J. Lipid Res. 2017, 58, 420–432. [Google Scholar] [CrossRef] [Green Version]
  48. Díaz, E.; Pfeffer, S.R. TIP47: A cargo selection device for mannose 6-phosphate receptor trafficking. Cell 1998, 93, 433–443. [Google Scholar] [CrossRef] [Green Version]
  49. Naslavsky, N.; McKenzie, J.; Altan-Bonnet, N.; Sheff, D.; Caplan, S. EHD3 regulates early-endosome-to-Golgi transport and preserves Golgi morphology. J. Cell Sci. 2009, 122, 389–400. [Google Scholar] [CrossRef] [Green Version]
  50. Tan, K.; Fujimoto, M.; Takii, R.; Takaki, E.; Hayashida, N.; Nakai, A. Mitochondrial SSBP1 protects cells from proteotoxic stresses by potentiating stress-induced HSF1 transcriptional activity. Nat. Commun. 2015, 6, 6580. [Google Scholar] [CrossRef] [Green Version]
  51. Becker, J.S.; McCarthy, R.L.; Sidoli, S.; Donahue, G.; Kaeding, K.E.; He, Z.; Lin, S.; Garcia, B.A.; Zaret, K.S. Genomic and proteomic resolution of heterochromatin and its restriction of alternate fate genes. Mol. Cell 2017, 68, 1023–1037.e15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Prieto, C.; Nguyen, D.T.T.; Liu, Z.; Wheat, J.; Perez, A.; Gourkanti, S.; Chou, T.; Barin, E.; Velleca, A.; Rohwetter, T.; et al. Transcriptional control of CBX5 by the RNA-binding proteins RBMX and RBMXL1 maintains chromatin state in myeloid leukemia. Nat. Cancer 2021, 2, 741–757. [Google Scholar] [CrossRef] [PubMed]
  53. Lande-Diner, L.; Boyault, C.; Kim, J.Y.; Weitz, C.J. A positive feedback loop links circadian clock factor CLOCK-BMAL1 to the basic transcriptional machinery. Proc. Natl. Acad. Sci. USA 2013, 110, 16021–16026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Beli, P.; Lukashchuk, N.; Wagner, S.A.; Weinert, B.T.; Olsen, J.V.; Baskcomb, L.; Mann, M.; Jackson, S.P.; Choudhary, C. Proteomic investigations reveal a role for RNA processing factor THRAP3 in the DNA damage response. Mol. Cell Mol. Cell 2012, 46, 212–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Pérez, V.I.; Lew, C.M.; Cortez, L.A.; Webb, C.R.; Rodriguez, M.; Liu, Y.; Qi, W.; Li, Y.; Chaudhuri, A.; Van Remmen, H.; et al. Thioredoxin 2 haploinsufficiency in mice results in impaired mitochondrial function and increased oxidative stress. Free Radic. Biol. Med. 2008, 44, 882–892. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, X.; Xing, Y.; Tang, Z.; Tang, Y.; Shen, J.; Zhang, F. Thioredoxin-2 impacts the inflammatory response via suppression of NF-κB and MAPK signaling in sepsis shock. Biochem. Biophys. Res. Commun. 2020, 524, 876–882. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, W.T.; Jing, Y.Y.; Han, Z.P.; Li, X.N.; Liu, Y.; Lai, F.B.; Li, R.; Zhao, Q.D.; Wu, M.C.; Wei, L.X. The injured liver induces hyperimmunoglobulinemia by failing to dispose of antigens and endotoxins in the portal system. PLoS ONE 2015, 10, e0122739. [Google Scholar] [CrossRef]
  58. Fang, X.; Zaman, M.H.; Guo, X.; Ding, H.; Xie, C.; Zhang, X.; Deng, G.M. Role of hepatic deposited immunoglobulin G in the pathogenesis of liver damage in systemic Lupus Erythematosus. Front. Immunol. 2018, 9, 1457. [Google Scholar] [CrossRef]
  59. Luo, X.Y.; Takahara, T.; Kawai, K.; Fujino, M.; Sugiyama, T.; Tsuneyama, K.; Tsukada, K.; Nakae, S.; Zhong, L.; Li, X.K. IFN-γ deficiency attenuates hepatic inflammation and fibrosis in a steatohepatitis model induced by a methionine- and choline-deficient high-fat diet. Am. J. Physiol. Gastrointest. Liver Physiol. 2013, 305, G891–G899. [Google Scholar] [CrossRef] [Green Version]
  60. Shimada, H.; Hashimoto, R.; Aoki, A.; Yamada, S.; Oba, K.I.; Kawase, A.; Nakanishi, T.; Iwaki, M. The regulatory mechanism involved in the prostaglandin E2 disposition in carbon tetrachloride-induced liver injury. Prostaglandins Leukot. Essent. Fat. Acids 2020, 155, 102081. [Google Scholar] [CrossRef]
  61. Gromovsky, A.D.; Schugar, R.C.; Brown, A.L.; Helsley, R.N.; Burrows, A.C.; Ferguson, D.; Zhang, R.; Sansbury, B.E.; Lee, R.G.; Morton, R.E.; et al. Δ-5 Fatty acid desaturase FADS1 impacts metabolic disease by balancing proinflammatory and proresolving lipid mediators. Arterioscler. Thromb. Vasc. Biol. 2018, 38, 218–231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zhang, K.; Kim, H.; Fu, Z.; Qiu, Y.; Yang, Z.; Wang, J.; Zhang, D.; Tong, X.; Yin, L.; Li, J.; et al. Deficiency of the mitochondrial NAD kinase causes stress-induced hepatic steatosis in mice. Gastroenterology 2018, 154, 224–237. [Google Scholar] [CrossRef] [PubMed]
  63. Garcia, D.; Hellberg, K.; Chaix, A.; Wallace, M.; Herzig, S.; Badur, M.G.; Lin, T.; Shokhirev, M.N.; Pinto, A.F.M.; Ross, D.S.; et al. Genetic liver-specific AMPK activation protects against diet-induced obesity and NAFLD. Cell Rep. 2019, 26, 192–208.e6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Bowman, C.E.; Rodriguez, S.; Selen Alpergin, E.S.; Acoba, M.G.; Zhao, L.; Hartung, T.; Claypool, S.M.; Watkins, P.A.; Wolfgang, M.J. The mammalian malonyl-CoA synthetase ACSF3 is required for mitochondrial protein malonylation and metabolic efficiency. Cell Chem. Biol. 2017, 24, 673–684.e4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Abu-Elheiga, L.; Brinkley, W.R.; Zhong, L.; Chirala, S.S.; Woldegiorgis, G.; Wakil, S.J. The subcellular localization of acetyl-CoA carboxylase 2. Proc. Natl. Acad. Sci. USA 2000, 97, 1444–1449. [Google Scholar] [CrossRef] [Green Version]
  66. Kim, C.W.; Addy, C.; Kusunoki, J.; Anderson, N.N.; Deja, S.; Fu, X.; Burgess, S.C.; Li, C.; Ruddy, M.; Chakravarthy, M.; et al. Acetyl CoA carboxylase inhibition reduces hepatic steatosis but elevates plasma triglycerides in mice and humans: A bedside to bench investigation. Cell Metab. 2017, 26, 394–406.e6. [Google Scholar] [CrossRef]
  67. Storch, J.; McDermott, L. Structural and functional analysis of fatty acid-binding proteins. J. Lipid Res. 2009, 50, S126–S131. [Google Scholar] [CrossRef] [Green Version]
  68. Zheng, Q.W.; Ding, X.F.; Cao, H.J.; Ni, Q.Z.; Zhu, B.; Ma, N.; Zhang, F.K.; Wang, Y.K.; Xu, S.; Chen, T.W.; et al. Ochratoxin A induces steatosis via PPARγ-CD36 axis. Toxins 2021, 13, 802. [Google Scholar] [CrossRef]
  69. Wang, G.; Li, M.; Yu, S.; Guan, M.; Ma, S.; Zhong, Z.; Guo, Y.; Leng, X.; Huang, H. Tandem mass tag-based proteomics analysis of type 2 diabetes mellitus with non-alcoholic fatty liver disease in mice treated with acupuncture. Biosci. Rep. 2022, 42, BSR20212248. [Google Scholar] [CrossRef]
  70. Ipsen, D.H.; Lykkesfeldt, J.; Tveden-Nyborg, P. Molecular mechanisms of hepatic lipid accumulation in non-alcoholic fatty liver disease. Cell Mol. Life Sci. 2018, 75, 3313–3327. [Google Scholar] [CrossRef] [Green Version]
  71. Zhu, Y.; Zhao, S.; Deng, Y.; Gordillo, R.; Ghaben, A.L.; Shao, M.; Zhang, F.; Xu, P.; Li, Y.; Cao, H.; et al. Hepatic GALE regulates whole-body glucose homeostasis by modulating Tff3 expression. Diabetes 2017, 66, 2789–2799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Wolosker, H.; Kline, D.; Bian, Y.; Blackshaw, S.; Cameron, A.M.; Fralich, T.J.; Schnaar, R.L.; Snyder, S.H. Molecularly cloned mammalian glucosamine-6-phosphate deaminase localizes to transporting epithelium and lacks oscillin activity. FASEB J. 1998, 12, 91–99. [Google Scholar] [PubMed]
  73. Rohacs, T.; Thyagarajan, B.; Lukacs, V. Phospholipase C mediated modulation of TRPV1 channels. Mol. Neurobiol. 2008, 37, 153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Zygmunt, P.M.; Ermund, A.; Movahed, P.; Andersson, D.A.; Simonsen, C.; Jönsson, B.A.; Blomgren, A.; Birnir, B.; Bevan, S.; Eschalier, A.; et al. Monoacylglycerols activate TRPV1–A link between phospholipase C and TRPV1. PLoS ONE 2013, 8, e81618. [Google Scholar] [CrossRef] [Green Version]
  75. Pecze, L.; Blum, W.; Henzi, T.; Schwaller, B. Endogenous TRPV1 stimulation leads to the activation of the inositol phospholipid pathway necessary for sustained Ca2+ oscillations. Biochim. Biophys. Acta 2016, 1863, 2905–2915. [Google Scholar] [CrossRef]
  76. Liu, L.; Yudin, Y.; Rohacs, T. Diacylglycerol kinases regulate TRPV1 channel activity. J. Biol. Chem. 2020, 295, 8174–8185. [Google Scholar] [CrossRef]
  77. Bartrons, R.; Hue, L.; Van Schaftingen, E.; Hers, H.G. Hormonal control of fructose 2,6-bisphosphate concentration in isolated rat hepatocytes. Biochem. J. 1983, 214, 829–837. [Google Scholar] [CrossRef] [Green Version]
  78. Rider, M.H.; Bertrand, L.; Vertommen, D.; Michels, P.A.; Rousseau, G.G.; Hue, L. 6-Phosphofructo-2-kinase/fructose-2,6-bisphosphatase: Head-to-head with a bifunctional enzyme that controls glycolysis. Biochem. J. 2004, 381, 561–579. [Google Scholar] [CrossRef] [Green Version]
  79. Marsin, A.S.; Bertrand, L.; Rider, M.H.; Deprez, J.; Beauloye, C.; Vincent, M.F.; Van den Berghe, G.; Carling, D.; Hue, L. Phosphorylation and activation of heart PFK-2 by AMPK has a role in the stimulation of glycolysis during ischaemia. Curr. Biol. 2000, 10, 1247–1255. [Google Scholar] [CrossRef] [Green Version]
  80. Hardie, D.G. Metformin-acting through Cyclic AMP as well as AMP? Cell Metab. 2013, 17, 313–314. [Google Scholar] [CrossRef] [Green Version]
  81. Rakus, D.; Gizak, A.; Kasprzak, A.A.; Zarzycki, M.; Maciaszczyk-Dziubinska, E.; Dzugaj, A. The mechanism of calcium-induced inhibition of muscle fructose 1,6-bisphosphatase and destabilization of glyconeogenic complex. PLoS ONE 2013, 8, e76669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Halls, M.L.; Cooper, D.M.F. Regulation by Ca2+-signaling pathways of adenylyl cyclases. Cold Spring Harb. Perspect. Biol. 2011, 3, a004143. [Google Scholar] [CrossRef] [PubMed]
  83. Niswender, C.M.; Willis, B.S.; Wallen, A.; Sweet, I.R.; Jetton, T.L.; Thompson, B.R.; Wu, C.; Lange, A.J.; McKnight, G.S. Cre recombinase-dependent expression of a constitutively active mutant allele of the catalytic subunit of protein kinase A. Genesis 2005, 43, 109–119. [Google Scholar] [CrossRef]
  84. Kotani, T.; Iemura, S.; Natsume, T.; Kawakami, K.; Yamashita, M. Mys protein regulates protein kinase A activity by interacting with regulatory type Ialpha subunit during vertebrate development. J. Biol. Chem. 2010, 285, 5106–5116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Li, H.; Adamik, R.; Pacheco-Rodriguez, G.; Moss, J.; Vaughan, M. Protein kinase A-anchoring (AKAP) domains in brefeldin A-inhibited guanine nucleotide-exchange protein 2 (BIG2). Proc. Natl. Acad. Sci. USA 2003, 100, 1627–1632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Ahn, J.H.; McAvoy, T.; Rakhilin, S.V.; Nishi, A.; Greengard, P.; Nairn, A.C. Protein kinase A activates protein phosphatase 2A by phosphorylation of the B56delta subunit. Proc. Natl. Acad. Sci. USA 2007, 104, 2979–2984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Guo, F.; Stanevich, V.; Wlodarchak, N.; Sengupta, R.; Jiang, L.; Satyshur, K.A.; Xing, Y. Structural basis of PP2A activation by PTPA, an ATP-dependent activation chaperone. Cell Res. 2014, 24, 190–203. [Google Scholar] [CrossRef] [Green Version]
  88. Abiria, S.A.; Krapivinsky, G.; Sah, R.; Santa-Cruz, A.G.; Chaudhuri, D.; Zhang, J.; Adstamongkonkul, P.; DeCaen, P.G.; Clapham, D.E. TRPM7 senses oxidative stress to release Zn2+ from unique intracellular vesicles. Proc. Natl. Acad. Sci. USA 2017, 114, E6079–E6088. [Google Scholar] [CrossRef] [Green Version]
  89. Thapa, C.; Roivas, P.; Haataja, T.; Permi, P.; Pentikäinen, U. Interaction mechanism of endogenous PP2A inhibitor protein ENSA with PP2A. FEBS J. 2022, 289, 519–534. [Google Scholar] [CrossRef]
  90. Yadav, H.; Devalaraja, S.; Chung, S.T.; Rane, S.G. TGF-β1/Smad3 pathway targets PP2A-AMPK-FoxO1 signaling to regulate hepatic gluconeogenesis. J. Biol. Chem. 2017, 292, 3420–3432. [Google Scholar] [CrossRef] [Green Version]
  91. Hernández, F.; Langa, E.; Cuadros, R.; Avila, J.; Villanueva, N. Regulation of GSK3 isoforms by phosphatases PP1 and PP2A. Mol. Cell. Biochem. 2010, 344, 211–215. [Google Scholar] [CrossRef] [PubMed]
  92. Zhang, T.; Wang, S.; Lin, Y.; Xu, W.; Ye, D.; Xiong, Y.; Zhao, S.; Guan, K.L. Acetylation negatively regulates glycogen phosphorylase by recruiting protein phosphatase 1. Cell Metab. 2012, 15, 75–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  94. Munhoz, A.C.; Vilas-Boas, E.A.; Panveloski-Costa, A.C.; Leite, J.S.M.; Lucena, C.F.; Riva, P.; Emilio, H.; Carpinelli, A.R. Intermittent fasting for twelve weeks leads to increases in fat mass and hyperinsulinemia in young female Wistar rats. Nutrients 2020, 12, 1029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Wiśniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359–362. [Google Scholar] [CrossRef]
  96. Wiśniewski, J.R.; Zielinska, D.F.; Mann, M. Comparison of ultrafiltration units for proteomic and N-glycoproteomic analysis by the filter-aided sample preparation method. Anal. Biochem. 2011, 410, 307–309. [Google Scholar] [CrossRef]
  97. Wiśniewski, J.R. Quantitative evaluation of filter aided sample preparation (FASP) and multienzyme digestion FASP protocols. Anal. Chem. 2016, 88, 5438–5443. [Google Scholar] [CrossRef] [Green Version]
  98. Ma, B.; Johnson, R. De novo sequencing and homology searching. Mol. Cell. Proteom. 2012, 11, O111.014902. [Google Scholar] [CrossRef] [Green Version]
  99. Zhang, J.; Xin, L.; Shan, B.; Chen, W.; Xie, M.; Yuen, D.; Zhang, W.; Zhang, Z.; Lajoie, G.A.; Ma, B. PEAKS DB: De novo sequencing assisted database search for sensitive and accurate peptide identification. Mol. Cell. Proteom. 2012, 11, M111.010587. [Google Scholar] [CrossRef] [Green Version]
  100. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  101. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  102. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  103. Carlin, D.E.; Demchak, B.; Pratt, D.; Sage, E.; Ideker, T. Network propagation in the cytoscape cyberinfrastructure. PLoS Comput. Biol. 2017, 13, e1005598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Cowen, L.; Ideker, T.; Raphael, B.J.; Sharan, R. Network propagation: A universal amplifier of genetic associations. Nat. Rev. Genet. 2017, 18, 551–562. [Google Scholar] [CrossRef]
  105. Martin-Sancho, L.; Lewinski, M.K.; Pache, L.; Stoneham, C.A.; Yin, X.; Becker, M.E.; Pratt, D.; Churas, C.; Rosenthal, S.B.; Liu, S.; et al. Functional landscape of SARS-CoV-2 cellular restriction. Mol. Cell 2021, 81, P2656–P2668.e8. [Google Scholar] [CrossRef]
  106. Su, G.; Kuchinsky, A.; Morris, J.H.; States, D.J.; Meng, F. GLay: Community structure analysis of biological networks. Bioinformatics 2010, 26, 3135–3137. [Google Scholar] [CrossRef] [Green Version]
  107. Morris, J.H.; Apeltsin, L.; Newman, A.M.; Baumbach, J.; Wittkop, T.; Su, G.; Bader, G.D.; Ferrin, T.E. ClusterMaker: A multi-algorithm clustering plugin for Cytoscape. BMC Bioinform. 2011, 12, 1–14. [Google Scholar] [CrossRef] [Green Version]
  108. Merico, D.; Isserlin, R.; Stueker, O.; Emili, A.; Bader, G.D. Enrichment map: A network-based method for gene-set enrichment visualization and interpretation. PLoS ONE 2010, 5, e13984. [Google Scholar] [CrossRef]
Ijms 23 07014 g001 550
Figure 1. Gene expression in the liver of C57BL/6J WT and Trpv1 KO mice. (A–E) Gene expression in the liver at ZT8 (8 h after lights on) and ZT20 (8 h after lights off). Results are shown as median, quantiles, and maximum and minimum expression values of genes of interest normalized by 18S ribosomal RNA and expressed relative to the minimal value of WT mice. + represents the mean (n = 4–5). In (A), statistical analysis was performed by Student’s t-test. (B–E) Significant differences between ZTs and between genotypes were demonstrated by two-Way ANOVA, followed by the Bonferroni post-test.
Figure 1. Gene expression in the liver of C57BL/6J WT and Trpv1 KO mice. (A–E) Gene expression in the liver at ZT8 (8 h after lights on) and ZT20 (8 h after lights off). Results are shown as median, quantiles, and maximum and minimum expression values of genes of interest normalized by 18S ribosomal RNA and expressed relative to the minimal value of WT mice. + represents the mean (n = 4–5). In (A), statistical analysis was performed by Student’s t-test. (B–E) Significant differences between ZTs and between genotypes were demonstrated by two-Way ANOVA, followed by the Bonferroni post-test.
Ijms 23 07014 g001
Ijms 23 07014 g002 550
Figure 2. Hepatic glucose metabolism in the absence of TRPV1 channels. (A–D,F) Relative transcripts of hepatic genes related to glucose metabolism in C57BL/6J WT and Trpv1 KO mice at ZT8 (physiological fasting) and ZT20 (eating phase). Results of mRNA quantification are shown as median, quantiles, and maximum and minimum expression values of gene transcripts normalized by 18S ribosomal RNA and expressed relative to the minimum mean value of WT mice. ₊ represents the mean (n = 4–5). Significant differences between ZTs and between genotypes were demonstrated by two-way ANOVA, followed by the Bonferroni post-test. (E) Hepatic glycogen quantification is expressed as mass (mg/g of tissue).
Figure 2. Hepatic glucose metabolism in the absence of TRPV1 channels. (A–D,F) Relative transcripts of hepatic genes related to glucose metabolism in C57BL/6J WT and Trpv1 KO mice at ZT8 (physiological fasting) and ZT20 (eating phase). Results of mRNA quantification are shown as median, quantiles, and maximum and minimum expression values of gene transcripts normalized by 18S ribosomal RNA and expressed relative to the minimum mean value of WT mice. ₊ represents the mean (n = 4–5). Significant differences between ZTs and between genotypes were demonstrated by two-way ANOVA, followed by the Bonferroni post-test. (E) Hepatic glycogen quantification is expressed as mass (mg/g of tissue).
Ijms 23 07014 g002
Ijms 23 07014 g003 550
Figure 3. The lack of a TRPV1 channel affects hepatic functionality, as revealed by the integration data of liver proteomes. (A) The UpSet plot presents the overlapping of proteins found in the liver proteome of Trpv1 KO and WT mice. (B) Scatter plots show the distribution of proteins shared between liver proteomes, indicating downregulated and upregulated proteins in Trpv1 KO mice. LFQ difference denotes log2 (fold change) among groups. (C) Gene ontology biological process (GO-BP) enrichment analyses of differentially identified proteins in the liver proteomes using the annotation of the DAVID database (p < 0.05) are highlighted for lipid metabolism, RNA processing, oxidative stress, and inflammation functional processes that were exclusively expressed in WT or Trpv1 KO animals.
Figure 3. The lack of a TRPV1 channel affects hepatic functionality, as revealed by the integration data of liver proteomes. (A) The UpSet plot presents the overlapping of proteins found in the liver proteome of Trpv1 KO and WT mice. (B) Scatter plots show the distribution of proteins shared between liver proteomes, indicating downregulated and upregulated proteins in Trpv1 KO mice. LFQ difference denotes log2 (fold change) among groups. (C) Gene ontology biological process (GO-BP) enrichment analyses of differentially identified proteins in the liver proteomes using the annotation of the DAVID database (p < 0.05) are highlighted for lipid metabolism, RNA processing, oxidative stress, and inflammation functional processes that were exclusively expressed in WT or Trpv1 KO animals.
Ijms 23 07014 g003
Ijms 23 07014 g004 550
Figure 4. Network propagation highlights important liver clustering pathways and biological networks in the absence of TRPV1 channels. Network propagation of the TRPV1 channel in the liver of WT (A) and Trpv1 KO mice (B). TRPV1 protein input under each network served as a driver for network propagation in Trpv1 KO and WT. Separately propagated subnetworks with 141 proteins found to be significant in the propagation were clustered into smaller subnetworks and annotated using integrated enrichment analysis (GO-BP and KEGG pathway) to identify functional networks that are different between two propagated subnetworks (EnrichmentMap FDR < 0.01). Node names represent the proteins exclusively found in each propagated subnetwork; black for proteins shared between Trpv1 KO and WT proteomes; blue for proteins only found or upregulated in the WT proteome; red for proteins only found or upregulated in the Trpv1 KO proteome. Cluster title indicates the top five most significant biological function(s) for each cluster, and (#) refers to the KEGG pathway/GO-BP annotations exclusively enriched in each propagated subnetwork. (C) Framework of propagated network integration for enriched glucagon signaling pathway from KEGG in the WT proteome. (D) Framework shows proposed modulation of FBP-2 by PKA and AMPK kinases in Trpv1 KO mice resulting in more active gluconeogenesis than in WT mice. Blue nodes for proteins only found in the WT proteome, dual colors (blue/red) for proteins shared among proteomes, and black for proteins not found in the proteomic data. (+p) refers to phosphorylation and (−p) to dephosphorylation.
Figure 4. Network propagation highlights important liver clustering pathways and biological networks in the absence of TRPV1 channels. Network propagation of the TRPV1 channel in the liver of WT (A) and Trpv1 KO mice (B). TRPV1 protein input under each network served as a driver for network propagation in Trpv1 KO and WT. Separately propagated subnetworks with 141 proteins found to be significant in the propagation were clustered into smaller subnetworks and annotated using integrated enrichment analysis (GO-BP and KEGG pathway) to identify functional networks that are different between two propagated subnetworks (EnrichmentMap FDR < 0.01). Node names represent the proteins exclusively found in each propagated subnetwork; black for proteins shared between Trpv1 KO and WT proteomes; blue for proteins only found or upregulated in the WT proteome; red for proteins only found or upregulated in the Trpv1 KO proteome. Cluster title indicates the top five most significant biological function(s) for each cluster, and (#) refers to the KEGG pathway/GO-BP annotations exclusively enriched in each propagated subnetwork. (C) Framework of propagated network integration for enriched glucagon signaling pathway from KEGG in the WT proteome. (D) Framework shows proposed modulation of FBP-2 by PKA and AMPK kinases in Trpv1 KO mice resulting in more active gluconeogenesis than in WT mice. Blue nodes for proteins only found in the WT proteome, dual colors (blue/red) for proteins shared among proteomes, and black for proteins not found in the proteomic data. (+p) refers to phosphorylation and (−p) to dephosphorylation.
Ijms 23 07014 g004
Ijms 23 07014 g005 550
Figure 5. Network integration between kinases and phosphatases for modulation of hepatic glycogenesis/glycogenolysis in the absence of TRPV1 channels. (A) Mapping of kinases, phosphatases, and their regulatory proteins exclusively identified in the liver of Trpv1 KO or WT mice or those that showed significant differences in the relative quantification by the LFQ approach. (B) Differentially identified kinases, phosphatases, and regulatory proteins of each proteome were integrated into the same network to identify mechanisms that were different between Trpv1 KO and WT livers. Blue color for proteins only found in the WT proteome, dual colors (blue/red) for proteins shared among proteomes, and black for proteins not found in the proteomic data. (+p) refers to phosphorylation and (−p) to dephosphorylation.
Figure 5. Network integration between kinases and phosphatases for modulation of hepatic glycogenesis/glycogenolysis in the absence of TRPV1 channels. (A) Mapping of kinases, phosphatases, and their regulatory proteins exclusively identified in the liver of Trpv1 KO or WT mice or those that showed significant differences in the relative quantification by the LFQ approach. (B) Differentially identified kinases, phosphatases, and regulatory proteins of each proteome were integrated into the same network to identify mechanisms that were different between Trpv1 KO and WT livers. Blue color for proteins only found in the WT proteome, dual colors (blue/red) for proteins shared among proteomes, and black for proteins not found in the proteomic data. (+p) refers to phosphorylation and (−p) to dephosphorylation.
Ijms 23 07014 g005
Table
Table 1. ANOVA parameters related to gene expression in the liver.
Table 1. ANOVA parameters related to gene expression in the liver.
Time F (DFn, DFd) and p-ValueGenotype F (DFn, DFd) and p-ValueInteraction F (DFn, DFd) and p-Value
Akt(1, 14) = 3.787p = 0.0720(1, 14) = 2.713p = 0.1218(1, 14) = 8.441p = 0.0115 *
Bmal(1, 13) = 34.92p < 0.0001 *(1, 13) = 6.57p = 0.0236 *(1, 12) = 0.318p = 0.5824
G6pase(1, 14) = 14.95p = 0.0017 *(1, 14) = 0.4033p = 0.5356(1, 14) = 0.1015p = 0.7547
Glut2(1, 14) = 1.429p = 0.2518(1, 14) = 6.557p = 0.0226 *(1, 14) = 6.231p = 0.0257 *
Gsk3b(1, 13) = 0.00001p = 0.9970(1, 13) = 5.553p = 0.0348 *(1, 13) = 3.349p = 0.0903
Pepck(1, 14) = 2.574p = 0.1310(1, 14) = 7.956p = 0.0136 *(1, 14) = 4.8750.0444 *
Per1(1, 13) = 0.7069p = 0.4157(1, 13) = 12.65p = 0.0035 *(1, 13) = 0.2005p = 0.6617
Pparα(1, 15) = 4.982p = 0.0413 *(1, 15) = 6.201p = 0.0250 *(1, 15) = 7.609p = 0.0146 *
Pparγ(1, 15) = 3.433p = 0.0837(1, 15) = 10.9p = 0.0048 *(1, 15) = 6.232p = 0.0247 *
* Represents statistically significant differences. DF = degree of freedom. DFn is calculated as the total group number—1 and DFd as the total number of mice—1 of each group. F distribution is based on a ratio of DFn and DFd.
Table
Table 2. List of 11 exclusively identified kinases, phosphatases, and associated regulator proteins reported to be upstream effectors in each proteome and their correspondingly mapped targets, followed by downstream signaling involved in glycogenesis and glycogenolysis.
Table 2. List of 11 exclusively identified kinases, phosphatases, and associated regulator proteins reported to be upstream effectors in each proteome and their correspondingly mapped targets, followed by downstream signaling involved in glycogenesis and glycogenolysis.
Proteins Identified in the Proteomes Downstream Target Mapped
GenotypeUpstream EffectorClassResearch SourceFunctionTargetClassResearch SourceFunctionDownstream Signaling
Trpv1 KODUSP3Phosphatasemap04010InhibitionMAPK14KinasePSP_454221InhibitionGSK3B
Trpv1 KOMAPK14KinasePSP_454221InhibitionGSK3BKinasemap04910InhibitionGYS2
Trpv1 KOMAPK14Kinasemap04068ActivationFOXOTFmap04910ActivationPEPCK, G6PC1, FBP1, FBP2
Trpv1 KOPTPARegulatorGuo et al., 2014ActivationPP2A complexPhosphataseHernández et al., 2010;
Yadav et al., 2017		Activation, Activation, InhibitionGSK3B,
FOXO,
AMPK
Trpv1 KOENSARegulatorThapa et al., 2020InhibitionPP2A complexPhosphataseHernández et al., 2010;
Yadav et al., 2017		Activation, Activation, InhibitionGSK3B,
FOXO,
AMPK
Trpv1 KOPRRC1RegulatorKotani et al., 2009ActivationPKAKinasePSP_449111, map04922,
Ahn et al., 2007		Inhibition, Inhibition, Activation, ActivationGSK3B,
GYS2,
PHKB,
PP2A complex
Trpv1 KOARFGEF2RegulatorLi et al., 2003ActivationPKAKinasePSP_449111, map04922,
Ahn et al., 2007		Inhibition, Inhibition, Activation, ActivationGSK3B,
GYS2,
PHKB,
PP2A complex
Trpv1 KOSGK2Kinasemap04068InhibitionFOXOTranscription factormap04910ActivationPEPCK, G6PC1, FBP1, FBP2
WTMAPK1Kinasemap04068InhibitionFOXOTranscription factormap04910ActivationPEPCK, G6PC1,FBP1, FBP2
WTATMKinasePSP_454221InhibitionGSK3BKinasemap04910InhibitionGYS2
WTPPP2R1BPhosphatasemap04152InhibitionAMPKKinasemap04152InhibitionGYS2
WTPPP1CBPhosphatasePSP_448786InhibitionAMPKKinasemap04152InhibitionGYS2
WTPPP1CBPhosphatasemap04910InhibitionPHKBKinasemap04922ActivationPYGL
WTPPP1CBPhosphatasemap04910ActivationGYS2Kinasemap04922Glycogen synthaseGlycogen
The research source included the KEGG database (ID: map), PhosphoSitePlus (ID: PSP), and literature data.
Table
Table 3. Gene access numbers and sequences of primers and probes.
Table 3. Gene access numbers and sequences of primers and probes.
Template
(Access Number)		Primers and Probes
18S rRNAForward: 5′-CGGCTACCACATCCAAGGAA-3′
Reverse: 5′-GCTGGAATTACCGCGGCT-3′
Akt
(NM_001165894.1)		Forward: 5′-CCGTGTGACCATGAACGAGT-3′
Reverse: 5′-GGTCGTGGGTCTGGAATGAG-3′
Bmal1
(NM_001243048)		Forward: 5′-AGCTTCTGCACAATCCACAGCAC-3′
Reverse: 5′-TGTCTGGCTCATTGTCTTCGTCCA-3′
Probe:5′/5HEX/AAAGCTGGCCACCCACGAAGATGGG/BHQ_1/-3′
G6pase
(NM_008061.4)		Forward: 5′-CAGTGGTCGGAGACTGGTTC-3′
Reverse: 5′-TATAGGCACGGAGCTGTTGC-3′
Glut2
(NM_0311979.2)		Forward: 5′-TGTTGGGGCCATCAACATGA-3′
Reverse: 5′-GGCGAATTTATCCAGCAGCAC-3′
Gsk3β
(NM_001347232.1)		Forward: 5′-TGGACAAAGGACTCACCAGG-3′
Reverse: 5′-AAGAGTGCAGGTGTGTCTCG-3′
Pepck
(NM_011044.3)		Forward: 5′-CGATGACATCGCCTGGATGA-3′
Reverse: 5′-TCTTGCCCTTGTGTTCTGCA-3′
Per1
(NM_0011065.3)		Forward: 5′-AGCAGGTTCAGGCTAACCAGGAAT-3′
Reverse: 5′-AGGTGTCCTGGTTTCGAAGTGTGT-3′
Probe:5’/5FAM/AGCTTGTGCCATGGACTGTCTACT/BHQ_1/-3′
Pparα
(NM_011144.6)		Forward: 5′-ACGTTTGTGGCTGGTCAAGT-3′
Reverse: 5′-TGGAGAGAGGGTGTCTGTGAT-3′
Pparγ
(NM_0011273302)		Forward: 5′-TGTGGGGATAAAGCATCAGGC-3′
Reverse: 5′-CCGGCAGTTAAGATCACACCTAT-3′
Trpv1
(NM_1001445.1)		Forward: 5′-CAGAGACCTGTGTCGGTTTATG-3′
Reverse: 5′-CATGTTGAGCAGGAGGATGTAG-3′
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lacerda, J.T.; Gomes, P.R.L.; Zanetti, G.; Mezzalira, N.; Lima, O.G.; de Assis, L.V.M.; Guler, A.; Castrucci, A.M.; Moraes, M.N. Lack of TRPV1 Channel Modulates Mouse Gene Expression and Liver Proteome with Glucose Metabolism Changes. Int. J. Mol. Sci. 2022, 23, 7014. https://doi.org/10.3390/ijms23137014

AMA Style

Lacerda JT, Gomes PRL, Zanetti G, Mezzalira N, Lima OG, de Assis LVM, Guler A, Castrucci AM, Moraes MN. Lack of TRPV1 Channel Modulates Mouse Gene Expression and Liver Proteome with Glucose Metabolism Changes. International Journal of Molecular Sciences. 2022; 23(13):7014. https://doi.org/10.3390/ijms23137014

Chicago/Turabian Style

Lacerda, José Thalles, Patrícia R. L. Gomes, Giovanna Zanetti, Nathana Mezzalira, Otoniel G. Lima, Leonardo V. M. de Assis, Ali Guler, Ana Maria Castrucci, and Maria Nathália Moraes. 2022. "Lack of TRPV1 Channel Modulates Mouse Gene Expression and Liver Proteome with Glucose Metabolism Changes" International Journal of Molecular Sciences 23, no. 13: 7014. https://doi.org/10.3390/ijms23137014

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop