Electrophysiological characterisation of iPSC-derived human ß-like cells and an SLC30A8 disease model.

iPSC-derived human ß-like cells (BLC) hold promise for both therapy and disease modelling, but their generation remains challenging and their functional analyses beyond transcriptomic and morphological assessments remain limited. Here, we validate an approach using multicellular and single cell electrophysiological tools to evaluate function of BLCs from pioneer protocols that can be easily adapted to more differentiated BLCs. The Multi-Electrode Arrays (MEAs) measuring the extracellular electrical activity revealed that BLCs are electrically coupled, produce slow potential (SP) signals like primary ß-cells that are closely linked to insulin secretion. We also used high-resolution single-cell patch-clamp measurements to capture the exocytotic properties, and characterise voltage-gated sodium and calcium currents and found that they were comparable to those in primary ß and EndoC-ßH1 cells. The KATP channel conductance is greater than in human primary ß-cells which may account for the limited glucose responsiveness observed with MEA. We used MEAs to study the impact of the type 2 diabetes protective SLC30A8 allele (p.Lys34Serfs*50) and found that BLCs with this allele have stronger electrical coupling activity. Our data suggest that BLCs can be used to evaluate the functional impact of genetic variants on ß-cell function and coupling.

Electrophysiological characterisation of iPSC-derived human ß-like cells and an SLC30A8 disease model.

iPSC-derived human ß-like cells (BLC) hold promise for both therapy and disease modelling, but their generation remains challenging and their functional analyses beyond transcriptomic and morphological assessments remain limited. Here, we validate an approach using multicellular and single cell electrophysiological tools to evaluate BLCs functions. The Multi-Electrode Arrays (MEAs) measuring the extracellular electrical activity revealed that BLCs are electrically coupled, produce slow potential (SP) signals like primary ß-cells that are closely linked to insulin secretion. We also used high-resolution single-cell patch-clamp measurements to capture the exocytotic properties, and characterize voltage-gated sodium and calcium currents. These were comparable to those in primary ß and EndoC-ßH1 cells. The KATP channel conductance is greater than in human primary ß cells which may account for the limited glucose responsiveness observed with MEA. We used MEAs to study the impact of the type 2 diabetes protective SLC30A8 allele (p.Lys34Serfs*50) and found that BLCs with this allele have stronger electrical coupling. Our data suggest that with an adapted approach BLCs from pioneer protocol can be used to evaluate the functional impact of genetic variants on ß-cell function and coupling.

WFS1 protein expression correlates with clinical progression of optic atrophy in patients with Wolfram syndrome.

BACKGROUND: Wolfram syndrome (WFS) is a rare disorder characterised by childhood-onset diabetes mellitus and progressive optic atrophy. Most patients have variants in the WFS1 gene. We undertook functional studies of WFS1 variants and correlated these with WFS1 protein expression and phenotype. METHODS: 9 patients with a clinical diagnosis of WFS were studied with quantitative PCR for markers of endoplasmic reticulum (ER) stress and immunoblotting of fibroblast protein extracts for WFS1 protein expression. Luciferase reporter assay was used to assess ATF-6 dependent unfolded protein response (UPR) activation. RESULTS: 6 patients with compound heterozygous nonsense mutations in WFS1 had no detectable WFS1 protein expression; 3 patients with missense variants had 4%, 45% and 48% WFS1 protein expression. One of these also had an OPA1 mutation and was reclassified as autosomal dominant optic atrophy-plus syndrome. There were no correlations between ER stress marker mRNA and WFS1 protein expression. ERSE-luciferase reporter indicated activation of the ATF6 branch of UPR in two patients tested. Patients with partial WFS1 expression showed milder visual acuity impairment (asymptomatic or colour blind only), compared with those with absent expression (registered severe vision impaired) (p=0.04). These differences remained after adjusting for duration of optic atrophy. CONCLUSIONS: Patients with WFS who have partial WFS1 protein expression present with milder visual impairment. This suggests a protective effect of partial WFS1 protein expression on the severity and perhaps progression of vision impairment and that therapies to increase residual WFS1 protein expression may be beneficial.

Dimethyl fumarate reduces hepatocyte senescence following paracetamol exposure.

Liver disease is a major cause of premature death. Oxidative stress in the liver represents a key disease driver. Compounds, such as dimethyl fumarate (DMF), can activate the antioxidant response and are used clinically to treat disease. In this study, we tested the protective properties of DMF before or after paracetamol exposure. Following DMF administration, Nrf2 nuclear translocation was tracked at the single-cell level and target gene transactivation confirmed. Next, the protective properties of DMF were examined following paracetamol exposure. Transcriptomic and biochemical analysis revealed that DMF rescue was underpinned by reduced Nf-kB and TGF-ß signaling and cell senescence. Following on from these studies, we employed a Zebrafish model to study paracetamol exposure in vivo. We combined a genetically modified Zebrafish model, expressing green fluorescent protein exclusively in the liver, with automated microscopy. Pre-treatment with DMF, prior to paracetamol exposure, led to reduced liver damage in Zebrafish demonstrating protective properties.

