LncRNA Snhg3 aggravates hepatic steatosis via PPARγ signaling
eLife assessment
This study provides useful evidence substantiating a role for long noncoding RNAs in liver metabolism and organismal physiology. Using murine knockout and knock-in models, the authors invoke a previously unidentified role for the lncRNA Snhg3 in fatty liver. The revised manuscript has improved and most studies are backed by solid evidence but the study was found to be incomplete and will require future studies to substantiate some of the claims.
https://doi.org/10.7554/eLife.96988.4.sa0Useful: Findings that have focused importance and scope
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Solid: Methods, data and analyses broadly support the claims with only minor weaknesses
Incomplete: Main claims are only partially supported
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
LncRNAs are involved in modulating the individual risk and the severity of progression in metabolic dysfunction-associated fatty liver disease (MASLD), but their precise roles remain largely unknown. This study aimed to investigate the role of lncRNA Snhg3 in the development and progression of MASLD, along with the underlying mechanisms. The result showed that Snhg3 was significantly downregulated in the liver of high-fat diet-induced obesity (DIO) mice. Notably, palmitic acid promoted the expression of Snhg3 and overexpression of Snhg3 increased lipid accumulation in primary hepatocytes. Furthermore, hepatocyte-specific Snhg3 deficiency decreased body and liver weight, alleviated hepatic steatosis and promoted hepatic fatty acid metabolism in DIO mice, whereas overexpression induced the opposite effect. Mechanistically, Snhg3 promoted the expression, stability and nuclear localization of SND1 protein via interacting with SND1, thereby inducing K63-linked ubiquitination modification of SND1. Moreover, Snhg3 decreased the H3K27me3 level and induced SND1-mediated chromatin loose remodeling, thus reducing H3K27me3 enrichment at the Pparg promoter and enhancing PPARγ expression. The administration of PPARγ antagonist T0070907 improved Snhg3-aggravated hepatic steatosis. Our study revealed a new signaling pathway, Snhg3/SND1/H3K27me3/PPARγ, responsible for mice MASLD and indicates that lncRNA-mediated epigenetic modification has a crucial role in the pathology of MASLD.
Introduction
Non-alcohol fatty liver disease (NAFLD) is characterized by excess liver fat in the absence of significant alcohol consumption. It can progress from simple steatosis to nonalcoholic steatohepatitis (NASH) and fibrosis and eventually to chronic progressive diseases such as cirrhosis, end-stage liver failure, and hepatocellular carcinoma (Loomba et al., 2021). In 2020, an international panel of experts led a consensus-driven process to develop a more appropriate term for the disease utilizing a two-stage Delphi consensus, that is, ‘metabolic dysfunction-associated fatty liver disease (MASLD)’ related to systemic metabolic dysregulation (Gofton et al., 2023; Rinella et al., 2023). The pathogenesis of MASLD has not been entirely elucidated. Multifarious factors such as genetic and epigenetic factors, nutritional factors, insulin resistance, lipotoxicity, microbiome, fibrogenesis and hormones secreted from the adipose tissue, are recognized to be involved in the development and progression of MASLD (Buzzetti et al., 2016; Friedman et al., 2018; Lee et al., 2017; Rada et al., 2020; Sakurai et al., 2021). Free fatty acids (FFAs), which are central to the pathogenesis of MASLD, originate from the periphery, mainly via lipolysis of triglyceride in the adipose tissue, or from increased hepatic de novo lipogenesis (DNL). Fatty acids in hepatocytes undergo mitochondrial β-oxidation and re-esterification to form triglyceride (TG), which are then exported into the blood as very low-density lipoproteins or stored in lipid droplets. Hepatic lipotoxicity occurs when the disposal of fatty acids through β-oxidation or the formation of TG is overwhelmed, which leads to endoplasmic reticulum (ER) stress, oxidative stress and inflammasome activation (Friedman et al., 2018). A cluster of differentiation 36/fatty acid translocase (CD36) and cell death-inducing DFF45-like effector proteins A/C (CIDEA/C) are critical for MASLD progression (Koonen et al., 2007; Matsusue et al., 2008; Sans et al., 2019). CD36 can increase FFAs uptake in the liver and drive hepatosteatosis onset. Overexpression of CD36 in hepatocytes increased FFAs uptake and TG storage; conversely, its deletion ameliorated hepatic steatosis and insulin resistance in DIO mice (Rada et al., 2020). Additionally, CIDEA/C can also regulate various aspects of lipid homeostasis, including lipid storage, lipolysis, and lipid secretion (Xu et al., 2024). As a transcription regulator of Cd36 and Cidea/c, peroxisome proliferator-activated receptor gamma (PPARγ) plays a crucial role in MASLD progression (Lee et al., 2023b; Lee et al., 2021; Lee et al., 2018; Matsusue et al., 2008; Puri et al., 2008; Skat-Rørdam et al., 2019).
Epigenetics, an inheritable phenomenon occurring without altering the DNA sequence, can regulate gene expression through different forms, including DNA methylation, histone modifications, chromatin remodeling, transcriptional control, and non-coding RNAs (Mann, 2014). Histone modifications, including acetylation, methylation, phosphorylation, ubiquitination, ribosylation, and ubiquitin-like protein modification (SUMO), are important epigenetic determinants of chromatin tightness and accessibility (Chen and Pikaard, 1997). Histone methylation is associated with chromatin-specific transcriptional activity states; for example, methylation of lysine 4 of histone H3 (H3K4), H3K36 and H3K79 are linked with a transcriptional activation state, and H3K9, H3K27, and H4K20 with transcriptional repression state (Pirola and Sookoian, 2022). Previous studies have illustrated that epigenetics factors including histone modification play key role in lipid metabolism (Bayoumi et al., 2020; Byun et al., 2020; Jun et al., 2012).
Long non-coding RNAs (lncRNAs) are non-coding RNAs with more than 200 bases in length, can be transcribed by RNA polymerase II, and are comparable to mRNA but lack the crucial open reading framework required for translation (Ng et al., 2013; Ulitsky and Bartel, 2013). LncRNAs are involved in epigenetic regulation of gene expression at different levels and through different molecular mechanisms such as chromatin remodeling, transcriptional regulation and post-transcriptional processing. Previous studies have indicated that lncRNAs are involved in the pathological progress of MASLD (Bayoumi et al., 2020; Sommerauer and Kutter, 2022). Although histone modification and lncRNAs influence the susceptibility to MASLD, their roles in MASLD remain largely unknown.
Small nucleolar RNA host genes (SNHG) family, a type of lncRNA, serve as host genes for producing intronic small nucleolar RNAs (snoRNAs) and are mainly related to tumor pathophysiology by regulating proliferation, apoptosis, invasion, and migration (Sen et al., 2020; Zimta et al., 2020). The family of mouse Snhg genes has 19 members including Snhg1-18, Snhg20 and Snhg7os. Here, we found that the expression of hepatic Snhg3 was decreased in high-fat diet (HFD)-induced obesity (DIO) mice. Experiments conducted using in vivo and in vitro models indicated that Snhg3 was involved in fatty acid metabolism and hepatic steatosis. Mechanistically, Snhg3 interacted with staphylococcal nuclease and Tudor domain containing 1 (SND1), a well-understood Tudor protein that participates in lipid metabolism and tumoral behavior by modulating cholesterol and glycerophospholipid metabolism and acylglyceride storage in lipid droplets (Navarro-Imaz et al., 2020). Furthermore, Snhg3 increased the expression of SND1 by promoting the stability of SND1 mediated by K63-linked ubiquitination and induced nuclear localization of SND1 protein, thereby reducing tri-methylation at H3K27 (H3K27me3) enrichment and boosting chromatin loose remodeling at Pparg promoter, eventually enhancing Pparg, Cd36 and Cidea/c expressions. Our result indicated that SND1/H3K27me3/PPARγ is partially responsible for Sngh3-induced hepatic steatosis.
Results
LncRNA-Snhg3 is downregulated in DIO mice
Firstly, we analyzed the lncRNAs expression profiles in the livers of DIO mice and normal chow-fed mice (control) by RNA-Seq, and found 18072 hepatic lncRNAs, including 338 differentially expressed lncRNAs (q-value ≤0.05, Figure 1A). Of all Snhgs, Snhg3 had the most prominent expression and exhibited more noticeable downregulation in the liver of the DIO mice compared to the control mice (Figure 1B), thus, it was selected for further study. The downregulation of Snhg3 was confirmed by RT-qPCR (Figure 1C). Additionally, the Coding Potential Calculator indicated that Snhg3 has a coding probability of 0.020757, classifying it as a noncoding sequence (Kang et al., 2017). Localization of Snhg3 was primarily observed in nuclei with a probability score of 0.451138, as predicted using software prediction (http://lin-group.cn/server/iLoc-LncRNA/predictor.php). The exact nuclear localization of Snhg3 was further confirmed by nuclear/cytoplasmic fractionation (Figure 1D). Interestingly, the expression of Snhg3 was induced by palmitic acid (PA) in primary hepatocytes (Figure 1E). Furthermore, overexpression of Snhg3 increased lipid accumulation in primary hepatocytes with PA treatment (Figure 1F and G).
Hepatocyte-specific Snhg3 knock-out alleviates hepatic steatosis in DIO mice
Given Snhg3 was associated with hepatic nutrition change, the role of Snhg3 was further confirmed by constructing hepatocyte-specific Snhg3 knock-out (Snhg3-HKO) mice that were then induced obesity with a high-fat diet (Figure 2A and Figure 2—figure supplement 1A, B). The result indicated that body weight was mildly decreased in Snhg3-HKO mice compared with the control Snhg3flox/flox (Snhg3-Flox) mice (Figure 2B). The energy consumption is mainly reflected as the sum of energy utilization during internal heat production using comprehensive laboratory animal monitoring system (CLAMS). Heat production showed an increasing trend but was not statistically significant in Snhg3-HKO mice (Figure 2—figure supplement 1C). Moreover, there were no obvious differences in total oxygen consumption, carbon dioxide production or respiratory exchange ratio (RER) between Snhg3-HKO and control mice (Figure 2—figure supplement 1C). Furthermore, insulin sensitivity, not glucose tolerance, was improved in Snhg3-HKO mice (Figure 2C). The Snhg3-HKO mice had a decrease in liver weight and the ratio of liver weight/body weight, and improved hepatic steatosis, including decreasing lipid accumulations and the ballooning degeneration of liver cells (Figure 2D–F). However, the hepatic fibrosis phenotype showed no difference (Figure 2—figure supplement 1D). Serum alanine transaminase (ALT) and aspartate transaminase (AST) levels were significantly decreased in Snhg3-HKO mice (Figure 2G). Moreover, serum FFAs, TG and TC were also reduced in Snhg3-HKO mice (Figure 2H). The Snhg3-HKO mice exhibited a decrease in inguinal white adipose tissue (iWAT) weight and weight/body weight ratio, while brown adipose (BAT) weight and weight/body weight ratio remained unaltered (Figure 2—figure supplement 1E). Additionally, there was no difference in serum insulin between Snhg3-HKO mice and control mice (Figure 2—figure supplement 1F). These results suggested that hepatic knockout of Snhg3 improves hepatic steatosis in mice.
Hepatocyte-specific Snhg3 knock-in aggravates hepatic steatosis in DIO mice
Furthermore, the hepatocyte-specific Snhg3 knock-in (Snhg3-HKI) mice were also constructed and subsequently induced obesity with a high-fat diet to detect the function of Snhg3 in the liver (Figure 3A and Figure 3—figure supplement 1A). The Snhg3-HKI mice showed greater weight gains than the control wild type (WT) mice (Figure 3B). Insulin sensitivity was also impaired in Snhg3-HKI mice (Figure 3C). The liver weight and the ratio of liver weight/body weight of Snhg3-HKI mice were markedly increased (Figure 3D). Also, Snhg3-HKI mice exhibited severe hepatic steatosis (Figure 3E and F) and higher serum ALT and AST levels (Figure 3G). Both serum TC and iWAT weight were increased in Snhg3-HKI mice (Figure 3H and Figure 3—figure supplement 1B). Similar to Snhg3-HKO mice, there was also no differences in heat production, total oxygen consumption, carbon dioxide production, RER, hepatic fibrosis phenotype, and serum insulin between Snhg3-HKI mice and WT mice (Figure 3—figure supplement 1C–E). These findings indicated that upregulation of Snhg3 could promote hepatic steatosis.
Snhg3 promotes hepatic steatosis by regulating chromatin remodeling
To clarify the molecular mechanism of Snhg3 in hepatic steatosis, we investigated the hepatic differentially expressed genes (DEGs) using RNA-Seq. There were 1393 DEGs between the Snhg3-HKI and control WT mice, with 1028 genes being upregulated and 365 genes downregulated (log2FC ≥1, q-value <0.001) in the liver of Snhg3-HKI mice (Figure 4A). A gene set enrichment analysis (GSEA) of DEGs revealed that Snhg3 exerts a global effect on the expression of genes involved in fatty acid metabolism and the PPAR signaling pathway (Figure 4B). RT-qPCR analysis confirmed that the hepatic expression levels of some genes involved in fatty acid metabolism, including Cd36, Cidea/c, and stearoyl-CoA desaturase (Scd1/2), the key enzymes involved in the biosynthesis of unsaturated fatty acids (Ntambi and Miyazaki, 2003), were upregulated in Snhg3-HKO mice and were downregulated in Snhg3-HKI mice compared to the controls (Figure 4C). Additionally, deficiency and overexpression of Snhg3 respectively decreased and increased the expression of profibrotic genes, such as collagen type I alpha 1/2 (Col1a1 and Col1a2), but had no effects on the pro-inflammatory factors, including transforming growth factor β1 (Tgfb1), tumor necrosis factor a (Tnfa), interleukin 6 and 1b (Il6 and Il1b; Figure 4—figure supplement 1A, B). LncRNAs in the nucleus can affect gene expression in multiple ways, such as chromatin remodeling, transcriptional regulation, and post-transcriptional processing (Morey and Avner, 2004; Thomson and Dinger, 2016). Since Snhg3 was mainly localized in the nuclei of hepatocytes, we next checked the genome-wide chromatin accessibility (log2FC >2, p-value <0.001) in the liver of Snhg3-HKI and WT mice using ATAC-Seq. We discovered that in all 6810 differentially accessible regions (DARs), 4305 (>63.2%) were more accessible in Snhg3-HKI mice and only 2505 (>36.8%) of peaks were more accessible in control mice, indicating that the chromatin states were ‘hyper-accessible’ in the liver of Snhg3-HKI mice. Moreover, DARs were with relatively few promoter-proximal (Up2k) and exon regions in both the control and Snhg3-HKI groups (Figure 4D), supporting the idea that gene activation depends on multiple regulatory regions, is not limited to its promoter and exon regions (Ackermann et al., 2016). Furthermore, 3966 genes were associated specifically with the accessible regions in the Snhg3-HKI group and only 2451 genes in the WT group (log2FC >2, p-value <0.001). Additionally, PPARg was identified as a potential transcription factor associated with hyper-accessible regions in the liver of the Snhg3-HKI group by HOMER and CREMA (Figure 4E and F).
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig4-data1-v1.xls
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig4-data2-v1.xls
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig4-data3-v1.xls
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig4-data4-v1.xls
To determine whether open chromatin regions were correlated with gene expression, we integrated ATAC-Seq data (genes associated with DARs, log2FC >2, p-value <0.001) with RNA-Seq data (DEGs in DIO Snhg3-HKI and control mice, log2FC >1, q-value <0.001). Overall, 233 upregulated genes shown in quadrant 2, including Cd36 and Cidea/c, had at least one associated open chromatin region, which accounted for >22.67% (total 1028) of DEGs mapped to ATAC-Seq peaks in the liver of Snhg3-HKI mice (Figure 4G and H). Meanwhile, at least one open chromatin region was associated with 65 downregulated genes in the quadrant 3, which accounted for >17.81% (total 365) of DEGs mapped to ATAC-Seq peaks in the liver of WT mice (Figure 4G).
Snhg3 induces SND1 expression by interacting with SND1 and enhancing the stability of SND1 protein
To further elucidate the molecular mechanism of Snhg3 in hepatic steatosis, an RNA pull-down followed by mass spectrometry (MS) assay was performed in primary hepatocytes. The result identified 234 specific Snhg3-associated proteins, involved in multiple signaling pathways, including PPAR, NAFLD and fatty acid degradation pathways (Figure 5A and B). Of these proteins, a well-understood Tudor protein SND1 was also predicted to interact with three fragments of Snhg3 by bioinformatic method (RBPsuite; Figure 5C and D). Snhg3 coprecipitation with SND1 was confirmed by RNA pull-down coupled with western blotting (Figure 5E), which was consistent with the RNA immunoprecipitation (RIP) assay results (Figure 5F). Meanwhile, Snhg3 regulated the protein, not mRNA, expression of SND1 in vivo and in vitro by mildly promoting the stability of SND1 protein (Figure 5G–J). Furthermore, we tested the effect of Snhg3 on the ubiquitin-modification of SND1 and found that Snhg3 enhanced SND1 ubiquitination in vivo and in vitro (Figure 5K and L). Previous studies indicated that K48-linked polyubiquitination aids in proteasome-mediated recognition and degradation and that K63-linked polyubiquitination participates in signaling assemblies and protein stability (Sun et al., 2020). As predicted, Snhg3 overexpression increased K63-linked ubiquitination modification in endogenous and exogenous SND1 protein, not K48- or K33-linked (Figure 5M and N). Additionally, Snhg3 overexpression enhanced the nuclear localization of SND1 in Hepa1-6 cells with PA treatment (Figure 5O). Collectively, these results suggested that Snhg3 promoted the K63-linked ubiquitination and stability of SND1 protein through interacting with SND1, thus resulting in SND1 protein increase and nuclear localization.
-
Figure 5—source data 1
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data1-v1.xls
-
Figure 5—source data 2
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data2-v1.xls
-
Figure 5—source data 3
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data3-v1.pdf
-
Figure 5—source data 4
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data4-v1.zip
-
Figure 5—source data 5
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data5-v1.pdf
-
Figure 5—source data 6
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data6-v1.zip
-
Figure 5—source data 7
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data7-v1.pdf
-
Figure 5—source data 8
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data8-v1.zip
-
Figure 5—source data 9
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data9-v1.pdf
-
Figure 5—source data 10
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data10-v1.zip
-
Figure 5—source data 11
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data11-v1.pdf
-
Figure 5—source data 12
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data12-v1.zip
-
Figure 5—source data 13
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data13-v1.pdf
-
Figure 5—source data 14
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data14-v1.zip
-
Figure 5—source data 15
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data15-v1.pdf
-
Figure 5—source data 16
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data16-v1.zip
-
Figure 5—source data 17
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data17-v1.pdf
-
Figure 5—source data 18
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data18-v1.zip
-
Figure 5—source data 19
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data19-v1.pdf
-
Figure 5—source data 20
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig5-data20-v1.zip
Snhg3 promotes PPARγ expression by decreasing H3K27me3 enrichment at the Pparg promoter
SND1, initially named as p100, is a highly conserved and ubiquitously expressed multifunctional Tudor domain-containing protein that participates in pivotal biological processes like double-stranded RNA editing, pre-mRNA splicing, microRNA-mediating gene silencing and piRNA biogenesis in germlines (Ying and Chen, 2012). Previous studies indicated that Tudor proteins participate in epigenetic regulation by binding to methyl-arginine⁄lysine residues (Ying and Chen, 2012). However, whether SND1 influences histone modification remains unclear. It is well known that histone modification dynamically regulates specific gene expression by altering the organization and function of chromatin and is involved in the pathophysiology of some diseases, such as histone H3 methylation modification, which may contribute to MASLD pathogenesis (Byun et al., 2020; Jun et al., 2012; Tessarz and Kouzarides, 2014). Considering that H3K27me3, a repressive chromatin mark, plays a role in autophagy-mediated lipid degradation (Byun et al., 2020), we tested the effect of SND1 on H3K27me3. The results revealed that both SND1 and Snhg3 overexpression reduced the H3K27me3 level (Figure 6A). Furthermore, disrupting SND1 expression increased the H3K27me3 level and reversed the Snhg3-induced H3K27me3 decrease (Figure 6B and C). Moreover, the hepatic H3K27me3 level was upregulated in Snhg3-HKO mice but downregulated in Snhg3-HKI mice (Figure 6D). The results indicated that Snhg3 negatively regulated the H3K27me3 level through SND1.
-
Figure 6—source data 1
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data1-v1.pdf
-
Figure 6—source data 2
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data2-v1.zip
-
Figure 6—source data 3
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data3-v1.pdf
-
Figure 6—source data 4
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data4-v1.zip
-
Figure 6—source data 5
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data5-v1.pdf
-
Figure 6—source data 6
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data6-v1.zip
-
Figure 6—source data 7
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data7-v1.pdf
-
Figure 6—source data 8
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data8-v1.zip
-
Figure 6—source data 9
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig6-data9-v1.xls
To further clarify whether Snhg3-induced H3K27me3 decrease is involved in hepatic steatosis, we examined the H3K27me3 enrichment in the liver of Snhg3-HKO mice using the CUT&Tag assay and detected 10915 peaks. The genomic locations of these peaks were divided into eight categories, and the H3K27me3 signals were predominantly enriched (about 54%) at the 2 kb promoter, 5’-untranslated region (5’-UTR), and exon categories. Meanwhile, very few signals (about 14%) were enriched in the 2 kb downstream and intergenic categories in the liver of Snhg3-HKO mice (Figure 6E). Moreover, the exon, upstream 2 k, 5’-UTR and intron regions of Pparg were enriched with the H3K27me3 mark (fold_enrichment = 4.15697) in the liver of Snhg3-HKO mice. Subsequently, ChIP assay revealed that hepatic H3K27me3 enrichment at the Pparg promoter was increased in Snhg3-HKO mice but decreased in Snhg3-HKI mice (Figure 6F). Snhg3-overexpression in Hepa1-6 cells yielded similar results (Figure 6G).
SND1 mediates Snhg3-induced PPARγ upregulation
PPARγ has been reported to influence MASLD progression by regulating target genes such as Cd36 and Cidea/c (Lee et al., 2023b; Lee et al., 2021; Lee et al., 2018; Matsusue et al., 2008; Skat-Rørdam et al., 2019). In this study, the mRNA and protein expression levels of hepatic PPARγ were decreased in Snhg3-HKO mice and increased in Snhg3-HKI mice (Figure 7A–C). Additionally, CD36 protein level was decreased in Snhg3-HKO mice and increased in Snhg3-HKI mice (Figure 7B and C). The upregulation of Snhg3 and SND1 also increased the expression of Pparg and Cd36 in vitro (Figure 7D–F). Meanwhile, disruption of SND1 expression alleviated Snhg3-induced PPARγ increase and lipid accumulation (Figure 7G–I). Collectively, these results demonstrated that SND1 mediated Snhg3-induced PPARγ and CD36 expression.
-
Figure 7—source data 1
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data1-v1.pdf
-
Figure 7—source data 2
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data2-v1.zip
-
Figure 7—source data 3
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data3-v1.pdf
-
Figure 7—source data 4
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data4-v1.zip
-
Figure 7—source data 5
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data5-v1.pdf
-
Figure 7—source data 6
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data6-v1.zip
-
Figure 7—source data 7
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data7-v1.pdf
-
Figure 7—source data 8
- https://cdn.elifesciences.org/articles/96988/elife-96988-fig7-data8-v1.zip
In addition, Snhg3 serves as host gene for producing intronic U17 snoRNAs, the H/ACA snoRNA. A previous study found that cholesterol trafficking phenotype was not due to reduced Snhg3 expression, but rather to haploinsufficiency of U17 snoRNA. Upregulation of hypoxia-upregulated mitochondrial movement regulator (HUMMR) in U17 snoRNA-deficient cells promoted the formation of ER-mitochondrial contacts, resulting in decreasing cholesterol esterification and facilitating cholesterol trafficking to mitochondria (Jinn et al., 2015). Additionally, disruption of U17 snoRNA caused resistance to lipid-induced cell death and general oxidative stress in cultured cells. Furthermore, knockdown of U17 snoRNA in vivo protected against hepatic steatosis and lipid-induced oxidative stress and inflammation (Sletten et al., 2021). In this study, the expression of U17 snoRNA decreased in the liver of Snhg3-HKO mice and unchanged in the liver of Snhg3-HKI mice, but overexpression of U17 snoRNA had no effect on the expression of SND1 and PPARγ (Figure 7—figure supplement 1A–C), indicating that Sngh3 induced hepatic steatosis was independent on U17 snoRNA.
PPARγ mediates Snhg3-induced hepatic steatosis
Hepatocyte-specific depletion of PPARγ is known to protect mice against NASH and boost the therapeutic efficacy of rosiglitazone, a synthetic PPARγ agonist, in the liver (Lee et al., 2021). Furthermore, PPARγ is an inducer of adipocyte differentiation and a reservoir for excess FFAs, thereby potentially preventing lipotoxicity in other tissues and organs (Medina-Gomez et al., 2007). To this end, we tested the effect of T0070907, a selective PPARγ antagonist, on Snhg3-induced hepatic steatosis in mice. The result showed that T0070907 treatment for 8 weeks had no effects on body weight, liver and iWAT weight, and serum FFAs, TG and TC in Snhg3-HKI mice, but improved Snhg3-induced hepatic steatosis in Snhg3-HKI mice (Figure 8A–D and Figure 8—figure supplement 1). Moreover, T0070907 mitigated the hepatic Cd36 and Cidea/c increase in Snhg3-HKI mice (Figure 8E). Additionally, Snhg3- and SND1-induced Cd36 increase also were abolished by T0070907 in hepa1-6 cells (Figure 8F). Collectively, these results suggested that PPARγ-mediated Snhg3-induced hepatic steatosis.
Discussion
Liver steatosis is common in various metabolic diseases and related disorders, including MASLD. Although lncRNAs are implicated in regulating numerous mechanisms related to liver steatosis and MASLD, their exact function remains to be determined. In this study, lncRNA-Snhg3 is downregulated in DIO mice and hepatocyte-specific Snhg3 deficiency improved hepatic steatosis and insulin resistance, while overexpression aggravated hepatic steatosis and insulin resistance in mice. Our results showed that the expression of Snhg3 was decreased in DIO mice which led us to speculate that the downregulation of Snhg3 might be a stress protective reaction to high nutritional state, but the specific details need to be clarified. This is probably similar to fibroblast growth factor 21 (FGF21) and growth differentiation factor 15 (GDF15), whose endogenous expression and circulating levels are elevated in obese humans and mice despite their beneficial effects on obesity and related metabolic complications (Keipert and Ost, 2021). Although FGF21 can be induced by oxidative stress and be activated in obese mice and in NASH patients, elevated FGF21 paradoxically protects against oxidative stress and reduces hepatic steatosis (Tillman and Rolph, 2020).
Excessive hepatic lipid deposition owing to increased FFAs uptake and hepatic DNL impairs autophagy and promotes ER stress and oxidative stress, insulin resistance, inflammation, and liver tissue damage, ultimately aggravating MASLD progression (Rada et al., 2020). In this study, Snhg3 induced the expression of fatty acid metabolism related genes such as Cd36, Cidea/c and Scd1/2. Under physiological conditions, CD36 expression in hepatocytes was found to be minimal; however, lipid overload or activation of nuclear receptors including PPARα/γ and liver X receptor (LXR), could significantly increase it (Rada et al., 2020). As a transcription regulator of Cd36 and Cidea/c, it is well known that PPARγ plays major adipogenic and lipogenic roles in adipose tissue. Although the expression of PPARγ in the liver is very low under healthy conditions, induced expression of PPARγ in both hepatocytes and non-parenchymal cells (Kupffer cells, immune cells, and hepatic stellate cells [HSCs]) in the liver has a crucial role in the pathophysiology of MASLD (Chen et al., 2023; Gross et al., 2017; Lee et al., 2023b). The activation of PPARγ in the liver induces the adipogenic program to store fatty acids in lipid droplets as observed in adipocytes (Lee et al., 2018). Moreover, the inactivation of liver PPARγ abolished rosiglitazone-induced an increase in hepatic TG and improved hepatic steatosis in lipoatrophic AZIP mice (Gavrilova et al., 2003). Furthermore, there is a strong correlation between the onset of hepatic steatosis and hepatocyte-specific PPARγ expression. Clinical trials have also indicated that increased insulin resistance and hepatic PPARγ expressions were associated with NASH scores in some obese patients (Lee et al., 2023a; Mukherjee et al., 2022). Even though PPARγ’s primary function is in adipose tissue, patients with MASLD have much higher hepatic expression levels of PPARγ, reflecting the fact that PPARγ plays different roles in different tissues and cell types (Mukherjee et al., 2022). As these studies mentioned above, our result also hinted at the importance of PPARγ in the pathophysiology of MASLD. Snhg3 deficiency or overexpression respectively induced the decrease or increase in hepatic PPARγ. Moreover, administration of PPARγ antagonist T0070907 mitigated the hepatic Cd36 and Cidea/c increase and improved Snhg3-induced hepatic steatosis. However, conflicting findings suggest that the expression of hepatic PPARγ is not increased as steatosis develops in humans and in clinical studies and that PPARγ agonists administration did not aggravate liver steatosis (Gross et al., 2017). Thus, understanding how the hepatic PPARγ expression is regulated may provide a new avenue to prevent and treat the MASLD (Lee et al., 2018).
Hepatotoxicity accelerates the development of progressive inflammation, oxidative stress, and fibrosis (Roehlen et al., 2020). Chronic liver injury including MASLD can progress to liver fibrosis with the formation of a fibrous scar. Injured hepatocytes can secrete fibrogenic factors or exosomes containing miRNAs that activate HSCs, the major source of the fibrous scar in liver fibrosis (Kisseleva and Brenner, 2021). Apart from promoting lipogenesis, PPARγ has also a crucial function in improving inflammation and fibrosis (Chen et al., 2023). In this study, no hepatic fibrosis phenotype was seen in Snhg3-HKO and Snhg3-HKI mice. Moreover, the expression levels of profibrotic genes including Col1a1 and Col1a2 were decreased in Snhg3-HKO mice and increased in Snhg3-HKI mice, but the pro-inflammatory factors including Tgfb1, Tnfa, Il6, and Il1b had no changes. Inflammation is an absolute requirement for fibrosis because factors from injured hepatocytes alone are not sufficient to directly activate HSCs and lead to fibrosis (Kisseleva and Brenner, 2021). Additionally, previous studies indicated that exposure to HFD for more 24 weeks causes less severe fibrosis (Alshawsh et al., 2022). In future, the effect of Snhg3 on hepatic fibrosis in mice need to be elucidated by prolonged high-fat diet feeding or adopting methionine- and choline deficient diet (MCD) feeding.
Epigenetics plays a crucial role in many physiological and pathological situations (Peixoto et al., 2020). Epigenetic regulation induces phenotypic changes that may respond to environmental cues through DNA methylation and histone modification, chromatin remodeling, and noncoding RNAs (Mann, 2014). Epigenetic changes interact with inherited risk factors to modulate the individual risk of MASLD development and the severity of progression. Epigenetic modifications, including DNA methylation, miRNAs, and histone modifications, have been associated with MASLD (Baffy, 2015; Eslam et al., 2018; Jonas and Schürmann, 2021). To date, there is no approved pharmacologic therapy for MASLD, and the mainstay of management remains lifestyle changes with exercise and dietary modifications (Bayoumi et al., 2020). Therefore, understanding the epigenetic modifications in MASLD pathogenesis might prove a rational strategy to prevent the disease and develop novel therapeutic interventions (Sodum et al., 2021).
LncRNAs, being abundant in the genome participate in regulating the expression of coding genes through various molecular mechanisms, including: (1) transcriptional regulation at the promoter of target genes; (2) inhibiting RNA polymerase II or mediating chromatin remodeling and histone modification; (3) interfering with the splicing and processing of mRNA or producing endogenous siRNA; (4) regulating the activity or cellular localization of the target protein; (5) acting as competitive endogenous RNAs; and (6) riboregulation by forming nucleic acid-protein complex as structural component (Morey and Avner, 2004; Sommerauer and Kutter, 2022; Thomson and Dinger, 2016). However, compared to the large number of lncRNAs, only few have been functionally well-characterized. Collective literature has shown that lncRNAs play a crucial role in MASLD (Sommerauer and Kutter, 2022). This study demonstrated that lncRNA-Snhg3 participated in the pathology of MASLD by epigenetic modification; that is, Snhg3 inhibited the H3K27me3 level, and promoted chromatin relaxation at the Pparg promoter and eventually increased PPARγ expression. The results from Ruan et al. demonstrated that more than a third of dynamically expressed lncRNAs were deregulated in a human MASLD cohort and the lncRNA human lncRNA metabolic regulator 1 (hLMR1) positively regulated transcription of genes involved in cholesterol metabolism (Ruan et al., 2021). Previous studies have also demonstrated that several lncRNAs, including FLRL2/3/6/7/8, H19, and MALAT-1, were associated with lipogenesis via proteins in the PPAR signaling pathway (Mukherjee et al., 2022). Recently, a murine long noncoding single-cell transcriptome analysis elucidated liver lncRNA cell-type specificities, spatial zonation patterns, associated regulatory networks, and temporal patterns of dysregulation during hepatic disease progression. Moreover, a subset of the liver disease-associated regulatory lncRNAs identified have human orthologs (Karri and Waxman, 2023). Based on the aforementioned information, lncRNAs emerge as promising candidates for biomarkers and therapeutic targets for MASLD.
Tudor proteins play vital roles in normal cell viability and growth by diverse epigenetic functions, including methylation dependent chromatin-remodeling, histone-binding, pre-RNA-processing, RNA-silencing, and transposon silencing in ligands (Ying and Chen, 2012). Tudor proteins are divided into four groups: Group1 Tudor proteins bind to methyl-lysine⁄arginine of histone tails, including Tdrd3, PHF1, PHF20, the Jumonji domain-containing protein (JMJD) family and TP53BP1; Group 2 Tudor proteins bind to methyl-arginine of ligands and representative members include SMN and SMNDC1; Group 3 is represented by SND1; and Group 4, contains many Tudor proteins, including Tdrd1-9 and Tdrd11, that have been identified in methylation-dependent association with PIWI proteins Ago3, Aub, and Piwi. In this study, Snhg3 induced the protein level of SND1 by promoting K63-linked ubiquitination of SND1 and increasing its protein stability. Additionally, Snhg3-induced SND1 protein stability seemed subtle, indicating there may be other way for Snhg3 promotion SND1, such as riboregulation. Some studies suggested that SND1 plays important roles in cancer by interacting with other transcription factors, including PPARγ, signal transducer and activator of transcription 5/6 (STAT6/5) and myeloblastosis oncogene (c-Myb) (Duan et al., 2014; Navarro-Imaz et al., 2020). SND1 could induce adipogenesis and promote the formation of lipid droplets in adipocytes through working as a co-activator of PPARγ and regulating H3 acetylation (Duan et al., 2014). Our study showed that both Snhg3 and SND1 decreased the H3K27me3 level and promoted the expression of PPARγ. SND1 could interact with Snhg3 and mediate the Snhg3-induced decrease in H3K27me3 and increase in PPARγ expression. Furthermore, inhibition of PPARγ with T0070907 alleviated Snhg3- and SND1-induced Cd36 and Cidea/c increase and improved Snhg3-aggravated hepatic steatosis. In lncRNA riboregulation, the actions of noncoding RNAs mostly rely on interactions with proteins, including canonical or noncanonical RNA-binding proteins (RBPs). Canonical RBPs, such as heterogeneous nuclear ribonucleoproteins (hnRNPs), polypyrimidine tract binding protein 1 (PTBP1) and human antigen R (HUR), are often involved in posttranscriptional regulation, including pre-mRNA processing, RNA stability, RNA decay, or nuclear export (Briata and Gherzi, 2020). Previous studies have demonstrated that some lncRNAs, for example LINC01018, MEG3, APOA4-AS, hLMR1, Blnc1, and LncARSR interact with canonical RBPs to govern the progression of MASLD (Sommerauer and Kutter, 2022).
In summary, our study demonstrates that lncRNA-Snhg3 influenced fatty acid metabolism and aggravated hepatic steatosis under DIO status. Furthermore, Snhg3 increased the expression, stability, and nuclear localization of SND1 protein by interacting with SND1, thus enhancing the expression of PPARγ via reducing H3K27me3 enrichment and boosting chromatin loose remodeling at the Pparg promoter, indicating that SND1/H3K27me3/PPARγ is partially responsible for Snhg3-induced hepatic steatosis. This study reveals a new signaling pathway, Snhg3/SND1/H3K27me3/PPARγ, responsible for hepatic steatosis and provides evidence of lncRNA-mediated epigenetics in the pathophysiology of MASLD (Figure 8G).
However, there are still some limitations to this study that require further investigation. Notably, the expression change of H3K27me3, a global repressive histone mark, may affect multiple downstream target genes, including Pparg; therefore, more target genes involved in MASLD need to be elucidated. Moreover, the precise mechanism by which SND1 regulates H3K27me3 is still unclear and hence requires further investigation. It is crucial to ascertain whether SND1 itself functions as a new demethylase or if it influences other demethylases, such as JMJD3, enhancer of zeste homolog 2 (EZH2), and ubiquitously transcribed tetratricopeptide repeat on chromosome X (UTX). SND1 has multiple roles through associating with different types of RNA molecules, including mRNA, miRNA, circRNA, dsRNA, and lncRNA. SND1 could bind negative-sense SARS-CoV-2 RNA and to promote viral RNA synthesis (Schmidt et al., 2023). SND1 is also involved in hypoxia by negatively regulating hypoxia‐related miRNAs (Saarikettu et al., 2023). Furthermore, a recent study revealed that lncRNA SNAI3-AS1 can competitively bind to SND1 and perturb the m6A-dependent recognition of Nrf2 mRNA 3'UTR by SND1, thereby reducing the mRNA stability of Nrf2 (Zheng et al., 2023). Huang et al. also reported that circMETTL9 can directly bind to and increase the expression of SND1 in astrocytes, leading to enhanced neuroinflammation (Huang et al., 2023). However, whether there is an independent-histone methylation role of SND1/lncRNA-Snhg3 involved in lipid metabolism in the liver needs to be further investigated.
Materials and methods
Animals and treatments
Request a detailed protocolC57BL/6 Snhg3flox/flox (Snhg3-Flox) mice and hepatocyte-specific Snhg3 knock-in (Snhg3-HKI) mice were created using the CRISPR-Cas9 system at Cyagen Biosciences. To engineer the targeting vector for Snhg3-Flox mice, the exon 3 of Snhg3 was selected as the conditional knockout region, and homology arms and the cKO region were generated by PCR using BAC clone as a template. Cas9 and gRNA were co-injected into fertilized eggs with a targeting vector for mice production. The obtained mice were identified by PCR followed by sequence analysis. Hepatocyte-specific Snhg3 knock-out (Snhg3-HKO) mice were generated by crossing Snhg3-Flox mice with C57BL/6-Alb-Cre mice. For Snhg3-HKI mice, the ‘Alb promoter-mouse Snhg3 cDNA-polyA’ cassette was inserted into an H11 locus (~0.7 kb 5' of the Eif4enif1 gene and ~4.5 kb 3' of the Drg1 gene), and homology arms were generated by PCR using BAC clone as template to engineer the targeting vector. Cas9 and gRNA were co-injected into fertilized eggs with targeting vector for mice production. The obtained mice were identified by PCR followed by sequence analysis. All mice were housed in the pathogen-free conditions (SPF) facility and maintained on a 12 hr light-dark cycle and a regular unrestricted diet. All mice were fed either a normal chow diet (9% fat; Lab Diet) or HFD (60% fat, Research Diets) for inducing obesity and libitum with free access to water. Unless otherwise noted, 6~8-week-old male mice were used for all experiments. 8-week-old mice fed on HFD were injected intraperitoneally (i.p.) with 1 mg/kg of T0070907 dissolved in DMSO for 5 days per week for 2 months. Liver tissue samples were analyzed by the High Fatty Sample Total Cholesterol (TC) Content Assay Kit (APPLYGEN, Cat#E1025-105) and the High Fatty Sample Triglyceride (TG) Content Assay Kit (APPLYGEN, Cat#E1026-105), respectively. Serum concentrations of ALT, AST, FFAs, TG and TC were determined using an automated Monarch device (Peking Union Medical College Hospital, Beijing, China). Serum insulin was detected using a mouse insulin ELISA kit (JM-02862M1, Beijing BioDee Biotechnology Co., Ltd.). All animal experiments were conducted under protocols approved by the Animal Research Committee of the Institute of Laboratory Animals, Institute of Basic Medical Sciences Chinese Academy of Medical Sciences & School of Basic Medicine Peking Union Medical College (ACUC-A01-2022-010).
Cell culture
Request a detailed protocolPrimary mouse hepatocytes were isolated from 8-week-old male C57BL/6 J mice and cultured in RPMI 1640 medium with 10% FBS as previously described (Matsumoto et al., 2002). Hepa1-6 cells (ATCC, Cat#CRL-1830) were cultured at 37 °C, 5% CO2 in DMEM (Gibco, Carlsbad, USA) supplemented with 10% FBS, and 1% penicillin-streptomycin. The species origin of Hepa1-6 cell was confirmed with PCR and the identity of Hepa1-6 cell was authenticated with STR profiling. Hepa1-6 cell line was checked negative for mycoplasma by PCR. After attachment, the cells were transfected with indicated plasmids or siSnd1 by Lipofectamine 3000 Transfection Kit (Invitrogen). Cells were treated with and without 0.25 mM PA for 12h - 24h before collection. The sequences of siSnd1 were seen key resources table.
Plasmid construction
Request a detailed protocolSnhg3 was amplified from liver cDNA and was then constructed into pcDNA3.1 using Kpn I and EcoR I. The primers were seen key resources table.
Real-time quantitative PCR (RT-qPCR)
Request a detailed protocolTotal RNA was extracted from mouse tissues using a Trizol-based method. Approximately 2 μg of total RNA was reverse-transcribed into a first-strand cDNA pool using reverse transcriptase and random primers, according to the manufacturer’s instructions. RT-qPCR was performed using SYBR Green PCR Master Mix (A6002, Promega) with the gene-specific primers (key resources table). All gene expression data were normalized to β-Actin expression levels.
Western blotting
Request a detailed protocolProtein was extracted from frozen tissue samples in cell lysis buffer. In total, protein was loaded onto a 10% SDS-polyacrylamide gel, and separated proteins were transferred to PVDF membranes. Western blot assays were performed using indicated specific antibodies (key resources table). The proteins were quantified by ImageJ software.
Coding potential prediction
Request a detailed protocolThe coding potential of Snhg3 was evaluated by the Coding Potential Calculator at CPC2@CBI,PKU(gao-lab.org) (Kang et al., 2017).
Histopathologic analysis
Request a detailed protocolLiver tissue sections were fixed in 4% paraformaldehyde, then embedded in paraffin and stained with H&E to visualize the general morphological and structural characteristics of tissues. Lipid droplet accumulation in the liver was visualized using Oil red O staining of frozen liver sections that were prepared in optimum cutting temperature (O.C.T.) compound. Liver fibrosis was visualized using Picro Sirius Red Stain.
Subcellular fractionation
Request a detailed protocolA Cytoplasmic & Nuclear fraction Kit (Beyotime Biotechnology, China) was used to detect Snhg3 expression in cytoplasmic and nuclear fractions. RNA was extracted from the cytoplasmic and nuclear fractions using a Trizol-based method and subjected to qPCR. Gapdh was used as a cytoplasmic marker, and Neat1 and Xist were used as a nuclear marker. The percentage of the transcript abundance was calculated using the following formula, Nucleus %=2^Ct(Nucleus)/(2^Ct(Cytoplasm)+2^Ct(Nucleus)), Cytoplasm %=1- Nucleus %.
Mouse calorimetry
Request a detailed protocolMale mice were housed individually in metabolic chambers of an Oxymax system (Columbus Instruments). The first readings were taken after a 24 hr acclimation period. Heat production, total carbon dioxide production and oxygen consumption, and RER were determined by Comprehensive laboratory animal monitoring system (CLAMS). The data were analyzed with CalR (Mina et al., 2018).
Insulin tolerance test (ITT) and glucose tolerance test (GTT)
Request a detailed protocolFor ITT, male mice fasted for 6 hr and received an intraperitoneal injection of human insulin (0.75 IU/kg). For GTT, male mice fasted for 6 hr or 16 hr received an intraperitoneal injection of glucose (1 g/kg). A series of blood glucose concentrations were measured from tail blood at the indicated times using a One-Touch Ultra glucometer (LifeScan Inc, Milpitas, CA).
Chromatin Immunoprecipitation (ChIP)
Request a detailed protocolThe ChIP assay was performed using Sonication ChIP Kit (Abclonal, Cat#RK20258). Briefly, the liver tissues or primary hepatocytes were collected and cross-linking fixed. Cross-linked chromatin fragments were precipitated with Rabbit control IgG (Abclonal, Cat#AC005) or anti-H3K27me3 antibody (Abclonal, Cat# A16199) for subsequent PCR analysis using the amplification primers for mouse Pparg promoter (+101 ~+ 420 bp; key resources table).
RNA immunoprecipitation (RIP)
Request a detailed protocolThe RIP assay was performed using RIP Assay Kit (BersinBio, Cat#Bes5101). Briefly, Hepa1-6 cells were transfected by indicated plasmids, respectively. The cells were collected, cross-linking fixed and precipitated with Mouse Control IgG (Abclonal, Cat#AC011) or anti-FLAG antibody (Abclonal, Cat#AE005) for subsequent RT-qPCR analysis using the amplification primers for Snhg3.
RNA sequencing (RNA-Seq)
Request a detailed protocolThe RNA-Seq was performed according to the manufacturer’s protocol (BGI-Shenzhen, https://www.yuque.com/yangyulan-ayaeq/oupzan/fuoao4). Briefly, total RNA was extracted from liver of three male DIO (27 weeks) mice and three male control mice for RNA-Seq to screen the differentially expressed lncRNAs using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. rRNA in total RNA removed using RNase H kit, was subsequently to construct library and perform sequencing analysis. The data were mapped to mouse genome (GRCm39) by using Bowtie2. The data for the differentially expressed lncRNAs had been deposited to National Genomics Data Center, China National Center for Bioinformation (NGDC-CNCB) (https://ngdc.cncb.ac.cn/) with the dataset identifier CRA009822. Additionally, total RNA was extracted from livers of three male DIO Snhg3-HKI mice and three male DIO WT mice for RNA-Seq to screen the differentially expressed mRNAs using Trizol. Total RNA was enriched by oligo (dT)-attached magnetic beads, followed by library construction and sequencing analysis. The data were mapped to mouse genome (GRCm39) by using Bowtie2. The data for the differentially expressed mRNAs have been deposited to the Sequence Read Archive (SRA) (submit.ncbi.nlm.nih.gov) with the dataset identifier SRR22368163, SRR22368164, SRR22368165, SRR22368166, SRR22368167, and SRR22368168.
Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq)
Request a detailed protocolThe ATAC-seq was performed according to manufacturer’s protocols (BGI_Shenzhen, https://www.yuque.com/yangyulan-ayaeq/oupzan/lllmzg). Briefly, fresh liver tissue samples from three male DIO Snhg3-HKI mice and three male DIO control mice were flash frozen by liquid nitrogen and then ground completely. The transposition reactions were initiated by adding transposase. The PCR reaction system was configured to initiate PCR amplification of the transposition products. The corresponding library quality control protocol would be selected depending on product requirements. Single-stranded PCR products were produced via denaturation. The reaction system and program for circularization were subsequently configured and set up. Single-stranded cyclized products were produced, while uncyclized linear DNA molecules were digested. Single-stranded circle DNA molecules were replicated via rolling cycle amplification, and a DNA nanoball (DNB) which contain multiple copies of DNA was generated. Sufficient quality DNBs were then loaded into patterned nanoarrays using high-intensity DNA nanochip technique and sequenced through combinatorial Probe-Anchor Synthesis (cPAS). Data were filtered by removing adaptor sequences, contamination and low-quality reads from raw reads. Bowtie2 was used to do genome alignment after evaluating its performance.
Peak Calling was performed by MACS (Model-based Analysis for ChIP-Seq). The candidate Peak region was extended to be long enough for modeling. Dynamic Possion Distribution was used to calculate the p-value of the specific region based on the unique mapped reads. The region would be defined as a Peak when the p-value <1e-05. MACS works well for the identification of the sharp peaks of most sequence-specific transcription factors.
Peak Distribution on Gene Elements. Peaks were classified based on the location (UCSC annotation data) and showed in the following genome regions: promoter (≤1 kb), promoter (1–2 kb), promoter (2–3 kb), 5'-UTR, 3'-UTR, intergenic, introns, and exons.
Differential Peaks were identified using MAnorm. First, the true intensities of the most common peaks were assumed to be the same between two ATAC-Seq samples. Second, the observed differences in sequence read density in common peaks were presumed to reflect the scaling relationship of ATAC-Seq signals between two samples, which could thus be applied to all peaks. Based on these hypotheses, the log2 ratio of read density between two samples M was plotted against the average log2 read density A for all peaks, and robust linear regression was applied to fit the global dependence between the M-A values of common peaks. Then the derived linear model was used as a reference for normalization and extrapolated to all peaks. Finally, the p-value for each Peak was calculated based on the Bayesian model, the significant regions were picked up if |M=log2FC|>2 and p-value <0.001. To identify DARs, the count matrix was input into MAnorm with the cutoff of abs(log2FC)>2 and p-value <0.001. Genomic features of DARs were annotated by R package ChIPseeker (v1.30.3).
Differential motifs analysis. After extracting the corresponding peak sequence, Hypergeometric Optimization of Motif EnRichment (HOMER, (ucsd.edu)) and Cis-Regulatory Element Motif Activities (CREMA (unibas.ch)) were used for motif analysis. Genomic regions with differential ATAC peaks were shown using IGV software.
The data of ATAC-seq has been deposited to National Genomics Data Center, China National Center for Bioinformation (NGDC-CNCB) (https://ngdc.cncb.ac.cn/) with the dataset identifier CRA009511.
Integrated analysis ATAC-Seq data with RNA-Seq data
Request a detailed protocolThe common number and unique number of genes associated with DARs and DEGs were counted using Wayne comparative analysis. The values of each quadrant satisfying the log2FC condition were selected from DARs-associated genes and DEGs to draw a nine-quadrant plot, respectively. Pearson correlation was calculated for both sets of data and p-value was calculated using the Z-test.
Cleavage Under Targets and Tagmentation sequencing (CUT&Tag-Seq)
Request a detailed protocolCUT&Tag was performed according to the Hyperactive Universal CUT&Tag Assay Kit for Illumina (Vazyme, Nanjing, Cat#TD903-01). Briefly, the mixed liver tissues from three DIO Snhg3-HKO mice were used for the CUT&Tag experiment. pA-Tn5 transposase was used to cut the genome and add a special adaptor sequence to build a library. The single-stranded PCR products were sequenced on illumina/DNBSEQ-T7 platform PE150 (Annoroad Gene Technology Co.Itd). The data were mapped to mouse genome (GRCm38) by using Bowtie2. Genomic features of DARs were annotated by R package ChIPseeker (v1.30.3). The data of CUT&Tag has been deposited to National Genomics Data Center, China National Center for Bioinformation (NGDC-CNCB) (https://ngdc.cncb.ac.cn/) with the dataset identifier CRA009582.
Biotin-RNA pull-down and mass spectrometry assay
Request a detailed protocolBiotin-RNA pull-down assay was performed as described in a previous study (Guo et al., 2018). Briefly, Snhg3 DNA was amplified from mouse liver cDNA using the primers listed in the key resources table and lacZ DNA fragment were constructed into pGEM-T easy vector. The pGEM-T-Snhg3 and pGEM-T-lacZ vectors were linearized by restriction enzyme digestion, then transcribed to Snhg3 and lacZ fragments. Biotinylated RNAs were transcribed in vitro with Biotin-RNA Labeling Mix (Roche, Indianapolis, IN) and purified with quick spin RNA columns (Roche, Indianapolis, IN). Biotin-labeled RNA or unbiotinylated RNAs was dissolved in RNA structure buffer (10 mM Tris, pH 7.0, 0.1 M KCl, 10 mM MgCl2) to allow formation of the secondary structure. Primary hepatocytes lysates were added to biotin-labeled RNA or unbiotinylated RNA. Streptavidin agarose beads (GE Healthcare, Little Chalfont, UK) were mixed with a pull-down reaction and then rotated constantly. RNA affinity captures were subjected to 12% SDS-PAGE followed by coomassie blue staining or Western blotting. The various bands that were visualized by coomassie blue staining were excised and subjected to mass spectrometry analyses (LC-MS/MS, A TripleTOF , ABsciex, Concord, ON). The data of RNA pull-down for Snhg3, control or lacZ were deposited to the iProX (https://www.iprox.cn/) with the dataset identifier PXD039526.
Ubiquitination assays
Request a detailed protocolFor endogenous ubiquitination assays, Hepa1-6 cells were transfected with the indicated combinations of plasmids, including HA-ubiquitin and Snhg3 plasmids. For exogenous ubiquitination assays, Hepa1-6 cells were transfected with the indicated combinations of plasmids, including HA-ubiquitin, HA-K33-ubiquitin, HA-K63-ubiquitin, HA-K48-ubiquitin, Flag-SND1 and Snhg3 plasmids. Cells were treated with 20 µM MG132 proteasome inhibitor (M1902, AbMole) for 6 hr prior to lyse in lysis buffer 200 mM NaCl, 20 mM Tris-HCl (pH 7.4), 2.5 mM MgCl2, 0.5% Triton X-100, 1 mM PMSF, and protease inhibitor cocktail and then were sonicated. After centrifugation at 14,000 × g, the cleared lysates were subjected to immunoprecipitation with anti-SND1 antibody (sc-166676, Santa Cruz) for endogenous ubiquitination assays or with anti-DDDDK-tag magnetic beads (M185-10R, MBL) for exogenous ubiquitination assays. The immunocomplexes were collected and subjected to estern blotting with the indicated antibodies.
Statistical analysis
Request a detailed protocolData analyses were performed with SPSS (Version 17.0, SPSS, Inc). The curves of body weight and ITT were analyzed using a repeated measure two-way ANOVA. For the other statistical analysis, the Kolmogorov-Smirnov test was firstly used for normality test. For the data conforming to the normal distribution, experiments were analyzed using Independent-Samples T-test or one-way ANOVA. All data were presented as the mean ± SD or the mean ± SEM. * p<0.05 was considered statistically significant.
Appendix 1
Data availability
The lncRNA-Seq data had been deposited to National Genomics Data Center, China National Center for Bioinformation (NGDC-CNCB) (https://ngdc.cncb.ac.cn/) with the dataset identifier CRA009822. The data of RNA pull-down for Snhg3, control or lacZ have been deposited to the iProX (https://www.iprox.cn/) with the dataset identifier IPX0005781000 (ProteomeXchange identifier: PXD039526). The data of RNA-seq have been deposited to the Sequence Read Archive (SRA) with the dataset identifier SRR22368163, SRR22368164, SRR22368165, SRR22368166, SRR22368167 and SRR22368168. The data of ATAC-seq and Cut&Tag have been deposited to National Genomics Data Center, China National Center for Bioinformation (NGDC-CNCB) (https://ngdc.cncb.ac.cn/) with the dataset identifier CRA009511 and CRA009582, respectively.
-
Genome Sequence ArchiveID CRA009822. The lncRNA-seq of mice liver.
-
Integrated Proteome ResourcesID IPX0005781000. RNA-Pull down and MS for Snhg3.
-
Genome Sequence ArchiveID CRA009511. The ATAC-seq of diet induced obesity Snhg3-LKI mice liver.
-
Genome Sequence ArchiveID CRA009582. The CUT&TAG-seq of Snhg3-LKI mice liver.
-
NCBI Sequence Read ArchiveID SRR22368163. lncRNA_Snhg3 regulates hepatic lipid metabolism.
-
NCBI Sequence Read ArchiveID SRR22368164. lncRNA_Snhg3 regulates hepatic lipid metabolism.
-
NCBI Sequence Read ArchiveID SRR22368165. lncRNA_Snhg3 regulates hepatic lipid metabolism.
-
NCBI Sequence Read ArchiveID SRR22368166. lncRNA_Snhg3 regulates hepatic lipid metabolism.
-
NCBI Sequence Read ArchiveID SRR22368167. lncRNA_Snhg3 regulates hepatic lipid metabolism.
-
NCBI Sequence Read ArchiveID SRR22368168. lncRNA_Snhg3 regulates hepatic lipid metabolism.
-
ProteomeXchangeID PXD039526. RNA-Pull down and MS for Snhg3.
References
-
MicroRNAs in nonalcoholic fatty liver diseaseJournal of Clinical Medicine 4:1977–1988.https://doi.org/10.3390/jcm4121953
-
PPAR-γ signaling in nonalcoholic fatty liver disease: Pathogenesis and therapeutic targetsPharmacology & Therapeutics 245:108391.https://doi.org/10.1016/j.pharmthera.2023.108391
-
Tudor-SN, a novel coactivator of peroxisome proliferator-activated receptor γ protein, is essential for adipogenesisThe Journal of Biological Chemistry 289:8364–8374.https://doi.org/10.1074/jbc.M113.523456
-
Genetics and epigenetics of NAFLD and NASH: Clinical impactJournal of Hepatology 68:268–279.https://doi.org/10.1016/j.jhep.2017.09.003
-
Mechanisms of NAFLD development and therapeutic strategiesNature Medicine 24:908–922.https://doi.org/10.1038/s41591-018-0104-9
-
Liver peroxisome proliferator-activated receptor gamma contributes to hepatic steatosis, triglyceride clearance, and regulation of body fat massThe Journal of Biological Chemistry 278:34268–34276.https://doi.org/10.1074/jbc.M300043200
-
MAFLD: How is it different from NAFLD?Clinical and Molecular Hepatology 29:S17–S31.https://doi.org/10.3350/cmh.2022.0367
-
PPARs in obesity-induced T2DM, dyslipidaemia and NAFLDNature Reviews. Endocrinology 13:36–49.https://doi.org/10.1038/nrendo.2016.135
-
snoRNA U17 regulates cellular cholesterol traffickingCell Metabolism 21:855–867.https://doi.org/10.1016/j.cmet.2015.04.010
-
Genetic and epigenetic factors determining NAFLD riskMolecular Metabolism 50:101111.https://doi.org/10.1016/j.molmet.2020.101111
-
CPC2: a fast and accurate coding potential calculator based on sequence intrinsic featuresNucleic Acids Research 45:W12–W16.https://doi.org/10.1093/nar/gkx428
-
Stress-induced FGF21 and GDF15 in obesity and obesity resistanceTrends in Endocrinology and Metabolism 32:904–915.https://doi.org/10.1016/j.tem.2021.08.008
-
Molecular and cellular mechanisms of liver fibrosis and its regressionNature Reviews. Gastroenterology & Hepatology 18:151–166.https://doi.org/10.1038/s41575-020-00372-7
-
Epigenetics in non-alcoholic fatty liver diseaseMolecular Aspects of Medicine 54:78–88.https://doi.org/10.1016/j.mam.2016.11.008
-
Hepatocyte-specific loss of PPARγ protects mice from NASH and increases the therapeutic effects of rosiglitazone in the liverCellular and Molecular Gastroenterology and Hepatology 11:1291–1311.https://doi.org/10.1016/j.jcmgh.2021.01.003
-
Hepatocyte PPARγ contributes to the progression of non-alcoholic steatohepatitis in male and female obese miceCellular and Molecular Life Sciences 80:39.https://doi.org/10.1007/s00018-022-04629-z
-
Role of hepatic peroxisome proliferator-activated receptor γ in non-alcoholic fatty liver diseaseThe Journal of Endocrinology 257:e220155.https://doi.org/10.1530/JOE-22-0155
-
Molecular and cellular insights into the role of SND1 in lipid metabolismBiochimica et Biophysica Acta. Molecular and Cell Biology of Lipids 1865:158589.https://doi.org/10.1016/j.bbalip.2019.158589
-
Recent insights into stearoyl-CoA desaturase-1Current Opinion in Lipidology 14:255–261.https://doi.org/10.1097/00041433-200306000-00005
-
From 1957 to Nowadays: A Brief History of EpigeneticsInternational Journal of Molecular Sciences 21:20.https://doi.org/10.3390/ijms21207571
-
Epigenetics factors in nonalcoholic fatty liver diseaseExpert Review of Gastroenterology & Hepatology 16:521–536.https://doi.org/10.1080/17474124.2020.1765772
-
Identification of human long noncoding RNAs associated with nonalcoholic fatty liver disease and metabolic homeostasisThe Journal of Clinical Investigation 131:e136336.https://doi.org/10.1172/JCI136336
-
Role of Insulin Resistance in MAFLDInternational Journal of Molecular Sciences 22:4156.https://doi.org/10.3390/ijms22084156
-
Are spliced ncRNA host genes distinct classes of lncRNAs?Theory in Biosciences = Theorie in Den Biowissenschaften 139:349–359.https://doi.org/10.1007/s12064-020-00330-6
-
A role of peroxisome proliferator-activated receptor γ in non-alcoholic fatty liver diseaseBasic & Clinical Pharmacology & Toxicology 124:528–537.https://doi.org/10.1111/bcpt.13190
-
Epigenetics in NAFLD/NASH: Targets and therapyPharmacological Research 167:105484.https://doi.org/10.1016/j.phrs.2021.105484
-
Noncoding RNAs and RNA-binding proteins: emerging governors of liver physiology and metabolic diseasesAmerican Journal of Physiology. Cell Physiology 323:C1003–C1017.https://doi.org/10.1152/ajpcell.00232.2022
-
Histone core modifications regulating nucleosome structure and dynamicsNature Reviews. Molecular Cell Biology 15:703–708.https://doi.org/10.1038/nrm3890
-
Endogenous microRNA sponges: evidence and controversyNature Reviews. Genetics 17:272–283.https://doi.org/10.1038/nrg.2016.20
-
FGF21: An emerging therapeutic target for non-alcoholic steatohepatitis and related metabolic diseasesFrontiers in Endocrinology 11:601290.https://doi.org/10.3389/fendo.2020.601290
-
Tudor domain-containing proteins of Drosophila melanogasterDevelopment, Growth & Differentiation 54:32–43.https://doi.org/10.1111/j.1440-169x.2011.01308.x
-
Epigenetically silenced lncRNA SNAI3-AS1 promotes ferroptosis in glioma via perturbing the m6A-dependent recognition of Nrf2 mRNA mediated by SND1Journal of Experimental & Clinical Cancer Research 42:127.https://doi.org/10.1186/s13046-023-02684-3
Article and author information
Author details
Funding
Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences (CIFMS2021-I2M-1-016)
- Xiaojun Liu
Beijing Natural Science Foundation (7242094)
- Xiaojun Liu
National Key R&D Program of China (2022YFC2504002)
- Li Yan
National Key R&D Program of China (2022YFC2504003)
- Xiaojun Liu
the National Natural Science Foundation of China (82270925)
- Aijun Qiao
High-level New R&D Institute of Department of Science and Technology of Guangdong Province (2019B090904008)
- Aijun Qiao
High-level Innovative Research Institute of Department of Science and Technology of Guangdong Province (2021B0909050003)
- Aijun Qiao
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences (CIFMS2021-I2M-1–016 to XL), National Key R&D Program of China (2022YFC2504002 to LY, 2022YFC2504003 to XL), Beijing Natural Science Foundation (7242094 to XL), the National Natural Science Foundation of China (82270925 to AQ), High-level New R&D Institute of Department of Science and Technology of Guangdong Province (2019B090904008 to AQ), and High-level Innovative Research Institute of Department of Science and Technology of Guangdong Province (2021B0909050003 to AQ). We thank Dr. Yi Li (Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College) for providing the plasmid of HA-Ub, HA-Ub-K63 and HA-Ub-K48. We thank Cyagen Biosciences for the collaborative efforts in the creation of Snhg3-loxP and liver-specific knock-in Snhg3 mice. We thank BGI-Shenzhen for the collaborative efforts in RNA-Seq and ATAC-Seq technology. We thank Anoroad gene technology (Beijing) for the collaborative efforts in CUT&Tag-Seq and Jingjie PTM BioLab (Hangzhou) Co. Ltd for the collaborative efforts in mass spectrometry technology resources.
Ethics
All animal experiments were conducted under protocols approved by the Animal Research Committee of the Institute of Laboratory Animals, Institute of Basic Medical Sciences Chinese Academy of Medical Sciences & School of Basic Medicine Peking Union Medical College (ACUC-A01-2022-010).
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.96988. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Xie et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 659
- views
-
- 70
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Chromosomes and Gene Expression
- Evolutionary Biology
Gene regulation is essential for life and controlled by regulatory DNA. Mutations can modify the activity of regulatory DNA, and also create new regulatory DNA, a process called regulatory emergence. Non-regulatory and regulatory DNA contain motifs to which transcription factors may bind. In prokaryotes, gene expression requires a stretch of DNA called a promoter, which contains two motifs called –10 and –35 boxes. However, these motifs may occur in both promoters and non-promoter DNA in multiple copies. They have been implicated in some studies to improve promoter activity, and in others to repress it. Here, we ask whether the presence of such motifs in different genetic sequences influences promoter evolution and emergence. To understand whether and how promoter motifs influence promoter emergence and evolution, we start from 50 ‘promoter islands’, DNA sequences enriched with –10 and –35 boxes. We mutagenize these starting ‘parent’ sequences, and measure gene expression driven by 240,000 of the resulting mutants. We find that the probability that mutations create an active promoter varies more than 200-fold, and is not correlated with the number of promoter motifs. For parent sequences without promoter activity, mutations created over 1500 new –10 and –35 boxes at unique positions in the library, but only ~0.3% of these resulted in de-novo promoter activity. Only ~13% of all –10 and –35 boxes contribute to de-novo promoter activity. For parent sequences with promoter activity, mutations created new –10 and –35 boxes in 11 specific positions that partially overlap with preexisting ones to modulate expression. We also find that –10 and –35 boxes do not repress promoter activity. Overall, our work demonstrates how promoter motifs influence promoter emergence and evolution. It has implications for predicting and understanding regulatory evolution, de novo genes, and phenotypic evolution.
-
- Chromosomes and Gene Expression
- Developmental Biology
The male-specific lethal complex (MSL), which consists of five proteins and two non-coding roX RNAs, is involved in the transcriptional enhancement of X-linked genes to compensate for the sex chromosome monosomy in Drosophila XY males compared with XX females. The MSL1 and MSL2 proteins form the heterotetrameric core of the MSL complex and are critical for the specific recruitment of the complex to the high-affinity ‘entry’ sites (HAS) on the X chromosome. In this study, we demonstrated that the N-terminal region of MSL1 is critical for stability and functions of MSL1. Amino acid deletions and substitutions in the N-terminal region of MSL1 strongly affect both the interaction with roX2 RNA and the MSL complex binding to HAS on the X chromosome. In particular, substitution of the conserved N-terminal amino-acids 3–7 in MSL1 (MSL1GS) affects male viability similar to the inactivation of genes encoding roX RNAs. In addition, MSL1GS binds to promoters such as MSL1WT but does not co-bind with MSL2 and MSL3 to X chromosomal HAS. However, overexpression of MSL2 partially restores the dosage compensation. Thus, the interaction of MSL1 with roX RNA is critical for the efficient assembly of the MSL complex on HAS of the male X chromosome.