Pancreatic islet enhancer clusters enriched in type 2 diabetes risk-associated variants

Long Non-coding RNAs as Local Regulators of Pancreatic Islet Transcription Factor Genes

The transcriptional programs of differentiated cells are tightly regulated by interactions between cell type-specific transcription factors and cis-regulatory elements. Long non-coding RNAs (lncRNAs) have emerged as additional regulators of gene transcription. Current evidence indicates that lncRNAs are a very heterogeneous group of molecules. For example, selected lncRNAs have been shown to regulate gene expression in cis or trans, although in most cases the precise underlying molecular mechanisms is unknown. Recent studies have uncovered a large number of lncRNAs that are selectively expressed in pancreatic islet cells, some of which were shown to regulate β cell transcriptional programs. A subset of such islet lncRNAs appears to control the expression of β cell-specific transcription factor (TF) genes by local cis-regulation. In this review, we discuss current knowledge of molecular mechanisms underlying cis-regulatory lncRNAs and discuss challenges involved in using genetic perturbations to define their function. We then discuss known examples of pancreatic islet lncRNAs that appear to exert cis-regulation of TF genes. We propose that cis-regulatory lncRNAs could represent a molecular target for modulation of diabetes-relevant genes.

Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps

We expanded GWAS discovery for type 2 diabetes (T2D) by combining data from 898,130 European-descent individuals (9% cases), after imputation to high-density reference panels. With these data, we (i) extend the inventory of T2D-risk variants (243 loci, 135 newly implicated in T2D predisposition, comprising 403 distinct association signals); (ii) enrich discovery of lower-frequency risk alleles (80 index variants with minor allele frequency <5%, 14 with estimated allelic odds ratio >2); (iii) substantially improve fine-mapping of causal variants (at 51 signals, one variant accounted for >80% posterior probability of association (PPA)); (iv) extend fine-mapping through integration of tissue-specific epigenomic information (islet regulatory annotations extend the number of variants with PPA >80% to 73); (v) highlight validated therapeutic targets (18 genes with associations attributable to coding variants); and (vi) demonstrate enhanced potential for clinical translation (genome-wide chip heritability explains 18% of T2D risk; individuals in the extremes of a T2D polygenic risk score differ more than ninefold in prevalence).

Type 2 Diabetes-Associated Genetic Variants Regulate Chromatin Accessibility in Human Islets

Type 2 diabetes (T2D) is a complex disorder in which both genetic and environmental risk factors contribute to islet dysfunction and failure. Genome-wide association studies (GWAS) have linked single nucleotide polymorphisms (SNPs), most of which are noncoding, in >200 loci to islet dysfunction and T2D. Identification of the putative causal variants and their target genes and whether they lead to gain or loss of function remains challenging. Here, we profiled chromatin accessibility in pancreatic islet samples from 19 genotyped individuals and identified 2,949 SNPs associated with in vivo cis-regulatory element use (i.e., chromatin accessibility quantitative trait loci [caQTL]). Among the caQTLs tested (n = 13) using luciferase reporter assays in MIN6 β-cells, more than half exhibited effects on enhancer activity that were consistent with in vivo chromatin accessibility changes. Importantly, islet caQTL analysis nominated putative causal SNPs in 13 T2D-associated GWAS loci, linking 7 and 6 T2D risk alleles, respectively, to gain or loss of in vivo chromatin accessibility. By investigating the effect of genetic variants on chromatin accessibility in islets, this study is an important step forward in translating T2D-associated GWAS SNP into functional molecular consequences.

Long noncoding RNA Bmncr regulates mesenchymal stem cell fate during skeletal aging

Bone marrow mesenchymal stem cells (BMSCs) exhibit an age-related lineage switch between osteogenic and adipogenic fates, which contributes to bone loss and adiposity. Here we identified a long noncoding RNA, Bmncr, which regulated the fate of BMSCs during aging. Mice depleted of Bmncr (Bmncr-KO) showed decreased bone mass and increased bone marrow adiposity, whereas transgenic overexpression of Bmncr (Bmncr-Tg) alleviated bone loss and bone marrow fat accumulation. Bmncr regulated the osteogenic niche of BMSCs by maintaining extracellular matrix protein fibromodulin (FMOD) and activation of the BMP2 pathway. Bmncr affected local 3D chromatin structure and transcription of Fmod. The absence of Fmod modified the bone phenotype of Bmncr-Tg mice. Further analysis revealed that Bmncr would serve as a scaffold to facilitate the interaction of TAZ and ABL, and thus facilitate the assembly of the TAZ and RUNX2/PPARG transcriptional complex, promoting osteogenesis and inhibiting adipogenesis. Adeno-associated viral-mediated overexpression of Taz in osteoprogenitors alleviated bone loss and marrow fat accumulation in Bmncr-KO mice. Furthermore, restoring BMNCR levels in human BMSCs reversed the age-related switch between osteoblast and adipocyte differentiation. Our findings indicate that Bmncr is a key regulator of the age-related osteogenic niche alteration and cell fate switch of BMSCs.

Super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers

Super-enhancers and stretch enhancers represent classes of transcriptional enhancers that have been shown to control the expression of cell identity genes and carry disease- and trait-associated variants. Specifically, super-enhancers are clusters of enhancers defined based on the binding occupancy of master transcription factors, chromatin regulators, or chromatin marks, while stretch enhancers are large chromatin-defined regulatory regions of at least 3,000 base pairs. Several studies have characterized these regulatory regions in numerous cell types and tissues to decipher their functional importance. However, the differences and similarities between these regulatory regions have not been fully assessed. We integrated genomic, epigenomic, and transcriptomic data from ten human cell types to perform a comparative analysis of super and stretch enhancers with respect to their chromatin profiles, cell type-specificity, and ability to control gene expression. We found that stretch enhancers are more abundant, more distal to transcription start sites, cover twice as much the genome, and are significantly less conserved than super-enhancers. In contrast, super-enhancers are significantly more enriched for active chromatin marks and cohesin complex, and more transcriptionally active than stretch enhancers. Importantly, a vast majority of super-enhancers (85%) overlap with only a small subset of stretch enhancers (13%), which are enriched for cell type-specific biological functions, and control cell identity genes. These results suggest that super-enhancers are transcriptionally more active and cell type-specific than stretch enhancers, and importantly, most of the stretch enhancers that are distinct from super-enhancers do not show an association with cell identity genes, are less active, and more likely to be poised enhancers.

Chronic d-serine supplementation impairs insulin secretion

OBJECTIVE

The metabolic role of d-serine, a non-proteinogenic NMDA receptor co-agonist, is poorly understood. Conversely, inhibition of pancreatic NMDA receptors as well as loss of the d-serine producing enzyme serine racemase have been shown to modulate insulin secretion. Thus, we aim to study the impact of chronic and acute d-serine supplementation on insulin secretion and other parameters of glucose homeostasis.

METHODS

We apply MALDI FT-ICR mass spectrometry imaging, NMR based metabolomics, 16s rRNA gene sequencing of gut microbiota in combination with a detailed physiological characterization to unravel the metabolic action of d-serine in mice acutely and chronically treated with 1% d-serine in drinking water in combination with either chow or high fat diet feeding. Moreover, we identify SNPs in SRR, the enzyme converting L-to d-serine and two subunits of the NMDA receptor to associate with insulin secretion in humans, based on the analysis of 2760 non-diabetic Caucasian individuals.

RESULTS

