Nikolay Oskolkov

TCF7L2 is a master regulator of insulin production and processing

Genome-wide association studies have revealed >60 loci associated with type 2 diabetes (T2D), but the underlying causal variants and functional mechanisms remain largely elusive. Although variants in TCF7L2 confer the strongest risk of T2D among common variants by presumed effects on islet function, the molecular mechanisms are not yet well understood. Using RNA-sequencing, we have identified a TCF7L2-regulated transcriptional network responsible for its effect on insulin secretion in rodent and human pancreatic islets. ISL1 is a primary target of TCF7L2 and regulates proinsulin production and processing via MAFA, PDX1, NKX6.1, PCSK1, PCSK2 and SLC30A8, thereby providing evidence for a coordinated regulation of insulin production and processing. The risk T-allele of rs7903146 was associated with increased TCF7L2 expression, and decreased insulin content and secretion. Using gene expression profiles of 66 human pancreatic islets donors’, we also show that the identified TCF7L2-ISL1 transcriptional network is regulated in a genotype-dependent manner. Taken together, these results demonstrate that not only synthesis of proinsulin is regulated by TCF7L2 but also processing and possibly clearance of proinsulin and insulin. These multiple targets in key pathways may explain why TCF7L2 has emerged as the gene showing one of the strongest associations with T2D.

A Syntenic Cross Species Aneuploidy Genetic Screen Links RCAN1 Expression to β-Cell Mitochondrial Dysfunction in Type 2 Diabetes

Type 2 diabetes (T2D) is a complex metabolic disease associated with obesity, insulin resistance and hypoinsulinemia due to pancreatic β-cell dysfunction. Reduced mitochondrial function is thought to be central to β-cell dysfunction. Mitochondrial dysfunction and reduced insulin secretion are also observed in β-cells of humans with the most common human genetic disorder, Down syndrome (DS, Trisomy 21). To identify regions of chromosome 21 that may be associated with perturbed glucose homeostasis we profiled the glycaemic status of different DS mouse models. The Ts65Dn and Dp16 DS mouse lines were hyperglycemic, while Tc1 and Ts1Rhr mice were not, providing us with a region of chromosome 21 containing genes that cause hyperglycemia. We then examined whether any of these genes were upregulated in a set of ~5,000 gene expression changes we had identified in a large gene expression analysis of human T2D β-cells. This approach produced a single gene, RCAN1, as a candidate gene linking hyperglycemia and functional changes in T2D β-cells. Further investigations demonstrated that RCAN1 methylation is reduced in human T2D islets at multiple sites, correlating with increased expression. RCAN1 protein expression was also increased in db/db mouse islets and in human and mouse islets exposed to high glucose. Mice overexpressing RCAN1 had reduced in vivo glucose-stimulated insulin secretion and their β-cells displayed mitochondrial dysfunction including hyperpolarised membrane potential, reduced oxidative phosphorylation and low ATP production. This lack of β-cell ATP had functional consequences by negatively affecting both glucose-stimulated membrane depolarisation and ATP-dependent insulin granule exocytosis. Thus, from amongst the myriad of gene expression changes occurring in T2D β-cells where we had little knowledge of which changes cause β-cell dysfunction, we applied a trisomy 21 screening approach which linked RCAN1 to β-cell mitochondrial dysfunction in T2D.

Genetic predisposition for beta cell fragility underlies type 1 and type 2 diabetes

Type 1 (T1D) and type 2 (T2D) diabetes share pathophysiological characteristics, yet mechanistic links have remained elusive. T1D results from autoimmune destruction of pancreatic beta cells, whereas beta cell failure in T2D is delayed and progressive. Here we find a new genetic component of diabetes susceptibility in T1D non-obese diabetic (NOD) mice, identifying immune-independent beta cell fragility. Genetic variation in Xrcc4 and Glis3 alters the response of NOD beta cells to unfolded protein stress, enhancing the apoptotic and senescent fates. The same transcriptional relationships were observed in human islets, demonstrating the role of beta cell fragility in genetic predisposition to diabetes.

