Dissecting the brown adipogenic regulatory network using integrative genomics

Vestigial like 4 regulates the adipogenesis of classical brown adipose tissue

Brown adipose tissue (BAT) is mammals' primary non-shivering thermogenesis organ, and the molecular mechanisms regulating BAT growth and adipogenesis are largely unknown. The Hippo-YAP pathway has been well-known for controlling organ size, and Vestigial like 4 (VGLL4) is a transcriptional regulator that modulates the Hippo-YAP pathway by competing against YAP for binding to TEAD proteins. In this study, we dissected the function of VGLL4 in regulating BAT development. We generated a conventional Vgll4 mutant mouse line, in which the two Tondu (TDU) domains of VGLL4 were disrupted. We found that deletion of the TDU domains of VGLL4 resulted in perinatal lethality and paucity of the interscapular BAT. Histological and magnetic resonance imaging studies confirmed that the adipogenesis of BAT was impaired in Vgll4 mutants. Adeno-associated virus (AAV) mediated, brown adipocyte-specific overexpression of VGLL4 increased BAT volume and protected the adult male mice from acute cold stress. Genomic studies suggest that VGLL4/TEAD1 complex directly regulates the myogenic and adipogenic gene expression programs of BAT. In conclusion, our data identify VGLL4 as a previously unrecognized adipogenesis factor that regulates classical BAT development.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

An insight into brown/beige adipose tissue whitening, a metabolic complication of obesity with the multifactorial origin

Brown adipose tissue (BAT), a thermoregulatory organ known to promote energy expenditure, has been extensively studied as a potential avenue to combat obesity. Although BAT is the opposite of white adipose tissue (WAT) which is responsible for energy storage, BAT shares thermogenic capacity with beige adipose tissue that emerges from WAT depots. This is unsurprising as both BAT and beige adipose tissue display a huge difference from WAT in terms of their secretory profile and physiological role. In obesity, the content of BAT and beige adipose tissue declines as these tissues acquire the WAT characteristics via the process called "whitening". This process has been rarely explored for its implication in obesity, whether it contributes to or exacerbates obesity. Emerging research has demonstrated that BAT/beige adipose tissue whitening is a sophisticated metabolic complication of obesity that is linked to multiple factors. The current review provides clarification on the influence of various factors such as diet, age, genetics, thermoneutrality, and chemical exposure on BAT/beige adipose tissue whitening. Moreover, the defects and mechanisms that underpin the whitening are described. Notably, the BAT/beige adipose tissue whitening can be marked by the accumulation of large unilocular lipid droplets, mitochondrial degeneration, and collapsed thermogenic capacity, by the virtue of mitochondrial dysfunction, devascularization, autophagy, and inflammation.

Live-seq enables temporal transcriptomic recording of single cells

Single-cell transcriptomics (scRNA-seq) has greatly advanced our ability to characterize cellular heterogeneity1. However, scRNA-seq requires lysing cells, which impedes further molecular or functional analyses on the same cells. Here, we established Live-seq, a single-cell transcriptome profiling approach that preserves cell viability during RNA extraction using fluidic force microscopy2,3, thus allowing to couple a cell's ground-state transcriptome to its downstream molecular or phenotypic behaviour. To benchmark Live-seq, we used cell growth, functional responses and whole-cell transcriptome read-outs to demonstrate that Live-seq can accurately stratify diverse cell types and states without inducing major cellular perturbations. As a proof of concept, we show that Live-seq can be used to directly map a cell's trajectory by sequentially profiling the transcriptomes of individual macrophages before and after lipopolysaccharide (LPS) stimulation, and of adipose stromal cells pre- and post-differentiation. In addition, we demonstrate that Live-seq can function as a transcriptomic recorder by preregistering the transcriptomes of individual macrophages that were subsequently monitored by time-lapse imaging after LPS exposure. This enabled the unsupervised, genome-wide ranking of genes on the basis of their ability to affect macrophage LPS response heterogeneity, revealing basal Nfkbia expression level and cell cycle state as important phenotypic determinants, which we experimentally validated. Thus, Live-seq can address a broad range of biological questions by transforming scRNA-seq from an end-point to a temporal analysis approach.

