3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

BACKGROUND

Genome-wide association studies (GWAS) have revealed a multitude of candidate genetic variants affecting the risk of developing complex traits and diseases. However, the highlighted regions are typically in the non-coding genome, and uncovering the functional causative single nucleotide variants (SNVs) is challenging. Prioritization of variants is commonly based on genomic annotation with markers of active regulatory elements, but current approaches still poorly predict functional variants. To address this, we systematically analyze six markers of active regulatory elements for their ability to identify functional variants.

RESULTS

We benchmark against molecular quantitative trait loci (molQTL) from assays of regulatory element activity that identify allelic effects on DNA-binding factor occupancy, reporter assay expression, and chromatin accessibility. We identify the combination of DNase footprints and divergent enhancer RNA (eRNA) as markers for functional variants. This signature provides high precision, but with a trade-off of low recall, thus substantially reducing candidate variant sets to prioritize variants for functional validation. We present this as a framework called FINDER-Functional SNV IdeNtification using DNase footprints and eRNA.

CONCLUSIONS

We demonstrate the utility to prioritize variants using leukocyte count trait and analyze variants in linkage disequilibrium with a lead variant to predict a functional variant in asthma. Our findings have implications for prioritizing variants from GWAS, in development of predictive scoring algorithms, and for functionally informed fine mapping approaches.

Untangling the genetics of beta cell dysfunction and death in type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is a complex multi-system disease which arises from both environmental and genetic factors, resulting in the destruction of insulin-producing pancreatic beta cells. Over the past two decades, human genetic studies have provided new insight into the etiology of T1D, including an appreciation for the role of beta cells in their own demise.

SCOPE OF REVIEW

Here, we outline models supported by human genetic data for the role of beta cell dysfunction and death in T1D. We highlight the importance of strong evidence linking T1D genetic associations to bona fide candidate genes for mechanistic and therapeutic consideration. To guide rigorous interpretation of genetic associations, we describe molecular profiling approaches, genomic resources, and disease models that may be used to construct variant-to-gene links and to investigate candidate genes and their role in T1D.

MAJOR CONCLUSIONS

We profile advances in understanding the genetic causes of beta cell dysfunction and death at individual T1D risk loci. We discuss how genetic risk prediction models can be used to address disease heterogeneity. Further, we present areas where investment will be critical for the future use of genetics to address open questions in the development of new treatment and prevention strategies for T1D.

Shared and unique 3D genomic features of substance use disorders across multiple cell types

Recent genome-wide association studies (GWAS) have revealed shared genetic components among alcohol, opioid, tobacco and cannabis use disorders. However, the extent of the underlying shared causal variants and effector genes, along with their cellular context, remain unclear. We leveraged our existing 3D genomic datasets comprising high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq and RNA-seq across >50 diverse human cell types to focus on genomic regions that coincide with GWAS loci. Using stratified LD regression, we determined the proportion of genomewide SNP heritability attributable to the features assayed across our cell types by integrating recent GWAS summary statistics for the relevant traits: alcohol use disorder (AUD), tobacco use disorder (TUD), opioid use disorder (OUD) and cannabis use disorder (CanUD). Statistically significant enrichments (P<0.05) were observed in 14 specific cell types, with heritability reaching 9.2-fold for iPSC-derived cortical neurons and neural progenitors, confirming that they are crucial cell types for further functional exploration. Additionally, several pancreatic cell types, notably pancreatic beta cells, showed enrichment for TUD, with heritability enrichments up to 4.8-fold, suggesting genomic overlap with metabolic processes. Further investigation revealed significant positive genetic correlations between T2D with both TUD and CanUD (FDR<0.05) and a significant negative genetic correlation with AUD. Interestingly, after partitioning the heritability for each cell type's cis-regulatory elements, the correlation between T2D and TUD for pancreatic beta cells was greater (r=0.2) than the global genetic correlation value. Our study provides new genomic insights into substance use disorders and implicates cell types where functional follow-up studies could reveal causal variant-gene mechanisms underpinning these disorders.

Primary osteoarthritis chondrocyte map of chromatin conformation reveals novel candidate effector genes

OBJECTIVES