The Genetic Landscape of Renal Complications in Type 1 Diabetes

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10-5) and the risk of type 2 diabetes (P=6.1×10-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10-6), and pentose and glucuronate interconversions (P=3.0×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.

Myocardin Family Members Drive Formation of Caveolae.

Caveolae are membrane organelles that play roles in glucose and lipid metabolism and in vascular function. Formation of caveolae requires caveolins and cavins. The make-up of caveolae and their density is considered to reflect cell-specific transcriptional control mechanisms for caveolins and cavins, but knowledge regarding regulation of caveolae genes is incomplete. Myocardin (MYOCD) and its relative MRTF-A (MKL1) are transcriptional coactivators that control genes which promote smooth muscle differentiation. MRTF-A communicates changes in actin polymerization to nuclear gene transcription. Here we tested if myocardin family proteins control biogenesis of caveolae via activation of caveolin and cavin transcription. Using human coronary artery smooth muscle cells we found that jasplakinolide and latrunculin B (LatB), substances that promote and inhibit actin polymerization, increased and decreased protein levels of caveolins and cavins, respectively. The effect of LatB was associated with reduced mRNA levels for these genes and this was replicated by the MRTF inhibitor CCG-1423 which was non-additive with LatB. Overexpression of myocardin and MRTF-A caused 5-10-fold induction of caveolins whereas cavin-1 and cavin-2 were induced 2-3-fold. PACSIN2 also increased, establishing positive regulation of caveolae genes from three families. Full regulation of CAV1 was retained in its proximal promoter. Knock down of the serum response factor (SRF), which mediates many of the effects of myocardin, decreased cavin-1 but increased caveolin-1 and -2 mRNAs. Viral transduction of myocardin increased the density of caveolae 5-fold in vitro. A decrease of CAV1 was observed concomitant with a decrease of the smooth muscle marker calponin in aortic aneurysms from mice (C57Bl/6) infused with angiotensin II. Human expression data disclosed correlations of MYOCD with CAV1 in a majority of human tissues and in the heart, correlation with MKL2 (MRTF-B) was observed. The myocardin family of transcriptional coactivators therefore drives formation of caveolae and this effect is largely independent of SRF.

RNA sequencing: current and prospective uses in metabolic research

Previous global RNA analysis was restricted to known transcripts in species with a defined transcriptome. Next generation sequencing has transformed transcriptomics by making it possible to analyse expressed genes with an exon level resolution from any tissue in any species without any a priori knowledge of which genes that are being expressed, splice patterns or their nucleotide sequence. In addition, RNA sequencing is a more sensitive technique compared with microarrays with a larger dynamic range, and it also allows for investigation of imprinting and allele-specific expression. This can be done for a cost that is able to compete with that of a microarray, making RNA sequencing a technique available to most researchers. Therefore RNA sequencing has recently become the state of the art with regards to large-scale RNA investigations and has to a large extent replaced microarrays. The only drawback is the large data amounts produced, which together with the complexity of the data can make a researcher spend far more time on analysis than performing the actual experiment.

Excess maternal transmission of variants in the THADA gene to offspring with type 2 diabetes

Aims/hypothesisGenome-wide association studies (GWAS) have identified more than 65 genetic loci associated with risk of type 2 diabetes. However, the contribution of distorted parental transmission of alleles to risk of type 2 diabetes has been mostly unexplored. Our goal was therefore to search for parent-of-origin effects (POE) among type 2 diabetes loci in families.MethodsFamilies from the Botnia study (n = 4,211, 1,083 families) were genotyped for 72 single-nucleotide polymorphisms (SNPs) associated with type 2 diabetes and assessed for POE on type 2 diabetes. The family-based Hungarian Transdanubian Biobank (HTB) (n = 1,463, >135 families) was used to replicate SNPs showing POE. Association of type 2 diabetes loci within families was also tested.ResultsThree loci showed nominal POE, including the previously reported variants in KCNQ1, for type 2 diabetes in families from Botnia (rs2237895: pPOE = 0.037), which can be considered positive controls. The strongest POE was seen for rs7578597 SNP in the THADA gene, showing excess transmission of the maternal risk allele T to diabetic offspring (Botnia: pPOE = 0.01; HTB pPOE = 0.045). These data are consistent with previous evidence of allelic imbalance for expression in islets, suggesting that the THADA gene can be imprinted in a POE-specific fashion. Five CpG sites, including those flanking rs7578597, showed differential methylation between diabetic and non-diabetic donor islets.Conclusions/interpretationTaken together, the data emphasise the need for genetic studies to consider from which parent an offspring has inherited a susceptibility allele.

LoFtool: a gene intolerance score based on loss-of-function variants in 60 706 individuals

Motivation: Depletion of loss‐of‐function (LoF) mutations may provide a rank of genic functional intolerance and consequently susceptibility to disease. Results: Here we have studied LoF mutations in 60 706 unrelated individuals and show that the most intolerant quartile of ranked genes is enriched in rare and early onset diseases and explains 87% of de novo haploinsufficient OMIM mutations, 17% more than any other gene scoring tool. We detected particular enrichment in expression of the depleted LoF genes in brain (odds ratio = 1.5; P‐value = 4.2e‐07). By searching for de novo haploinsufficient mutations putatively associated with neurodevelopmental disorders in four recent studies, we were able to explain 81% of them. Taken together, this study provides a novel gene intolerance ranking system, called LoFtool, which may help in ranking genes of interest based on their LoF intolerance and tissue expression. Availability and implementation: The LoFtool gene scores are available in the Supplementary data. Contact: joaofadista@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Glucose-induced Changes in Gene Expression in Human Pancreatic Islets - Causes or Consequences of Chronic Hyperglycemia

Dysregulation of gene expression in islets from patients with type 2 diabetes (T2D) might be causally involved in the development of hyperglycemia, or it could develop as a consequence of hyperglycemia (i.e., glucotoxicity). To separate the genes that could be causally involved in pathogenesis from those likely to be secondary to hyperglycemia, we exposed islets from human donors to normal or high glucose concentrations for 24 h and analyzed gene expression. We compared these findings with gene expression in islets from donors with normal glucose tolerance and hyperglycemia (including T2D). The genes whose expression changed in the same direction after short-term glucose exposure, as in T2D, were considered most likely to be a consequence of hyperglycemia. Genes whose expression changed in hyperglycemia but not after short-term glucose exposure, particularly those that also correlated with insulin secretion, were considered the strongest candidates for causal involvement in T2D. For example, ERO1LB, DOCK10, IGSF11, and PRR14L were downregulated in donors with hyperglycemia and correlated positively with insulin secretion, suggesting a protective role, whereas TMEM132C was upregulated in hyperglycemia and correlated negatively with insulin secretion, suggesting a potential pathogenic role. This study provides a catalog of gene expression changes in human pancreatic islets after exposure to glucose.

N 1 -methylnicotinamide is a signalling molecule produced in skeletal muscle coordinating energy metabolism

Obesity is a major health problem, and although caloric restriction and exercise are successful strategies to lose adipose tissue in obese individuals, a simultaneous decrease in skeletal muscle mass, negatively effects metabolism and muscle function. To deeper understand molecular events occurring in muscle during weight-loss, we measured the expressional change in human skeletal muscle following a combination of severe caloric restriction and exercise over 4 days in 15 Swedish men. Key metabolic genes were regulated after the intervention, indicating a shift from carbohydrate to fat metabolism. Nicotinamide N-methyltransferase (NNMT) was the most consistently upregulated gene following the energy-deficit exercise. Circulating levels of N1-methylnicotinamide (MNA), the product of NNMT activity, were doubled after the intervention. The fasting-fed state was an important determinant of plasma MNA levels, peaking at ~18 h of fasting and being lowest ~3 h after a meal. In culture, MNA was secreted by isolated human myotubes and stimulated lipolysis directly, with no effect on glucagon or insulin secretion. We propose that MNA is a novel myokine that enhances the utilization of energy stores in response to low muscle energy availability. Future research should focus on applying MNA as a biomarker to identify individuals with metabolic disturbances at an early stage.