Development of a functional beige fat cell line uncovers independent subclasses of cells expressing UCP1 and the futile creatine cycle.

Although uncoupling protein 1 (UCP1) is established as a major contributor to adipose thermogenesis, recent data have illustrated an important role for alternative pathways, particularly the futile creatine cycle (FCC). How these pathways co-exist in cells and tissues has not been explored. Beige cell adipogenesis occurs in vivo but has been difficult to model in vitro; here, we describe the development of a murine beige cell line that executes a robust respiratory response, including uncoupled respiration and the FCC. The key FCC enzyme, tissue-nonspecific alkaline phosphatase (TNAP), is localized almost exclusively to mitochondria in these cells. Surprisingly, single-cell cloning from this cell line shows that cells with the highest levels of UCP1 express little TNAP, and cells with the highest expression of TNAP express little UCP1. Immunofluorescence analysis of subcutaneous fat from cold-exposed mice confirms that the highest levels of these critical thermogenic components are expressed in distinct fat cell populations.

A spatially-resolved transcriptional atlas of the murine dorsal pons at single-cell resolution.

The "dorsal pons", or "dorsal pontine tegmentum" (dPnTg), is part of the brainstem. It is a complex, densely packed region whose nuclei are involved in regulating many vital functions. Notable among them are the parabrachial nucleus, the Kölliker Fuse, the Barrington nucleus, the locus coeruleus, and the dorsal, laterodorsal, and ventral tegmental nuclei. In this study, we applied single-nucleus RNA-seq (snRNA-seq) to resolve neuronal subtypes based on their unique transcriptional profiles and then used multiplexed error robust fluorescence in situ hybridization (MERFISH) to map them spatially. We sampled ~1 million cells across the dPnTg and defined the spatial distribution of over 120 neuronal subtypes. Our analysis identified an unpredicted high transcriptional diversity in this region and pinpointed the unique marker genes of many neuronal subtypes. We also demonstrated that many neuronal subtypes are transcriptionally similar between humans and mice, enhancing this study's translational value. Finally, we developed a freely accessible, GPU and CPU-powered dashboard ( http://harvard.heavy.ai:6273/ ) that combines interactive visual analytics and hardware-accelerated SQL into a data science framework to allow the scientific community to query and gain insights into the data.

Protocol for flow cytometry-assisted single-nucleus RNA sequencing of human and mouse adipose tissue with sample multiplexing.

Adipocyte size and fragility and commercial kit costs impose significant limitations on single-cell RNA sequencing of adipose tissue. Accordingly, we developed a workflow to isolate and sample-barcode nuclei from individual adipose tissue samples, integrating flow cytometry for quality control, counting, and precise nuclei pooling for direct loading onto the popular 10× Chromium controller. This approach can eliminate batch confounding, and significantly reduces poor-quality nuclei, ambient RNA contamination, and droplet loading-associated reagent waste, resulting in pronounced improvements in information content and cost efficiency.

A spatially-resolved transcriptional atlas of the murine dorsal pons at single-cell resolution.

The "dorsal pons", or "dorsal pontine tegmentum" (dPnTg), is part of the brainstem. It is a complex, densely packed region whose nuclei are involved in regulating many vital functions. Notable among them are the parabrachial nucleus, the Kölliker Fuse, the Barrington nucleus, the locus coeruleus, and the dorsal, laterodorsal, and ventral tegmental nuclei. In this study, we applied single-nucleus RNA-seq (snRNA-seq) to resolve neuronal subtypes based on their unique transcriptional profiles and then used multiplexed error robust fluorescence in situ hybridization (MERFISH) to map them spatially. We sampled ~1 million cells across the dPnTg and defined the spatial distribution of over 120 neuronal subtypes. Our analysis identified an unpredicted high transcriptional diversity in this region and pinpointed many neuronal subtypes' unique marker genes. We also demonstrated that many neuronal subtypes are transcriptionally similar between humans and mice, enhancing this study's translational value. Finally, we developed a freely accessible, GPU and CPU-powered dashboard (http://harvard.heavy.ai:6273/) that combines interactive visual analytics and hardware-accelerated SQL into a data science framework to allow the scientific community to query and gain insights into the data.

Genetic variants in genes involved in creatine biosynthesis in patients with severe obesity or anorexia nervosa.

Increased thermogenesis in brown adipose tissue might have an obesity-reducing effect in humans. In transgenic mice, depletion of genes involved in creatine metabolism results in disrupted thermogenic capacity and altered effects of high-fat feeding on body weight. Data analyses of a sex-stratified genome-wide association study (GWAS) for body mass index (BMI) within the genomic regions of genes of this pathway (CKB, CKMT1B, and GATM) revealed one sex-dimorphic BMI-associated SNP in CKB (rs1136165). The effect size was larger in females than in males. A mutation screen of the coding regions of these three candidate genes in a screening group (192 children and adolescents with severe obesity, 192 female patients with anorexia nervosa, and 192 healthy-lean controls) identified five variants in each, CKB and GATM, and nine variants in the coding sequence of CKMT1B. Non-synonymous variants identified in CKB and CKMT1B were genotyped in an independent confirmation study group (781 families with severe obesity (trios), 320 children and adolescents with severe obesity, and 253 healthy-lean controls). In silico tools predicted mainly benign yet protein-destabilizing potentials. A transmission disequilibrium test in trios with severe obesity indicated an obesity-protective effect of the infrequent allele at rs149544188 located in CKMT1B. Subsequent correlation analyses in 1,479 individuals of the Leipzig Obesity BioBank revealed distinct correlations of CKB with the other two genes in omental visceral adipose tissue (VAT) and abdominal subcutaneous adipose tissue (SAT). Furthermore, between-subject comparisons of gene expression levels showed generally higher expressions of all three genes of interest in VAT than in SAT. Future in vitro analyses are needed to assess the functional implications of these findings.

HGFAC is a ChREBP-regulated hepatokine that enhances glucose and lipid homeostasis.

Carbohydrate response element-binding protein (ChREBP) is a carbohydrate-sensing transcription factor that regulates both adaptive and maladaptive genomic responses in coordination of systemic fuel homeostasis. Genetic variants in the ChREBP locus associate with diverse metabolic traits in humans, including circulating lipids. To identify novel ChREBP-regulated hepatokines that contribute to its systemic metabolic effects, we integrated ChREBP ChIP-Seq analysis in mouse liver with human genetic and genomic data for lipid traits and identified hepatocyte growth factor activator (HGFAC) as a promising ChREBP-regulated candidate in mice and humans. HGFAC is a protease that activates the pleiotropic hormone hepatocyte growth factor. We demonstrate that HGFAC-KO mice had phenotypes concordant with putative loss-of-function variants in human HGFAC. Moreover, in gain- and loss-of-function genetic mouse models, we demonstrate that HGFAC enhanced lipid and glucose homeostasis, which may be mediated in part through actions to activate hepatic PPARγ activity. Together, our studies show that ChREBP mediated an adaptive response to overnutrition via activation of HGFAC in the liver to preserve glucose and lipid homeostasis.

Architecture of the outbred brown fat proteome defines regulators of metabolic physiology.

Brown adipose tissue (BAT) regulates metabolic physiology. However, nearly all mechanistic studies of BAT protein function occur in a single inbred mouse strain, which has limited the understanding of generalizable mechanisms of BAT regulation over physiology. Here, we perform deep quantitative proteomics of BAT across a cohort of 163 genetically defined diversity outbred mice, a model that parallels the genetic and phenotypic variation found in humans. We leverage this diversity to define the functional architecture of the outbred BAT proteome, comprising 10,479 proteins. We assign co-operative functions to 2,578 proteins, enabling systematic discovery of regulators of BAT. We also identify 638 proteins that correlate with protection from, or sensitivity to, at least one parameter of metabolic disease. We use these findings to uncover SFXN5, LETMD1, and ATP1A2 as modulators of BAT thermogenesis or adiposity, and provide OPABAT as a resource for understanding the conserved mechanisms of BAT regulation over metabolic physiology.

A single-cell atlas of human and mouse white adipose tissue.

White adipose tissue, once regarded as morphologically and functionally bland, is now recognized to be dynamic, plastic and heterogenous, and is involved in a wide array of biological processes including energy homeostasis, glucose and lipid handling, blood pressure control and host defence1. High-fat feeding and other metabolic stressors cause marked changes in adipose morphology, physiology and cellular composition1, and alterations in adiposity are associated with insulin resistance, dyslipidemia and type 2 diabetes2. Here we provide detailed cellular atlases of human and mouse subcutaneous and visceral white fat at single-cell resolution across a range of body weight. We identify subpopulations of adipocytes, adipose stem and progenitor cells, vascular and immune cells and demonstrate commonalities and differences across species and dietary conditions. We link specific cell types to increased risk of metabolic disease and provide an initial blueprint for a comprehensive set of interactions between individual cell types in the adipose niche in leanness and obesity. These data comprise an extensive resource for the exploration of genes, traits and cell types in the function of white adipose tissue across species, depots and nutritional conditions.

Highly selective brain-to-gut communication via genetically defined vagus neurons.

The vagus nerve innervates many organs, and most, if not all, of its motor fibers are cholinergic. However, no one knows its organizing principles-whether or not there are dedicated neurons with restricted targets that act as "labeled lines" to perform certain functions, including two opposing ones (gastric contraction versus relaxation). By performing unbiased transcriptional profiling of DMV cholinergic neurons, we discovered seven molecularly distinct subtypes of motor neurons. Then, by using subtype-specific Cre driver mice, we show that two of these subtypes exclusively innervate the glandular domain of the stomach where, remarkably, they contact different enteric neurons releasing functionally opposing neurotransmitters (acetylcholine versus nitric oxide). Thus, the vagus motor nerve communicates via genetically defined labeled lines to control functionally unique enteric neurons within discrete subregions of the gastrointestinal tract. This discovery reveals that the parasympathetic nervous system utilizes a striking division of labor to control autonomic function.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63+ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2.

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Obesity-Linked PPARγ S273 Phosphorylation Promotes Insulin Resistance through Growth Differentiation Factor 3.

The thiazolidinediones (TZDs) are ligands of PPARγ that improve insulin sensitivity, but their use is limited by significant side effects. Recently, we demonstrated a mechanism wherein TZDs improve insulin sensitivity distinct from receptor agonism and adipogenesis: reversal of obesity-linked phosphorylation of PPARγ at serine 273. However, the role of this modification hasn't been tested genetically. Here we demonstrate that mice encoding an allele of PPARγ that cannot be phosphorylated at S273 are protected from insulin resistance, without exhibiting differences in body weight or TZD-associated side effects. Indeed, hyperinsulinemic-euglycemic clamp experiments confirm insulin sensitivity. RNA-seq in these mice reveals reduced expression of Gdf3, a BMP family member. Ectopic expression of Gdf3 is sufficient to induce insulin resistance in lean, healthy mice. We find Gdf3 inhibits BMP signaling and insulin signaling in vitro. Together, these results highlight the diabetogenic role of PPARγ S273 phosphorylation and focus attention on a putative target, Gdf3.

Adipocytes fail to maintain cellular identity during obesity due to reduced PPARγ activity and elevated TGFß-SMAD signaling.

OBJECTIVE: Obesity due to overnutrition causes adipose tissue dysfunction, which is a critical pathological step on the road to type 2 diabetes (T2D) and other metabolic disorders. In this study, we conducted an unbiased investigation into the fundamental molecular mechanisms by which adipocytes transition to an unhealthy state during obesity. METHODS: We used nuclear tagging and translating ribosome affinity purification (NuTRAP) reporter mice crossed with Adipoq-Cre mice to determine adipocyte-specific 1) transcriptional profiles (RNA-seq), 2) promoter and enhancer activity (H3K27ac ChIP-seq), 3) and PPARγ cistrome (ChIP-seq) profiles in mice fed chow or a high-fat diet (HFD) for 10 weeks. We also assessed the impact of the PPARγ agonist rosiglitazone (Rosi) on gene expression and cellular state of adipocytes from the HFD-fed mice. We integrated these data to determine the transcription factors underlying adipocyte responses to HFD and conducted functional studies using shRNA-mediated loss-of-function approaches in 3T3-L1 adipocytes. RESULTS: Adipocytes from the HFD-fed mice exhibited reduced expression of adipocyte markers and metabolic genes and enhanced expression of myofibroblast marker genes involved in cytoskeletal organization, accompanied by the formation of actin filament structures within the cell. PPARγ binding was globally reduced in adipocytes after HFD feeding, and Rosi restored the molecular and cellular phenotypes of adipocytes associated with HFD feeding. We identified the TGFß1 effector protein SMAD to be enriched at HFD-induced promoters and enhancers and associated with myofibroblast signature genes. TGFß1 treatment of mature 3T3-L1 adipocytes induced gene expression and cellular changes similar to those seen after HFD in vivo, and knockdown of Smad3 blunted the effects of TGFß1. CONCLUSIONS: Our data demonstrate that adipocytes fail to maintain cellular identity after HFD feeding, acquiring characteristics of a myofibroblast-like cell type through reduced PPARγ activity and elevated TGFß-SMAD signaling. This cellular identity crisis may be a fundamental mechanism that drives functional decline of adipose tissues during obesity.

Warming Induces Significant Reprogramming of Beige, but Not Brown, Adipocyte Cellular Identity.

Beige and brown adipocytes generate heat in response to reductions in ambient temperature. When warmed, both beige and brown adipocytes exhibit morphological "whitening," but it is unknown whether or to what extent this represents a true shift in cellular identity. Using cell-type-specific profiling in vivo, we uncover a unique paradigm of temperature-dependent epigenomic plasticity of beige, but not brown, adipocytes, with conversion from a brown to a white chromatin state. Despite this profound shift in cellular identity, warm whitened beige adipocytes retain an epigenomic memory of prior cold exposure defined by an array of poised enhancers that prime thermogenic genes for rapid response during a second bout of cold exposure. We further show that a transcriptional cascade involving glucocorticoid receptor and Zfp423 can drive warm-induced whitening of beige adipocytes. These studies identify the epigenomic and transcriptional bases of an extraordinary example of cellular plasticity in response to environmental signals.

Aldosterone-Sensing Neurons in the NTS Exhibit State-Dependent Pacemaker Activity and Drive Sodium Appetite via Synergy with Angiotensin II Signaling.

Sodium deficiency increases angiotensin II (ATII) and aldosterone, which synergistically stimulate sodium retention and consumption. Recently, ATII-responsive neurons in the subfornical organ (SFO) and aldosterone-sensitive neurons in the nucleus of the solitary tract (NTSHSD2 neurons) were shown to drive sodium appetite. Here we investigate the basis for NTSHSD2 neuron activation, identify the circuit by which NTSHSD2 neurons drive appetite, and uncover an interaction between the NTSHSD2 circuit and ATII signaling. NTSHSD2 neurons respond to sodium deficiency with spontaneous pacemaker-like activity-the consequence of "cardiac" HCN and Nav1.5 channels. Remarkably, NTSHSD2 neurons are necessary for sodium appetite, and with concurrent ATII signaling their activity is sufficient to produce rapid consumption. Importantly, NTSHSD2 neurons stimulate appetite via projections to the vlBNST, which is also the effector site for ATII-responsive SFO neurons. The interaction between angiotensin signaling and NTSHSD2 neurons provides a neuronal context for the long-standing "synergy hypothesis" of sodium appetite regulation.

A molecular census of arcuate hypothalamus and median eminence cell types.

The hypothalamic arcuate-median eminence complex (Arc-ME) controls energy balance, fertility and growth through molecularly distinct cell types, many of which remain unknown. To catalog cell types in an unbiased way, we profiled gene expression in 20,921 individual cells in and around the adult mouse Arc-ME using Drop-seq. We identify 50 transcriptionally distinct Arc-ME cell populations, including a rare tanycyte population at the Arc-ME diffusion barrier, a new leptin-sensing neuron population, multiple agouti-related peptide (AgRP) and pro-opiomelanocortin (POMC) subtypes, and an orexigenic somatostatin neuron population. We extended Drop-seq to detect dynamic expression changes across relevant physiological perturbations, revealing cell type-specific responses to energy status, including distinct responses in AgRP and POMC neuron subtypes. Finally, integrating our data with human genome-wide association study data implicates two previously unknown neuron populations in the genetic control of obesity. This resource will accelerate biological discovery by providing insights into molecular and cell type diversity from which function can be inferred.

Simultaneous Transcriptional and Epigenomic Profiling from Specific Cell Types within Heterogeneous Tissues In Vivo.

Epigenomic mechanisms direct distinct gene expression programs for different cell types. Various in vivo tissues have been subjected to epigenomic analysis; however, these studies have been limited by cellular heterogeneity, resulting in composite gene expression and epigenomic profiles. Here, we introduce "NuTRAP," a transgenic mouse that allows simultaneous isolation of cell-type-specific translating mRNA and chromatin from complex tissues. Using NuTRAP, we successfully characterize gene expression and epigenomic states of various adipocyte populations in vivo, revealing significant differences compared to either whole adipose tissue or in vitro adipocyte cell lines. We find that chromatin immunoprecipitation sequencing (ChIP-seq) using NuTRAP is highly efficient, scalable, and robust with even limited cell input. We further demonstrate the general utility of NuTRAP by analyzing hepatocyte-specific epigenomic states. The NuTRAP mouse is a resource that provides a powerful system for cell-type-specific gene expression and epigenomic profiling.

Nuclear Mechanisms of Insulin Resistance.

Insulin resistance is a sine qua non of type 2 diabetes and is associated with many other clinical conditions. Decades of research into mechanisms underlying insulin resistance have mostly focused on problems in insulin signal transduction and other mitochondrial and cytosolic pathways. By contrast, relatively little attention has been focused on transcriptional and epigenetic contributors to insulin resistance, despite strong evidence that such nuclear mechanisms play a major role in the etiopathogenesis of this condition. In this review, we summarize the evidence for nuclear mechanisms of insulin resistance, focusing on three transcription factors with a major impact on insulin action in liver, muscle, and fat.

Identification of nuclear hormone receptor pathways causing insulin resistance by transcriptional and epigenomic analysis.

Insulin resistance is a cardinal feature of Type 2 diabetes (T2D) and a frequent complication of multiple clinical conditions, including obesity, ageing and steroid use, among others. How such a panoply of insults can result in a common phenotype is incompletely understood. Furthermore, very little is known about the transcriptional and epigenetic basis of this disorder, despite evidence that such pathways are likely to play a fundamental role. Here, we compare cell autonomous models of insulin resistance induced by the cytokine tumour necrosis factor-α or by the steroid dexamethasone to construct detailed transcriptional and epigenomic maps associated with cellular insulin resistance. These data predict that the glucocorticoid receptor and vitamin D receptor are common mediators of insulin resistance, which we validate using gain- and loss-of-function studies. These studies define a common transcriptional and epigenomic signature in cellular insulin resistance enabling the identification of pathogenic mechanisms.