Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D.

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

AKR1D1 is a novel regulator of metabolic phenotype in human hepatocytes and is dysregulated in non-alcoholic fatty liver disease.

OBJECTIVE: Non-alcoholic fatty liver disease (NAFLD) is the hepatic manifestation of metabolic syndrome. Steroid hormones and bile acids are potent regulators of hepatic carbohydrate and lipid metabolism. Steroid 5ß-reductase (AKR1D1) is highly expressed in human liver where it inactivates steroid hormones and catalyzes a fundamental step in bile acid synthesis. METHODS: Human liver biopsies were obtained from 34 obese patients and AKR1D1 mRNA expression levels were measured using qPCR. Genetic manipulation of AKR1D1 was performed in human HepG2 and Huh7 liver cell lines. Metabolic assessments were made using transcriptome analysis, western blotting, mass spectrometry, clinical biochemistry, and enzyme immunoassays. RESULTS: In human liver biopsies, AKR1D1 expression decreased with advancing steatosis, fibrosis and inflammation. Expression was decreased in patients with type 2 diabetes. In human liver cell lines, AKR1D1 knockdown decreased primary bile acid biosynthesis and steroid hormone clearance. RNA-sequencing identified disruption of key metabolic pathways, including insulin action and fatty acid metabolism. AKR1D1 knockdown increased hepatocyte triglyceride accumulation, insulin sensitivity, and glycogen synthesis, through increased de novo lipogenesis and decreased ß-oxidation, fueling hepatocyte inflammation. Pharmacological manipulation of bile acid receptor activation prevented the induction of lipogenic and carbohydrate genes, suggesting that the observed metabolic phenotype is driven through bile acid rather than steroid hormone availability. CONCLUSIONS: Genetic manipulation of AKR1D1 regulates the metabolic phenotype of human hepatoma cell lines, driving steatosis and inflammation. Taken together, the observation that AKR1D1 mRNA is down-regulated with advancing NAFLD suggests that it may have a crucial role in the pathogenesis and progression of the disease.

Diabetes causes marked inhibition of mitochondrial metabolism in pancreatic ß-cells.

Diabetes is a global health problem caused primarily by the inability of pancreatic ß-cells to secrete adequate levels of insulin. The molecular mechanisms underlying the progressive failure of ß-cells to respond to glucose in type-2 diabetes remain unresolved. Using a combination of transcriptomics and proteomics, we find significant dysregulation of major metabolic pathways in islets of diabetic ßV59M mice, a non-obese, eulipidaemic diabetes model. Multiple genes/proteins involved in glycolysis/gluconeogenesis are upregulated, whereas those involved in oxidative phosphorylation are downregulated. In isolated islets, glucose-induced increases in NADH and ATP are impaired and both oxidative and glycolytic glucose metabolism are reduced. INS-1 ß-cells cultured chronically at high glucose show similar changes in protein expression and reduced glucose-stimulated oxygen consumption: targeted metabolomics reveals impaired metabolism. These data indicate hyperglycaemia induces metabolic changes in ß-cells that markedly reduce mitochondrial metabolism and ATP synthesis. We propose this underlies the progressive failure of ß-cells in diabetes.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk.

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Systematic Functional Characterization of Candidate Causal Genes for Type 2 Diabetes Risk Variants.

Most genetic association signals for type 2 diabetes risk are located in noncoding regions of the genome, hindering translation into molecular mechanisms. Physiological studies have shown a majority of disease-associated variants to exert their effects through pancreatic islet dysfunction. Systematically characterizing the role of regional transcripts in ß-cell function could identify the underlying disease-causing genes, but large-scale studies in human cellular models have previously been impractical. We developed a robust and scalable strategy based on arrayed gene silencing in the human ß-cell line EndoC-ßH1. In a screen of 300 positional candidates selected from 75 type 2 diabetes regions, each gene was assayed for effects on multiple disease-relevant phenotypes, including insulin secretion and cellular proliferation. We identified a total of 45 genes involved in ß-cell function, pointing to possible causal mechanisms at 37 disease-associated loci. The results showed a strong enrichment for genes implicated in monogenic diabetes. Selected effects were validated in a follow-up study, including several genes (ARL15, ZMIZ1, and THADA) with previously unknown or poorly described roles in ß-cell biology. We have demonstrated the feasibility of systematic functional screening in a human ß-cell model and successfully prioritized plausible disease-causing genes at more than half of the regions investigated.

