Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes

OBJECTIVE

To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D).

RESEARCH DESIGN AND METHODS

We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD.

RESULTS

Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs11444867 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16).

CONCLUSIONS

The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Silencing ANGPTL8 reduces mouse preadipocyte differentiation and insulin signaling

ANGPTL8, expressed mainly in the liver and adipose tissue, regulates the activity of lipoprotein lipase (LPL) present in the extracellular space and triglyceride (TG) metabolism through its interaction with ANGPTL3 and ANGPTL4. Whether intracellular ANGPTL8 can also exert effects in tissues where it is expressed is uncertain. ANGPTL8 expression was low in preadipocytes and much increased during differentiation. To better understand the role of intracellular ANGPTL8 in adipocytes and assess whether it may play a role in adipocyte differentiation, we knocked down its expression in normal mouse subcutaneous preadipocytes. ANGPTL8 knockdown reduced adipocyte differentiation, cellular TG accumulation and also isoproterenol-stimulated lipolysis at day 7 of differentiation. RNA-Seq analysis of ANGPTL8 siRNA or control siRNA transfected SC preadipocytes on days 0, 2, 4 and 7 of differentiation showed that ANGPTL8 knockdown impeded the early (day 2) expression of adipogenic and insulin signaling genes, PPARγ, as well as genes related to extracellular matrix and NF-κB signaling. Insulin mediated Akt phosphorylation was reduced at an early stage during adipocyte differentiation. This study based on normal primary cells shows that ANGPTL8 has intracellular actions in addition to effects in the extracellular space, like modulating LPL activity. Preadipocyte ANGPTL8 expression modulates their differentiation possibly via changes in insulin signaling gene expression.

Differential CpG methylation at Nnat in the early establishment of beta cell heterogeneity

AIMS/HYPOTHESIS

Beta cells within the pancreatic islet represent a heterogenous population wherein individual sub-groups of cells make distinct contributions to the overall control of insulin secretion. These include a subpopulation of highly connected 'hub' cells, important for the propagation of intercellular Ca2+ waves. Functional subpopulations have also been demonstrated in human beta cells, with an altered subtype distribution apparent in type 2 diabetes. At present, the molecular mechanisms through which beta cell hierarchy is established are poorly understood. Changes at the level of the epigenome provide one such possibility, which we explore here by focusing on the imprinted gene Nnat (encoding neuronatin [NNAT]), which is required for normal insulin synthesis and secretion.

METHODS

Single-cell RNA-seq datasets were examined using Seurat 4.0 and ClusterProfiler running under R. Transgenic mice expressing enhanced GFP under the control of the Nnat enhancer/promoter regions were generated for FACS of beta cells and downstream analysis of CpG methylation by bisulphite sequencing and RNA-seq, respectively. Animals deleted for the de novo methyltransferase DNA methyltransferase 3 alpha (DNMT3A) from the pancreatic progenitor stage were used to explore control of promoter methylation. Proteomics was performed using affinity purification mass spectrometry and Ca2+ dynamics explored by rapid confocal imaging of Cal-520 AM and Cal-590 AM. Insulin secretion was measured using homogeneous time-resolved fluorescence imaging.

RESULTS

Nnat mRNA was differentially expressed in a discrete beta cell population in a developmental stage- and DNA methylation (DNMT3A)-dependent manner. Thus, pseudo-time analysis of embryonic datasets demonstrated the early establishment of Nnat-positive and -negative subpopulations during embryogenesis. NNAT expression is also restricted to a subset of beta cells across the human islet that is maintained throughout adult life. NNAT+ beta cells also displayed a discrete transcriptome at adult stages, representing a subpopulation specialised for insulin production, and were diminished in db/db mice. 'Hub' cells were less abundant in the NNAT+ population, consistent with epigenetic control of this functional specialisation.

CONCLUSIONS/INTERPRETATION

These findings demonstrate that differential DNA methylation at Nnat represents a novel means through which beta cell heterogeneity is established during development. We therefore hypothesise that changes in methylation at this locus may contribute to a loss of beta cell hierarchy and connectivity, potentially contributing to defective insulin secretion in some forms of diabetes.

DATA AVAILABILITY

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD048465.

Differential CpG methylation at Nnat in the early establishment of beta cell heterogeneity

Aims/hypothesis

Beta cells within the pancreatic islet represent a heterogenous population wherein individual sub-groups of cells make distinct contributions to the overall control of insulin secretion. These include a subpopulation of highly-connected 'hub' cells, important for the propagation of intercellular Ca2+ waves. Functional subpopulations have also been demonstrated in human beta cells, with an altered subtype distribution apparent in type 2 diabetes. At present, the molecular mechanisms through which beta cell hierarchy is established are poorly understood. Changes at the level of the epigenome provide one such possibility which we explore here by focussing on the imprinted gene neuronatin (Nnat), which is required for normal insulin synthesis and secretion.

Methods

Single cell RNA-seq datasets were examined using Seurat 4.0 and ClusterProfiler running under R. Transgenic mice expressing eGFP under the control of the Nnat enhancer/promoter regions were generated for fluorescence-activated cell (FAC) sorting of beta cells and downstream analysis of CpG methylation by bisulphite and RNA sequencing, respectively. Animals deleted for the de novo methyltransferase, DNMT3A from the pancreatic progenitor stage were used to explore control of promoter methylation. Proteomics was performed using affinity purification mass spectrometry and Ca2+ dynamics explored by rapid confocal imaging of Cal-520 and Cal-590. Insulin secretion was measured using Homogeneous Time Resolved Fluorescence Imaging.

Results

Nnat mRNA was differentially expressed in a discrete beta cell population in a developmental stage- and DNA methylation (DNMT3A)-dependent manner. Thus, pseudo-time analysis of embryonic data sets demonstrated the early establishment of Nnat-positive and negative subpopulations during embryogenesis. NNAT expression is also restricted to a subset of beta cells across the human islet that is maintained throughout adult life. NNAT+ beta cells also displayed a discrete transcriptome at adult stages, representing a sub-population specialised for insulin production, reminiscent of recently-described "βHI" cells and were diminished in db/db mice. 'Hub' cells were less abundant in the NNAT+ population, consistent with epigenetic control of this functional specialization.

Conclusions/interpretation

These findings demonstrate that differential DNA methylation at Nnat represents a novel means through which beta cell heterogeneity is established during development. We therefore hypothesise that changes in methylation at this locus may thus contribute to a loss of beta cell hierarchy and connectivity, potentially contributing to defective insulin secretion in some forms of diabetes.

Comparative proximity biotinylation implicates the small GTPase RAB18 in sterol mobilization and biosynthesis