Machine Learning based histology phenotyping to investigate the epidemiologic and genetic basis of adipocyte morphology and cardiometabolic traits.

Genetic studies have recently highlighted the importance of fat distribution, as well as overall adiposity, in the pathogenesis of obesity-associated diseases. Using a large study (n = 1,288) from 4 independent cohorts, we aimed to investigate the relationship between mean adipocyte area and obesity-related traits, and identify genetic factors associated with adipocyte cell size. To perform the first large-scale study of automatic adipocyte phenotyping using both histological and genetic data, we developed a deep learning-based method, the Adipocyte U-Net, to rapidly derive mean adipocyte area estimates from histology images. We validate our method using three state-of-the-art approaches; CellProfiler, Adiposoft and floating adipocytes fractions, all run blindly on two external cohorts. We observe high concordance between our method and the state-of-the-art approaches (Adipocyte U-net vs. CellProfiler: R2visceral = 0.94, P < 2.2 × 10-16, R2subcutaneous = 0.91, P < 2.2 × 10-16), and faster run times (10,000 images: 6mins vs 3.5hrs). We applied the Adipocyte U-Net to 4 cohorts with histology, genetic, and phenotypic data (total N = 820). After meta-analysis, we found that mean adipocyte area positively correlated with body mass index (BMI) (Psubq = 8.13 × 10-69, ßsubq = 0.45; Pvisc = 2.5 × 10-55, ßvisc = 0.49; average R2 across cohorts = 0.49) and that adipocytes in subcutaneous depots are larger than their visceral counterparts (Pmeta = 9.8 × 10-7). Lastly, we performed the largest GWAS and subsequent meta-analysis of mean adipocyte area and intra-individual adipocyte variation (N = 820). Despite having twice the number of samples than any similar study, we found no genome-wide significant associations, suggesting that larger sample sizes and a homogenous collection of adipose tissue are likely needed to identify robust genetic associations.

Analysis of Differentiation Protocols Defines a Common Pancreatic Progenitor Molecular Signature and Guides Refinement of Endocrine Differentiation.

Several distinct differentiation protocols for deriving pancreatic progenitors (PPs) from human pluripotent stem cells have been described, but it remains to be shown how similar the PPs are across protocols and how well they resemble their in vivo counterparts. Here, we evaluated three differentiation protocols, performed RNA and assay for transposase-accessible chromatin using sequencing on isolated PPs derived with these, and compared them with fetal human pancreas populations. This enabled us to define a shared transcriptional and epigenomic signature of the PPs, including several genes not previously implicated in pancreas development. Furthermore, we identified a significant and previously unappreciated cross-protocol variation of the PPs through multi-omics analysis and demonstrate how such information can be applied to refine differentiation protocols for derivation of insulin-producing beta-like cells. Together, our study highlights the importance of a detailed characterization of defined cell populations derived from distinct differentiation protocols and provides a valuable resource for exploring human pancreatic development.

Loss of ZnT8 function protects against diabetes by enhanced insulin secretion.

A rare loss-of-function allele p.Arg138* in SLC30A8 encoding the zinc transporter 8 (ZnT8), which is enriched in Western Finland, protects against type 2 diabetes (T2D). We recruited relatives of the identified carriers and showed that protection was associated with better insulin secretion due to enhanced glucose responsiveness and proinsulin conversion, particularly when compared with individuals matched for the genotype of a common T2D-risk allele in SLC30A8, p.Arg325. In genome-edited human induced pluripotent stem cell (iPSC)-derived ß-like cells, we establish that the p.Arg138* allele results in reduced SLC30A8 expression due to haploinsufficiency. In human ß cells, loss of SLC30A8 leads to increased glucose responsiveness and reduced KATP channel function similar to isolated islets from carriers of the T2D-protective allele p.Trp325. These data position ZnT8 as an appealing target for treatment aimed at maintaining insulin secretion capacity in T2D.

Understanding human fetal pancreas development using subpopulation sorting, RNA sequencing and single-cell profiling.

