Open Chromatin Profiling in Adipose Tissue Marks Genomic Regions with Functional Roles in Cardiometabolic Traits

Inhibiting sorting nexin 10 promotes mucosal healing through SREBP2-mediated stemness restoration of intestinal stem cells

Intestinal stem cell (ISC) is a promising therapeutic target for inflammatory bowel disease. Cholesterol availability is critical for ISC stemness. Low plasma cholesterol is a typical feature of Crohn's disease (CD); however, its impact on mucosal healing remains unclear. Here, we identified an essential role of sorting nexin 10 (SNX10) in maintaining the stemness of ISCs. SNX10 expression in intestinal tissues positively correlates with the severity of human CD and mouse colitis. Conditional SNX10 knockout in intestinal epithelial cells or ISCs promotes intestinal mucosal repair by maintaining the ISC population associated with increased intracellular cholesterol synthesis. Disassociation of ERLIN2 with SCAP by SNX10 deletion enhances the activation of SREBP2, resulting in increased cholesterol biosynthesis. DC-SX029, a small-molecule inhibitor of SNX10, was used to verify the druggable potential of SNX10 for the treatment of patients with CD. Our study provides a strategy for mucosal healing through SREBP2-mediated stemness restoration of ISCs.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population

BACKGROUND

Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource.

RESULTS

We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue.

CONCLUSIONS

By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Genetics of sexually dimorphic adipose distribution in humans

Obesity-associated morbidity is exacerbated by abdominal obesity, which can be measured as the waist-to-hip ratio adjusted for the body mass index (WHRadjBMI). Here we identify genes associated with obesity and WHRadjBMI and characterize allele-sensitive enhancers that are predicted to regulate WHRadjBMI genes in women. We found that several waist-to-hip ratio-associated variants map within primate-specific Alu retrotransposons harboring a DNA motif associated with adipocyte differentiation. This suggests that a genetic component of adipose distribution in humans may involve co-option of retrotransposons as adipose enhancers. We evaluated the role of the strongest female WHRadjBMI-associated gene, SNX10, in adipose biology. We determined that it is required for human adipocyte differentiation and function and participates in diet-induced adipose expansion in female mice, but not males. Our data identify genes and regulatory mechanisms that underlie female-specific adipose distribution and mediate metabolic dysfunction in women.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Optimization of the Omni-ATAC protocol to chromatin accessibility profiling in snap-frozen rat adipose and muscle tissues

ATAC-seq is a fast and sensitive method for the epigenomic profiling of open chromatin and for mapping of transcription factor binding sites [1]. Despite the development of the Omni-ATAC protocol for the profiling of chromatin accessibility in frozen tissues [2], studies in adipose tissue have been restricted due to technical challenges including the high lipid content of adipocytes and reproducibility issues between replicates. Here, we provide a modified Omni-ATAC protocol that achieves high data reproducibility in various tissue types from rat, including adipose and muscle tissues [3].•This protocol describes a methodology that enables chromatin accessibility profiling from snap-frozen rat adipose and muscle tissues.•The technique comprises an optimized bead-based tissue homogenization process that substitutes to Dounce homogenization, reduces variability in the experimental procedure, and is adaptable to various tissue types.•In comparison with the Omni-ATAC protocol, the method described here results in improved ATAC-seq data quality that complies with ENCODE quality standards.

Human and rat skeletal muscle single-nuclei multi-omic integrative analyses nominate causal cell types, regulatory elements, and SNPs for complex traits