Loss of functional RAB18 causes the autosomal recessive condition Warburg Micro syndrome. To better understand this disease, we used proximity biotinylation to generate an inventory of potential RAB18 effectors. A restricted set of 28 RAB18 interactions were dependent on the binary RAB3GAP1-RAB3GAP2 RAB18-guanine nucleotide exchange factor complex. Twelve of these 28 interactions are supported by prior reports, and we have directly validated novel interactions with SEC22A, TMCO4, and INPP5B. Consistent with a role for RAB18 in regulating membrane contact sites, interactors included groups of microtubule/membrane-remodeling proteins, membrane-tethering and docking proteins, and lipid-modifying/transporting proteins. Two of the putative interactors, EBP and OSBPL2/ORP2, have sterol substrates. EBP is a Δ8-Δ7 sterol isomerase, and ORP2 is a lipid transport protein. This prompted us to investigate a role for RAB18 in cholesterol biosynthesis. We found that the cholesterol precursor and EBP-product lathosterol accumulates in both RAB18-null HeLa cells and RAB3GAP1-null fibroblasts derived from an affected individual. Furthermore, de novo cholesterol biosynthesis is impaired in cells in which RAB18 is absent or dysregulated or in which ORP2 expression is disrupted. Our data demonstrate that guanine nucleotide exchange factor-dependent Rab interactions are highly amenable to interrogation by proximity biotinylation and may suggest that Micro syndrome is a cholesterol biosynthesis disorder.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus

BACKGROUND

Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD.

METHODS

From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D.

RESULTS

A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase (P=1.0×10-16).

CONCLUSIONS

The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Molecular phenotyping of single pancreatic islet leader beta cells by "Flash-Seq"

AIMS

Spatially-organized increases in cytosolic Ca²⁺ within pancreatic beta cells in the pancreatic islet underlie the stimulation of insulin secretion by high glucose. Recent data have revealed the existence of subpopulations of beta cells including "leaders" which initiate Ca²⁺ waves. Whether leader cells possess unique molecular features, or localisation, is unknown.

MAIN METHODS

High speed confocal Ca²⁺ imaging was used to identify leader cells and connectivity analysis, running under MATLAB and Python, to identify highly connected "hub" cells. To explore transcriptomic differences between beta cell sub-groups, individual leaders or followers were labelled by photo-activation of the cryptic fluorescent protein PA-mCherry and subjected to single cell RNA sequencing ("Flash-Seq").

KEY FINDINGS

Distinct Ca²⁺ wave types were identified in individual islets, with leader cells present in 73 % (28 of 38 islets imaged). Scale-free, power law-adherent behaviour was also observed in 29 % of islets, though "hub" cells in these islets did not overlap with leaders. Transcripts differentially expressed (295; padj < 0.05) between leader and follower cells included genes involved in cilium biogenesis and transcriptional regulation. Providing some support for these findings, ADCY6 immunoreactivity tended to be higher in leader than follower cells, whereas cilia number and length tended to be lower in the former. Finally, leader cells were located significantly closer to delta, but not alpha, cells in Euclidian space than were follower cells.

SIGNIFICANCE

The existence of both a discrete transcriptome and unique localisation implies a role for these features in defining the specialized function of leaders. These data also raise the possibility that localised signalling between delta and leader cells contributes to the initiation and propagation of islet Ca²⁺ waves.

Single-Cell RNA Sequencing Reveals a Role for Reactive Oxygen Species and Peroxiredoxins in Fatty Acid-Induced Rat β-Cell Proliferation

The functional mass of insulin-secreting pancreatic β-cells expands to maintain glucose homeostasis in the face of nutrient excess, in part via replication of existing β-cells. Type 2 diabetes appears when these compensatory mechanisms fail. Nutrients including glucose and fatty acids are important contributors to the β-cell compensatory response, but their underlying mechanisms of action remain poorly understood. We investigated the transcriptional mechanisms of β-cell proliferation in response to fatty acids. Isolated rat islets were exposed to 16.7 mmol/L glucose with or without 0.5 mmol/L oleate (C18:1) or palmitate (C16:0) for 48 h. The islet transcriptome was assessed by single-cell RNA sequencing. β-Cell proliferation was measured by flow cytometry. Unsupervised clustering of pooled β-cells identified different subclusters, including proliferating β-cells. β-Cell proliferation increased in response to oleate but not palmitate. Both fatty acids enhanced the expression of genes involved in energy metabolism and mitochondrial activity. Comparison of proliferating versus nonproliferating β-cells and pseudotime ordering suggested the involvement of reactive oxygen species (ROS) and peroxiredoxin signaling. Accordingly, N-acetyl cysteine and the peroxiredoxin inhibitor conoidin A both blocked oleate-induced β-cell proliferation. Our study reveals a key role for ROS signaling through peroxiredoxin activation in oleate-induced β-cell proliferation.

Hepatic glycerol shunt and glycerol-3-phosphate phosphatase control liver metabolism and glucodetoxification under hyperglycemia

OBJECTIVE

Glycerol-3-phosphate (Gro3P) phosphatase (G3PP) hydrolyzes Gro3P to glycerol that exits the cell, thereby operating a "glycerol shunt", a metabolic pathway that we identified recently in mammalian cells. We have investigated the role of G3PP and the glycerol shunt in the regulation of glucose metabolism and lipogenesis in mouse liver.

METHODS

We generated hepatocyte-specific G3PP-KO mice (LKO), by injecting AAV8-TBG-iCre to male G3PP^fl/fl mice. Controls received AAV8-TBG-eGFP. Both groups were fed chow diet for 10 weeks. Hyperglycemia (16-20 mM) was induced by glucose infusion for 55 h. Hepatocytes were isolated from normoglycemic mice for ex vivo studies and targeted metabolomics were measured in mice liver after glucose infusion.

RESULTS

LKO mice showed no change in body weight, food intake, fed and fasted glycemia but had increased fed plasma triglycerides. Hepatic glucose production from glycerol was increased in fasted LKO mice. LKO mouse hepatocytes displayed reduced glycerol production, elevated triglyceride and lactate production at high glucose concentration. Hyperglycemia in LKO mice led to increased liver weight and accumulation of triglycerides, glycogen and cholesterol together with elevated levels of Gro3P, dihydroxyacetone phosphate, acetyl-CoA and some Krebs cycle intermediates in liver. Hyperglycemic LKO mouse liver showed elevated expression of proinflammatory cytokines and M1-macrophage markers accompanied by increased plasma triglycerides, LDL/VLDL, urea and uric acid and myocardial triglycerides.

CONCLUSIONS