To decipher the populations of cells present in the human fetal pancreas and their lineage relationships, we developed strategies to isolate pancreatic progenitors, endocrine progenitors and endocrine cells. Transcriptome analysis of the individual populations revealed a large degree of conservation among vertebrates in the drivers of gene expression changes that occur at different steps of differentiation, although notably, sometimes, different members of the same gene family are expressed. The transcriptome analysis establishes a resource to identify novel genes and pathways involved in human pancreas development. Single-cell profiling further captured intermediate stages of differentiation and enabled us to decipher the sequence of transcriptional events occurring during human endocrine differentiation. Furthermore, we evaluate how well individual pancreatic cells derived in vitro from human pluripotent stem cells mirror the natural process occurring in human fetuses. This comparison uncovers a few differences at the progenitor steps, a convergence at the steps of endocrine induction, and the current inability to fully resolve endocrine cell subtypes in vitro.

Type 2 diabetes risk alleles in PAM impact insulin release from human pancreatic ß-cells.

The molecular mechanisms underpinning susceptibility loci for type 2 diabetes (T2D) remain poorly understood. Coding variants in peptidylglycine α-amidating monooxygenase (PAM) are associated with both T2D risk and insulinogenic index. Here, we demonstrate that the T2D risk alleles impact negatively on overall PAM activity via defects in expression and catalytic function. PAM deficiency results in reduced insulin content and altered dynamics of insulin secretion in a human ß-cell model and primary islets from cadaveric donors. Thus, our results demonstrate a role for PAM in ß-cell function, and establish molecular mechanisms for T2D risk alleles at this locus.

Derivation and molecular characterization of pancreatic differentiated MODY1-iPSCs.

Maturity onset diabetes of the young (MODY) is a hereditary form of diabetes mellitus presenting at childhood or adolescence, which eventually leads to pancreatic ß-cells dysfunction. The underlying genetic basis of MODY disorders is haploinsufficiency, where loss-of-function mutations in a single allele cause the diabetic phenotype in heterozygous patients. MODY1 is a type of MODY disorder resulting from a mutation in the transcription factor hepatocyte nuclear factor 4 alpha (HNF4α). In order to establish a human based model to study MODY1, we generated patient-derived induced pluripotent stem cells (iPSCs). Differentiation of these pluripotent cells towards the pancreatic lineage enabled to evaluate the effects of the MODY1 mutation and its impact on endodermal and pancreatic cells. Analyzing the gene expression profiles of differentiated MODY1 cells, revealed the outcome of HNF4α haploinsufficiency on its targets. This molecular analysis suggests that the differential expression of HNF4α target genes in MODY1 is affected by the number of HNF4α binding sites, their distance from the transcription start site, and the number of other transcription factor binding sites. These features may help explain the molecular manifestations of haploinsufficiency in MODY1 disease.

NKX6.1 induced pluripotent stem cell reporter lines for isolation and analysis of functionally relevant neuronal and pancreas populations.

Recent studies have reported significant advances in the differentiation of human pluripotent stem cells to clinically relevant cell types such as the insulin producing beta-like cells and motor neurons. However, many of the current differentiation protocols lead to heterogeneous cell cultures containing cell types other than the targeted cell fate. Genetically modified human pluripotent stem cells reporting the expression of specific genes are of great value for differentiation protocol optimization and for the purification of relevant cell populations from heterogeneous cell cultures. Here we present the generation of human induced pluripotent stem cell (iPSC) lines with a GFP reporter inserted in the endogenous NKX6.1 locus. Characterization of the reporter lines demonstrated faithful GFP labelling of NKX6.1 expression during pancreas and motor neuron differentiation. Cell sorting and gene expression profiling by RNA sequencing revealed that NKX6.1-positive cells from pancreatic differentiations closely resemble human beta cells. Furthermore, functional characterization of the isolated cells demonstrated that glucose-stimulated insulin secretion is mainly confined to the NKX6.1-positive cells. We expect that the NKX6.1-GFP iPSC lines and the results presented here will contribute to the further refinement of differentiation protocols and characterization of hPSC-derived beta cells and motor neurons for disease modelling and cell replacement therapies.

Patterns of differential gene expression in a cellular model of human islet development, and relationship to type 2 diabetes predisposition.