Skeletal muscle accounts for the largest proportion of human body mass, on average, and is a key tissue in complex diseases and mobility. It is composed of several different cell and muscle fiber types. Here, we optimize single-nucleus ATAC-seq (snATAC-seq) to map skeletal muscle cell-specific chromatin accessibility landscapes in frozen human and rat samples, and single-nucleus RNA-seq (snRNA-seq) to map cell-specific transcriptomes in human. We additionally perform multi-omics profiling (gene expression and chromatin accessibility) on human and rat muscle samples. We capture type I and type II muscle fiber signatures, which are generally missed by existing single-cell RNA-seq methods. We perform cross-modality and cross-species integrative analyses on 33,862 nuclei and identify seven cell types ranging in abundance from 59.6% to 1.0% of all nuclei. We introduce a regression-based approach to infer cell types by comparing transcription start site-distal ATAC-seq peaks to reference enhancer maps and show consistency with RNA-based marker gene cell type assignments. We find heterogeneity in enrichment of genetic variants linked to complex phenotypes from the UK Biobank and diabetes genome-wide association studies in cell-specific ATAC-seq peaks, with the most striking enrichment patterns in muscle mesenchymal stem cells (∼3.5% of nuclei). Finally, we overlay these chromatin accessibility maps on GWAS data to nominate causal cell types, SNPs, transcription factor motifs, and target genes for type 2 diabetes signals. These chromatin accessibility profiles for human and rat skeletal muscle cell types are a useful resource for nominating causal GWAS SNPs and cell types.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Addiction-Associated Genetic Variants Implicate Brain Cell Type- and Region-Specific Cis-Regulatory Elements in Addiction Neurobiology

Recent large genome-wide association studies have identified multiple confident risk loci linked to addiction-associated behavioral traits. Most genetic variants linked to addiction-associated traits lie in noncoding regions of the genome, likely disrupting cis-regulatory element (CRE) function. CREs tend to be highly cell type-specific and may contribute to the functional development of the neural circuits underlying addiction. Yet, a systematic approach for predicting the impact of risk variants on the CREs of specific cell populations is lacking. To dissect the cell types and brain regions underlying addiction-associated traits, we applied stratified linkage disequilibrium score regression to compare genome-wide association studies to genomic regions collected from human and mouse assays for open chromatin, which is associated with CRE activity. We found enrichment of addiction-associated variants in putative CREs marked by open chromatin in neuronal (NeuN+) nuclei collected from multiple prefrontal cortical areas and striatal regions known to play major roles in reward and addiction. To further dissect the cell type-specific basis of addiction-associated traits, we also identified enrichments in human orthologs of open chromatin regions of female and male mouse neuronal subtypes: cortical excitatory, D1, D2, and PV. Last, we developed machine learning models to predict mouse cell type-specific open chromatin, enabling us to further categorize human NeuN+ open chromatin regions into cortical excitatory or striatal D1 and D2 neurons and predict the functional impact of addiction-associated genetic variants. Our results suggest that different neuronal subtypes within the reward system play distinct roles in the variety of traits that contribute to addiction.SIGNIFICANCE STATEMENT We combine statistical genetic and machine learning techniques to find that the predisposition to for nicotine, alcohol, and cannabis use behaviors can be partially explained by genetic variants in conserved regulatory elements within specific brain regions and neuronal subtypes of the reward system. Our computational framework can flexibly integrate open chromatin data across species to screen for putative causal variants in a cell type- and tissue-specific manner for numerous complex traits.

Chromatin accessibility and gene expression during adipocyte differentiation identify context-dependent effects at cardiometabolic GWAS loci

Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and mechanisms at genome-wide association study (GWAS) loci. To identify regulatory elements that display differential activity across adipocyte differentiation, we performed ATAC-seq and RNA-seq in a human cell model of preadipocytes and adipocytes at days 4 and 14 of differentiation. For comparison, we created a consensus map of ATAC-seq peaks in 11 human subcutaneous adipose tissue samples. We identified 58,387 context-dependent chromatin accessibility peaks and 3,090 context-dependent genes between all timepoint comparisons (log2 fold change>1, FDR<5%) with 15,919 adipocyte- and 18,244 preadipocyte-dependent peaks. Adipocyte-dependent peaks showed increased overlap (60.1%) with Roadmap Epigenomics adipocyte nuclei enhancers compared to preadipocyte-dependent peaks (11.5%). We linked context-dependent peaks to genes based on adipocyte promoter capture Hi-C data, overlap with adipose eQTL variants, and context-dependent gene expression. Of 16,167 context-dependent peaks linked to a gene, 5,145 were linked by two or more strategies to 1,670 genes. Among GWAS loci for cardiometabolic traits, adipocyte-dependent peaks, but not preadipocyte-dependent peaks, showed significant enrichment (LD score regression P<0.005) for waist-to-hip ratio and modest enrichment (P < 0.05) for HDL-cholesterol. We identified 659 peaks linked to 503 genes by two or more approaches and overlapping a GWAS signal, suggesting a regulatory mechanism at these loci. To identify variants that may alter chromatin accessibility between timepoints, we identified 582 variants in 454 context-dependent peaks that demonstrated allelic imbalance in accessibility (FDR<5%), of which 55 peaks also overlapped GWAS variants. At one GWAS locus for palmitoleic acid, rs603424 was located in an adipocyte-dependent peak linked to SCD and exhibited allelic differences in transcriptional activity in adipocytes (P = 0.003) but not preadipocytes (P = 0.09). These results demonstrate that context-dependent peaks and genes can guide discovery of regulatory variants at GWAS loci and aid identification of regulatory mechanisms.

