Long-range chromosomal interactions increase and mark repressed gene expression during adipogenesis

Obesity perturbs central functions of human adipose tissue, centred on differentiation of preadipocytes to adipocytes, i.e., adipogenesis. The large environmental component of obesity makes it important to elucidate epigenetic regulatory factors impacting adipogenesis. Promoter Capture Hi-C (pCHi-C) has been used to identify chromosomal interactions between promoters and associated regulatory elements. However, long range interactions (LRIs) greater than 1 Mb are often filtered out of pCHi-C datasets, due to technical challenges and their low prevalence. To elucidate the unknown role of LRIs in adipogenesis, we investigated preadipocyte differentiation to adipocytes using pCHi-C and bulk and single nucleus RNA-seq data. We first show that LRIs are reproducible between biological replicates, and they increase >2-fold in frequency across adipogenesis. We further demonstrate that genomic loci containing LRIs are more epigenetically repressed than regions without LRIs, corresponding to lower gene expression in the LRI regions. Accordingly, as preadipocytes differentiate into adipocytes, LRI regions are more likely to contain repressed preadipocyte marker genes; whereas these same LRI regions are depleted of actively expressed adipocyte marker genes. Finally, we show that LRIs can be used to restrict multiple testing of the long-range cis-eQTL analysis to identify variants that regulate genes via LRIs. We exemplify this by identifying a putative long range cis regulatory mechanism at the LYPLAL1/TGFB2 obesity locus. In summary, we identify LRIs that mark repressed regions of the genome, and these interactions increase across adipogenesis, pinpointing developmental regions that need to be repressed in a cell-type specific way for adipogenesis to proceed.Peer reviewe

Identification of TBX15 as an adipose master trans regulator of abdominal obesity genes

Background: Obesity predisposes individuals to multiple cardiometabolic disorders, including type 2 diabetes (T2D). As body mass index (BMI) cannot reliably differentiate fat from lean mass, the metabolically detrimental abdominal obesity has been estimated using waist-hip ratio (WHR). Waist-hip ratio adjusted for body mass index (WHRadjBMI) in turn is a well-established sex-specific marker for abdominal fat and adiposity, and a predictor of adverse metabolic outcomes, such as T2D. However, the underlying genes and regulatory mechanisms orchestrating the sex differences in obesity and body fat distribution in humans are not well understood. Methods: We searched for genetic master regulators of WHRadjBMI by employing integrative genomics approaches on human subcutaneous adipose RNA sequencing (RNA-seq) data (n similar to 1400) and WHRadjBMI GWAS data (n similar to 700,000) from the WHRadjBMI GWAS cohorts and the UK Biobank (UKB), using co-expression network, transcriptome-wide association study (TWAS), and polygenic risk score (PRS) approaches. Finally, we functionally verified our genomic results using gene knockdown experiments in a human primary cell type that is critical for adipose tissue function. Results: Here, we identified an adipose gene co-expression network that contains 35 obesity GWAS genes and explains a significant amount of polygenic risk for abdominal obesity and T2D in the UKB (n = 392,551) in a sex-dependent way. We showed that this network is preserved in the adipose tissue data from the Finnish Kuopio Obesity Study and Mexican Obesity Study. The network is controlled by a novel adipose master transcription factor (TF), TBX15, a WHRadjBMI GWAS gene that regulates the network in trans. Knockdown of TBX15 in human primary preadipocytes resulted in changes in expression of 130 network genes, including the key adipose TFs, PPARG and KLF15, which were significantly impacted (FDR < 0.05), thus functionally verifying the trans regulatory effect of TBX15 on the WHRadjBMI co-expression network. Conclusions: Our study discovers a novel key function for the TBX15 TF in trans regulating an adipose co-expression network of 347 adipose, mitochondrial, and metabolically important genes, including PPARG, KLF15, PPARA, ADIPOQ, and 35 obesity GWAS genes. Thus, based on our converging genomic, transcriptional, and functional evidence, we interpret the role of TBX15 to be a main transcriptional regulator in the adipose tissue and discover its importance in human abdominal obesity.Peer reviewe

Reverse gene-environment interaction approach to identify variants influencing body-mass index in humans.

Identifying gene-environment interactions (GxEs) contributing to human cardiometabolic disorders is challenging. Here we apply a reverse GxE candidate search by deriving candidate variants from promoter-enhancer interactions that respond to dietary fatty acid challenge through altered chromatin accessibility in human primary adipocytes. We then test all variants residing in the lipid-responsive open chromatin sites within adipocyte promoter-enhancer contacts for interaction effects between the genotype and dietary saturated fat intake on body mass index (BMI) in the UK Biobank. We discover 14 novel GxE variants in 12 lipid-responsive promoters, including well-known lipid genes (LIPE, CARM1, and PLIN2) and novel genes, such as LDB3, for which we provide further functional and integrative genomics evidence. We further identify 24 GxE variants in enhancers, totaling 38 new GxE variants for BMI in the UK Biobank, demonstrating that molecular genomics data produced in physiologically relevant contexts can discover new functional GxE mechanisms in humans

Increased body mass index is linked to systemic inflammation through altered chromatin co-accessibility in human preadipocytes

Obesity-induced adipose tissue dysfunction can cause low-grade inflammation and downstream obesity comorbidities. Although preadipocytes may contribute to this pro-inflammatory environment, the underlying mechanisms are unclear. We used human primary preadipocytes from body mass index (BMI)-discordant monozygotic (MZ) twin pairs to generate epigenetic (ATAC-sequence) and transcriptomic (RNA-sequence) data for testing whether increased BMI alters the subnuclear compartmentalization of open chromatin in the twins' preadipocytes, causing downstream inflammation. Here we show that the co-accessibility of open chromatin, i.e. compartmentalization of chromatin activity, is altered in the higher vs lower BMI MZ siblings for a large subset (similar to 88.5 Mb) of the active subnuclear compartments. Using the UK Biobank we show that variants within these regions contribute to systemic inflammation through interactions with BMI on C-reactive protein. In summary, open chromatin co-accessibility in human preadipocytes is disrupted among the higher BMI siblings, suggesting a mechanism how obesity may lead to inflammation via gene-environment interactions.Peer reviewe

Reverse gene-environment interaction approach to identify variants influencing body-mass index in humans.

Identifying gene-environment interactions (GxEs) contributing to human cardiometabolic disorders is challenging. Here we apply a reverse GxE candidate search by deriving candidate variants from promoter-enhancer interactions that respond to dietary fatty acid challenge through altered chromatin accessibility in human primary adipocytes. We then test all variants residing in the lipid-responsive open chromatin sites within adipocyte promoter-enhancer contacts for interaction effects between the genotype and dietary saturated fat intake on body mass index (BMI) in the UK Biobank. We discover 14 novel GxE variants in 12 lipid-responsive promoters, including well-known lipid genes (LIPE, CARM1, and PLIN2) and novel genes, such as LDB3, for which we provide further functional and integrative genomics evidence. We further identify 24 GxE variants in enhancers, totaling 38 new GxE variants for BMI in the UK Biobank, demonstrating that molecular genomics data produced in physiologically relevant contexts can discover new functional GxE mechanisms in humans