We show that chronic elevation of d-serine results in reduced high fat diet intake. In addition, d-serine leads to diet-independent hyperglycemia due to blunted insulin secretion from pancreatic beta cells. Inhibition of alpha 2-adrenergic receptors rapidly restores glycemia and glucose tolerance in d-serine supplemented mice. Moreover, we show that single nucleotide polymorphisms (SNPs) in SRR as well as in individual NMDAR subunits are associated with insulin secretion in humans.

CONCLUSION

Thus, we identify a novel role of d-serine in regulating systemic glucose metabolism through modulating insulin secretion.

Epigenetics and Epigenomics: Implications for Diabetes and Obesity

The American Diabetes Association convened a research symposium, "Epigenetics and Epigenomics: Implications for Diabetes and Obesity" on 17-19 November 2017. International experts in genetics, epigenetics, computational biology, and physiology discussed the current state of understanding of the relationships between genetics, epigenetics, and environment in diabetes and examined existing evidence for the role of epigenetic factors in regulating metabolism and the risk of diabetes and its complications. The authors summarize the presentations, which highlight how the complex interactions between genes and environment may in part be mediated through epigenetic changes and how information about nutritional and other environmental stimuli can be transmitted to the next generation. In addition, the authors present expert consensus on knowledge gaps and research recommendations for the field.

Positional integration of lung adenocarcinoma susceptibility loci with primary human alveolar epithelial cell epigenomes

AIM

To identify functional lung adenocarcinoma (LUAD) risk SNPs.

MATERIALS & METHODS

Eighteen validated LUAD risk SNPs (p ≤ 5 × 10^-8) and 930 SNPs in high linkage disequilibrium (r² > 0.5) were integrated with epigenomic information from primary human alveolar epithelial cells. Enhancer-associated SNPs likely affecting transcription factor-binding sites were predicted. Three SNPs were functionally investigated using luciferase assays, expression quantitative trait loci and cancer-specific expression.

RESULTS

Forty-seven SNPs mapped to putative enhancers; 11 located to open chromatin. Of these, seven altered predicted transcription factor-binding motifs. Rs6942067 showed allele-specific luciferase expression and expression quantitative trait loci analysis indicates that it influences expression of DCBLD1, a gene that encodes an unknown membrane protein and is overexpressed in LUAD.

CONCLUSION

Integration of candidate LUAD risk SNPS with epigenomic marks from normal alveolar epithelium identified numerous candidate functional LUAD risk SNPs including rs6942067, which appears to affect DCBLD1 expression. Data deposition: Data are provided in GEO record GSE84273.

The Dysregulation of the DLK1-MEG3 Locus in Islets From Patients With Type 2 Diabetes Is Mimicked by Targeted Epimutation of Its Promoter With TALE-DNMT Constructs

Type 2 diabetes mellitus (T2DM) is characterized by the inability of the insulin-producing β-cells to overcome insulin resistance. We previously identified an imprinted region on chromosome 14, the DLK1-MEG3 locus, as being downregulated in islets from humans with T2DM. In this study, using targeted epigenetic modifiers, we prove that increased methylation at the promoter of Meg3 in mouse βTC6 β-cells results in decreased transcription of the maternal transcripts associated with this locus. As a result, the sensitivity of β-cells to cytokine-mediated oxidative stress was increased. Additionally, we demonstrate that an evolutionarily conserved intronic region at the MEG3 locus can function as an enhancer in βTC6 β-cells. Using circular chromosome conformation capture followed by high-throughput sequencing, we demonstrate that the promoter of MEG3 physically interacts with this novel enhancer and other putative regulatory elements in this imprinted region in human islets. Remarkably, this enhancer is bound in an allele-specific manner by the transcription factors FOXA2, PDX1, and NKX2.2. Overall, these data suggest that the intronic MEG3 enhancer plays an important role in the regulation of allele-specific expression at the imprinted DLK1-MEG3 locus in human β-cells, which in turn impacts the sensitivity of β-cells to cytokine-mediated oxidative stress.

A Chromatin Basis for Cell Lineage and Disease Risk in the Human Pancreas

Understanding the genomic logic that underlies cellular diversity and developmental potential in the human pancreas will accelerate the growth of cell replacement therapies and reveal genetic risk mechanisms in diabetes. Here, we identified and characterized thousands of chromatin regions governing cell-specific gene regulation in human pancreatic endocrine and exocrine lineages, including islet β cells, α cells, duct, and acinar cells. Our findings have captured cellular ontogenies at the chromatin level, identified lineage-specific regulators potentially acting on these sites, and uncovered hallmarks of regulatory plasticity between cell types that suggest mechanisms to regenerate β cells from pancreatic endocrine or exocrine cells. Our work shows that disease risk variants related to pancreas are significantly enriched in these regulatory regions and reveals previously unrecognized links between endocrine and exocrine pancreas in diabetes risk.

Islet Long Noncoding RNAs: A Playbook for Discovery and Characterization

Diabetes is a complex group of metabolic disorders that can be accompanied by several comorbidities, including increased risk of early death. Decades of diabetes research have elucidated many genetic drivers of normal islet function and dysfunction; however, a lack of suitable treatment options suggests our knowledge about the disease remains incomplete. The establishment of long noncoding RNAs (lncRNAs), once dismissed as "junk" DNA, as essential gene regulators in many biological processes has redefined the central role for RNA in cells. Studies showing that misregulation of lncRNAs can lead to disease have contributed to the emergence of lncRNAs as attractive candidates for drug targeting. These findings underscore the need to reexamine islet biology in the context of a regulatory role for RNA. This review will 1) highlight what is known about lncRNAs in the context of diabetes, 2) summarize the strategies used in lncRNA discovery pipelines, and 3) discuss future directions and the potential impact of studying the role of lncRNAs in diabetes.

Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes

Type 2 diabetes (T2D) is a very common disease in humans. Here we conduct a meta-analysis of genome-wide association studies (GWAS) with ~16 million genetic variants in 62,892 T2D cases and 596,424 controls of European ancestry. We identify 139 common and 4 rare variants associated with T2D, 42 of which (39 common and 3 rare variants) are independent of the known variants. Integration of the gene expression data from blood (n = 14,115 and 2765) with the GWAS results identifies 33 putative functional genes for T2D, 3 of which were targeted by approved drugs. A further integration of DNA methylation (n = 1980) and epigenomic annotation data highlight 3 genes (CAMK1D, TP53INP1, and ATP5G1) with plausible regulatory mechanisms, whereby a genetic variant exerts an effect on T2D through epigenetic regulation of gene expression. Our study uncovers additional loci, proposes putative genetic regulatory mechanisms for T2D, and provides evidence of purifying selection for T2D-associated variants.

Genomic insights into the causes of type 2 diabetes

Genome-wide association studies have implicated around 250 genomic regions in predisposition to type 2 diabetes, with evidence for causal variants and genes emerging for several of these regions. Understanding of the underlying mechanisms, including the interplay between β-cell failure, insulin sensitivity, appetite regulation, and adipose storage has been facilitated by the integration of multidimensional data for diabetes-related intermediate phenotypes, detailed genomic annotations, functional experiments, and now multiomic molecular features. Studies in diverse ethnic groups and examples from population isolates have shown the value and need for a broad genomic approach to this global disease. Transethnic discovery efforts and large-scale biobanks in diverse populations and ancestries could help to address some of the Eurocentric bias. Despite rapid progress in the discovery of the highly polygenic architecture of type 2 diabetes, dominated by common alleles with small, cumulative effects on disease risk, these insights have been of little clinical use in terms of disease prediction or prevention, and have made only small contributions to subtype classification or stratified approaches to treatment. Successful development of academia-industry partnerships for exome or genome sequencing in large biobanks could help to deliver economies of scale, with implications for the future of genomics-focused research.

