Perturbed rhythmic activation of signaling pathways in mice deficient for Sterol Carrier Protein 2-dependent diurnal lipid transport and metabolism

Transcriptomics-driven metabolic pathway analysis reveals similar alterations in lipid metabolism in mouse MASH model and human

BACKGROUND

Metabolic dysfunction-associated steatotic liver disease (MASLD) is a prevalent chronic liver disease worldwide, and can rapidly progress to metabolic dysfunction-associated steatohepatitis (MASH). Accurate preclinical models and methodologies are needed to understand underlying metabolic mechanisms and develop treatment strategies. Through meta-analysis of currently proposed mouse models, we hypothesized that a diet- and chemical-induced MASH model closely resembles the observed lipid metabolism alterations in humans.

METHODS

We developed transcriptomics-driven metabolic pathway analysis (TDMPA), a method to aid in the evaluation of metabolic resemblance. TDMPA uses genome-scale metabolic models to calculate enzymatic reaction perturbations from gene expression data. We performed TDMPA to score and compare metabolic pathway alterations in MASH mouse models to human MASH signatures. We used an already-established WD+CCl4-induced MASH model and performed functional assays and lipidomics to confirm TDMPA findings.

RESULTS

Both human MASH and mouse models exhibit numerous altered metabolic pathways, including triglyceride biosynthesis, fatty acid beta-oxidation, bile acid biosynthesis, cholesterol metabolism, and oxidative phosphorylation. We confirm a significant reduction in mitochondrial functions and bioenergetics, as well as in acylcarnitines for the mouse model. We identify a wide range of lipid species within the most perturbed pathways predicted by TDMPA. Triglycerides, phospholipids, and bile acids are increased significantly in mouse MASH liver, confirming our initial observations.

CONCLUSIONS

We introduce TDMPA, a methodology for evaluating metabolic pathway alterations in metabolic disorders. By comparing metabolic signatures that typify human MASH, we show a good metabolic resemblance of the WD+CCl4 mouse model. Our presented approach provides a valuable tool for defining metabolic space to aid experimental design for assessing metabolism.

Characterisation of hepatic lipid signature distributed across the liver zonation using mass spectrometry imaging

Background & Aims

Lipid metabolism plays an important role in liver pathophysiology. The liver lobule asymmetrically distributes oxygen and nutrition, resulting in heterogeneous metabolic functions. Periportal and pericentral hepatocytes have different metabolic functions, which lead to generating liver zonation. We developed spatial metabolic imaging using desorption electrospray ionisation mass spectrometry to investigate lipid distribution across liver zonation with high reproducibility and accuracy.

Methods

Fresh frozen livers from healthy mice with control diet were analysed using desorption electrospray ionisation mass spectrometry imaging. Imaging was performed at 50 μm × 50 μm pixel size. Regions of interest (ROIs) were manually created by co-registering with histological data to determine the spatial hepatic lipids across liver zonation. The ROIs were confirmed by double immunofluorescence. The mass list of specific ROIs was automatically created, and univariate and multivariate statistical analysis were performed to identify statistically significant lipids across liver zonation.

Results

A wide range of lipid species was identified, including fatty acids, phospholipids, triacylglycerols, diacylglycerols, ceramides, and sphingolipids. We characterised hepatic lipid signatures in three different liver zones (periportal zone, midzone, and pericentral zone) and validated the reproducibility of our method for measuring a wide range of lipids. Fatty acids were predominantly detected in the periportal region, whereas phospholipids were distributed in both the periportal and pericentral zones. Interestingly, phosphatidylinositols, PI(36:2), PI(36:3), PI(36:4), PI(38:5), and PI(40:6) were located predominantly in the midzone (zone 2). Triacylglycerols and diacylglycerols were detected mainly in the pericentral region. De novo triacylglycerol biosynthesis appeared to be the most influenced pathway across the three zones.

Conclusions

The ability to accurately assess zone-specific hepatic lipid distribution in the liver could lead to a better understanding of lipid metabolism during the progression of liver disease.

Impact and Implications