The glycerol shunt orchestrated by G3PP acts as a glucose excess detoxification pathway in hepatocytes by preventing metabolic disturbances that contribute to enhanced liver fat, glycogen storage, inflammation and lipid build-up in the heart. We propose G3PP as a novel therapeutic target for hepatic disorders linked to nutrient excess.

Science fiction authors' perspectives on human genetic engineering

Participants in the human gene editing debate often consider examples from science fiction but have rarely engaged directly with the science fiction community as stakeholders. To understand how science fiction authors develop and spread their views on gene editing, we created an online questionnaire that was answered by 78 authors, including 71 who had previously written about genetic engineering. When asked which ethical issues science fiction should explore, respondents most frequently mentioned affordability, new social divisions, consent and unforeseen safety risks. They rarely advocated exploring psychological effects or religious objections. When asked which works of fiction had influenced their perceptions of gene editing, the most frequent responses were the film Gattaca, the Star Trek franchise and the novels The Island of Doctor Moreau and Brave New World Unlike other stakeholders, they rarely cited Frankenstein as an influence. This article examines several differences between bioethicists, the general public and science fiction authors, and discusses how this community's involvement might benefit proponents and opponents of gene editing. It also provides an overview of works mentioned by our respondents that might serve as useful references in the debate.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

Genome-wide landscape establishes novel association signals for metabolic traits in the Arab population

While the Arabian population has a high prevalence of metabolic disorders, it has not been included in global studies that identify genetic risk loci for metabolic traits. Determining the transferability of such largely Euro-centric established risk loci is essential to transfer the research tools/resources, and drug targets generated by global studies to a broad range of ethnic populations. Further, consideration of populations such as Arabs, that are characterized by consanguinity and a high level of inbreeding, can lead to identification of novel risk loci. We imputed published GWAS data from two Kuwaiti Arab cohorts (n = 1434 and 1298) to the 1000 Genomes Project haplotypes and performed meta-analysis for associations with 13 metabolic traits. We compared the observed association signals with those established for metabolic traits. Our study highlighted 70 variants from 9 different genes, some of which have established links to metabolic disorders. By relaxing the genome-wide significance threshold, we identified ‘novel’ risk variants from 11 genes for metabolic traits. Many novel risk variant association signals were observed at or borderline to genome-wide significance. Furthermore, 349 previously established variants from 187 genes were validated in our study. Pleiotropic effect of risk variants on multiple metabolic traits were observed. Fine-mapping illuminated rs7838666/CSMD1 rs1864163/CETP and rs112861901/[INTS10,LPL] as candidate causal variants influencing fasting plasma glucose and high-density lipoprotein levels. Computational functional analysis identified a variety of gene regulatory signals around several variants. This study enlarges the population ancestry diversity of available GWAS and elucidates new variants in an ethnic group burdened with metabolic disorders.

A Long Non-coding RNA, LOC157273, Is an Effector Transcript at the Chromosome 8p23.1-PPP1R3B Metabolic Traits and Type 2 Diabetes Risk Locus

Aims

Causal transcripts at genomic loci associated with type 2 diabetes (T2D) are mostly unknown. The chr8p23.1 variant rs4841132, associated with an insulin-resistant diabetes risk phenotype, lies in the second exon of a long non-coding RNA (lncRNA) gene, LOC157273, located 175 kilobases from PPP1R3B, which encodes a key protein regulating insulin-mediated hepatic glycogen storage in humans. We hypothesized that LOC157273 regulates expression of PPP1R3B in human hepatocytes.

Methods

We tested our hypothesis using Stellaris fluorescent in situ hybridization to assess subcellular localization of LOC157273; small interfering RNA (siRNA) knockdown of LOC157273, followed by RT-PCR to quantify LOC157273 and PPP1R3B expression; RNA-seq to quantify the whole-transcriptome gene expression response to LOC157273 knockdown; and an insulin-stimulated assay to measure hepatocyte glycogen deposition before and after knockdown.

Results

We found that siRNA knockdown decreased LOC157273 transcript levels by approximately 80%, increased PPP1R3B mRNA levels by 1.7-fold, and increased glycogen deposition by >50% in primary human hepatocytes. An A/G heterozygous carrier (vs. three G/G carriers) had reduced LOC157273 abundance due to reduced transcription of the A allele and increased PPP1R3B expression and glycogen deposition.

Conclusion

We show that the lncRNA LOC157273 is a negative regulator of PPP1R3B expression and glycogen deposition in human hepatocytes and a causal transcript at an insulin-resistant T2D risk locus.

Sex-specific effects of a microsatellite polymorphism on human growth hormone receptor gene expression

Growth hormone (GH) binds to its specific receptor (GHR) at the surface of target cells activating multiple signaling pathways implicated in growth and metabolism. Dysregulation of GHRs leads to pathophysiological states that most commonly affect stature. We previously showed the association of a polymorphic (n = 15-37) GT microsatellite in the human GHR gene promoter with short stature in a sex-specific manner. In the present study we evaluated the functional relevance of this polymorphism in regulating GHR expression. Using luciferase reporter assays, we found that the GT repeat had a significant cis regulatory effect in response to HIF1α and a potential repressor role following C/EBPβ stimulation. Using a digital PCR application to measure allelic imbalance (AI), we showed a high prevalence of AI (∼76%) at the GHR locus in lymphoblastoid cell lines (LCLs), with a significantly higher degree of imbalance in LCLs derived from males. Examination of expression of GHR as well as other members of the GH-IGF1 axis in the LCLs revealed significant associations of GHR, IGF1 and BCL2 expression with GT genotype in a sex-specific manner. Our results suggest that this GT microsatellite exerts both cis and trans effects in a sex-specific context, revealing a new mechanism by which GHR gene expression is regulated.

Characterization of a FOXG1:TLE1 transcriptional network in glioblastoma-initiating cells

Glioblastoma (GBM) is the most common and deadly malignant brain cancer of glial cell origin, with a median patient survival of less than 20 months. Transcription factors FOXG1 and TLE1 promote GBM propagation by supporting maintenance of brain tumour-initiating cells (BTICs) with stem-like properties. Here, we characterize FOXG1 and TLE1 target genes in GBM patient-derived BTICs using ChIP-Seq and RNA-Seq approaches. These studies identify 150 direct FOXG1 targets, several of which are also TLE1 targets, involved in cell proliferation, differentiation, survival, chemotaxis and angiogenesis. Negative regulators of NOTCH signalling, including CHAC1, are among the transcriptional repression targets of FOXG1:TLE1 complexes, suggesting a crosstalk between FOXG1:TLE1 and NOTCH-mediated pathways in GBM. These results provide previously unavailable insight into the transcriptional programs underlying the tumour-promoting functions of FOXG1:TLE1 in GBM.

Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis

BACKGROUND

Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes.

METHODS & FINDINGS

Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants.

CONCLUSIONS

As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa

OBJECTIVE

The authors conducted a genome-wide association study of anorexia nervosa and calculated genetic correlations with a series of psychiatric, educational, and metabolic phenotypes.

METHOD

Following uniform quality control and imputation procedures using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, the authors performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression was used to calculate genome-wide common variant heritability (single-nucleotide polymorphism [SNP]-based heritability [h²_SNP]), partitioned heritability, and genetic correlations (r_g) between anorexia nervosa and 159 other phenotypes.

RESULTS

Results were obtained for 10,641,224 SNPs and insertion-deletion variants with minor allele frequencies >1% and imputation quality scores >0.6. The h²_SNP of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. The authors identified one genome-wide significant locus on chromosome 12 (rs4622308) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high-density lipoprotein cholesterol, and significant negative genetic correlations were observed between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes.

CONCLUSIONS

Anorexia nervosa is a complex heritable phenotype for which this study has uncovered the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. The study results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

Independent test assessment using the extreme value distribution theory

The new generation of whole genome sequencing platforms offers great possibilities and challenges for dissecting the genetic basis of complex traits. With a very high number of sequence variants, a naïve multiple hypothesis threshold correction hinders the identification of reliable associations by the overreduction of statistical power. In this report, we examine 2 alternative approaches to improve the statistical power of a whole genome association study to detect reliable genetic associations. The approaches were tested using the Genetic Analysis Workshop 19 (GAW19) whole genome sequencing data. The first tested method estimates the real number of effective independent tests actually being performed in whole genome association project by the use of an extreme value distribution and a set of phenotype simulations. Given the familiar nature of the GAW19 data and the finite number of pedigree founders in the sample, the number of correlations between genotypes is greater than in a set of unrelated samples. Using our procedure, we estimate that the effective number represents only 15 % of the total number of independent tests performed. However, even using this corrected significance threshold, no genome-wide significant association could be detected for systolic and diastolic blood pressure traits.

The second approach implements a biological relevance-driven hypothesis tested by exploiting prior computational predictions on the effect of nonsynonymous genetic variants detected in a whole genome sequencing association study. This guided testing approach was able to identify 2 promising single-nucleotide polymorphisms (SNPs), 1 for each trait, targeting biologically relevant genes that could help shed light on the genesis of the human hypertension. The first gene, PFH14, associated with systolic blood pressure, interacts directly with genes involved in calcium-channel formation and the second gene, MAP4, encodes a microtubule-associated protein and had already been detected by previous genome-wide association study experiments conducted in an Asian population. Our results highlight the necessity of the development of alternative approached to improve the efficiency on the detection of reasonable candidate associations in whole genome sequencing studies.

Alu repeats as transcriptional regulatory platforms in macrophage responses to M. tuberculosis infection

To understand the epigenetic regulation of transcriptional response of macrophages during early-stage M. tuberculosis (Mtb) infection, we performed ChIPseq analysis of H3K4 monomethylation (H3K4me1), a marker of poised or active enhancers. De novo H3K4me1 peaks in infected cells were associated with genes implicated in host defenses and apoptosis. Our analysis revealed that 40% of de novo regions contained human/primate-specific Alu transposable elements, enriched in the AluJ and S subtypes. These contained several transcription factor binding sites, including those for members of the MEF2 and ATF families, and LXR and RAR nuclear receptors, all of which have been implicated in macrophage differentiation, survival, and responses to stress and infection. Combining bioinformatics, molecular genetics, and biochemical approaches, we linked genes adjacent to H3K4me1-associated Alu repeats to macrophage metabolic responses against Mtb infection. In particular, we show that LXRα signaling, which reduced Mtb viability 18-fold by altering cholesterol metabolism and enhancing macrophage apoptosis, can be initiated at response elements present in Alu repeats. These studies decipher the mechanism of early macrophage transcriptional responses to Mtb, highlighting the role of Alu element transposition in shaping human transcription programs during innate immunity.

Pancreatic β-Cell Dysfunction in Diet-Induced Obese Mice: Roles of AMP-Kinase, Protein Kinase Cε, Mitochondrial and Cholesterol Metabolism, and Alterations in Gene Expression

Diet induced obese (DIO) mice can be stratified according to their weight gain in response to high fat diet as low responders (LDR) and high responders (HDR). This allows the study of β-cell failure and the transitions to prediabetes (LDR) and early diabetes (HDR). C57BL/6N mice were fed for 8 weeks with a normal chow diet (ND) or a high fat diet and stratified as LDR and HDR. Freshly isolated islets from ND, LDR and HDR mice were studied ex-vivo for mitochondrial metabolism, AMPK activity and signalling, the expression and activity of key enzymes of energy metabolism, cholesterol synthesis, and mRNA profiling. Severely compromised glucose-induced insulin secretion in HDR islets, as compared to ND and LDR islets, was associated with suppressed AMP-kinase activity. HDR islets also showed reduced acetyl-CoA carboxylase activity and enhanced activity of 3-hydroxy-3-methylglutaryl-CoA reductase, which led respectively to elevated fatty acid oxidation and increased cholesterol biosynthesis. HDR islets also displayed mitochondrial membrane hyperpolarization and reduced ATP turnover in the presence of elevated glucose. Expression of protein kinase Cε, which reduces both lipolysis and production of signals for insulin secretion, was elevated in DIO islets. Genes whose expression increased or decreased by more than 1.2-fold were minor between LDR and ND islets (17 differentially expressed), but were prominent between HDR and ND islets (1508 differentially expressed). In HDR islets, particularly affected genes were related to cell cycle and proliferation, AMPK signaling, mitochondrial metabolism and cholesterol metabolism. In conclusion, chronically reduced AMPK activity, mitochondrial dysfunction, elevated cholesterol biosynthesis in islets, and substantial alterations in gene expression accompany β-cell failure in HDR islets. The β-cell compensation process in the prediabetic state (LDR) is largely independent of transcriptional adaptive changes, whereas the transition to early diabetes (HDR) is associated with major alterations in gene expression.

Development of a platform for single cell genomics using convex lens-induced confinement