Editing the Epigenome: Reshaping the Genomic Landscape

The eukaryotic epigenome has an instrumental role in determining and maintaining cell identity and function. Epigenetic components such as DNA methylation, histone tail modifications, chromatin accessibility, and DNA architecture are tightly correlated with central cellular processes, while their dysregulation manifests in aberrant gene expression and disease. The ability to specifically edit the epigenome holds the promise of enhancing understanding of how epigenetic modifications function and enabling manipulation of cell phenotype for research or therapeutic purposes. Genome engineering technologies use highly specific DNA-targeting tools to precisely deposit epigenetic changes in a locus-specific manner, creating diverse epigenome editing platforms. This review summarizes these technologies and insights from recent studies, describes the complex relationship between epigenetic components and gene regulation, and highlights caveats and promises of the emerging field of epigenome editing, including applications for translational purposes, such as epigenetic therapy and regenerative medicine.

CDKN2A/B T2D Genome-Wide Association Study Risk SNPs Impact Locus Gene Expression and Proliferation in Human Islets

Genome-wide association studies link the CDKN2A/B locus with type 2 diabetes (T2D) risk, but mechanisms increasing risk remain unknown. The CDKN2A/B locus encodes cell cycle inhibitors p14, p15, and p16; MTAP; and ANRIL, a long noncoding RNA. The goal of this study was to determine whether CDKN2A/B T2D risk SNPs impact locus gene expression, insulin secretion, or β-cell proliferation in human islets. Islets from donors without diabetes (n = 95) were tested for SNP genotype (rs10811661, rs2383208, rs564398, and rs10757283), gene expression (p14, p15, p16, MTAP, ANRIL, PCNA, KI67, and CCND2), insulin secretion (n = 61), and β-cell proliferation (n = 47). Intriguingly, locus genes were coregulated in islets in two physically overlapping cassettes: p14-p16-ANRIL, which increased with age, and MTAP-p15, which did not. Risk alleles at rs10811661 and rs2383208 were differentially associated with expression of ANRIL, but not p14, p15, p16, or MTAP, in age-dependent fashion, such that younger homozygous risk donors had higher ANRIL expression, equivalent to older donor levels. We identified several risk SNP combinations that may impact locus gene expression, suggesting possible mechanisms by which SNPs impact locus biology. Risk allele carriers at ANRIL coding SNP rs564398 had reduced β-cell proliferation index. In conclusion, CDKN2A/B locus SNPs may impact T2D risk by modulating islet gene expression and β-cell proliferation.

Stretching the limits: from homeostasis to stem cell plasticity in wound healing and cancer

Stem cells (SCs) govern tissue homeostasis and wound repair. They reside within niches, the special microenvironments within tissues that control SC lineage outputs. Upon injury or stress, new signals emanating from damaged tissue can divert nearby cells into adopting behaviours that are not part of their homeostatic repertoire. This behaviour, known as SC plasticity, typically resolves as wounds heal. However, in cancer, it can endure. Recent studies have yielded insights into the orchestrators of maintenance and lineage commitment for SCs belonging to three mammalian tissues: the haematopoietic system, the skin epithelium and the intestinal epithelium. We delineate the multifactorial determinants and general principles underlying the remarkable facets of SC plasticity, which lend promise for regenerative medicine and cancer therapeutics.

Insulin promoter in human pancreatic β cells contacts diabetes susceptibility loci and regulates genes affecting insulin metabolism

Both type 1 and type 2 diabetes involve a complex interplay between genetic, epigenetic, and environmental factors. Our laboratory has been interested in the physical interactions, in nuclei of human pancreatic β cells, between the insulin (INS) gene and other genes that are involved in insulin metabolism. We have identified, using Circularized Chromosome Conformation Capture (4C), many physical contacts in a human pancreatic β cell line between the INS promoter on chromosome 11 and sites on most other chromosomes. Many of these contacts are associated with type 1 or type 2 diabetes susceptibility loci. To determine whether physical contact is correlated with an ability of the INS locus to affect expression of these genes, we knock down INS expression by targeting the promoter; 259 genes are either up or down-regulated. Of these, 46 make physical contact with INS We analyze a subset of the contacted genes and show that all are associated with acetylation of histone H3 lysine 27, a marker of actively expressed genes. To demonstrate the usefulness of this approach in revealing regulatory pathways, we identify from among the contacted sites the previously uncharacterized gene SSTR5-AS1 and show that it plays an important role in controlling the effect of somatostatin-28 on insulin secretion. These results are consistent with models in which clustering of genes supports transcriptional activity. This may be a particularly important mechanism in pancreatic β cells and in other cells where a small subset of genes is expressed at high levels.

A Common Type 2 Diabetes Risk Variant Potentiates Activity of an Evolutionarily Conserved Islet Stretch Enhancer and Increases C2CD4A and C2CD4B Expression

Genome-wide association studies (GWASs) and functional genomics approaches implicate enhancer disruption in islet dysfunction and type 2 diabetes (T2D) risk. We applied genetic fine-mapping and functional (epi)genomic approaches to a T2D- and proinsulin-associated 15q22.2 locus to identify a most likely causal variant, determine its direction of effect, and elucidate plausible target genes. Fine-mapping and conditional analyses of proinsulin levels of 8,635 non-diabetic individuals from the METSIM study support a single association signal represented by a cluster of 16 strongly associated (p < 10^-17) variants in high linkage disequilibrium (r² > 0.8) with the GWAS index SNP rs7172432. These variants reside in an evolutionarily and functionally conserved islet and β cell stretch or super enhancer; the most strongly associated variant (rs7163757, p = 3 × 10^-19) overlaps a conserved islet open chromatin site. DNA sequence containing the rs7163757 risk allele displayed 2-fold higher enhancer activity than the non-risk allele in reporter assays (p < 0.01) and was differentially bound by β cell nuclear extract proteins. Transcription factor NFAT specifically potentiated risk-allele enhancer activity and altered patterns of nuclear protein binding to the risk allele in vitro, suggesting that it could be a factor mediating risk-allele effects. Finally, the rs7163757 proinsulin-raising and T2D risk allele (C) was associated with increased expression of C2CD4B, and possibly C2CD4A, both of which were induced by inflammatory cytokines, in human islets. Together, these data suggest that rs7163757 contributes to genetic risk of islet dysfunction and T2D by increasing NFAT-mediated islet enhancer activity and modulating C2CD4B, and possibly C2CD4A, expression in (patho)physiologic states.

Genomic Effects of the Vitamin D Receptor: Potentially the Link between Vitamin D, Immune Cells, and Multiple Sclerosis