Isocitrate-to-SENP1 signaling amplifies insulin secretion and rescues dysfunctional ß cells.

Insulin secretion from ß cells of the pancreatic islets of Langerhans controls metabolic homeostasis and is impaired in individuals with type 2 diabetes (T2D). Increases in blood glucose trigger insulin release by closing ATP-sensitive K+ channels, depolarizing ß cells, and opening voltage-dependent Ca2+ channels to elicit insulin exocytosis. However, one or more additional pathway(s) amplify the secretory response, likely at the distal exocytotic site. The mitochondrial export of isocitrate and engagement with cytosolic isocitrate dehydrogenase (ICDc) may be one key pathway, but the mechanism linking this to insulin secretion and its role in T2D have not been defined. Here, we show that the ICDc-dependent generation of NADPH and subsequent glutathione (GSH) reduction contribute to the amplification of insulin exocytosis via sentrin/SUMO-specific protease-1 (SENP1). In human T2D and an in vitro model of human islet dysfunction, the glucose-dependent amplification of exocytosis was impaired and could be rescued by introduction of signaling intermediates from this pathway. Moreover, islet-specific Senp1 deletion in mice caused impaired glucose tolerance by reducing the amplification of insulin exocytosis. Together, our results identify a pathway that links glucose metabolism to the amplification of insulin secretion and demonstrate that restoration of this axis rescues ß cell function in T2D.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease.

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.

PTEN mutations as a cause of constitutive insulin sensitivity and obesity.

BACKGROUND: Epidemiologic and genetic evidence links type 2 diabetes, obesity, and cancer. The tumor-suppressor phosphatase and tensin homologue (PTEN) has roles in both cellular growth and metabolic signaling. Germline PTEN mutations cause a cancer-predisposition syndrome, providing an opportunity to study the effect of PTEN haploinsufficiency in humans. METHODS: We measured insulin sensitivity and beta-cell function in 15 PTEN mutation carriers and 15 matched controls. Insulin signaling was measured in muscle and adipose-tissue biopsy specimens from 5 mutation carriers and 5 well-matched controls. We also assessed the effect of PTEN haploinsufficiency on obesity by comparing anthropometric indexes between the 15 patients and 2097 controls from a population-based study of healthy adults. Body composition was evaluated by means of dual-emission x-ray absorptiometry and skinfold thickness. RESULTS: Measures of insulin resistance were lower in the patients with a PTEN mutation than in controls (e.g., mean fasting plasma insulin level, 29 pmol per liter [range, 9 to 99] vs. 74 pmol per liter [range, 22 to 185]; P=0.001). This finding was confirmed with the use of hyperinsulinemic euglycemic clamping, showing a glucose infusion rate among carriers 2 times that among controls (P=0.009). The patients' insulin sensitivity could be explained by the presence of enhanced insulin signaling through the PI3K-AKT pathway, as evidenced by increased AKT phosphorylation. The PTEN mutation carriers were obese as compared with population-based controls (mean body-mass index [the weight in kilograms divided by the square of the height in meters], 32 [range, 23 to 42] vs. 26 [range, 15 to 48]; P<0.001). This increased body mass in the patients was due to augmented adiposity without corresponding changes in fat distribution. CONCLUSIONS: PTEN haploinsufficiency is a monogenic cause of profound constitutive insulin sensitization that is apparently obesogenic. We demonstrate an apparently divergent effect of PTEN mutations: increased risks of obesity and cancer but a decreased risk of type 2 diabetes owing to enhanced insulin sensitivity. (Funded by the Wellcome Trust and others.).