We demonstrate a lab-on-a-chip that combines micro/nano-fabricated features with a Convex Lens-Induced Confinement (CLIC) device for the in situ analysis of single cells. A complete cycle of single cell analysis was achieved that includes: cell trapping, cell isolation, lysis, protein digestion, genomic DNA extraction and on-chip genomic DNA linearization. The ability to dynamically alter the flow-cell dimensions using the CLIC method was coupled with a flow-control mechanism for achieving efficient cell trapping, buffer exchange, and loading of long DNA molecules into nanofluidic arrays. Finite element simulation of fluid flow gives rise to optimized design parameters for overcoming the high hydraulic resistance present in the micro/nano-confinement region. By tuning design parameters such as the pressure gradient and CLIC confinement, an efficient on-chip single cell analysis protocol can be obtained. We demonstrate that we can extract Mbp long genomic DNA molecules from a single human lybphoblastoid cell and stretch these molecules in the nanochannels for optical interrogation.

A Mendelian randomization study of the effect of type-2 diabetes on coronary heart disease

In observational studies, type-2 diabetes (T2D) is associated with an increased risk of coronary heart disease (CHD), yet interventional trials have shown no clear effect of glucose-lowering on CHD. Confounding may have therefore influenced these observational estimates. Here we use Mendelian randomization to obtain unconfounded estimates of the influence of T2D and fasting glucose (FG) on CHD risk. Using multiple genetic variants associated with T2D and FG, we find that risk of T2D increases CHD risk (odds ratio (OR)=1.11 (1.05–1.17), per unit increase in odds of T2D, P=8.8 × 10⁻⁵; using data from 34,840/114,981 T2D cases/controls and 63,746/130,681 CHD cases/controls). FG in non-diabetic individuals tends to increase CHD risk (OR=1.15 (1.00–1.32), per mmol·per l, P=0.05; 133,010 non-diabetic individuals and 63,746/130,681 CHD cases/controls). These findings provide evidence supporting a causal relationship between T2D and CHD and suggest that long-term trials may be required to discern the effects of T2D therapies on CHD risk.

In order to effectively design interventions, it is useful to understand the complex interplay between multiple syndromes. Here, Ahmad et al. use genome-wide association study data and Mendelian randomisation to examine the influence of Type 2 diabetes and fasting glucose levels on coronary heart disease.

Association of levels of fasting glucose and insulin with rare variants at the chromosome 11p11.2-MADD locus: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study

BACKGROUND

Common variation at the 11p11.2 locus, encompassing MADD, ACP2, NR1H3, MYBPC3, and SPI1, has been associated in genome-wide association studies with fasting glucose and insulin (FI). In the Cohorts for Heart and Aging Research in Genomic Epidemiology Targeted Sequencing Study, we sequenced 5 gene regions at 11p11.2 to identify rare, potentially functional variants influencing fasting glucose or FI levels.

METHODS AND RESULTS

Sequencing (mean depth, 38×) across 16.1 kb in 3566 individuals without diabetes mellitus identified 653 variants, 79.9% of which were rare (minor allele frequency <1%) and novel. We analyzed rare variants in 5 gene regions with FI or fasting glucose using the sequence kernel association test. At NR1H3, 53 rare variants were jointly associated with FI (P=2.73×10(-3)); of these, 7 were predicted to have regulatory function and showed association with FI (P=1.28×10(-3)). Conditioning on 2 previously associated variants at MADD (rs7944584, rs10838687) did not attenuate this association, suggesting that there are >2 independent signals at 11p11.2. One predicted regulatory variant, chr11:47227430 (hg18; minor allele frequency=0.00068), contributed 20.6% to the overall sequence kernel association test score at NR1H3, lies in intron 2 of NR1H3, and is a predicted binding site for forkhead box A1 (FOXA1), a transcription factor associated with insulin regulation. In human HepG2 hepatoma cells, the rare chr11:47227430 A allele disrupted FOXA1 binding and reduced FOXA1-dependent transcriptional activity.

CONCLUSIONS

Sequencing at 11p11.2-NR1H3 identified rare variation associated with FI. One variant, chr11:47227430, seems to be functional, with the rare A allele reducing transcription factor FOXA1 binding and FOXA1-dependent transcriptional activity.

Increased subcutaneous adipose tissue expression of genes involved in glycerolipid-fatty acid cycling in obese insulin-resistant versus -sensitive individuals

CONTEXT

A subpopulation of obese individuals remains insulin sensitive (ISO). They represent a unique human model to investigate factors underlying insulin resistance (IR) without the confounding effect of major differences in weight/adiposity. Altered fatty-acid (FA) metabolism in sc adipose tissue (SAT) contributes to obesity-associated IR.

OBJECTIVE

To test the hypothesis that ISO and body mass index-matched insulin-resistant obese (IRO) patients demonstrate differential SAT expression profiles of genes involved in glycerolipid-FA metabolism and that weight loss-induced improvement of IR ameliorates these changes.

DESIGN AND SETTING

A cross-sectional and longitudinal study.

PATIENTS AND INTERVENTION

Thirty-eight nondiabetic obese women were stratified into ISO (n = 25) or IRO (n = 13) groups based on hyperinsulinemic-euglycemic clamp results. Subjects were studied before and after a 6-month hypocaloric diet intervention.

MAIN OUTCOME MEASURES

mRNA (quantitative RT-PCR) and protein (mass spectrometry and immunoblots) levels were measured in SAT biopsies.

RESULTS

Despite having age, body mass index, and fat mass similar to ISO individuals, IRO patients had lower insulin sensitivity and glucose tolerance (P < .05). Baseline SAT mRNA and protein levels of genes involved in both the synthesis and lipolysis of glycerolipid-FAs were higher in IRO individuals (P < .05), even when groups were matched for visceral adipose tissue content. The dietary intervention resulted in approximately 6% weight loss in both the IRO and ISO groups (P < .05) but only ameliorated insulin sensitivity in IRO individuals (P < .05). Likewise, the intervention reduced the expression of most glycerolipid-FA metabolism genes (P < .05), with expression levels in IRO individuals being restored to ISO levels.

CONCLUSIONS

Increased SAT expression of genes involved in both the synthesis and hydrolysis of glycerolipid-FAs is closely associated with IR in obese women. The results suggest that enhanced glycerolipid-FA cycling in SAT contributes to obesity-associated IR.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Mapping of a chromosome 12 region associated with airway hyperresponsiveness in a recombinant congenic mouse strain and selection of potential candidate genes by expression and sequence variation analyses