Vitamin D has a plethora of functions that are important for the maintenance of general health and in particular, the functional integrity of the immune system, such as promoting an anti-inflammatory cytokine profile and reducing the Treg/Th17 ratio. Multiple sclerosis (MS) is a chronic, inflammatory, and neurodegenerative central nervous system (CNS) disorder of probable autoimmune origin. MS is characterized by recurring or progressive demyelination and degeneration of the CNS due in part to a misguided immune response to as yet undefined (CNS) antigens, potentially including myelin basic protein and proteolipid protein. MS has also been shown to be associated significantly with environmental factors such as the lack of vitamin D. The role of vitamin D in the pathogenesis and progression of MS is complex. Recent genetic studies have shown that various common MS-associated risk-single-nucleotide polymorphisms (SNPs) are located within or in the vicinity of genes associated with the complex metabolism of vitamin D. The functional aspects of these genetic associations may be explained either by a direct SNP-associated loss- or gain-of-function in a vitamin D-associated gene or due to a change in the regulation of gene expression in certain immune cell types. The development of new genetic tools using next-generation sequencing: e.g., chromatin immunoprecipitation sequencing (ChIP-seq) and the accompanying rapid progress of epigenomics has made it possible to recognize that the association between vitamin D and MS could be based on the extensive and characteristic genomic binding of the vitamin D receptor (VDR). Therefore, it is important to analyze comprehensively the spatiotemporal VDR binding patterns that have been identified using ChIP-seq in multiple immune cell types to reveal an integral profile of genomic VDR interaction. In summary, the aim of this review is to connect genomic effects vitamin D has on immune cells with MS and thus, to contribute to a better understanding of the influence of vitamin D on the etiology and the pathogenesis of this complex autoimmune disease.

BRICHOS domain of Bri2 inhibits islet amyloid polypeptide (IAPP) fibril formation and toxicity in human beta cells

Accumulation of islet amyloid polypeptide (IAPP)-containing amyloid fibrils is the main pathological finding in pancreatic islets in type 2 diabetes. The formation of these IAPP amyloid fibrils is considered toxic and may constitute a major cause for the loss of insulin-producing beta cells. The protein domain BRICHOS is present in several different proproteins and possesses antiamyloid chaperone activity. This study demonstrates expression of the BRICHOS-containing protein Bri2 in human pancreatic beta cells and its colocalization with IAPP. The Bri2 BRICHOS domain effectively prevents IAPP from forming fibrils and protects cells from the toxicity associated with IAPP fibrillation. It is concluded that the Bri2 BRICHOS domain may act as an endogenous inhibitor of IAPP amyloid formation in pancreatic beta cells.

Aggregation of islet amyloid polypeptide (IAPP) into amyloid fibrils in islets of Langerhans is associated with type 2 diabetes, and formation of toxic IAPP species is believed to contribute to the loss of insulin-producing beta cells. The BRICHOS domain of integral membrane protein 2B (Bri2), a transmembrane protein expressed in several peripheral tissues and in the brain, has recently been shown to prevent fibril formation and toxicity of Aβ42, an amyloid-forming peptide in Alzheimer disease. In this study, we demonstrate expression of Bri2 in human islets and in the human beta-cell line EndoC-βH1. Bri2 colocalizes with IAPP intracellularly and is present in amyloid deposits in patients with type 2 diabetes. The BRICHOS domain of Bri2 effectively inhibits fibril formation in vitro and instead redirects IAPP into formation of amorphous aggregates. Reduction of endogenous Bri2 in EndoC-βH1 cells with siRNA increases sensitivity to metabolic stress leading to cell death while a concomitant overexpression of Bri2 BRICHOS is protective. Also, coexpression of IAPP and Bri2 BRICHOS in lateral ventral neurons of Drosophila melanogaster results in an increased cell survival. IAPP is considered to be the most amyloidogenic peptide known, and described findings identify Bri2, or in particular its BRICHOS domain, as an important potential endogenous inhibitor of IAPP aggregation and toxicity, with the potential to be a possible target for the treatment of type 2 diabetes.

Characterising cis-regulatory variation in the transcriptome of histologically normal and tumour-derived pancreatic tissues

OBJECTIVE

To elucidate the genetic architecture of gene expression in pancreatic tissues.

DESIGN

We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison.

RESULTS

We identified 38 615 cis-eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis-eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis-eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10-8) and tumour-derived (p=8.3×10-5) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the 'O' mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO 'O' mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL.

CONCLUSIONS

We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

A type 2 diabetes-associated SNP in KCNQ1 (rs163184) modulates the binding activity of the locus for Sp3 and Lsd1/Kdm1a, potentially affecting CDKN1C expression

Although genome-wide association studies have shown that potassium voltage-gated channel subfamily Q member 1 (KCNQ1) is one of the genes that is most significantly associated with type 2 diabetes mellitus (T2DM), functionally annotating disease-associated single nucleotide polymorphisms (SNPs) remains a challenge. Recently, our group described a novel strategy to identify proteins that bind to SNP-containing loci in an allele-specific manner. The present study successfully applied this strategy to investigate rs163184, a T2DM susceptibility SNP located in the intronic region of KCNQ1. Comparative analysis of DNA-binding proteins revealed that the binding activities for the genomic region containing SNP rs163184 differed between alleles for several proteins, including Sp3 and Lsd1/Kdm1a. Sp3 preferentially bound to the non-risk rs163184 allele and stimulated transcriptional activity in an artificial promoter containing this region. Lsd1/Kdm1a was identified to be preferentially recruited to the non-risk allele of the rs163184 region and reduced Sp3-dependent transcriptional activity in the artificial promoter. In addition, expression of the nearby cyclin‑dependent kinase inhibitor 1C (CDKN1C) gene was revealed to be upregulated after SP3 knockdown in cells that possessed non-risk alleles. This suggests that CDKN1C is potentially one of the functional targets of SNP rs163184, which modulates the binding activity of the locus for Sp3 and Lsd1/Kdm1a.

Genome-wide analysis of PDX1 target genes in human pancreatic progenitors

Objective

Homozygous loss-of-function mutations in the gene coding for the homeobox transcription factor (TF) PDX1 leads to pancreatic agenesis, whereas heterozygous mutations can cause Maturity-Onset Diabetes of the Young 4 (MODY4). Although the function of Pdx1 is well studied in pre-clinical models during insulin-producing β-cell development and homeostasis, it remains elusive how this TF controls human pancreas development by regulating a downstream transcriptional program. Also, comparative studies of PDX1 binding patterns in pancreatic progenitors and adult β-cells have not been conducted so far. Furthermore, many studies reported the association between single nucleotide polymorphisms (SNPs) and T2DM, and it has been shown that islet enhancers are enriched in T2DM-associated SNPs. Whether regions, harboring T2DM-associated SNPs are PDX1 bound and active at the pancreatic progenitor stage has not been reported so far.

Methods

In this study, we have generated a novel induced pluripotent stem cell (iPSC) line that efficiently differentiates into human pancreatic progenitors (PPs). Furthermore, PDX1 and H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify PDX1 transcriptional targets and active enhancer and promoter regions. To address potential differences in the function of PDX1 during development and adulthood, we compared PDX1 binding profiles from PPs and adult islets. Moreover, combining ChIP-seq and GWAS meta-analysis data we identified T2DM-associated SNPs in PDX1 binding sites and active chromatin regions.

Results

ChIP-seq for PDX1 revealed a total of 8088 PDX1-bound regions that map to 5664 genes in iPSC-derived PPs. The PDX1 target regions include important pancreatic TFs, such as PDX1 itself, RFX6, HNF1B, and MEIS1, which were activated during the differentiation process as revealed by the active chromatin mark H3K27ac and mRNA expression profiling, suggesting that auto-regulatory feedback regulation maintains PDX1 expression and initiates a pancreatic TF program. Remarkably, we identified several PDX1 target genes that have not been reported in the literature in human so far, including RFX3, required for ciliogenesis and endocrine differentiation in mouse, and the ligand of the Notch receptor DLL1, which is important for endocrine induction and tip-trunk patterning. The comparison of PDX1 profiles from PPs and adult human islets identified sets of stage-specific target genes, associated with early pancreas development and adult β-cell function, respectively. Furthermore, we found an enrichment of T2DM-associated SNPs in active chromatin regions from iPSC-derived PPs. Two of these SNPs fall into PDX1 occupied sites that are located in the intronic regions of TCF7L2 and HNF1B. Both of these genes are key transcriptional regulators of endocrine induction and mutations in cis-regulatory regions predispose to diabetes.