Zone-specific hepatic lipid metabolism could play an important role in lipid homoeostasis during disease progression. Herein, we defined the zone-specific references of hepatic lipid species in the three liver zones using molecular imaging. The de novo triacylglycerol biosynthesis was highlighted as the most influenced pathway across the three zones.

Steroidogenic Factor 1 Regulation of the Hypothalamic-Pituitary-Ovarian Axis of Adult Female Mice

The orphan nuclear receptor steroidogenic factor-1 (SF-1 or NR5A1) is an indispensable regulator of adrenal and gonadal formation, playing roles in sex determination, hypothalamic development, and pituitary function. This study aimed to identify the roles of SF-1 in postnatal female reproductive function. Using a progesterone receptor-driven Cre recombinase, we developed a novel murine model, characterized by conditional depletion of SF-1 [PR-Cre;Nr5a1f/f; conditional knockout (cKO)] in the hypothalamic-pituitary-gonadal axis. Mature female cKO were infertile due to the absence of ovulation. Reduced gonadotropin concentrations in the pituitary gland that were nevertheless sufficient to maintain regular estrous cycles were observed in mature cKO females. The cKO ovaries showed abnormal lipid accumulation in the stroma, associated with an irregular expression of cholesterol homeostatic genes such as Star, Scp2, and Acat1. The depletion of SF-1 in granulosa cells prevented appropriate cumulus oöphorus expansion, characterized by reduced expression of Areg, Ereg, and Ptgs2. Exogenous delivery of gonadotropins to cKO females to induce ovulation did not restore fertility and was associated with impaired formation and function of corpora lutea accompanied by reduced expression of the steroidogenic genes Cyp11a1 and Cyp19a1 and attenuated progesterone production. Surgical transplantation of cKO ovaries to ovariectomized control animals (Nr5a1f/f) resulted in 2 separate phenotypes, either sterility or apparently normal fertility. The deletion of SF-1 in the pituitary and in granulosa cells near the moment of ovulation demonstrated that this nuclear receptor functions across the pituitary-gonadal axis and plays essential roles in gonadotropin synthesis, cumulus expansion, and luteinization.

Systematic analysis of differential rhythmic liver gene expression mediated by the circadian clock and feeding rhythms

The circadian clock and feeding rhythms are both important regulators of rhythmic gene expression in the liver. To further dissect the respective contributions of feeding and the clock, we analyzed differential rhythmicity of liver tissue samples across several conditions. We developed a statistical method tailored to compare rhythmic liver messenger RNA (mRNA) expression in mouse knockout models of multiple clock genes, as well as PARbZip output transcription factors (Hlf/Dbp/Tef). Mice were exposed to ad libitum or night-restricted feeding under regular light-dark cycles. During ad libitum feeding, genetic ablation of the core clock attenuated rhythmic-feeding patterns, which could be restored by the night-restricted feeding regimen. High-amplitude mRNA expression rhythms in wild-type livers were driven by the circadian clock, but rhythmic feeding also contributed to rhythmic gene expression, albeit with significantly lower amplitudes. We observed that Bmal1 and Cry1/2 knockouts differed in their residual rhythmic gene expression. Differences in mean expression levels between wild types and knockouts correlated with rhythmic gene expression in wild type. Surprisingly, in PARbZip knockout mice, the mean expression levels of PARbZip targets were more strongly impacted than their rhythms, potentially due to the rhythmic activity of the D-box-repressor NFIL3. Genes that lost rhythmicity in PARbZip knockouts were identified to be indirect targets. Our findings provide insights into the diurnal transcriptome in mouse liver as we identified the differential contributions of several core clock regulators. In addition, we gained more insights on the specific effects of the feeding-fasting cycle.

Reproducibility across single-cell RNA-seq protocols for spatial ordering analysis