In a previous study we determined that BcA86 mice, a strain belonging to a panel of AcB/BcA recombinant congenic strains, have an airway responsiveness phenotype resembling mice from the airway hyperresponsive A/J strain. The majority of the BcA86 genome is however from the hyporesponsive C57BL/6J strain. The aim of this study was to identify candidate regions and genes associated with airway hyperresponsiveness (AHR) by quantitative trait locus (QTL) analysis using the BcA86 strain. Airway responsiveness of 205 F2 mice generated from backcrossing BcA86 strain to C57BL/6J strain was measured and used for QTL analysis to identify genomic regions in linkage with AHR. Consomic mice for the QTL containing chromosomes were phenotyped to study the contribution of each chromosome to lung responsiveness. Candidate genes within the QTL were selected based on expression differences in mRNA from whole lungs, and the presence of coding non-synonymous mutations that were predicted to have a functional effect by amino acid substitution prediction tools. One QTL for AHR was identified on Chromosome 12 with its 95% confidence interval ranging from 54.6 to 82.6 Mbp and a maximum LOD score of 5.11 (p = 3.68×10⁻³). We confirmed that the genotype of mouse Chromosome 12 is an important determinant of lung responsiveness using a Chromosome 12 substitution strain. Mice with an A/J Chromosome 12 on a C57BL/6J background have an AHR phenotype similar to hyperresponsive strains A/J and BcA86. Within the QTL, genes with deleterious coding variants, such as Foxa1, and genes with expression differences, such as Mettl21d and Snapc1, were selected as possible candidates for the AHR phenotype. Overall, through QTL analysis of a recombinant congenic strain, microarray analysis and coding variant analysis we identified Chromosome 12 and three potential candidate genes to be in linkage with airway responsiveness.

Fusion of TTYH1 with the C19MC microRNA cluster drives expression of a brain-specific DNMT3B isoform in the embryonal brain tumor ETMR

Embryonal tumors with multilayered rosettes (ETMRs) are rare, deadly pediatric brain tumors characterized by high-level amplification of the microRNA cluster C19MC. We performed integrated genetic and epigenetic analyses of 12 ETMR samples and identified, in all cases, C19MC fusions to TTYH1 driving expression of the microRNAs. ETMR tumors, cell lines and xenografts showed a specific DNA methylation pattern distinct from those of other tumors and normal tissues. We detected extreme overexpression of a previously uncharacterized isoform of DNMT3B originating at an alternative promoter that is active only in the first weeks of neural tube development. Transcriptional and immunohistochemical analyses suggest that C19MC-dependent DNMT3B deregulation is mediated by RBL2, a known repressor of DNMT3B. Transfection with individual C19MC microRNAs resulted in DNMT3B upregulation and RBL2 downregulation in cultured cells. Our data suggest a potential oncogenic re-engagement of an early developmental program in ETMR via epigenetic alteration mediated by an embryonic, brain-specific DNMT3B isoform.

An electrochemically controlled microcantilever biosensor

An oligonucleotide-based electrochemically controlled gold-coated microcantilever biosensor that can transduce specific biomolecular interactions is reported. The derivatized microcantilever exhibits characteristic surface stress time course patterns in response to an externally applied periodic square wave potential. Experiments demonstrate that control of the surface charge density with an electrode potential is essential to producing a sensor that exhibits large, reproducible surface stress changes. The time course of surface stress changes are proposed to be linked to an electrochemically mediated competition between the adsorption of solution-based ions and the single- or double-stranded oligonucleotides tethered to the gold surface. A similar potential-actuated change in surface stress also results from the interaction between an oligonucleotide aptamer and its cognate ligand, demonstrating the broad applicability of this methodology.

Denaturation mapping of Saccharomyces cerevisiae

Optical mapping of DNA provides large-scale genomic information that can be used to assemble contigs from next-generation sequencing, and to detect rearrangements between single cells. A recent optical mapping technique called denaturation mapping has the advantage of using physical principles rather than the action of enzymes to probe genomic structure. Denaturation mapping uses fluorescence microscopy to image the pattern of partial melting along a DNA molecule extended in a channel of cross-section 120 nm at the heart of a nanofluidic device. We used denaturation mapping to locate single DNA molecules on the yeast genome (12.1 Mbp) by comparing images to a computationally predicted map for the entire genome sequence. By locating 84 molecules we assembled an optical map of the yeast genome with > 50% coverage.

Mapping of clinical and expression quantitative trait loci in a sex-dependent effect of host susceptibility to mouse-adapted influenza H3N2/HK/1/68

Seasonal influenza outbreaks and recurrent influenza pandemics present major challenges to public health. By studying immunological responses to influenza in different host species, it may be possible to discover common mechanisms of susceptibility in response to various influenza strains. This could lead to novel therapeutic targets with wide clinical application. Using a mouse-adapted strain of influenza (A/HK/1/68-MA20 [H3N2]), we produced a mouse model of severe influenza that reproduces the hallmark high viral load and overexpression of cytokines associated with susceptibility to severe influenza in humans. We mapped genetic determinants of the host response using a panel of 29 closely related mouse strains (AcB/BcA panel of recombinant congenic strains) created from influenza-susceptible A/J and influenza-resistant C57BL/6J (B6) mice. Combined clinical quantitative trait loci (QTL) and lung expression QTL mapping identified candidate genes for two sex-specific QTL on chromosomes 2 and 17. The former includes the previously described Hc gene, a deficit of which is associated with the susceptibility phenotype in females. The latter includes the phospholipase gene Pla2g7 and Tnfrsf21, a member of the TNFR superfamily. Confirmation of the gene underlying the chromosome 17 QTL may reveal new strategies for influenza treatment.

An N-ethyl-N-nitrosourea (ENU)-induced dominant negative mutation in the JAK3 kinase protects against cerebral malaria

Cerebral malaria (CM) is a lethal neurological complication of malaria. We implemented a genome-wide screen in mutagenized mice to identify host proteins involved in CM pathogenesis and whose inhibition may be of therapeutic value. One pedigree (P48) segregated a resistance trait whose CM-protective effect was fully penetrant, mapped to chromosome 8, and identified by genome sequencing as homozygosity for a mis-sense mutation (W81R) in the FERM domain of Janus-associated kinase 3 (Jak3). The causative effect of Jak3(W81R) was verified by complementation testing in Jak3(W81R/-) double heterozygotes that were fully protected against CM. Jak3(W81R) homozygotes showed defects in thymic development with depletion of CD8(+) T cell, B cell, and NK cell compartments, and defective T cell-dependent production of IFN-γ. Adoptive transfer of normal splenocytes abrogates CM resistance in Jak3(W81R) homozygotes, an effect attributed to the CD8(+) T cells. Jak3(W81R) behaves as a dominant negative variant, with significant CM resistance of Jak3(W81R/+) heterozygotes, compared to CM-susceptible Jak3(+/+) and Jak3(+/-) controls. CM resistance in Jak3(W81R/+) heterozygotes occurs in presence of normal T, B and NK cell numbers. These findings highlight the pathological role of CD8(+) T cells and Jak3-dependent IFN-γ-mediated Th1 responses in CM pathogenesis.