AIMS/HYPOTHESIS: Most type 2 diabetes-associated genetic variants identified via genome-wide association studies (GWASs) appear to act via the pancreatic islet. Observed defects in insulin secretion could result from an impact of these variants on islet development and/or the function of mature islets. Most functional studies have focused on the latter, given limitations regarding access to human fetal islet tissue. Capitalising upon advances in in vitro differentiation, we characterised the transcriptomes of human induced pluripotent stem cell (iPSC) lines differentiated along the pancreatic endocrine lineage, and explored the contribution of altered islet development to the pathogenesis of type 2 diabetes. METHODS: We performed whole-transcriptome RNA sequencing of human iPSC lines from three independent donors, at baseline and at seven subsequent stages during in vitro islet differentiation. Differentially expressed genes (q < 0.01, log2 fold change [FC] > 1) were assigned to the stages at which they were most markedly upregulated. We used these data to characterise upstream transcription factors directing different stages of development, and to explore the relationship between RNA expression profiles and genes mapping to type 2 diabetes GWAS signals. RESULTS: We identified 9409 differentially expressed genes across all stages, including many known markers of islet development. Integration of differential expression data with information on transcription factor motifs highlighted the potential contribution of REST to islet development. Over 70% of genes mapping within type 2 diabetes-associated credible intervals showed peak differential expression during islet development, and type 2 diabetes GWAS loci of largest effect (including TCF7L2; log2FC = 1.2; q = 8.5 × 10-10) were notably enriched in genes differentially expressed at the posterior foregut stage (q = 0.002), as calculated by gene set enrichment analyses. In a complementary analysis of enrichment, genes differentially expressed in the final, beta-like cell stage of in vitro differentiation were significantly enriched (hypergeometric test, permuted p value <0.05) for genes within the credible intervals of type 2 diabetes GWAS loci. CONCLUSIONS/INTERPRETATION: The present study characterises RNA expression profiles during human islet differentiation, identifies potential transcriptional regulators of the differentiation process, and suggests that the inherited predisposition to type 2 diabetes is partly mediated through modulation of islet development. DATA AVAILABILITY: Sequence data for this study has been deposited at the European Genome-phenome Archive (EGA), under accession number EGAS00001002721.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls.

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

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

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

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

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

Genome-edited human stem cell-derived beta cells: a powerful tool for drilling down on type 2 diabetes GWAS biology.

Type 2 diabetes (T2D) is a disease of pandemic proportions, one defined by a complex aetiological mix of genetic, epigenetic, environmental, and lifestyle risk factors. Whilst the last decade of T2D genetic research has identified more than 100 loci showing strong statistical association with disease susceptibility, our inability to capitalise upon these signals reflects, in part, a lack of appropriate human cell models for study. This review discusses the impact of two complementary, state-of-the-art technologies on T2D genetic research: the generation of stem cell-derived, endocrine pancreas-lineage cells and the editing of their genomes. Such models facilitate investigation of diabetes-associated genomic perturbations in a physiologically representative cell context and allow the role of both developmental and adult islet dysfunction in T2D pathogenesis to be investigated. Accordingly, we interrogate the role that patient-derived induced pluripotent stem cell models are playing in understanding cellular dysfunction in monogenic diabetes, and how site-specific nucleases such as the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system are helping to confirm genes crucial to human endocrine pancreas development. We also highlight the novel biology gleaned in the absence of patient lines, including an ability to model the whole phenotypic spectrum of diabetes phenotypes occurring both in utero and in adult cells, interrogating the non-coding 'islet regulome' for disease-causing perturbations, and understanding the role of other islet cell types in aberrant glycaemia. This article aims to reinforce the importance of investigating T2D signals in cell models reflecting appropriate species, genomic context, developmental time point, and tissue type.

The genetic architecture of type 2 diabetes.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

Insights into islet development and biology through characterization of a human iPSC-derived endocrine pancreas model.

Directed differentiation of stem cells offers a scalable solution to the need for human cell models recapitulating islet biology and T2D pathogenesis. We profiled mRNA expression at 6 stages of an induced pluripotent stem cell (iPSC) model of endocrine pancreas development from 2 donors, and characterized the distinct transcriptomic profiles associated with each stage. Established regulators of endodermal lineage commitment, such as SOX17 (log2 fold change [FC] compared to iPSCs = 14.2, p-value = 4.9 × 10(-5)) and the pancreatic agenesis gene GATA6 (log2 FC = 12.1, p-value = 8.6 × 10(-5)), showed transcriptional variation consistent with their known developmental roles. However, these analyses highlighted many other genes with stage-specific expression patterns, some of which may be novel drivers or markers of islet development. For example, the leptin receptor gene, LEPR, was most highly expressed in published data from in vivo-matured cells compared to our endocrine pancreas-like cells (log2 FC = 5.5, p-value = 2.0 × 10(-12)), suggesting a role for the leptin pathway in the maturation process. Endocrine pancreas-like cells showed significant stage-selective expression of adult islet genes, including INS, ABCC8, and GLP1R, and enrichment of relevant GO-terms (e.g. "insulin secretion"; odds ratio = 4.2, p-value = 1.9 × 10(-3)): however, principal component analysis indicated that in vitro-differentiated cells were more immature than adult islets. Integration of the stage-specific expression information with genetic data from T2D genome-wide association studies revealed that 46 of 82 T2D-associated loci harbor genes present in at least one developmental stage, facilitating refinement of potential effector transcripts. Together, these data show that expression profiling in an iPSC islet development model can further understanding of islet biology and T2D pathogenesis.