Osteoarthritis is a complex disease with a huge public health burden. Genome-wide association studies (GWAS) have identified hundreds of osteoarthritis-associated sequence variants, but the effector genes underpinning these signals remain largely elusive. Understanding chromosome organisation in three-dimensional (3D) space is essential for identifying long-range contacts between distant genomic features (e.g., between genes and regulatory elements), in a tissue-specific manner. Here, we generate the first whole genome chromosome conformation analysis (Hi-C) map of primary osteoarthritis chondrocytes and identify novel candidate effector genes for the disease.

METHODS

Primary chondrocytes collected from 8 patients with knee osteoarthritis underwent Hi-C analysis to link chromosomal structure to genomic sequence. The identified loops were then combined with osteoarthritis GWAS results and epigenomic data from primary knee osteoarthritis chondrocytes to identify variants involved in gene regulation via enhancer-promoter interactions.

RESULTS

We identified 345 genetic variants residing within chromatin loop anchors that are associated with 77 osteoarthritis GWAS signals. Ten of these variants reside directly in enhancer regions of 10 newly described active enhancer-promoter loops, identified with multiomics analysis of publicly available chromatin immunoprecipitation sequencing (ChIP-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data from primary knee chondrocyte cells, pointing to two new candidate effector genes SPRY4 and PAPPA (pregnancy-associated plasma protein A) as well as further support for the gene SLC44A2 known to be involved in osteoarthritis. For example, PAPPA is directly associated with the turnover of insulin-like growth factor 1 (IGF-1) proteins, and IGF-1 is an important factor in the repair of damaged chondrocytes.

CONCLUSIONS

We have constructed the first Hi-C map of primary human chondrocytes and have made it available as a resource for the scientific community. By integrating 3D genomics with large-scale genetic association and epigenetic data, we identify novel candidate effector genes for osteoarthritis, which enhance our understanding of disease and can serve as putative high-value novel drug targets.

Unbiasedly decoding the tumor microenvironment with single-cell multiomics analysis in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive malignancy with a poor prognosis and limited therapeutic options. Research on the tumor microenvironment (TME) of PDAC has propelled the development of immunotherapeutic and targeted therapeutic strategies with a promising future. The emergence of single-cell sequencing and mass spectrometry technologies, coupled with spatial omics, has collectively revealed the heterogeneity of the TME from a multiomics perspective, outlined the development trajectories of cell lineages, and revealed important functions of previously underrated myeloid cells and tumor stroma cells. Concurrently, these findings necessitated more refined annotations of biological functions at the cell cluster or single-cell level. Precise identification of all cell clusters is urgently needed to determine whether they have been investigated adequately and to identify target cell clusters with antitumor potential, design compatible treatment strategies, and determine treatment resistance. Here, we summarize recent research on the PDAC TME at the single-cell multiomics level, with an unbiased focus on the functions and potential classification bases of every cellular component within the TME, and look forward to the prospects of integrating single-cell multiomics data and retrospectively reusing bulk sequencing data, hoping to provide new insights into the PDAC TME.

Gluco-regulation & type 2 diabetes: entrenched misconceptions updated to new governing principles for gold standard management

Our understanding of type 2 diabetes (T2D) has evolved dramatically. Advances have upended entrenched dogmas pertaining to the onset and progression of T2D, beliefs that have prevailed from the early era of diabetes research-and continue to populate our medical textbooks and continuing medical education materials. This review article highlights key insights that lend new governing principles for gold standard management of T2D. From the historical context upon which old beliefs arose to new findings, this article outlines evidence and perspectives on beta cell function, the underlying defects in glucoregulation, the remediable nature of T2D, and, the rationale supporting the shift to complication-centric prescribing. Practical approaches translate this rectified understanding of T2D into strategies that fill gaps in current management practices of prediabetes through late type 2 diabetes.

LONP1 regulation of mitochondrial protein folding provides insight into beta cell failure in type 2 diabetes

Proteotoxicity is a contributor to the development of type 2 diabetes (T2D), but it is unknown whether protein misfolding in T2D is generalized or has special features. Here, we report a robust accumulation of misfolded proteins within the mitochondria of human pancreatic islets in T2D and elucidate its impact on β cell viability. Surprisingly, quantitative proteomics studies of protein aggregates reveal that human islets from donors with T2D have a signature more closely resembling mitochondrial rather than ER protein misfolding. The matrix protease LonP1 and its chaperone partner mtHSP70 were among the proteins enriched in protein aggregates. Deletion of LONP1 in mice yields mitochondrial protein misfolding and reduced respiratory function, ultimately leading to β cell apoptosis and hyperglycemia. Intriguingly, LONP1 gain of function ameliorates mitochondrial protein misfolding and restores human β cell survival following glucolipotoxicity via a protease-independent effect requiring LONP1-mtHSP70 chaperone activity. Thus, LONP1 promotes β cell survival and prevents hyperglycemia by facilitating mitochondrial protein folding. These observations may open novel insights into the nature of impaired proteostasis on β cell loss in the pathogenesis of T2D that could be considered as future therapeutic targets.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

Glucagon-producing α-cell transcriptional identity and reprogramming towards insulin production

β-Cell replacement by in situ reprogramming of non-β-cells is a promising diabetes therapy. Following the observation that near-total β-cell ablation in adult mice triggers the reprogramming of pancreatic α-, δ-, and γ-cells into insulin (INS)-producing cells, recent studies are delving deep into the mechanisms controlling adult α-cell identity. Systematic analyses of the α-cell transcriptome and epigenome have started to pinpoint features that could be crucial for maintaining α-cell identity. Using different transgenic and chemical approaches, significant advances have been made in reprogramming α-cells in vivo into INS-secreting cells in mice. The recent reprogramming of human α-cells in vitro is an important step forward that must now be complemented with a comprehensive molecular dissection of the mechanisms controlling α-cell identity.

Variant to gene mapping for carpal tunnel syndrome risk loci implicates skeletal muscle regulatory elements

BACKGROUND

Carpal tunnel syndrome (CTS) is a common disorder caused by compression of the median nerve in the wrist, resulting in pain and numbness throughout the hand and forearm. While multiple behavioural and physiological factors influence CTS risk, a growing body of evidence supports a strong genetic contribution. Recent genome-wide association study (GWAS) efforts have reported 53 independent signals associated with CTS. While GWAS can identify genetic loci conferring risk, it does not determine which cell types drive the genetic aetiology of the trait, which variants are "causal" at a given signal, and which effector genes correspond to these non-coding variants. These obstacles limit interpretation of potential disease mechanisms.

METHODS

We analysed CTS GWAS findings in the context of chromatin conformation between gene promoters and accessible chromatin regions across cellular models of bone, skeletal muscle, adipocytes and neurons. We identified proxy variants in high LD with the lead CTS sentinel SNPs residing in promoter connected open chromatin in the skeletal muscle and bone contexts.

FINDINGS

We detected significant enrichment for heritability in skeletal muscle myotubes, as well as a weaker correlation in human mesenchymal stem cell-derived osteoblasts. In myotubes, our approach implicated 117 genes contacting 60 proxy variants corresponding to 20 of the 53 GWAS signals. In the osteoblast context we implicated 30 genes contacting 24 proxy variants coinciding with 12 signals, of which 19 genes shared. We subsequently prioritized BZW2 as a candidate effector gene in CTS and implicated it as novel gene that perturbs myocyte differentiation in vitro.

INTERPRETATION

Taken together our results suggest that the CTS genetic component influences the size, integrity, and organization of multiple tissues surrounding the carpal tunnel, in particular muscle and bone, to predispose the nerve to being compressed in this disease setting.

FUNDING

This work was supported by NIH Grant UM1 DK126194 (SFAG and WY), R01AG072705 (SFAG & KDH) and the Center for Spatial and Functional Genomics at CHOP (SFAG & ADW). SFAG is supported by the Daniel B. Burke Endowed Chair for Diabetes Research. WY is supported by the Perelman School of Medicine of the University of Pennsylvania.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Nutrigenomics and redox regulation: Concepts relating to the Special Issue on nutrigenomics

During our whole lifespan, from conception to death, the epigenomes of all tissues and cell types of our body integrate signals from the environment. This includes signals derived from our diet and the uptake of macro- and micronutrients. In most cases, this leads only to transient changes, but some effects of this epigenome programming process are persistent and can even be transferred to the next generation. Both epigenetic programming and redox processes are affected by the individual choice of diet and other lifestyle decisions like physical activity. The nutrient-gene communication pathways have adapted during human evolution and are essential for maintaining health. However, when they are maladaptive, such as in long-term obesity, they significantly contribute to diseases like type 2 diabetes and cancer. The field of nutrigenomics investigates nutrition-related signal transduction pathways and their effect on gene expression involving interactions both with the genome and the epigenomes. Several of these diet-(epi)genome interactions and the involved signal transduction cascades are redox-regulated. Examples include the effects of the NAD+/NADH ratio, vitamin C levels and secondary metabolites of dietary molecules from plants on the acetylation and methylation state of the epigenome as well as on gene expression through redox-sensitive pathways via the transcription factors NFE2L2 and FOXO. In this review, we summarize and extend on these topics as well as those discussed in the articles of this Special Issue and take them into the context of redox biology.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic human embryonic stem cell-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional roles of many loci remain unexplored. Here, we engineered isogenic knockout human embryonic stem cell lines for 20 genes associated with T2D risk. We examined the impacts of each knockout on β cell differentiation, functions, and survival. We generated gene expression and chromatin accessibility profiles on β cells derived from each knockout line. Analyses of T2D-association signals overlapping HNF4A-dependent ATAC peaks identified a likely causal variant at the FAIM2 T2D-association signal. Additionally, the integrative association analyses identified four genes (CP, RNASE1, PCSK1N, and GSTA2) associated with insulin production, and two genes (TAGLN3 and DHRS2) associated with β cell sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental line and identified a single likely functional variant at each of 23 T2D-association signals.

Multi-Scale Imaging of the Dynamic Organization of Chromatin

Chromatin is now regarded as a heterogeneous and dynamic structure occupying a non-random position within the cell nucleus, where it plays a key role in regulating various functions of the genome. This current view of chromatin has emerged thanks to high spatiotemporal resolution imaging, among other new technologies developed in the last decade. In addition to challenging early assumptions of chromatin being regular and static, high spatiotemporal resolution imaging made it possible to visualize and characterize different chromatin structures such as clutches, domains and compartments. More specifically, super-resolution microscopy facilitates the study of different cellular processes at a nucleosome scale, providing a multi-scale view of chromatin behavior within the nucleus in different environments. In this review, we describe recent imaging techniques to study the dynamic organization of chromatin at high spatiotemporal resolution. We also discuss recent findings, elucidated by these techniques, on the chromatin landscape during different cellular processes, with an emphasis on the DNA damage response.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.

Inferring causal genes at type 2 diabetes GWAS loci through chromosome interactions in islet cells

Background: Resolving causal genes for type 2 diabetes at loci implicated by genome-wide association studies (GWAS) requires integrating functional genomic data from relevant cell types. Chromatin features in endocrine cells of the pancreatic islet are particularly informative and recent studies leveraging chromosome conformation capture (3C) with Hi-C based methods have elucidated regulatory mechanisms in human islets. However, these genome-wide approaches are less sensitive and afford lower resolution than methods that target specific loci. Methods: To gauge the extent to which targeted 3C further resolves chromatin-mediated regulatory mechanisms at GWAS loci, we generated interaction profiles at 23 loci using next-generation (NG) capture-C in a human beta cell model (EndoC-βH1) and contrasted these maps with Hi-C maps in EndoC-βH1 cells and human islets and a promoter capture Hi-C map in human islets. Results: We found improvements in assay sensitivity of up to 33-fold and resolved ~3.6X more chromatin interactions. At a subset of 18 loci with 25 co-localised GWAS and eQTL signals, NG Capture-C interactions implicated effector transcripts at five additional genetic signals relative to promoter capture Hi-C through physical contact with gene promoters. Conclusions: High resolution chromatin interaction profiles at selectively targeted loci can complement genome- and promoter-wide maps.

Interpreting non-coding disease-associated human variants using single-cell epigenomics

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

ChromaFold predicts the 3D contact map from single-cell chromatin accessibility

The identification of cell-type-specific 3D chromatin interactions between regulatory elements can help to decipher gene regulation and to interpret the function of disease-associated non-coding variants. However, current chromosome conformation capture (3C) technologies are unable to resolve interactions at this resolution when only small numbers of cells are available as input. We therefore present ChromaFold, a deep learning model that predicts 3D contact maps and regulatory interactions from single-cell ATAC sequencing (scATAC-seq) data alone. ChromaFold uses pseudobulk chromatin accessibility, co-accessibility profiles across metacells, and predicted CTCF motif tracks as input features and employs a lightweight architecture to enable training on standard GPUs. Once trained on paired scATAC-seq and Hi-C data in human cell lines and tissues, ChromaFold can accurately predict both the 3D contact map and peak-level interactions across diverse human and mouse test cell types. In benchmarking against a recent deep learning method that uses bulk ATAC-seq, DNA sequence, and CTCF ChIP-seq to make cell-type-specific predictions, ChromaFold yields superior prediction performance when including CTCF ChIP-seq data as an input and comparable performance without. Finally, fine-tuning ChromaFold on paired scATAC-seq and Hi-C in a complex tissue enables deconvolution of chromatin interactions across cell subpopulations. ChromaFold thus achieves state-of-the-art prediction of 3D contact maps and regulatory interactions using scATAC-seq alone as input data, enabling accurate inference of cell-type-specific interactions in settings where 3C-based assays are infeasible.

An mTOR/RNA pol I axis shapes chromatin architecture in response to fasting

Chromatin architecture is a fundamental mediator of genome function. Fasting is a major environmental cue across the animal kingdom. Yet, how it impacts on 3D genome organization is unknown. Here, we show that fasting induces a reversible and large-scale spatial reorganization of chromatin in C. elegans . This fasting-induced 3D genome reorganization requires inhibition of the nutrient-sensing mTOR pathway, a major regulator of ribosome biogenesis. Remarkably, loss of transcription by RNA Pol I, but not RNA Pol II nor Pol III, induces a similar 3D genome reorganization in fed animals, and prevents the restoration of the fed-state architecture upon restoring nutrients to fasted animals. Our work documents the first large-scale chromatin reorganization triggered by fasting and reveals that mTOR and RNA Pol I shape genome architecture in response to nutrients.

New discoveries in the field of metabolism by applying single-cell and spatial omics

Single-cell multi-Omics (SCM-Omics) and spatial multi-Omics (SM-Omics) technologies provide state-of-the-art methods for exploring the composition and function of cell types in tissues/organs. Since its emergence in 2009, single-cell RNA sequencing (scRNA-seq) has yielded many groundbreaking new discoveries. The combination of this method with the emergence and development of SM-Omics techniques has been a pioneering strategy in neuroscience, developmental biology, and cancer research, especially for assessing tumor heterogeneity and T-cell infiltration. In recent years, the application of these methods in the study of metabolic diseases has also increased. The emerging SCM-Omics and SM-Omics approaches allow the molecular and spatial analysis of cells to explore regulatory states and determine cell fate, and thus provide promising tools for unraveling heterogeneous metabolic processes and making them amenable to intervention. Here, we review the evolution of SCM-Omics and SM-Omics technologies, and describe the progress in the application of SCM-Omics and SM-Omics in metabolism-related diseases, including obesity, diabetes, nonalcoholic fatty liver disease (NAFLD) and cardiovascular disease (CVD). We also conclude that the application of SCM-Omics and SM-Omics approaches can help resolve the molecular mechanisms underlying the pathogenesis of metabolic diseases in the body and facilitate therapeutic measures for metabolism-related diseases. This review concludes with an overview of the current status of this emerging field and the outlook for its future.

Single-nucleus multi-omics of human stem cell-derived islets identifies deficiencies in lineage specification

Insulin-producing β cells created from human pluripotent stem cells have potential as a therapy for insulin-dependent diabetes, but human pluripotent stem cell-derived islets (SC-islets) still differ from their in vivo counterparts. To better understand the state of cell types within SC-islets and identify lineage specification deficiencies, we used single-nucleus multi-omic sequencing to analyse chromatin accessibility and transcriptional profiles of SC-islets and primary human islets. Here we provide an analysis that enabled the derivation of gene lists and activity for identifying each SC-islet cell type compared with primary islets. Within SC-islets, we found that the difference between β cells and awry enterochromaffin-like cells is a gradient of cell states rather than a stark difference in identity. Furthermore, transplantation of SC-islets in vivo improved cellular identities overtime, while long-term in vitro culture did not. Collectively, our results highlight the importance of chromatin and transcriptional landscapes during islet cell specification and maturation.

Prioritization of genes associated with type 2 diabetes mellitus for functional studies

Existing therapies for type 2 diabetes mellitus (T2DM) show limited efficacy or have adverse effects. Numerous genetic variants associated with T2DM have been identified, but progress in translating these findings into potential drug targets has been limited. Here, we describe the tools and platforms available to identify effector genes from T2DM-associated coding and non-coding variants and prioritize them for functional studies. We discuss QSER1 and SLC12A8 as examples of genes that have been identified as possible T2DM candidate genes using these tools and platforms. We suggest further approaches, including the use of sequencing data with increased sample size and ethnic diversity, single-cell omics data for analyses, glycaemic trait associations to predict gene function and, potentially, human induced pluripotent stem cell 'village' cultures, to strengthen current gene functionalization workflows. Effective prioritization of T2DM-associated genes for experimental validation could expedite our understanding of the genetic mechanisms responsible for T2DM to facilitate the use of precision medicine in its treatment.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic hESC-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined β-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-β cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.

Vitamin D: A master example of nutrigenomics

Nutrigenomics attempts to characterize and integrate the relation between dietary molecules and gene expression on a genome-wide level. One of the biologically active nutritional compounds is vitamin D₃, which activates via its metabolite 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) the nuclear receptor VDR (vitamin D receptor). Vitamin D₃ can be synthesized endogenously in our skin, but since we spend long times indoors and often live at higher latitudes where for many winter months UV-B radiation is too low, it became a true vitamin. The ligand-inducible transcription factor VDR is expressed in the majority of human tissues and cell types, where it modulates the epigenome at thousands of genomic sites. In a tissue-specific fashion this results in the up- and downregulation of primary vitamin D target genes, some of which are involved in attenuating oxidative stress. Vitamin D affects a wide range of physiological functions including the control of metabolism, bone formation and immunity. In this review, we will discuss how the epigenome- and transcriptome-wide effects of 1,25(OH)₂D₃ and its receptor VDR serve as a master example in nutrigenomics. In this context, we will outline the basis of a mechanistic understanding for personalized nutrition with vitamin D₃.