Rare MTNR1B variants impairing melatonin receptor 1B function contribute to type 2 diabetes

Genome-wide association studies have revealed that common noncoding variants in MTNR1B (encoding melatonin receptor 1B, also known as MT(2)) increase type 2 diabetes (T2D) risk(1,2). Although the strongest association signal was highly significant (P < 1 × 10(-20)), its contribution to T2D risk was modest (odds ratio (OR) of ∼1.10-1.15)(1-3). We performed large-scale exon resequencing in 7,632 Europeans, including 2,186 individuals with T2D, and identified 40 nonsynonymous variants, including 36 very rare variants (minor allele frequency (MAF) <0.1%), associated with T2D (OR = 3.31, 95% confidence interval (CI) = 1.78-6.18; P = 1.64 × 10(-4)). A four-tiered functional investigation of all 40 mutants revealed that 14 were non-functional and rare (MAF < 1%), and 4 were very rare with complete loss of melatonin binding and signaling capabilities. Among the very rare variants, the partial- or total-loss-of-function variants but not the neutral ones contributed to T2D (OR = 5.67, CI = 2.17-14.82; P = 4.09 × 10(-4)). Genotyping the four complete loss-of-function variants in 11,854 additional individuals revealed their association with T2D risk (8,153 individuals with T2D and 10,100 controls; OR = 3.88, CI = 1.49-10.07; P = 5.37 × 10(-3)). This study establishes a firm functional link between MTNR1B and T2D risk.

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

Common variants at 10 genomic loci influence hemoglobin A₁(C) levels via glycemic and nonglycemic pathways

OBJECTIVE

Glycated hemoglobin (HbA₁(c)), used to monitor and diagnose diabetes, is influenced by average glycemia over a 2- to 3-month period. Genetic factors affecting expression, turnover, and abnormal glycation of hemoglobin could also be associated with increased levels of HbA₁(c). We aimed to identify such genetic factors and investigate the extent to which they influence diabetes classification based on HbA₁(c) levels.

RESEARCH DESIGN AND METHODS

We studied associations with HbA₁(c) in up to 46,368 nondiabetic adults of European descent from 23 genome-wide association studies (GWAS) and 8 cohorts with de novo genotyped single nucleotide polymorphisms (SNPs). We combined studies using inverse-variance meta-analysis and tested mediation by glycemia using conditional analyses. We estimated the global effect of HbA₁(c) loci using a multilocus risk score, and used net reclassification to estimate genetic effects on diabetes screening.

RESULTS

Ten loci reached genome-wide significant association with HbA(1c), including six new loci near FN3K (lead SNP/P value, rs1046896/P = 1.6 × 10⁻²⁶), HFE (rs1800562/P = 2.6 × 10⁻²⁰), TMPRSS6 (rs855791/P = 2.7 × 10⁻¹⁴), ANK1 (rs4737009/P = 6.1 × 10⁻¹²), SPTA1 (rs2779116/P = 2.8 × 10⁻⁹) and ATP11A/TUBGCP3 (rs7998202/P = 5.2 × 10⁻⁹), and four known HbA₁(c) loci: HK1 (rs16926246/P = 3.1 × 10⁻⁵⁴), MTNR1B (rs1387153/P = 4.0 × 10⁻¹¹), GCK (rs1799884/P = 1.5 × 10⁻²⁰) and G6PC2/ABCB11 (rs552976/P = 8.2 × 10⁻¹⁸). We show that associations with HbA₁(c) are partly a function of hyperglycemia associated with 3 of the 10 loci (GCK, G6PC2 and MTNR1B). The seven nonglycemic loci accounted for a 0.19 (% HbA₁(c)) difference between the extreme 10% tails of the risk score, and would reclassify ∼2% of a general white population screened for diabetes with HbA₁(c).

CONCLUSIONS

GWAS identified 10 genetic loci reproducibly associated with HbA₁(c). Six are novel and seven map to loci where rarer variants cause hereditary anemias and iron storage disorders. Common variants at these loci likely influence HbA₁(c) levels via erythrocyte biology, and confer a small but detectable reclassification of diabetes diagnosis by HbA₁(c).

Molecular diagnosis of neonatal diabetes mellitus using next-generation sequencing of the whole exome

Background

Accurate molecular diagnosis of monogenic non-autoimmune neonatal diabetes mellitus (NDM) is critical for patient care, as patients carrying a mutation in KCNJ11 or ABCC8 can be treated by oral sulfonylurea drugs instead of insulin therapy. This diagnosis is currently based on Sanger sequencing of at least 42 PCR fragments from the KCNJ11, ABCC8, and INS genes. Here, we assessed the feasibility of using the next-generation whole exome sequencing (WES) for the NDM molecular diagnosis.

Methodology/Principal Findings

We carried out WES for a patient presenting with permanent NDM, for whom mutations in KCNJ11, ABCC8 and INS and abnormalities in chromosome 6q24 had been previously excluded. A solution hybridization selection was performed to generate WES in 76 bp paired-end reads, by using two channels of the sequencing instrument. WES quality was assessed using a high-resolution oligonucleotide whole-genome genotyping array. From our WES with high-quality reads, we identified a novel non-synonymous mutation in ABCC8 (c.1455G>C/p.Q485H), despite a previous negative sequencing of this gene. This mutation, confirmed by Sanger sequencing, was not present in 348 controls and in the patient's mother, father and young brother, all of whom are normoglycemic.

Conclusions/Significance

WES identified a novel de novo ABCC8 mutation in a NDM patient. Compared to the current Sanger protocol, WES is a comprehensive, cost-efficient and rapid method to identify mutations in NDM patients. We suggest WES as a near future tool of choice for further molecular diagnosis of NDM cases, negative for chr6q24, KCNJ11 and INS abnormalities.

Characterizing the profile of obese patients who are metabolically healthy

The presence of obesity-related metabolic disturbances varies widely among obese individuals. Accordingly, a unique subset of obese individuals has been described in the medical literature, which seems to be protected or more resistant to the development of metabolic abnormalities associated with obesity. These individuals, now known as 'metabolically healthy but obese' (MHO), despite having excessive body fatness, display a favorable metabolic profile characterized by high levels of insulin sensitivity, no hypertension as well as a favorable lipid, inflammation, hormonal, liver enzyme and immune profile. However, recent studies have indicated that this healthier metabolic profile may not translate into a lower risk for mortality. Mechanisms that could explain the favorable metabolic profile of MHO individuals are poorly understood. However, preliminary evidence suggests that differences in visceral fat accumulation, birth weight, adipose cell size and gene expression-encoding markers of adipose cell differentiation may favor the development of the MHO phenotype. Despite the uncertainty regarding the exact degree of protection related to the MHO status, identification of underlying factors and mechanisms associated with this phenotype will eventually be invaluable in helping us understand factors that predispose, delay or protect obese individuals from metabolic disturbances. Collectively, a greater understanding of the MHO individual has important implications for therapeutic decision making, the characterization of subjects in research protocols and medical education.