As newer single-cell protocols generate increasingly more cells at reduced sequencing depths, the value of a higher read depth may be overlooked. Using data from three different single-cell RNA-seq protocols that lend themselves to having either higher read depth (Smart-seq) or many cells (MARS-seq and 10X), we evaluate their ability to recapitulate biological signals in the context of spatial reconstruction. Overall, we find gene expression profiles after spatial reconstruction analysis are highly reproducible between datasets despite being generated by different protocols and using different computational algorithms. While UMI-based protocols such as 10X and MARS-seq allow for capturing more cells, Smart-seq's higher sensitivity and read-depth allow for analysis of lower expressed genes and isoforms. Additionally, we evaluate trade-offs for each protocol by performing subsampling analyses and find that optimizing the balance between sequencing depth and number of cells within a protocol is necessary for efficient use of resources. Our analysis emphasizes the importance of selecting a protocol based on the biological questions and features of interest.

Atlas of Circadian Metabolism Reveals System-wide Coordination and Communication between Clocks

Metabolic diseases are often characterized by circadian misalignment in different tissues, yet how altered coordination and communication among tissue clocks relate to specific pathogenic mechanisms remains largely unknown. Applying an integrated systems biology approach, we performed 24-hr metabolomics profiling of eight mouse tissues simultaneously. We present a temporal and spatial atlas of circadian metabolism in the context of systemic energy balance and under chronic nutrient stress (high-fat diet [HFD]). Comparative analysis reveals how the repertoires of tissue metabolism are linked and gated to specific temporal windows and how this highly specialized communication and coherence among tissue clocks is rewired by nutrient challenge. Overall, we illustrate how dynamic metabolic relationships can be reconstructed across time and space and how integration of circadian metabolomics data from multiple tissues can improve our understanding of health and disease.

Transcriptional programming of lipid and amino acid metabolism by the skeletal muscle circadian clock

Circadian clocks are fundamental physiological regulators of energy homeostasis, but direct transcriptional targets of the muscle clock machinery are unknown. To understand how the muscle clock directs rhythmic metabolism, we determined genome-wide binding of the master clock regulators brain and muscle ARNT-like protein 1 (BMAL1) and REV-ERBα in murine muscles. Integrating occupancy with 24-hr gene expression and metabolomics after muscle-specific loss of BMAL1 and REV-ERBα, here we unravel novel molecular mechanisms connecting muscle clock function to daily cycles of lipid and protein metabolism. Validating BMAL1 and REV-ERBα targets using luciferase assays and in vivo rescue, we demonstrate how a major role of the muscle clock is to promote diurnal cycles of neutral lipid storage while coordinately inhibiting lipid and protein catabolism prior to awakening. This occurs by BMAL1-dependent activation of Dgat2 and REV-ERBα-dependent repression of major targets involved in lipid metabolism and protein turnover (MuRF-1, Atrogin-1). Accordingly, muscle-specific loss of BMAL1 is associated with metabolic inefficiency, impaired muscle triglyceride biosynthesis, and accumulation of bioactive lipids and amino acids. Taken together, our data provide a comprehensive overview of how genomic binding of BMAL1 and REV-ERBα is related to temporal changes in gene expression and metabolite fluctuations.

Circadian- and UPR-dependent control of CPEB4 mediates a translational response to counteract hepatic steatosis under ER stress

The cytoplasmic polyadenylation element-binding (CPEB) proteins regulate pre-mRNA processing and translation of CPE-containing mRNAs in early embryonic development and synaptic activity. However, specific functions in adult organisms are poorly understood. Here we show that CPEB4 is required for adaptation to high-fat-diet- and ageing-induced endoplasmic reticulum (ER) stress, and subsequent hepatosteatosis. Stress-activated liver CPEB4 expression is dual-mode regulated. First, Cpeb4 mRNA transcription is controlled by the circadian clock, and then its translation is regulated by the unfolded protein response (UPR) through upstream open reading frames within the 5'UTR. Thus, the CPEB4 protein is synthesized only following ER stress but the induction amplitude is circadian. In turn, CPEB4 activates a second wave of UPR translation required to maintain ER and mitochondrial homeostasis. Our results suggest that combined transcriptional and translational Cpeb4 regulation generates a 'circadian mediator', which coordinates hepatic UPR activity with periods of high ER-protein-folding demand. Accordingly, CPEB4 deficiency results in non-alcoholic fatty liver disease.