Conclusions

Our data provide stage-specific target genes of PDX1 during in vitro differentiation of stem cells into pancreatic progenitors that could be useful to identify pathways and molecular targets that predispose for diabetes. In addition, we show that T2DM-associated SNPs are enriched in active chromatin regions at the pancreatic progenitor stage, suggesting that the susceptibility to T2DM might originate from imperfect execution of a β-cell developmental program.

•

PDX1 ChIP-seq analysis revealed 5664 target genes in human pancreatic progenitors, including unreported target genes.

•

Comparison of PDX1 profiles from PPs and adult human islets identified stage-specific PDX1 target gene sets.

•

T2DM-associated SNPs are enriched in active chromatin regions from iPSC-derived PPs.

•

Three SNPs fall into PDX1 occupied sites, located in intronic regions of the developmental regulatory TFs TCF7L2 and HNF1B.

Re-analysis of public genetic data reveals a rare X-chromosomal variant associated with type 2 diabetes

The reanalysis of existing GWAS data represents a powerful and cost-effective opportunity to gain insights into the genetics of complex diseases. By reanalyzing publicly available type 2 diabetes (T2D) genome-wide association studies (GWAS) data for 70,127 subjects, we identify seven novel associated regions, five driven by common variants (LYPLAL1, NEUROG3, CAMKK2, ABO, and GIP genes), one by a low-frequency (EHMT2), and one driven by a rare variant in chromosome Xq23, rs146662075, associated with a twofold increased risk for T2D in males. rs146662075 is located within an active enhancer associated with the expression of Angiotensin II Receptor type 2 gene (AGTR2), a modulator of insulin sensitivity, and exhibits allelic specific activity in muscle cells. Beyond providing insights into the genetics and pathophysiology of T2D, these results also underscore the value of reanalyzing publicly available data using novel genetic resources and analytical approaches.

Genome-wide association studies have uncovered several loci associated with diabetes risk. Here, the authors reanalyse public type 2 diabetes GWAS data to fine map 50 known loci and identify seven new ones, including one near ATGR2 on the X-chromosome that doubles the risk of diabetes in men.

Identification of a functional variant for colorectal cancer risk mapping to chromosome 5q31.1

Genome-wide association studies (GWASs) have established chromosome 5q31.1 as a risk locus for colorectal cancer (CRC). We previously identified a potentially regulatory single nucleotide polymorphism (SNP) rs17716310 within 5q31.1. Now, we extended our study with another independent Chinese population, functional assays and analyses of TCGA (The Cancer Genome Atlas) data. Significant associations between rs17716310 and CRC risk were found in Present Study including 1075 CRC cases and 1999 controls (additive model: OR = 1.149, 95% CI = 1.027–1.286, P = 0.016), and in Combined Study including 1766 cases and 2708 controls (additive model: OR = 1.145, 95% CI = 1.045–1.254, P = 0.004). Dual luciferase reporter gene assays indicated that the variant C allele obviously increased transcriptional activity. Using TCGA datasets, we indicated rs17716310 as a cis expression quantitative trait locus (eQTL) for the gene SMAD5, whose expression was significantly higher in CRC tissues. These findings suggested that the functional polymorphism rs17716310 A > C might be a genetic modifier for CRC, promoting the expression of SMAD5 that belonged to the transforming growth factor beta (TGF-β) signaling pathway.

Integrative network analysis highlights biological processes underlying GLP-1 stimulated insulin secretion: A DIRECT study

Glucagon-like peptide 1 (GLP-1) stimulated insulin secretion has a considerable heritable component as estimated from twin studies, yet few genetic variants influencing this phenotype have been identified. We performed the first genome-wide association study (GWAS) of GLP-1 stimulated insulin secretion in non-diabetic individuals from the Netherlands Twin register (n = 126). This GWAS was enhanced using a tissue-specific protein-protein interaction network approach. We identified a beta-cell protein-protein interaction module that was significantly enriched for low gene scores based on the GWAS P-values and found support at the network level in an independent cohort from Tübingen, Germany (n = 100). Additionally, a polygenic risk score based on SNPs prioritized from the network was associated (P < 0.05) with glucose-stimulated insulin secretion phenotypes in up to 5,318 individuals in MAGIC cohorts. The network contains both known and novel genes in the context of insulin secretion and is enriched for members of the focal adhesion, extracellular-matrix receptor interaction, actin cytoskeleton regulation, Rap1 and PI3K-Akt signaling pathways. Adipose tissue is, like the beta-cell, one of the target tissues of GLP-1 and we thus hypothesized that similar networks might be functional in both tissues. In order to verify peripheral effects of GLP-1 stimulation, we compared the transcriptome profiling of ob/ob mice treated with liraglutide, a clinically used GLP-1 receptor agonist, versus baseline controls. Some of the upstream regulators of differentially expressed genes in the white adipose tissue of ob/ob mice were also detected in the human beta-cell network of genes associated with GLP-1 stimulated insulin secretion. The findings provide biological insight into the mechanisms through which the effects of GLP-1 may be modulated and highlight a potential role of the beta-cell expressed genes RYR2, GDI2, KIAA0232, COL4A1 and COL4A2 in GLP-1 stimulated insulin secretion.

Integration of human pancreatic islet genomic data refines regulatory mechanisms at Type 2 Diabetes susceptibility loci

Human genetic studies have emphasised the dominant contribution of pancreatic islet dysfunction to development of Type 2 Diabetes (T2D). However, limited annotation of the islet epigenome has constrained efforts to define the molecular mechanisms mediating the, largely regulatory, signals revealed by Genome-Wide Association Studies (GWAS). We characterised patterns of chromatin accessibility (ATAC-seq, n = 17) and DNA methylation (whole-genome bisulphite sequencing, n = 10) in human islets, generating high-resolution chromatin state maps through integration with established ChIP-seq marks. We found enrichment of GWAS signals for T2D and fasting glucose was concentrated in subsets of islet enhancers characterised by open chromatin and hypomethylation, with the former annotation predominant. At several loci (including CDC123, ADCY5, KLHDC5) the combination of fine-mapping genetic data and chromatin state enrichment maps, supplemented by allelic imbalance in chromatin accessibility pinpointed likely causal variants. The combination of increasingly-precise genetic and islet epigenomic information accelerates definition of causal mechanisms implicated in T2D pathogenesis.

Assay for Transposase Accessible Chromatin (ATAC-Seq) to Chart the Open Chromatin Landscape of Human Pancreatic Islets

The regulatory mechanisms that ensure an accurate control of gene transcription are central to cellular function, development and disease. Such mechanisms rely largely on noncoding regulatory sequences that allow the establishment and maintenance of cell identity and tissue-specific cellular functions.The study of chromatin structure and nucleosome positioning allowed revealing transcription factor accessible genomic sites with regulatory potential, facilitating the comprehension of tissue-specific cis-regulatory networks. Recently a new technique coupled with high-throughput sequencing named Assay for Transposase Accessible Chromatin (ATAC-seq) emerged as an efficient method to chart open chromatin genome wide. The application of such technique to different cell types allowed unmasking tissue-specific regulatory elements and characterizing cis-regulatory networks. Herein we describe the implementation of the ATAC-seq method to human pancreatic islets, a tissue playing a central role in the control of glucose metabolism.