Vav-1 expression correlates with NFκB activation and CD40-mediated cell death in diffuse large B-cell lymphoma cell lines

Diffuse large B-cell lymphoma (DLBCL) is an aggressive malignancy with a variable response to therapy. We have previously shown that DLBCL cell lines differ in their susceptibility to CD40-mediated cell death, and that resistance to CD40-targeted antibodies correlated with increased expression of markers of immature B-cell and absence of Vav-1 mRNA. We used gene expression profiling to investigate the mechanism of CD40 resistance in these cell lines, and found that resistance correlated with lack of Vav-1 and inability to activate NFκB upon CD40 ligation. Analysis of tissue microarrays of 213 DLBCL cases revealed that Vav-1 expression correlated with a higher proliferative index and the presence of the post-germinal centre marker Irf-4. Our results suggest that Vav-1 expression may be associated with activated B-cell DLBCL origin and higher proliferative activity, and indicate Vav-1 as a potential marker to identify tumours likely to respond to CD40-targeted therapies.

Glucolipotoxicity alters lipid partitioning and causes mitochondrial dysfunction, cholesterol, and ceramide deposition and reactive oxygen species production in INS832/13 ss-cells

Elevated glucose and saturated fatty acids synergize in inducing apoptosis in INS832/13 cells and in human islet cells. In order to gain insight into the molecular mechanism(s) of glucolipotoxicity (Gltox), gene profiling and metabolic analyses were performed in INS832/13 cells cultured at 5 or 20 mm glucose in the absence or presence of palmitate. Expression changes were observed for transcripts involved in mitochondrial, lipid, and glucose metabolism. At 24 h after Gltox, increased expression of lipid partitioning genes suggested a promotion of fatty acid esterification and reduced lipid oxidation/detoxification, whereas changes in the expression of energy metabolism genes suggested mitochondrial dysfunction. These changes were associated with decreased glucose-induced insulin secretion, total insulin content, ATP levels, AMP-kinase activity, mitochondrial membrane potential and fat oxidation, unchanged de novo fatty acid synthesis, and increased reactive oxygen species, cholesterol, ceramide, and triglyceride levels. However, the synergy between elevated glucose and palmitate to cause ss-cell toxicity in term of apoptosis and reduced glucose-induced insulin secretion only correlated with triglyceride and ceramide depositions. Overexpression of endoplasmic reticulum glycerol-3-phosphate acyl transferase to enhance lipid esterification amplified Gltox at intermediate glucose (11 mm), whereas reducing acetyl-coenzyme A carboxylase 1 expression by small interfering RNA to shift lipid partitioning to fat oxidation reduced Gltox. The results suggest that Gltox entails alterations in lipid partitioning, sterol and ceramide accumulation, mitochondrial dysfunction, and reactive oxygen species production, all contributing to altering ss-cell function. The data also suggest that the early promotion of lipid esterification processes is instrumental in the Gltox process.

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

A genome-wide association study reveals variants in ARL15 that influence adiponectin levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P< or =5x10(-8)). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P< or =0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2x10(-19) for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9x10(-8), n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5x10(-6), n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2x10(-3), n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Genetic variant in HK1 is associated with a proanemic state and A1C but not other glycemic control-related traits

OBJECTIVE

A1C is widely considered the gold standard for monitoring effective blood glucose levels. Recently, a genome-wide association study reported an association between A1C and rs7072268 within HK1 (encoding hexokinase 1), which catalyzes the first step of glycolysis. HK1 deficiency in erythrocytes (red blood cells [RBCs]) causes severe nonspherocytic hemolytic anemia in both humans and mice.

RESEARCH DESIGN AND METHODS

The contribution of rs7072268 to A1C and the RBC-related traits was assessed in 6,953 nondiabetic European participants. We additionally analyzed the association with hematologic traits in 5,229 nondiabetic European individuals (in whom A1C was not measured) and 1,924 diabetic patients. Glucose control-related markers other than A1C were analyzed in 18,694 nondiabetic European individuals. A type 2 diabetes case-control study included 7,447 French diabetic patients.

RESULTS

Our study confirms a strong association between the rs7072268-T allele and increased A1C (beta = 0.029%; P = 2.22 x 10(-7)). Surprisingly, despite adequate study power, rs7072268 showed no association with any other markers of glucose control (fasting- and 2-h post-OGTT-related parameters, n = 18,694). In contrast, rs7072268-T allele decreases hemoglobin levels (n = 13,416; beta = -0.054 g/dl; P = 3.74 x 10(-6)) and hematocrit (n = 11,492; beta = -0.13%; P = 2.26 x 10(-4)), suggesting a proanemic effect. The T allele also increases risk for anemia (836 cases; odds ratio 1.13; P = 0.018).

CONCLUSIONS

HK1 variation, although strongly associated with A1C, does not seem to be involved in blood glucose control. Since HK1 rs7072268 is associated with reduced hemoglobin levels and favors anemia, we propose that HK1 may influence A1C levels through its anemic effect or its effect on glucose metabolism in RBCs. These findings may have implications for type 2 diabetes diagnosis and clinical management because anemia is a frequent complication of the diabetes state.

Genetic variant near IRS1 is associated with type 2 diabetes, insulin resistance and hyperinsulinemia

Genome-wide association studies have identified common variants that only partially explain the genetic risk for type 2 diabetes (T2D). Using genome-wide association data from 1,376 French individuals, we identified 16,360 SNPs nominally associated with T2D and studied these SNPs in an independent sample of 4,977 French individuals. We then selected the 28 best hits for replication in 7,698 Danish subjects and identified 4 SNPs showing strong association with T2D, one of which (rs2943641, P = 9.3 x 10(-12), OR = 1.19) was located adjacent to the insulin receptor substrate 1 gene (IRS1). Unlike previously reported T2D risk loci, which predominantly associate with impaired beta cell function, the C allele of rs2943641 was associated with insulin resistance and hyperinsulinemia in 14,358 French, Danish and Finnish participants from population-based cohorts; this allele was also associated with reduced basal levels of IRS1 protein and decreased insulin induction of IRS1-associated phosphatidylinositol-3-OH kinase activity in human skeletal muscle biopsies.