Cardiovascular and metabolic diseases are widespread, and an increased understanding of genetic mechanisms behind these diseases could improve treatment. Chromatin accessibility and gene expression in relevant cell contexts can guide identification of regulatory elements and genetic mechanisms for disease traits. A relevant context for cardiovascular and metabolic disease traits is adipocyte differentiation. To identify regulatory elements and genes that display differences in activity during adipocyte differentiation, we profiled chromatin accessibility and gene expression in a human cell model of preadipocytes and adipocytes. We identified chromatin regions that change accessibility during differentiation and predicted genes they may affect. We also linked these chromatin regions to genetic variants associated with risk of disease. At one genomic region linked to fatty acids, a chromatin region more accessible in adipocytes linked to a fatty acid synthesis gene and exhibited allelic differences in transcriptional activity in adipocytes but not preadipocytes. These results demonstrate that chromatin regions and genes that change during cell context can guide discovery of regulatory variants and aid identification of disease mechanisms.

Aspirin inhibits prostaglandins to prevents colon tumor formation via down-regulating Wnt production

According to numerous epidemiological studies, aspirin is a non-steroidal anti-inflammatory drug (NSAID) that reduces the occurrence and mortality of colorectal cancer (CRC). However, the underlying mechanisms are not well identified. In an effort to fill these gaps, we administered aspirin on mice one day before induction in an azoxymethane (AOM)/dextran sulfate sodium (DSS) induced CRC model. In this study, we assessed the effects of aspirin on tumorigenesis and tumor cell proliferation. Multi-layer analyses were carried out to identify changes in cytokines, metabolites, level of gene expressions, and proteins associated with tumorigenesis and aspirin treatment. The results showed that aspirin-treated mice developed fewer colon tumors in response to AOM/DSS, and aspirin can actively block cyclooxygenase (COX) metabolism and reduce levels of pro-inflammatory cytokines. In addition, the transcriptomic and proteomic analyses both indicated that aspirin has an inhibitory effect on the Wnt pathway. The in vitro results further indicated that aspirin inhibits WNT6 production, possibly by suppressing its transcription factor NR4A2, which in turn is regulated by prostaglandin E2, thereby ultimately inhibiting the Wnt pathway. These findings improve our understanding of the mechanisms behind aspirin's chemoprevention effect on CRC.

The molecular structure and function of sorting nexin 10 in skeletal disorders, cancers, and other pathological conditions

SNX10 is a member of the phox homology domain-containing family of phosphoinositide-binding proteins. Intracellularly, SNX10 localizes to endosomes where it mediates intracellular trafficking, endosome organization, and protein localization to the centrosome and cilium. It is highly expressed in bone and the gut where it participates in bone mineral and calcium homeostasis through the regulation of osteoclastic bone resorption and gastric acid secretion, respectively. Not surprisingly, patients harboring mutations in SNX10 mutation manifest a phenotype of autosomal recessive osteopetrosis or malignant infantile osteopetrosis, which is clinically characterized by dense bones with increased cortical bone into the medullary space with bone marrow occlusion or depletion, bone marrow failure, and anemia. Accordingly, SNX10 mutant osteoclasts exhibit impaired bone resorptive capacity. Beyond the skeleton, there is emerging evidence implicating SNX10 in cancer development, metabolic disorders, inflammation, and chaperone-mediated autophagy. Understanding the structural basis through which SNX10 exerts its diverse biological functions in both cell and tissue-specific manners may therefore inform new therapeutic opportunities toward the treatment and management of SNX10-related diseases.