Circadian Transcription from Beta Cell Function to Diabetes Pathophysiology

The mammalian circadian clock plays a central role in the temporal coordination of physiology across the 24-h light-dark cycle. A major function of the clock is to maintain energy constancy in anticipation of alternating periods of fasting and feeding that correspond with sleep and wakefulness. While it has long been recognized that humans exhibit robust variation in glucose tolerance and insulin sensitivity across the sleep-wake cycle, experimental genetic analysis has now revealed that the clock transcription cycle plays an essential role in insulin secretion and metabolic function within pancreatic beta cells. This review addresses how studies of the beta cell clock may elucidate the etiology of subtypes of diabetes associated with circadian and sleep cycle disruption, in addition to more general forms of the disease.

Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium

AIMS/HYPOTHESIS

Elevated levels of fasting glucose and fasting insulin in non-diabetic individuals are markers of dysregulation of glucose metabolism and are strong risk factors for type 2 diabetes. Genome-wide association studies have discovered over 50 SNPs associated with these traits. Most of these loci were discovered in European populations and have not been tested in a well-powered multi-ethnic study. We hypothesised that a large, ancestrally diverse, fine-mapping genetic study of glycaemic traits would identify novel and population-specific associations that were previously undetectable by European-centric studies.

METHODS

A multiethnic study of up to 26,760 unrelated individuals without diabetes, of predominantly Hispanic/Latino and African ancestries, were genotyped using the Metabochip. Transethnic meta-analysis of racial/ethnic-specific linear regression analyses were performed for fasting glucose and fasting insulin. We attempted to replicate 39 fasting glucose and 17 fasting insulin loci. Genetic fine-mapping was performed through sequential conditional analyses in 15 regions that included both the initially reported SNP association(s) and denser coverage of SNP markers. In addition, Metabochip-wide analyses were performed to discover novel fasting glucose and fasting insulin loci. The most significant SNP associations were further examined using bioinformatic functional annotation.

RESULTS

Previously reported SNP associations were significantly replicated (p ≤ 0.05) in 31/39 fasting glucose loci and 14/17 fasting insulin loci. Eleven glycaemic trait loci were refined to a smaller list of potentially causal variants through transethnic meta-analysis. Stepwise conditional analysis identified two loci with independent secondary signals (G6PC2-rs477224 and GCK-rs2908290), which had not previously been reported. Population-specific conditional analyses identified an independent signal in G6PC2 tagged by the rare variant rs77719485 in African ancestry. Further Metabochip-wide analysis uncovered one novel fasting insulin locus at SLC17A2-rs75862513.

CONCLUSIONS/INTERPRETATION

These findings suggest that while glycaemic trait loci often have generalisable effects across the studied populations, transethnic genetic studies help to prioritise likely functional SNPs, identify novel associations that may be population-specific and in turn have the potential to influence screening efforts or therapeutic discoveries.

DATA AVAILABILITY

The summary statistics from each of the ancestry-specific and transethnic (combined ancestry) results can be found under the PAGE study on dbGaP here: https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000356.v1.p1.

Genes Associated with Pancreas Development and Function Maintain Open Chromatin in iPSCs Generated from Human Pancreatic Beta Cells

Current in vitro islet differentiation protocols suffer from heterogeneity and low efficiency. Induced pluripotent stem cells (iPSCs) derived from pancreatic beta cells (BiPSCs) preferentially differentiate toward endocrine pancreas-like cells versus those from fibroblasts (FiPSCs). We interrogated genome-wide open chromatin in BiPSCs and FiPSCs via ATAC-seq and identified ∼8.3k significant, differential open chromatin sites (DOCS) between the two iPSC subtypes (false discovery rate [FDR] < 0.05). DOCS where chromatin was more accessible in BiPSCs (Bi-DOCS) were significantly enriched for known regulators of endodermal development, including bivalent and weak enhancers, and FOXA2 binding sites (FDR < 0.05). Bi-DOCS were associated with genes related to pancreas development and beta-cell function, including transcription factors mutated in monogenic diabetes (PDX1, NKX2-2, HNF1A; FDR < 0.05). Moreover, Bi-DOCS correlated with enhanced gene expression in BiPSC-derived definitive endoderm and pancreatic progenitor cells. Bi-DOCS therefore highlight genes and pathways governing islet-lineage commitment, which can be exploited for differentiation protocol optimization, diabetes disease modeling, and therapeutic purposes.

An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans

To characterize type 2 diabetes (T2D)-associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci (P < 5 × 10^-8), including variants near the GLP2R, GIP, and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common single nucleotide variants. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion and in adipocytes, monocytes, and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology.

Defining a Novel Role for the Pdx1 Transcription Factor in Islet β-Cell Maturation and Proliferation During Weaning

The transcription factor encoded by the Pdx1 gene is a critical transcriptional regulator, as it has fundamental actions in the formation of all pancreatic cell types, islet β-cell development, and adult islet β-cell function. Transgenic- and cell line-based experiments have identified 5'-flanking conserved sequences that control pancreatic and β-cell type-specific transcription, which are found within areas I (bp -2694 to -2561), II (bp -2139 to -1958), III (bp -1879 to -1799), and IV (bp -6200 to -5670). Because of the presence in area IV of binding sites for transcription factors associated with pancreas development and islet cell function, we analyzed how an endogenous deletion mutant affected Pdx1 expression embryonically and postnatally. The most striking result was observed in male Pdx1ΔIV mutant mice after 3 weeks of birth (i.e., the onset of weaning), with only a small effect on pancreas organogenesis and no deficiencies in their female counterparts. Compromised Pdx1 mRNA and protein levels in weaned male mutant β-cells were tightly linked with hyperglycemia, decreased β-cell proliferation, reduced β-cell area, and altered expression of Pdx1-bound genes that are important in β-cell replication, endoplasmic reticulum function, and mitochondrial activity. We discuss the impact of these novel findings to Pdx1 gene regulation and islet β-cell maturation postnatally.

Disentangling the Role of Melatonin and its Receptor MTNR1B in Type 2 Diabetes: Still a Long Way to Go?

PURPOSE OF REVIEW

Type 2 diabetes (T2D) is a complex genetic metabolic disorder. T2D heritability has been estimated around 40-70%. In the last decade, exponential progress has been made in identifying T2D genetic determinants, through genome-wide association studies (GWAS). Among single-nucleotide polymorphisms mostly associated with T2D risk, rs10830963 is located in the MTNR1B gene, encoding one of the two receptors of melatonin, a neurohormone involved in circadian rhythms. Subsequent studies aiming to disentangle the role of MTNR1B in T2D pathophysiology led to controversies. In this review, we will tackle them and will try to give some directions to get a better view of MTNR1B contribution to T2D pathophysiology.

RECENT FINDINGS

Recent studies either based on genetic/genomic analyses, clinical/epidemiology data, functional analyses at rs10830963 locus, insulin secretion assays in response to melatonin (involving or not MTNR1B) or animal model analyses have led to strong controversies at each level of interpretation. We discuss possible caveats in these studies and present ways to go beyond these issues, towards a better understanding of T2D molecular mechanisms, keeping in mind that melatonin is a versatile hormone and regulates many functions via its primary role in the body clock.

GWAS on prolonged gestation (post-term birth): analysis of successive Finnish birth cohorts

Background

Gestation is a crucial timepoint in human development. Deviation from a term gestational age correlates with both acute and long-term adverse health effects for the child. Both being born preterm and post-term, that is, having short and long gestational ages, are heritable and influenced by the prenatal and perinatal environment. Despite the obvious heritable component, specific genetic influences underlying differences in gestational age are poorly understood.

Methods

We investigated the genetic architecture of gestational age in 9141 individuals, including 1167 born post-term, across two Northern Finland cohorts born in 1966 or 1986.

Results

Here we identify one globally significant intronic genetic variant within the ADAMTS13 gene that is associated with prolonged gestation (p=4.85×10−8). Additional variants that reached suggestive levels of significance were identified within introns at the ARGHAP42 and TKT genes, and in the upstream (5’) intergenic regions of the B3GALT5 and SSBP2 genes. The variants near the ADAMTS13, B3GALT5, SSBP2 and TKT loci are linked to alterations in gene expression levels (cis-eQTLs). Luciferase assays confirmed the allele specific enhancer activity for the BGALT5 and TKT loci.

Conclusions

Our findings provide the first evidence of a specific genetic influence associated with prolonged gestation. This study forms a foundation for a better understanding of the genetic and long-term health risks faced by induced and post-term individuals. The long-term risks for induced individuals who have a previously overlooked post-term potential may be a major issue for current health providers.

Genetic aetiology of glycaemic traits: approaches and insights

Glycaemic traits such as fasting and post-challenge glucose and insulin measures, as well as glycated haemoglobin (HbA1c), are used to diagnose and monitor diabetes. These traits are risk factors for cardiovascular disease even below the diabetic threshold, and their study can additionally yield insights into the pathophysiology of type 2 diabetes. To date, a diverse set of genetic approaches have led to the discovery of over 97 loci influencing glycaemic traits. In this review, we will focus on recent advances in the genetic aetiology of glycaemic traits, and the resulting biological insights. We will provide a brief overview of results ranging from common, to low- and rare-frequency variant-trait association studies, studies leveraging the diversity across populations, and studies harnessing the power of genetic and genomic approaches to gain insights into the biological underpinnings of these traits.

Diverse type 2 diabetes genetic risk factors functionally converge in a phenotype-focused gene network

Type 2 Diabetes (T2D) constitutes a global health burden. Efforts to uncover predisposing genetic variation have been considerable, yet detailed knowledge of the underlying pathogenesis remains poor. Here, we constructed a T2D phenotypic-linkage network (T2D-PLN), by integrating diverse gene functional information that highlight genes, which when disrupted in mice, elicit similar T2D-relevant phenotypes. Sensitising the network to T2D-relevant phenotypes enabled significant functional convergence to be detected between genes implicated in monogenic or syndromic diabetes and genes lying within genomic regions associated with T2D common risk. We extended these analyses to a recent multiethnic T2D case-control exome of 12,940 individuals that found no evidence of T2D risk association for rare frequency variants outside of previously known T2D risk loci. Examining associations involving protein-truncating variants (PTV), most at low population frequencies, the T2D-PLN was able to identify a convergent set of biological pathways that were perturbed within four of five independent T2D case/control ethnic sets of 2000 to 5000 exomes each. These same pathways were found to be over-represented among both known monogenic or syndromic diabetes genes and genes within T2D-associated common risk loci. Our study demonstrates convergent biology amongst variants representing different classes of T2D genetic risk. Although convergence was observed at the pathway level, few of the contributing genes were found in common between different cohorts or variant classes, most notably between the exome variant sets which suggests that future rare variant studies may be better focusing their power onto a single population of recent common ancestry.

Three-dimensional genome architecture and emerging technologies: looping in disease

Genome compaction is a universal feature of cells and has emerged as a global regulator of gene expression. Compaction is maintained by a multitude of architectural proteins, long non-coding RNAs (lncRNAs), and regulatory DNA. Each component comprises interlinked regulatory circuits that organize the genome in three-dimensional (3D) space to manage gene expression. In this review, we update the current state of 3D genome catalogues and focus on how recent technological advances in 3D genomics are leading to an enhanced understanding of disease mechanisms. We highlight the use of genome-wide chromatin conformation capture (Hi-C) coupled with oligonucleotide capture technology (capture Hi-C) to map interactions between gene promoters and distal regulatory elements such as enhancers that are enriched for disease variants from genome-wide association studies (GWASs). We discuss how aberrations in architectural units are associated with various pathological outcomes, and explore how recent advances in genome and epigenome editing show great promise for a systematic understanding of complex genetic disorders. Our growing understanding of 3D genome architecture—coupled with the ability to engineer changes in it—may create novel therapeutic opportunities.

Alpha TC1 and Beta-TC-6 genomic profiling uncovers both shared and distinct transcriptional regulatory features with their primary islet counterparts

Alpha TC1 (αTC1) and Beta-TC-6 (βTC6) mouse islet cell lines are cellular models of islet (dys)function and type 2 diabetes (T2D). However, genomic characteristics of these cells, and their similarities to primary islet alpha and beta cells, are undefined. Here, we report the epigenomic (ATAC-seq) and transcriptomic (RNA-seq) landscapes of αTC1 and βTC6 cells. Each cell type exhibits hallmarks of its primary islet cell counterpart including cell-specific expression of beta (e.g., Pdx1) and alpha (e.g., Arx) cell transcription factors (TFs), and enrichment of binding motifs for these TFs in αTC1/βTC6 cis-regulatory elements. αTC1/βTC6 transcriptomes overlap significantly with the transcriptomes of primary mouse/human alpha and beta cells. Our data further indicate that ATAC-seq detects cell-specific regulatory elements for cell types comprising ≥ 20% of a mixed cell population. We identified αTC1/βTC6 cis-regulatory elements orthologous to those containing type 2 diabetes (T2D)-associated SNPs in human islets for 33 loci, suggesting these cells’ utility to dissect T2D molecular genetics in these regions. Together, these maps provide important insights into the conserved regulatory architecture between αTC1/βTC6 and primary islet cells that can be leveraged in functional (epi)genomic approaches to dissect the genetic and molecular factors controlling islet cell identity and function.

Sasquatch: predicting the impact of regulatory SNPs on transcription factor binding from cell- and tissue-specific DNase footprints

In the era of genome-wide association studies (GWAS) and personalized medicine, predicting the impact of single nucleotide polymorphisms (SNPs) in regulatory elements is an important goal. Current approaches to determine the potential of regulatory SNPs depend on inadequate knowledge of cell-specific DNA binding motifs. Here, we present Sasquatch, a new computational approach that uses DNase footprint data to estimate and visualize the effects of noncoding variants on transcription factor binding. Sasquatch performs a comprehensive k-mer-based analysis of DNase footprints to determine any k-mer's potential for protein binding in a specific cell type and how this may be changed by sequence variants. Therefore, Sasquatch uses an unbiased approach, independent of known transcription factor binding sites and motifs. Sasquatch only requires a single DNase-seq data set per cell type, from any genotype, and produces consistent predictions from data generated by different experimental procedures and at different sequence depths. Here we demonstrate the effectiveness of Sasquatch using previously validated functional SNPs and benchmark its performance against existing approaches. Sasquatch is available as a versatile webtool incorporating publicly available data, including the human ENCODE collection. Thus, Sasquatch provides a powerful tool and repository for prioritizing likely regulatory SNPs in the noncoding genome.

Complex cis-regulatory landscape of the insulin receptor gene underlies the broad expression of a central signaling regulator