Nutrigenomics in the context of evolution

Nutrigenomics describes the interaction between nutrients and our genome. Since the origin of our species most of these nutrient-gene communication pathways have not changed. However, our genome experienced over the past 50,000 years a number of evolutionary pressures, which are based on the migration to new environments concerning geography and climate, the transition from hunter-gatherers to farmers including the zoonotic transfer of many pathogenic microbes and the rather recent change of societies to a preferentially sedentary lifestyle and the dominance of Western diet. Human populations responded to these challenges not only by specific anthropometric adaptations, such as skin color and body stature, but also through diversity in dietary intake and different resistance to complex diseases like the metabolic syndrome, cancer and immune disorders. The genetic basis of this adaptation process has been investigated by whole genome genotyping and sequencing including that of DNA extracted from ancient bones. In addition to genomic changes, also the programming of epigenomes in pre- and postnatal phases of life has an important contribution to the response to environmental changes. Thus, insight into the variation of our (epi)genome in the context of our individual's risk for developing complex diseases, helps to understand the evolutionary basis how and why we become ill. This review will discuss the relation of diet, modern environment and our (epi)genome including aspects of redox biology. This has numerous implications for the interpretation of the risks for disease and their prevention.

Nutrition and epigenetic programming

PURPOSE OF REVIEW

The aim of this study is to highlight the epigenomic programming properties of nutritional molecules and their metabolites in human tissues and cell types.

RECENT FINDINGS

Chromatin is the physical expression of the epigenome and has a memory function on the level of DNA methylation, histone modification and 3-dimensional (3D) organization. This epigenetic memory does not only affect transient gene expression but also represents long-lasting decisions on cellular fate. The memory is based on an epigenetic programming process, which is directed by extracellular and intracellular signals that are sensed by transcription factors and chromatin modifiers. Many dietary molecules and their intermediary metabolites serve as such signals, that is they contribute to epigenetic programming and memory. In this context, we will discuss about molecules of intermediary energy metabolism affecting chromatin modifier actions, nutrition-triggered epigenetic memory in pre- and postnatal phases of life; and epigenetic programming of immune cells by vitamin D. These mechanisms explain some of the susceptibility for complex diseases, such as the metabolic syndrome, cancer and immune disorders.

SUMMARY

The observation that nutritional molecules are able to modulate the epigenome initiated the new nutrigenomic subdiscipline nutritional epigenetics. The concept that epigenetic memory and programming is directed by our diet has numerous implications for the interpretation of disease risk including their prevention.