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.

Creatine kinase B controls futile creatine cycling in thermogenic fat.

Obesity increases the risk of mortality because of metabolic sequelae such as type 2 diabetes and cardiovascular disease1. Thermogenesis by adipocytes can counteract obesity and metabolic diseases2,3. In thermogenic fat, creatine liberates a molar excess of mitochondrial ADP-purportedly via a phosphorylation cycle4-to drive thermogenic respiration. However, the proteins that control this futile creatine cycle are unknown. Here we show that creatine kinase B (CKB) is indispensable for thermogenesis resulting from the futile creatine cycle, during which it traffics to mitochondria using an internal mitochondrial targeting sequence. CKB is powerfully induced by thermogenic stimuli in both mouse and human adipocytes. Adipocyte-selective inactivation of Ckb in mice diminishes thermogenic capacity, increases predisposition to obesity, and disrupts glucose homeostasis. CKB is therefore a key effector of the futile creatine cycle.

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.

NCoR1-independent mechanism plays a role in the action of the unliganded thyroid hormone receptor.

Nuclear receptor corepressor 1 (NCoR1) is considered to be the major corepressor that mediates ligand-independent actions of the thyroid hormone receptor (TR) during development and in hypothyroidism. We tested this by expressing a hypomorphic NCoR1 allele (NCoR1ΔID), which cannot interact with the TR, in Pax8-KO mice, which make no thyroid hormone. Surprisingly, abrogation of NCoR1 function did not reverse the ligand-independent action of the TR on many gene targets and did not fully rescue the high mortality rate due to congenital hypothyroidism in these mice. To further examine NCoR1's role in repression by the unliganded TR, we deleted NCoR1 in the livers of euthyroid and hypothyroid mice and examined the effects on gene expression and enhancer activity measured by histone 3 lysine 27 (H3K27) acetylation. Even in the absence of NCoR1 function, we observed strong repression of more than 43% of positive T3 (3,3',5-triiodothyronine) targets in hypothyroid mice. Regulation of approximately half of those genes correlated with decreased H3K27 acetylation, and nearly 80% of these regions with affected H3K27 acetylation contained a bona fide TRß1-binding site. Moreover, using liver-specific TRß1-KO mice, we demonstrate that hypothyroidism-associated changes in gene expression and histone acetylation require TRß1. Thus, many of the genomic changes mediated by the TR in hypothyroidism are independent of NCoR1, suggesting a role for additional signaling modulators in hypothyroidism.

Charting a dynamic DNA methylation landscape of the human genome.

DNA methylation is a defining feature of mammalian cellular identity and is essential for normal development. Most cell types, except germ cells and pre-implantation embryos, display relatively stable DNA methylation patterns, with 70-80% of all CpGs being methylated. Despite recent advances, we still have a limited understanding of when, where and how many CpGs participate in genomic regulation. Here we report the in-depth analysis of 42 whole-genome bisulphite sequencing data sets across 30 diverse human cell and tissue types. We observe dynamic regulation for only 21.8% of autosomal CpGs within a normal developmental context, most of which are distal to transcription start sites. These dynamic CpGs co-localize with gene regulatory elements, particularly enhancers and transcription-factor-binding sites, which allow identification of key lineage-specific regulators. In addition, differentially methylated regions (DMRs) often contain single nucleotide polymorphisms associated with cell-type-related diseases as determined by genome-wide association studies. The results also highlight the general inefficiency of whole-genome bisulphite sequencing, as 70-80% of the sequencing reads across these data sets provided little or no relevant information about CpG methylation. To demonstrate further the utility of our DMR set, we use it to classify unknown samples and identify representative signature regions that recapitulate major DNA methylation dynamics. In summary, although in theory every CpG can change its methylation state, our results suggest that only a fraction does so as part of coordinated regulatory programs. Therefore, our selected DMRs can serve as a starting point to guide new, more effective reduced representation approaches to capture the most informative fraction of CpGs, as well as further pinpoint putative regulatory elements.

Novel mechanism of positive versus negative regulation by thyroid hormone receptor ß1 (TRß1) identified by genome-wide profiling of binding sites in mouse liver.

Triiodothyronine (T3) regulates key metabolic processes in the liver through the thyroid hormone receptor, TRß1. However, the number of known target genes directly regulated by TRß1 is limited, and the mechanisms by which positive and especially negative transcriptional regulation occur are not well understood. To characterize the TRß1 cistrome in vivo, we expressed a biotinylated TRß1 in hypo- and hyperthyroid mouse livers, used ChIP-seq to identify genomic TRß1 targets, and correlated these data with gene expression changes. As with other nuclear receptors, the majority of TRß1 binding sites were not in proximal promoters but in the gene body of known genes. Remarkably, T3 can dictate changes in TRß1 binding, with strong correlation to T3-induced gene expression changes, suggesting that differential TRß1 binding regulates transcriptional outcome. Additionally, DR-4 and DR-0 motifs were significantly enriched at binding sites where T3 induced an increase or decrease in TRß1 binding, respectively, leading to either positive or negative regulation by T3. Taken together, the results of this study provide new insights into the mechanisms of transcriptional regulation by TRß1 in vivo.

Aging impairs cold-induced beige adipogenesis and adipocyte metabolic reprogramming.

The energy-burning capability of beige adipose tissue is a potential therapeutic tool for reducing obesity and metabolic disease, but this capacity is decreased by aging. Here, we evaluate the impact of aging on the profile and activity of adipocyte stem and progenitor cells (ASPCs) and adipocytes during the beiging process. We found that aging increases the expression of Cd9 and other fibro-inflammatory genes in fibroblastic ASPCs and blocks their differentiation into beige adipocytes. Fibroblastic ASPC populations from young and aged mice were equally competent for beige differentiation in vitro, suggesting that environmental factors suppress adipogenesis in vivo. Examination of adipocytes by single nucleus RNA-sequencing identified compositional and transcriptional differences in adipocyte populations with age and cold exposure. Notably, cold exposure induced an adipocyte population expressing high levels of de novo lipogenesis (DNL) genes, and this response was severely blunted in aged animals. We further identified natriuretic peptide clearance receptor Npr3, a beige fat repressor, as a marker gene for a subset of white adipocytes and an aging-upregulated gene in adipocytes. In summary, this study indicates that aging blocks beige adipogenesis and dysregulates adipocyte responses to cold exposure and provides a unique resource for identifying cold and aging-regulated pathways in adipose tissue.

Aging impairs cold-induced beige adipogenesis and adipocyte metabolic reprogramming.

The energy-burning capability of beige adipose tissue is a potential therapeutic tool for reducing obesity and metabolic disease, but this capacity is decreased by aging. Here, we evaluate the impact of aging on the profile and activity of adipocyte stem and progenitor cells (ASPCs) and adipocytes during the beiging process in mice. We found that aging increases the expression of Cd9 and other fibro-inflammatory genes in fibroblastic ASPCs and blocks their differentiation into beige adipocytes. Fibroblastic ASPC populations from young and aged mice were equally competent for beige differentiation in vitro, suggesting that environmental factors suppress adipogenesis in vivo. Examination of adipocytes by single nucleus RNA-sequencing identified compositional and transcriptional differences in adipocyte populations with aging and cold exposure. Notably, cold exposure induced an adipocyte population expressing high levels of de novo lipogenesis (DNL) genes, and this response was severely blunted in aged animals. We further identified Npr3, which encodes the natriuretic peptide clearance receptor, as a marker gene for a subset of white adipocytes and an aging-upregulated gene in adipocytes. In summary, this study indicates that aging blocks beige adipogenesis and dysregulates adipocyte responses to cold exposure and provides a resource for identifying cold and aging-regulated pathways in adipose tissue.

AgRP neuron cis-regulatory analysis across hunger states reveals that IRF3 mediates leptin's acute effects.

AgRP neurons in the arcuate nucleus of the hypothalamus (ARC) coordinate homeostatic changes in appetite associated with fluctuations in food availability and leptin signaling. Identifying the relevant transcriptional regulatory pathways in these neurons has been a priority, yet such attempts have been stymied due to their low abundance and the rich cellular diversity of the ARC. Here we generated AgRP neuron-specific transcriptomic and chromatin accessibility profiles from male mice during three distinct hunger states of satiety, fasting-induced hunger, and leptin-induced hunger suppression. Cis-regulatory analysis of these integrated datasets enabled the identification of 18 putative hunger-promoting and 29 putative hunger-suppressing transcriptional regulators in AgRP neurons, 16 of which were predicted to be transcriptional effectors of leptin. Within our dataset, Interferon regulatory factor 3 (IRF3) emerged as a leading candidate mediator of leptin-induced hunger-suppression. Measures of IRF3 activation in vitro and in vivo reveal an increase in IRF3 nuclear occupancy following leptin administration. Finally, gain- and loss-of-function experiments in vivo confirm the role of IRF3 in mediating the acute satiety-evoking effects of leptin in AgRP neurons. Thus, our findings identify IRF3 as a key mediator of the acute hunger-suppressing effects of leptin in AgRP neurons.

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.

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.

Hypoxia-inducible factor 1 signaling drives placental aging and can provoke preterm labor.

Most cases of preterm labor have unknown cause, and the burden of preterm birth is immense. Placental aging has been proposed to promote labor onset, but specific mechanisms remain elusive. We report findings stemming from unbiased transcriptomic analysis of mouse placenta, which revealed that hypoxia-inducible factor 1 (HIF-1) stabilization is a hallmark of advanced gestational timepoints, accompanied by mitochondrial dysregulation and cellular senescence; we detected similar effects in aging human placenta. In parallel in primary mouse trophoblasts and human choriocarcinoma cells, we modeled HIF-1 induction and demonstrated resultant mitochondrial dysfunction and cellular senescence. Transcriptomic analysis revealed that HIF-1 stabilization recapitulated gene signatures observed in aged placenta. Further, conditioned media from trophoblasts following HIF-1 induction promoted contractility in immortalized uterine myocytes, suggesting a mechanism by which the aging placenta may drive the transition from uterine quiescence to contractility at the onset of labor. Finally, pharmacological induction of HIF-1 via intraperitoneal administration of dimethyloxalyl glycine (DMOG) to pregnant mice caused preterm labor. These results provide clear evidence for placental aging in normal pregnancy, and demonstrate how HIF-1 signaling in late gestation may be a causal determinant of the mitochondrial dysfunction and senescence observed within the trophoblast as well as a trigger for uterine contraction.

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.

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.

Drosophila tao controls mushroom body development and ethanol-stimulated behavior through par-1.

In both mammalian and insect models of ethanol-induced behavior, low doses of ethanol stimulate locomotion. However, the mechanisms of the stimulant effects of ethanol on the CNS are mostly unknown. We have identified tao, encoding a serine-threonine kinase of the Ste20 family, as a gene necessary for ethanol-induced locomotor hyperactivity in Drosophila. Mutations in tao also affect behavioral responses to cocaine and nicotine, making flies resistant to the effects of both drugs. We show that tao function is required during the development of the adult nervous system and that tao mutations cause defects in the development of central brain structures, including the mushroom body. Silencing of a subset of mushroom body neurons is sufficient to reduce ethanol-induced hyperactivity, revealing the mushroom body as an important locus mediating the stimulant effects of ethanol. We also show that mutations in par-1 suppress both the mushroom body morphology and behavioral phenotypes of tao mutations and that the phosphorylation state of the microtubule-binding protein Tau can be altered by RNA interference knockdown of tao, suggesting that tao and par-1 act in a pathway to control microtubule dynamics during neural development.

Genetic analysis of dietary intake identifies new loci and functional links with metabolic traits.

Dietary intake is a major contributor to the global obesity epidemic and represents a complex behavioural phenotype that is partially affected by innate biological differences. Here, we present a multivariate genome-wide association analysis of overall variation in dietary intake to account for the correlation between dietary carbohydrate, fat and protein in 282,271 participants of European ancestry from the UK Biobank (n = 191,157) and Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n = 91,114), and identify 26 distinct genome-wide significant loci. Dietary intake signals map exclusively to specific brain regions and are enriched for genes expressed in specialized subtypes of GABAergic, dopaminergic and glutamatergic neurons. We identified two main clusters of genetic variants for overall variation in dietary intake that were differently associated with obesity and coronary artery disease. These results enhance the biological understanding of interindividual differences in dietary intake by highlighting neural mechanisms, supporting functional follow-up experiments and possibly providing new avenues for the prevention and treatment of prevalent complex metabolic diseases.