Insulin signaling plays key roles in development, growth and metabolism through dynamic control of glucose uptake, global protein translation and transcriptional regulation. Altered levels of insulin signaling are known to play key roles in development and disease, yet the molecular basis of such differential signaling remains obscure. Expression of the insulin receptor (InR) gene itself appears to play an important role, but the nature of the molecular wiring controlling InR transcription has not been elucidated. We characterized the regulatory elements driving Drosophila InR expression and found that the generally broad expression of this gene is belied by complex individual switch elements, the dynamic regulation of which reflects direct and indirect contributions of FOXO, EcR, Rbf and additional transcription factors through redundant elements dispersed throughout ∼40 kb of non-coding regions. The control of InR transcription in response to nutritional and tissue-specific inputs represents an integration of multiple cis-regulatory elements, the structure and function of which may have been sculpted by evolutionary selection to provide a highly tailored set of signaling responses on developmental and tissue-specific levels.

A Type 2 Diabetes-Associated Functional Regulatory Variant in a Pancreatic Islet Enhancer at the ADCY5 Locus

Molecular mechanisms remain unknown for most type 2 diabetes genome-wide association study identified loci. Variants associated with type 2 diabetes and fasting glucose levels reside in introns of ADCY5, a gene that encodes adenylate cyclase 5. Adenylate cyclase 5 catalyzes the production of cyclic AMP, which is a second messenger molecule involved in cell signaling and pancreatic β-cell insulin secretion. We demonstrated that type 2 diabetes risk alleles are associated with decreased ADCY5 expression in human islets and examined candidate variants for regulatory function. rs11708067 overlaps a predicted enhancer region in pancreatic islets. The type 2 diabetes risk rs11708067-A allele showed fewer H3K27ac ChIP-seq reads in human islets, lower transcriptional activity in reporter assays in rodent β-cells (rat 832/13 and mouse MIN6), and increased nuclear protein binding compared with the rs11708067-G allele. Homozygous deletion of the orthologous enhancer region in 832/13 cells resulted in a 64% reduction in expression level of Adcy5, but not adjacent gene Sec22a, and a 39% reduction in insulin secretion. Together, these data suggest that rs11708067-A risk allele contributes to type 2 diabetes by disrupting an islet enhancer, which results in reduced ADCY5 expression and impaired insulin secretion.

Mechanisms of Type 2 Diabetes Risk Loci

PURPOSE OF REVIEW

Deciphering the mechanisms of type 2 diabetes (T2DM) risk loci can greatly inform on disease pathology. This review discusses current knowledge of mechanisms through which genetic variants influence T2DM risk and considerations for future studies.

RECENT FINDINGS

Over 100 T2DM risk loci to date have been identified. Candidate causal variants at risk loci map predominantly to non-coding sequence. Physiological, epigenomic and gene expression data suggest that variants at many known T2DM risk loci affect pancreatic islet regulation, although variants at other loci also affect protein function and regulatory processes in adipose, pre-adipose, liver, skeletal muscle and brain. The effects of T2DM variants on regulatory activity in these tissues appear largely, but not exclusively, due to altered transcription factor binding. Putative target genes of T2DM variants have been defined at an increasing number of loci and some, such as FTO, may entail several genes and multiple tissues. Gene networks in islets and adipocytes have been implicated in T2DM risk, although the molecular pathways of risk genes remain largely undefined. Efforts to fully define the mechanisms of T2DM risk loci are just beginning. Continued identification of risk mechanisms will benefit from combining genetic fine-mapping with detailed phenotypic association data, high-throughput epigenomics data from diabetes-relevant tissue, functional screening of candidate genes and genome editing of cellular and animal models.

Genome-wide association study identifies novel type II diabetes risk loci in Jordan subpopulations

The prevalence of Type II Diabetes (T2D) has been increasing and has become a disease of significant public health burden in Jordan. None of the previous genome-wide association studies (GWAS) have specifically investigated the Middle East populations. The Circassian and Chechen communities in Jordan represent unique populations that are genetically distinct from the Arab population and other populations in the Caucasus. Prevalence of T2D is very high in both the Circassian and Chechen communities in Jordan despite low obesity prevalence. We conducted GWAS on T2D in these two populations and further performed meta-analysis of the results. We identified a novel T2D locus at chr20p12.2 at genome-wide significance (rs6134031, P = 1.12 × 10⁻⁸) and we replicated the results in the Wellcome Trust Case Control Consortium (WTCCC) dataset. Another locus at chr12q24.31 is associated with T2D at suggestive significance level (top SNP rs4758690, P = 4.20 × 10⁻⁵) and it is a robust eQTL for the gene, MLXIP (P = 1.10 × 10⁻¹⁴), and is significantly associated with methylation level in MLXIP, the functions of which involves cellular glucose response. Therefore, in this first GWAS of T2D in Jordan subpopulations, we identified novel and unique susceptibility loci which may help inform the genetic underpinnings of T2D in other populations.

Epigenetic dysregulation in pancreatic islets and pathogenesis of type 2 diabetes

Most of type 2 diabetes (T2D) is thought to be the result of interaction between genetic and environmental factors. However, the genetic components discovered to date can explain only a small proportion of the observed heritability. The "missing heritability" may be accounted for by rare variants, gene-environment interactions, and epigenetics.

Long non-coding RNAs as novel players in β cell function and type 1 diabetes

Background

Long non-coding RNAs (lncRNAs) are a sub-class within non-coding RNA repertoire that have emerged as crucial regulators of the gene expression in various pathophysiological conditions. lncRNAs display remarkable versatility and wield their functions through interactions with RNA, DNA, or proteins. Accumulating body of evidence based on multitude studies has highlighted the role of lncRNAs in many autoimmune and inflammatory diseases, including type 1 diabetes (T1D).

Main body of abstract

This review highlights emerging roles of lncRNAs in immune and islet β cell function as well as some of the challenges and opportunities in understanding the pathogenesis of T1D and its complications.

Conclusion

We accentuate that the lncRNAs within T1D-loci regions in consort with regulatory variants and enhancer clusters orchestrate the chromatin remodeling in β cells and thereby act as cis/trans-regulatory determinants of islet cell transcriptional programs.

10 Years of GWAS Discovery: Biology, Function, and Translation

Application of the experimental design of genome-wide association studies (GWASs) is now 10 years old (young), and here we review the remarkable range of discoveries it has facilitated in population and complex-trait genetics, the biology of diseases, and translation toward new therapeutics. We predict the likely discoveries in the next 10 years, when GWASs will be based on millions of samples with array data imputed to a large fully sequenced reference panel and on hundreds of thousands of samples with whole-genome sequencing data.

GARLIC: a bioinformatic toolkit for aetiologically connecting diseases and cell type-specific regulatory maps

Genome-wide association studies (GWAS) have emerged as a powerful tool to uncover the genetic basis of human common diseases, which often show a complex, polygenic and multi-factorial aetiology. These studies have revealed that 70–90% of all single nucleotide polymorphisms (SNPs) associated with common complex diseases do not occur within genes (i.e. they are non-coding), making the discovery of disease-causative genetic variants and the elucidation of the underlying pathological mechanisms far from straightforward. Based on emerging evidences suggesting that disease-associated SNPs are frequently found within cell type-specific regulatory sequences, here we present GARLIC (GWAS-based Prediction Toolkit for Connecting Diseases and Cell Types), a user-friendly, multi-purpose software with an associated database and online viewer that, using global maps of cis-regulatory elements, can aetiologically connect human diseases with relevant cell types. Additionally, GARLIC can be used to retrieve potential disease-causative genetic variants overlapping regulatory sequences of interest. Overall, GARLIC can satisfy several important needs within the field of medical genetics, thus potentially assisting in the ultimate goal of uncovering the elusive and complex genetic basis of common human disorders.