Transcriptomic profiling of the telomerase transformed Mesenchymal stromal cells derived adipocytes in response to rosiglitazone

Background

Differentiation of Immortalized Human Bone Marrow Mesenchymal Stromal Cells - hTERT (iMSC3) into adipocytes is in vitro model of obesity. In our earlier study, rosiglitazone enhanced adipogenesis particularly the brown adipogenesis of iMSC3. In this study, the transcriptomic profiles of iMSC3 derived adipocytes with and without rosiglitazone were analyzed through mRNA sequencing.

Results

A total of 1508 genes were differentially expressed between iMSC3 and the derived adipocytes without rosiglitazone treatment. GO and KEGG enrichment analyses revealed that rosiglitazone regulates PPAR and PI3K-Akt pathways. The constant rosiglitazone treatment enhanced the expression of Fatty Acid Binding Protein 4 (FABP4) which enriched GO terms such as fatty acid binding, lipid droplet, as well as white and brown fat cell differentiation. Moreover, the constant treatment upregulated several lipid droplets (LDs) associated proteins such as PLIN1. Rosiglitazone also activated the receptor complex PTK2B that has essential roles in beige adipocytes thermogenic program. Several uniquely expressed novel regulators of brown adipogenesis were also expressed in adipocytes derived with rosiglitazone: PRDM16, ZBTB16, HOXA4, and KLF15 in addition to other uniquely expressed genes.

Conclusions

Rosiglitazone regulated several differentially regulated genes and non-coding RNAs that warrant further investigation about their roles in adipogenesis particularly brown adipogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01027-z.

Loss of cytoplasmic FMR1-interacting protein 2 (CYFIP2) induces browning in 3T3-L1 adipocytes via repression of GABA-BR and activation of mTORC1

Obesity and related metabolic disorders are epidemic diseases. Promoting thermogenesis and a functional increase in the browning of white adipocytes may counteract obesity. On the other hand, the molecular mechanism that regulates brown and beige fat-mediated thermogenesis is unclear. This article reports a molecular network led by cytoplasmic FMR1-interacting protein 2 (CYFIP2) that negatively regulates adipocyte browning in white adipocytes. Although the function of CYFIP2 in Fragile X Syndrome (FXS) and autism have been reported, its physiological roles in adipocytes remain elusive. Therefore, this study examined the physiological consequences of its deprivation in cultured 3T3-L1 white adipocytes using loss-of-function studies. Combined real-time quantitative reverse-transcription polymerase chain reaction and immunoblot analysis showed that the loss of CYFIP2 induces fat browning, as evidenced by the gene and protein expression levels of the brown fat-associated markers. A deficiency of CYFIP2 promoted mitochondrial biogenesis and significantly enhanced the expression of the core set beige fat-specific genes (Cd137, Cidea, Cited1, Tbx1, and Tmem26) and proteins (PGC-1α, PRDM16, and UCP1). In addition, a CYFIP2 deficiency promoted lipid catabolism and suppressed adipogenesis, lipogenesis, and autophagy. A mechanistic study showed that the loss of CYFIP2 induces browning in white adipocytes, independently via the activation of mTORC1 and suppression of the GABA-BR signaling pathway. The present data revealed a previously unidentified mechanism of CYFIP2 in the browning of white adipocytes and emphasized the potential of CYFIP2 as a pharmacotherapeutic target for treating obesity and other metabolic disorders.

A universal framework for detecting cis-regulatory diversity in DNA regulatory regions

High-throughput sequencing-based assays measure different biochemical activities pertaining to gene regulation, genome-wide. These activities include transcription factor (TF)–DNA binding, enhancer activity, open chromatin, and more. A major goal is to understand underlying sequence components, or motifs, that can explain the measured activity. It is usually not one motif but a combination of motifs bound by cooperatively acting proteins that confers activity to such regions. Furthermore, regions can be diverse, governed by different combinations of TFs/motifs. Current approaches do not take into account this issue of combinatorial diversity. We present a new statistical framework, cisDIVERSITY, which models regions as diverse modules characterized by combinations of motifs while simultaneously learning the motifs themselves. Because cisDIVERSITY does not rely on knowledge of motifs, modules, cell type, or organism, it is general enough to be applied to regions reported by most high-throughput assays. For example, in enhancer predictions resulting from different assays—GRO-cap, STARR-seq, and those measuring chromatin structure—cisDIVERSITY discovers distinct modules and combinations of TF binding sites, some specific to the assay. From protein–DNA binding data, cisDIVERSITY identifies potential cofactors of the profiled TF, whereas from ATAC-seq data, it identifies tissue-specific regulatory modules. Finally, analysis of single-cell ATAC-seq data suggests that regions open in one cell-state encode information about future states, with certain modules staying open and others closing down in the next time point.

Systems-Genetics-Based Inference of a Core Regulatory Network Underlying White Fat Browning

Recruitment of brite/beige cells, known as browning of white adipose tissue (WAT), is an efficient way to turn an energy-storing organ into an energy-dissipating one and may therefore be of therapeutic value in combating obesity. However, a comprehensive understanding of the regulatory mechanisms mediating WAT browning is still lacking. Here, we exploit the large natural variation in WAT browning propensity between inbred mouse strains to gain an inclusive view of the core regulatory network coordinating this cellular process. Combining comparative transcriptomics, perturbation-based validations, and gene network analyses, we present a comprehensive gene regulatory network of inguinal WAT browning, revealing up to four distinct regulatory modules with key roles for uncovered transcriptional factors, while also providing deep insights into the genetic architecture of brite adipogenesis. The presented findings therefore greatly increase our understanding of the molecular drivers mediating the intriguing cellular heterogeneity and plasticity of adipose tissue.

Homeobox A4 suppresses vascular remodeling by repressing YAP/TEAD transcriptional activity

The Hippo signaling pathway is involved in the pathophysiology of various cardiovascular diseases. Yes-associated protein (YAP) and transcriptional enhancer activator domain (TEAD) transcriptional factors, the main transcriptional complex of the Hippo pathway, were recently identified as modulators of phenotypic switching of vascular smooth muscle cells (VSMCs). However, the intrinsic regulator of YAP/TEAD-mediated gene expressions involved in vascular pathophysiology remains to be elucidated. Here, we identified Homeobox A4 (HOXA4) as a potent repressor of YAP/TEAD transcriptional activity using lentiviral shRNA screen. Mechanistically, HOXA4 interacts with TEADs and attenuates YAP/TEAD-mediated transcription by competing with YAP for TEAD binding. We also clarified that the expression of HOXA4 is relatively abundant in the vasculature, especially in VSMCs. In vitro experiments in human VSMCs showed HOXA4 maintains the differentiation state of VSMCs via inhibition of YAP/TEAD-induced phenotypic switching. We generated Hoxa4-deficient mice and confirmed the downregulation of smooth muscle-specific contractile genes and the exacerbation of vascular remodeling after carotid artery ligation in vivo. Our results demonstrate that HOXA4 is a repressor of VSMC phenotypic switching by inhibiting YAP/TEAD-mediated transcription.

Comparative Analyses of Chromatin Landscape in White Adipose Tissue Suggest Humans May Have Less Beigeing Potential than Other Primates

Humans carry a much larger percentage of body fat than other primates. Despite the central role of adipose tissue in metabolism, little is known about the evolution of white adipose tissue in primates. Phenotypic divergence is often caused by genetic divergence in cis-regulatory regions. We examined the cis-regulatory landscape of fat during human origins by performing comparative analyses of chromatin accessibility in human and chimpanzee adipose tissue using rhesus macaque as an outgroup. We find that many regions that have decreased accessibility in humans are enriched for promoter and enhancer sequences, are depleted for signatures of negative selection, are located near genes involved with lipid metabolism, and contain a short sequence motif involved in the beigeing of fat, the process in which lipid-storing white adipocytes are transdifferentiated into thermogenic beige adipocytes. The collective closing of many putative regulatory regions associated with beigeing of fat suggests a mechanism that increases body fat in humans.

BRB-seq: ultra-affordable high-throughput transcriptomics enabled by bulk RNA barcoding and sequencing

Despite its widespread use, RNA-seq is still too laborious and expensive to replace RT-qPCR as the default gene expression analysis method. We present a novel approach, BRB-seq, which uses early multiplexing to produce 3' cDNA libraries for dozens of samples, requiring just 2 hours of hands-on time. BRB-seq has a comparable performance to the standard TruSeq approach while showing greater tolerance for lower RNA quality and being up to 25 times cheaper. We anticipate that BRB-seq will transform basic laboratory practice given its capacity to generate genome-wide transcriptomic data at a similar cost as profiling four genes using RT-qPCR.

ZFP30 promotes adipogenesis through the KAP1-mediated activation of a retrotransposon-derived Pparg2 enhancer

Krüppel-associated box zinc finger proteins (KZFPs) constitute the largest family of mammalian transcription factors, but most remain completely uncharacterized. While initially proposed to primarily repress transposable elements, recent reports have revealed that KFZPs contribute to a wide variety of other biological processes. Using murine and human in vitro and in vivo models, we demonstrate here that one poorly studied KZFP, ZFP30, promotes adipogenesis by directly targeting and activating a retrotransposon-derived Pparg2 enhancer. Through mechanistic studies, we further show that ZFP30 recruits the co-regulator KRAB-associated protein 1 (KAP1), which, surprisingly, acts as a ZFP30 co-activator in this adipogenic context. Our findings provide an understanding of both adipogenic and KZFP-KAP1 complex-mediated gene regulation, showing that the KZFP-KAP1 axis can also function in a non-repressive manner.

Temporal enhancer profiling of parallel lineages identifies AHR and GLIS1 as regulators of mesenchymal multipotency

Temporal data on gene expression and context-specific open chromatin states can improve identification of key transcription factors (TFs) and the gene regulatory networks (GRNs) controlling cellular differentiation. However, their integration remains challenging. Here, we delineate a general approach for data-driven and unbiased identification of key TFs and dynamic GRNs, called EPIC-DREM. We generated time-series transcriptomic and epigenomic profiles during differentiation of mouse multipotent bone marrow stromal cell line (ST2) toward adipocytes and osteoblasts. Using our novel approach we constructed time-resolved GRNs for both lineages and identifed the shared TFs involved in both differentiation processes. To take an alternative approach to prioritize the identified shared regulators, we mapped dynamic super-enhancers in both lineages and associated them to target genes with correlated expression profiles. The combination of the two approaches identified aryl hydrocarbon receptor (AHR) and Glis family zinc finger 1 (GLIS1) as mesenchymal key TFs controlled by dynamic cell type-specific super-enhancers that become repressed in both lineages. AHR and GLIS1 control differentiation-induced genes and their overexpression can inhibit the lineage commitment of the multipotent bone marrow-derived ST2 cells.

Reversible De-differentiation of Mature White Adipocytes into Preadipocyte-like Precursors during Lactation

Adipose tissue in the mammary gland undergoes dramatic remodeling during reproduction. Adipocytes are replaced by mammary alveolar structures during pregnancy and lactation, then reappear upon weaning. The fate of the original adipocytes during lactation and the developmental origin of the re-appearing adipocyte post involution are unclear. Here, we reveal that adipocytes in the mammary gland de-differentiate into Pdgfrα⁺ preadipocyte- and fibroblast-like cells during pregnancy and remain de-differentiated during lactation. Upon weaning, de-differentiated fibroblasts proliferate and re-differentiate into adipocytes. This cycle occurs over multiple pregnancies. These observations reveal the potential of terminally differentiated adipocytes to undergo repeated cycles of de-differentiation and re-differentiation in a physiological setting.

A stromal cell population that inhibits adipogenesis in mammalian fat depots

Adipocyte development and differentiation have an important role in the aetiology of obesity and its co-morbidities^1,2. Although multiple studies have investigated the adipogenic stem and precursor cells that give rise to mature adipocytes^3-14, our understanding of their in vivo origin and properties is incomplete^2,15,16. This is partially due to the highly heterogeneous and unstructured nature of adipose tissue depots¹⁷, which has proven difficult to molecularly dissect using classical approaches such as fluorescence-activated cell sorting and Cre-lox lines based on candidate marker genes^16,18. Here, using the resolving power of single-cell transcriptomics¹⁹ in a mouse model, we reveal distinct subpopulations of adipose stem and precursor cells in the stromal vascular fraction of subcutaneous adipose tissue. We identify one of these subpopulations as CD142⁺ adipogenesis-regulatory cells, which can suppress adipocyte formation in vivo and in vitro in a paracrine manner. We show that adipogenesis-regulatory cells are refractory to adipogenesis and that they are functionally conserved in humans. Our findings point to a potentially critical role for adipogenesis-regulatory cells in modulating adipose tissue plasticity, which is linked to metabolic control, differential insulin sensitivity and type 2 diabetes.

Epigenomic Control of Thermogenic Adipocyte Differentiation and Function

Obesity and its associated metabolic disorders are spreading at a fast pace throughout the world; thus, effective therapeutic approaches are necessary to combat this epidemic. Obesity develops when there is a greater caloric intake than energy expenditure. Promoting energy expenditure has recently attracted much attention as a promising approach for the management of body weight. Thermogenic adipocytes are capable of burning fat to dissipate chemical energy into heat, thereby enhancing energy expenditure. After the recent re-discovery of thermogenic adipocytes in adult humans, much effort has focused on understanding the molecular mechanisms, especially the epigenetic mechanisms, which regulate thermogenic adipocyte development and function. A number of chromatin signatures, such as histone modifications, DNA methylation, chromatin accessibilities, and interactions, have been profiled at the genome level and analyzed in various murine and human thermogenic fat cell systems. Moreover, writers and erasers, as well as readers of the epigenome are also investigated using genomic tools in thermogenic adipocytes. In this review, we summarize and discuss the recent advance in these studies and highlight the insights gained into the epigenomic regulation of thermogenic program as well as the pathogenesis of human metabolic diseases.

NFIA co-localizes with PPARγ and transcriptionally controls the brown fat gene program

Brown fat dissipates energy as heat and protects against obesity. Here, we identified nuclear factor I-A (NFIA) as a transcriptional regulator of brown fat by a genome-wide open chromatin analysis of murine brown and white fat followed by motif analysis of brown-fat-specific open chromatin regions. NFIA and the master transcriptional regulator of adipogenesis, PPARγ, co-localize at the brown-fat-specific enhancers. Moreover, the binding of NFIA precedes and facilitates the binding of PPARγ, leading to increased chromatin accessibility and active transcription. Introduction of NFIA into myoblasts results in brown adipocyte differentiation. Conversely, the brown fat of NFIA-knockout mice displays impaired expression of the brown-fat-specific genes and reciprocal elevation of muscle genes. Finally, expression of NFIA and the brown-fat-specific genes is positively correlated in human brown fat. These results indicate that NFIA activates the cell-type-specific enhancers and facilitates the binding of PPARγ to control the brown fat gene program.