Sex differences in human adipose tissue gene expression and genetic regulation involve adipogenesis

Sex differences in adipose tissue distribution and function are associated with sex differences in cardiometabolic disease. While many studies have revealed sex differences in adipocyte cell signaling and physiology, there is a relative dearth of information regarding sex differences in transcript abundance and regulation. We investigated sex differences in subcutaneous adipose tissue transcriptional regulation using omic-scale data from ∼3000 geographically and ethnically diverse human samples. We identified 162 genes with robust sex differences in expression. Differentially expressed genes were implicated in oxidative phosphorylation and adipogenesis. We further determined that sex differences in gene expression levels could be related to sex differences in the genetics of gene expression regulation. Our analyses revealed sex-specific genetic associations, and this finding was replicated in a study of 98 inbred mouse strains. The genes under genetic regulation in human and mouse were enriched for oxidative phosphorylation and adipogenesis. Enrichment analysis showed that the associated genetic loci resided within binding motifs for adipogenic transcription factors (e.g., PPARG and EGR1). We demonstrated that sex differences in gene expression could be influenced by sex differences in genetic regulation for six genes (e.g., FADS1 and MAP1B). These genes exhibited dynamic expression patterns during adipogenesis and robust expression in mature human adipocytes. Our results support a role for adipogenesis-related genes in subcutaneous adipose tissue sex differences in the genetic and environmental regulation of gene expression.

Adiponectin GWAS loci harboring extensive allelic heterogeneity exhibit distinct molecular consequences

Loci identified in genome-wide association studies (GWAS) can include multiple distinct association signals. We sought to identify the molecular basis of multiple association signals for adiponectin, a hormone involved in glucose regulation secreted almost exclusively from adipose tissue, identified in the Metabolic Syndrome in Men (METSIM) study. With GWAS data for 9,262 men, four loci were significantly associated with adiponectin: ADIPOQ, CDH13, IRS1, and PBRM1. We performed stepwise conditional analyses to identify distinct association signals, a subset of which are also nearly independent (lead variant pairwise r²<0.01). Two loci exhibited allelic heterogeneity, ADIPOQ and CDH13. Of seven association signals at the ADIPOQ locus, two signals colocalized with adipose tissue expression quantitative trait loci (eQTLs) for three transcripts: trait-increasing alleles at one signal were associated with increased ADIPOQ and LINC02043, while trait-increasing alleles at the other signal were associated with decreased ADIPOQ-AS1. In reporter assays, adiponectin-increasing alleles at two signals showed corresponding directions of effect on transcriptional activity. Putative mechanisms for the seven ADIPOQ signals include a missense variant (ADIPOQ G90S), a splice variant, a promoter variant, and four enhancer variants. Of two association signals at the CDH13 locus, the first signal consisted of promoter variants, including the lead adipose tissue eQTL variant for CDH13, while a second signal included a distal intron 1 enhancer variant that showed ~2-fold allelic differences in transcriptional reporter activity. Fine-mapping and experimental validation demonstrated that multiple, distinct association signals at these loci can influence multiple transcripts through multiple molecular mechanisms.

Many DNA variants affect common human traits, and distinct variants can have different effects on the function or expression level of the same gene. We identified variants associated with levels of adiponectin, a hormone involved in glucose regulation. Among these variants, we specifically studied the sets of variants located near two genes, ADIPOQ and CDH13, to determine how the variants affect gene expression or function. We focused on sets of variants that can be inherited together but are not always inherited together. Of the variants associated with adiponectin and located near ADIPOQ, one set were also associated with higher expression levels of the protein-coding ADIPOQ gene and a nearby non-coding gene, a second set of variants were associated with lower levels of the ADIPOQ-AS1 antisense gene, and additional variants changed the amino acid sequence or size of the adiponectin protein. Our examples show the benefits of identifying multiple sets of trait-associated variants in the same DNA region. These variants explain more trait variation, help identify genes that affect the trait, and guide studies of gene regulation and biological processes.

Adipose Tissue Gene Expression Associations Reveal Hundreds of Candidate Genes for Cardiometabolic Traits

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits.