An in vivo cis-regulatory screen at the type 2 diabetes associated TCF7L2 locus identifies multiple tissue-specific enhancers

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Mechanistic basis and efficacy of targeting the β-catenin-TCF7L2-JMJD6-c-Myc axis to overcome resistance to BET inhibitors

The promising activity of BET protein inhibitors (BETi's) is compromised by adaptive or innate resistance in acute myeloid leukemia (AML). Here, modeling of BETi-persister/resistance (BETi-P/R) in human postmyeloproliferative neoplasm (post-MPN) secondary AML (sAML) cells demonstrated accessible and active chromatin in specific superenhancers/enhancers, which was associated with increased levels of nuclear β-catenin, TCF7L2, JMJD6, and c-Myc in BETi-P/R sAML cells. Following BETi treatment, c-Myc levels were rapidly restored in BETi-P/R sAML cells. CRISPR/Cas9-mediated knockout of TCF7L2 or JMJD6 reversed BETi-P/R, whereas ectopic overexpression conferred BETi-P/R in sAML cells, confirming the mechanistic role of the β-catenin-TCF7L2-JMJD6-c-Myc axis in BETi resistance. Patient-derived, post-MPN, CD34+ sAML blasts exhibiting relative resistance to BETi, as compared with sensitive sAML blasts, displayed higher messenger RNA and protein expression of TCF7L2, JMJD6, and c-Myc and following BETi washout exhibited rapid restoration of c-Myc and JMJD6. CRISPR/Cas9 knockout of TCF7L2 and JMJD6 depleted their levels, inducing loss of viability of the sAML blasts. Disruption of colocalization of nuclear β-catenin with TBL1 and TCF7L2 by the small-molecule inhibitor BC2059 combined with depletion of BRD4 by BET proteolysis-targeting chimera reduced c-Myc levels and exerted synergistic lethality in BETi-P/R sAML cells. This combination also reduced leukemia burden and improved survival of mice engrafted with BETi-P/R sAML cells or patient-derived AML blasts innately resistant to BETi. Therefore, multitargeted disruption of the β-catenin-TCF7L2-JMJD6-c-Myc axis overcomes adaptive and innate BETi resistance, exhibiting preclinical efficacy against human post-MPN sAML cells.

Admixture mapping and fine-mapping of type 2 diabetes susceptibility loci in African American women

African American women are disproportionately affected by type 2 diabetes. Genetic factors may explain part of the excess risk. More than 100 genetic variants have been associated with risk of type 2 diabetes, but most studies have been conducted in white populations. Two genome-wide association studies (GWAS) in African Americans have identified three novel genetic variants only. We conducted admixture mapping using 2,918 ancestral informative markers in 2,632 cases of type 2 diabetes and 2,596 controls nested in the ongoing Black Women’s Health Study cohort, with the goal of identifying genomic loci with local African ancestry associated with type 2 diabetes. In addition, we performed replication analysis of 71 previously identified index SNPs, and fine-mapped those genetic loci to identify better or new genetic variants associated with type 2 diabetes in African Americans. We found that individual African ancestry was associated with higher risk of type 2 diabetes. In addition, we identified two genomic regions, 3q26 and 12q23, with excess of African ancestry associated with higher risk of type 2 diabetes. Lastly, we replicated 8 out of 71 index SNPs from previous GWAS, including, for the first time in African Americans, the X-linked rs5945326 SNP near the DUSP9 gene. In addition, our fine-mapping efforts suggest independent signals at five loci. Our detailed analysis identified two genomic regions associated with risk of type 2 diabetes, and showed that many genetic risk variants are shared across ancestries.

A genetic risk score that includes common type 2 diabetes risk variants is associated with gestational diabetes

OBJECTIVE

Gestational diabetes (GDM) is characterized by maternal glucose intolerance that manifests during pregnancy. Because GDM resembles type 2 diabetes (T2DM), shared genetic predisposition is likely but has not been established. We tested the hypothesis that a genetic risk score (GRS) that included variants known to be associated with T2DM is associated with GDM.

STUDY DESIGN

We conducted a case-control study using the Vanderbilt Medical Center biobank (BioVU) and calculated a simple-count GRS using 34 variants previously associated with T2DM or fasting glucose in the general population, or with GDM or glucose intolerance in pregnancy. We assessed the association of the GRS with GDM adjusting for maternal age, parity, and body mass index (BMI) and calculated the area under the curve for the receiver-operating characteristic curve (c-statistic).

STUDY POPULATION

Among Caucasian women, we identified 458 cases of GDM and 1538 pregnant controls with normal glucose tolerance.

RESULTS

Cases of GDM had a higher number of risk alleles compared to controls (38.9±4.0 vs 37.4±4.0 risk alleles, P=1.6×10-11 ). The GRS was significantly associated with GDM; the adjusted odds ratio associated with each additional risk allele was 1.10 (95% CI: 1.07-1.13, P=6×10-11 ). Clinical variables predicted the risk of GDM (c-statistic 0.67, 95% CI: 0.64-0.70), and adding the GRS modestly improved prediction (0.70, 95% CI: 0.67-0.73).

CONCLUSIONS

Among Caucasian women, a GRS that included common T2DM genetic risk variants was associated with increased risk of GDM but showed limited utility in the identification of GDM cases.

Effects of mutations in Wnt/β-catenin, hedgehog, Notch and PI3K pathways on GSK-3 activity-Diverse effects on cell growth, metabolism and cancer

Glycogen synthase kinase-3 (GSK-3) is a serine/threonine kinase that participates in an array of critical cellular processes. GSK-3 was first characterized as an enzyme that phosphorylated and inactivated glycogen synthase. However, subsequent studies have revealed that this moon-lighting protein is involved in numerous signaling pathways that regulate not only metabolism but also have roles in: apoptosis, cell cycle progression, cell renewal, differentiation, embryogenesis, migration, regulation of gene transcription, stem cell biology and survival. In this review, we will discuss the roles that GSK-3 plays in various diseases as well as how this pivotal kinase interacts with multiple signaling pathways such as: PI3K/PTEN/Akt/mTOR, Ras/Raf/MEK/ERK, Wnt/beta-catenin, hedgehog, Notch and TP53. Mutations that occur in these and other pathways can alter the effects that natural GSK-3 activity has on regulating these signaling circuits that can lead to cancer as well as other diseases. The novel roles that microRNAs play in regulation of the effects of GSK-3 will also be evaluated. Targeting GSK-3 and these other pathways may improve therapy and overcome therapeutic resistance.

Living long and ageing well: is epigenomics the missing link between nature and nurture?

Human longevity is a complex trait and increasingly we understand that both genes and lifestyle interact in the longevity phenotype. Non-genetic factors, including diet, physical activity, health habits, and psychosocial factors contribute approximately 50% of the variability in human lifespan with another 25% explained by genetic differences. Family clusters of nonagenarian and centenarian siblings, who show both exceptional age-span and health-span, are likely to have inherited facilitatory gene groups, but also have nine decades of life experiences and behaviours which have interacted with their genetic profiles. Identification of their shared genes is just one small step in the link from genes to their physical and psychological profiles. Behavioural genomics is beginning to demonstrate links to biological mechanisms through regulation of gene expression, which directs the proteome and influences the personal phenotype. Epigenetics has been considered the missing link between nature and nurture. Although there is much that remains to be discovered, this article will discuss some of genetic and environmental factors which appear important in good quality longevity and link known epigenetic mechanisms to themes identified by nonagenarians themselves related to their longevity. Here we suggest that exceptional 90-year old siblings have adopted a range of behaviours and life-styles which have contributed to their ageing-well-phenotype and which link with important public health messages.

Dissecting the causal genetic mechanisms of coronary heart disease

Large-scale genome-wide association studies (GWAS) have identified 46 loci that are associated with coronary heart disease (CHD). Additionally, 104 independent candidate variants (false discovery rate of 5 %) have been identified (Schunkert H, Konig IR, Kathiresan S, Reilly MP, Assimes TL, Holm H et al. Nat Genet 43:333-8, 2011; Deloukas P, Kanoni S, Willenborg C, Farrall M, Assimes TL, Thompson JR et al. Nat Genet 45:25-33, 2012; C4D Genetics Consortium. Nat Genet 43:339-44, 2011). The majority of the causal genes in these loci function independently of conventional risk factors. It is postulated that a number of the CHD-associated genes regulate basic processes in the vascular cells involved in atherosclerosis, and that study of the signaling pathways that are modulated in this cell type by causal regulatory variation will provide critical new insights for targeting the initiation and progression of disease. In this review, we will discuss the types of experimental approaches and data that are critical to understanding the molecular processes that underlie the disease risk at 9p21.3, TCF21, SORT1, and other CHD-associated loci.

Evidence of non-pancreatic beta cell-dependent roles of Tcf7l2 in the regulation of glucose metabolism in mice

Non-coding variation within TCF7L2 remains the strongest genetic determinant of type 2 diabetes risk in humans. A considerable effort has been placed in understanding the functional roles of TCF7L2 in pancreatic beta cells, despite evidence of TCF7L2 expression in various peripheral tissues important in glucose homeostasis. Here, we use a humanized mouse model overexpressing Tcf7l2, resulting in glucose intolerance, to infer the contribution of Tcf7l2 overexpression in beta cells and in other tissues to the metabolic phenotypes displayed by these mice. Restoring Tcf7l2 expression specifically in beta cells to endogenous levels, in face of its overexpression elsewhere, results in impaired insulin secretion, reduced beta cell number and islet area, corroborating data obtained in humans showing similar phenotypes as a result of manipulations leading to Tcf7l2 loss of function. Interestingly, the persistent overexpression of Tcf7l2 in non-pancreatic tissues results in a significant worsening in glucose tolerance in vivo, indicating that Tcf7l2 overexpression in beta cells does not account for the glucose intolerance in the Tcf7l2 overexpression mouse model. Collectively, these data posit that Tcf7l2 plays key roles in glucose metabolism through actions beyond pancreatic beta cells, and further points to functionally opposing cell-type specific effects for Tcf7l2 on the maintenance of balanced glucose metabolism, thereby urging a careful examination of its role in non-pancreatic tissues as well as its composite metabolic effects across distinct tissues. Uncovering these roles may lead to new therapeutic targets for type 2 diabetes.

Selective disruption of Tcf7l2 in the pancreatic β cell impairs secretory function and lowers β cell mass

Type 2 diabetes (T2D) is characterized by β cell dysfunction and loss. Single nucleotide polymorphisms in the T-cell factor 7-like 2 (TCF7L2) gene, associated with T2D by genome-wide association studies, lead to impaired β cell function. While deletion of the homologous murine Tcf7l2 gene throughout the developing pancreas leads to impaired glucose tolerance, deletion in the β cell in adult mice reportedly has more modest effects. To inactivate Tcf7l2 highly selectively in β cells from the earliest expression of the Ins1 gene (∼E11.5) we have therefore used a Cre recombinase introduced at the Ins1 locus. Tcfl2^fl/fl::Ins1Cre mice display impaired oral and intraperitoneal glucose tolerance by 8 and 16 weeks, respectively, and defective responses to the GLP-1 analogue liraglutide at 8 weeks. Tcfl2^fl/fl::Ins1Cre islets displayed defective glucose- and GLP-1-stimulated insulin secretion and the expression of both the Ins2 (∼20%) and Glp1r (∼40%) genes were significantly reduced. Glucose- and GLP-1-induced intracellular free Ca²⁺ increases, and connectivity between individual β cells, were both lowered by Tcf7l2 deletion in islets from mice maintained on a high (60%) fat diet. Finally, analysis by optical projection tomography revealed ∼30% decrease in β cell mass in pancreata from Tcfl2^fl/fl::Ins1Cre mice. These data demonstrate that Tcf7l2 plays a cell autonomous role in the control of β cell function and mass, serving as an important regulator of gene expression and islet cell coordination. The possible relevance of these findings for the action of TCF7L2 polymorphisms associated with Type 2 diabetes in man is discussed.

Spatially conserved regulatory elements identified within human and mouse Cd247 gene using high-throughput sequencing data from the ENCODE project

The Cd247 gene encodes for a transmembrane protein important for the expression and assembly of TCR/CD3 complex on the surface of T lymphocytes. Down-regulation of CD247 has functional consequences in systemic autoimmunity and has been shown to be associated with Type 1 Diabetes in NOD mouse. In this study, we have utilized the wealth of high-throughput sequencing data produced during the Encyclopedia of DNA Elements (ENCODE) project to identify spatially conserved regulatory elements within the Cd247 gene from human and mouse. We show the presence of two transcription factor binding sites, supported by histone marks and ChIP-seq data, that specifically have features of an enhancer and a promoter, respectively. We also identified a putative long non-coding RNA from the characteristically long first intron of the Cd247 gene. The long non-coding RNA annotation is supported by manual annotations from the GENCODE project in human and our expression quantification analysis performed in NOD and B6 mice using qRT-PCR. Furthermore, 17 of the 23 SNPs already known to be implicated with T1D were observed within the long non-coding RNA region in mouse. The spatially conserved regulatory elements identified in this study have the potential to enrich our understanding of the role of Cd247 gene in autoimmune diabetes.

How does the Mediterranean diet promote cardiovascular health? Current progress toward molecular mechanisms: gene-diet interactions at the genomic, transcriptomic, and epigenomic levels provide novel insights into new mechanisms

Epidemiological evidence supports a health-promoting effect of the Mediterranean Diet (MedDiet), especially in the prevention of cardiovascular diseases. These cardiovascular benefits have been attributed to a number of components of the MedDiet such as monounsaturated fatty acids, antioxidant vitamins and phytochemicals. However, the underlying mechanisms remain unknown. Likewise, little is known about the genes that define inter-individual variation in response to the MedDiet, although the TCF7L2 gene is emerging as an illustrative candidate for determining relative risk of cardiovascular events in response to the MedDiet. Moreover, omics technologies are providing evidence supporting potential mechanisms, some of them implicating epigenetics (i.e. microRNAs, methylation), and certain data suggest that some traditional foods could contribute via microRNAs possibly acting as exogenous regulators of gene expression. Future research should aim at increasing and consolidating the nutrigenetic and nutrigenomic knowledge of the MedDiet in order to provide sound, personalized and optimized nutritional recommendations.

Exploring the effects of polymorphisms on cis-regulatory signal transduction response

cis-Regulatory sequences (CRSs) direct cell-specific and inducible gene expression in response to signal transduction networks, and it is becoming apparent that many cases of disease susceptibility and drug response stratification are due to polymorphisms that alter CRS responses in a context-dependent manner. In the current review, we describe successful methods for identifying CRSs and analyzing the effects of allelic variation on their responses to signal transduction. The technologies described build on the successes of ENCODE (ENCyclopedia Of DNA Elements) by exploring the effects of polymorphisms on CRS context dependency. This understanding is essential to uncover the genomic basis of disease susceptibility and will play a major role in delivering on the promise of personalized medicine.

In vitro scan for enhancers at the TCF7L2 locus

AIMS/HYPOTHESIS

Recent functional characterisations of genome-wide association study (GWAS) loci suggest that cis-regulatory variation may be a common paradigm for complex disease susceptibility. Several studies point to a similar mechanism at the transcription factor 7-like 2 (TCF7L2) GWAS locus for type 2 diabetes. To address this possibility, we carried out an in vitro scan of this diabetes-associated locus to fine-map cis-regulatory sequences within this genomic interval.

METHODS

A systematic cell-based enhancer strategy was employed to interrogate all sequences within the 92 kb type-2-diabetes-association interval for cis-regulatory activity in a panel of cell lines (HCT-116, Neuro-2a, C2C12, U2OS, MIN6 and HepG2). We further evaluated chromatin state at a subset of these regions in HCT-116 and U2OS cells and examined allelic-specific enhancer properties at the type-2-diabetes-associated single nucleotide polymorphism (SNP) rs7903146.

RESULTS

In total, we assigned cis-regulatory activity to approximately 30% (9/28) of constructs tested. Notably, a subset of enhancers was active across multiple cell lines and overlapped with key epigenetic markers suggestive of cis-regulatory sequences. We further replicated the allelic-specific properties for SNP rs7903146 in pancreatic beta cells and additionally demonstrate identical allelic-specific enhancer effects in other cell lines.

CONCLUSIONS

These results provide a detailed map of cis-regulatory elements within the TCF7L2 GWAS locus and support the hypothesis of cis-regulatory variation leading to type 2 diabetes predisposition. The detection of allelic-specific effects for SNP rs7903146 in multiple cell lines further alludes to the likelihood of a peripheral defect in disease aetiology.

From mice to humans

The genomes of many species have now been completely sequenced including human and mouse. Great progress has been made in understanding the complex genetics that underlie diabetes and obesity in human populations. One of the current challenges is the functional identification and characterization of the genes within loci that are being mapped. There are many approaches to this problem and this review outlines the valuable role that the mouse can play. We outline the mouse resources that are available to the research community, including knockouts with conditional potential for every gene, and the efforts of the International Mouse Phenotyping Consortium to attach phenotype information to these genes. We also briefly consider the potential of TALEN technology to tailor-make new mouse models of specific mutations discovered in humans. Finally, we consider the recent progress in characterizing the GWAS genes FTO, TCF7L2, CDKAL1, and SLC30A8 in engineered mouse models.