Epigenetic programming of adipose-derived stem cells in low birthweight individuals.

AIMS/HYPOTHESIS: Low birthweight (LBW) is associated with dysfunctions of adipose tissue and metabolic disease in adult life. We hypothesised that altered epigenetic and transcriptional regulation of adipose-derived stem cells (ADSCs) could play a role in programming adipose tissue dysfunction in LBW individuals. METHODS: ADSCs were isolated from the subcutaneous adipose tissue of 13 normal birthweight (NBW) and 13 LBW adult men. The adipocytes were cultured in vitro, and genome-wide differences in RNA expression and DNA methylation profiles were analysed in ADSCs and differentiated adipocytes. RESULTS: We demonstrated that ADSCs from LBW individuals exhibit multiple expression changes as well as genome-wide alterations in methylation pattern. Reduced expression of the transcription factor cyclin T2 encoded by CCNT2 may play a key role in orchestrating several of the gene expression changes in ADSCs from LBW individuals. Indeed, silencing of CCNT2 in human adipocytes decreased leptin secretion as well as the mRNA expression of several genes involved in adipogenesis, including MGLL, LIPE, PPARG, LEP and ADIPOQ. Only subtle genome-wide mRNA expression and DNA methylation changes were seen in mature cultured adipocytes from LBW individuals. CONCLUSIONS/INTERPRETATION: Epigenetic and transcriptional changes in LBW individuals are most pronounced in immature ADSCs that in turn may programme physiological characteristics of the mature adipocytes that influence the risk of metabolic diseases. Reduced expression of CCNT2 may play a key role in the developmental programming of adipose tissue.

A Genome-Wide mQTL Analysis in Human Adipose Tissue Identifies Genetic Variants Associated with DNA Methylation, Gene Expression and Metabolic Traits.

Little is known about the extent to which interactions between genetics and epigenetics may affect the risk of complex metabolic diseases and/or their intermediary phenotypes. We performed a genome-wide DNA methylation quantitative trait locus (mQTL) analysis in human adipose tissue of 119 men, where 592,794 single nucleotide polymorphisms (SNPs) were related to DNA methylation of 477,891 CpG sites, covering 99% of RefSeq genes. SNPs in significant mQTLs were further related to gene expression in adipose tissue and obesity related traits. We found 101,911 SNP-CpG pairs (mQTLs) in cis and 5,342 SNP-CpG pairs in trans showing significant associations between genotype and DNA methylation in adipose tissue after correction for multiple testing, where cis is defined as distance less than 500 kb between a SNP and CpG site. These mQTLs include reported obesity, lipid and type 2 diabetes loci, e.g. ADCY3/POMC, APOA5, CETP, FADS2, GCKR, SORT1 and LEPR. Significant mQTLs were overrepresented in intergenic regions meanwhile underrepresented in promoter regions and CpG islands. We further identified 635 SNPs in significant cis-mQTLs associated with expression of 86 genes in adipose tissue including CHRNA5, G6PC2, GPX7, RPL27A, THNSL2 and ZFP57. SNPs in significant mQTLs were also associated with body mass index (BMI), lipid traits and glucose and insulin levels in our study cohort and public available consortia data. Importantly, the Causal Inference Test (CIT) demonstrates how genetic variants mediate their effects on metabolic traits (e.g. BMI, cholesterol, high-density lipoprotein (HDL), hemoglobin A1c (HbA1c) and homeostatic model assessment of insulin resistance (HOMA-IR)) via altered DNA methylation in human adipose tissue. This study identifies genome-wide interactions between genetic and epigenetic variation in both cis and trans positions influencing gene expression in adipose tissue and in vivo (dys)metabolic traits associated with the development of obesity and diabetes.

Blood-based biomarkers of age-associated epigenetic changes in human islets associate with insulin secretion and diabetes.

Aging associates with impaired pancreatic islet function and increased type 2 diabetes (T2D) risk. Here we examine whether age-related epigenetic changes affect human islet function and if blood-based epigenetic biomarkers reflect these changes and associate with future T2D. We analyse DNA methylation genome-wide in islets from 87 non-diabetic donors, aged 26-74 years. Aging associates with increased DNA methylation of 241 sites. These sites cover loci previously associated with T2D, for example, KLF14. Blood-based epigenetic biomarkers reflect age-related methylation changes in 83 genes identified in human islets (for example, KLF14, FHL2, ZNF518B and FAM123C) and some associate with insulin secretion and T2D. DNA methylation correlates with islet expression of multiple genes, including FHL2, ZNF518B, GNPNAT1 and HLTF. Silencing these genes in ß-cells alter insulin secretion. Together, we demonstrate that blood-based epigenetic biomarkers reflect age-related DNA methylation changes in human islets, and associate with insulin secretion in vivo and T2D.

Genome-wide associations between genetic and epigenetic variation influence mRNA expression and insulin secretion in human pancreatic islets.

Genetic and epigenetic mechanisms may interact and together affect biological processes and disease development. However, most previous studies have investigated genetic and epigenetic mechanisms independently, and studies examining their interactions throughout the human genome are lacking. To identify genetic loci that interact with the epigenome, we performed the first genome-wide DNA methylation quantitative trait locus (mQTL) analysis in human pancreatic islets. We related 574,553 single nucleotide polymorphisms (SNPs) with genome-wide DNA methylation data of 468,787 CpG sites targeting 99% of RefSeq genes in islets from 89 donors. We identified 67,438 SNP-CpG pairs in cis, corresponding to 36,783 SNPs (6.4% of tested SNPs) and 11,735 CpG sites (2.5% of tested CpGs), and 2,562 significant SNP-CpG pairs in trans, corresponding to 1,465 SNPs (0.3% of tested SNPs) and 383 CpG sites (0.08% of tested CpGs), showing significant associations after correction for multiple testing. These include reported diabetes loci, e.g. ADCY5, KCNJ11, HLA-DQA1, INS, PDX1 and GRB10. CpGs of significant cis-mQTLs were overrepresented in the gene body and outside of CpG islands. Follow-up analyses further identified mQTLs associated with gene expression and insulin secretion in human islets. Causal inference test (CIT) identified SNP-CpG pairs where DNA methylation in human islets is the potential mediator of the genetic association with gene expression or insulin secretion. Functional analyses further demonstrated that identified candidate genes (GPX7, GSTT1 and SNX19) directly affect key biological processes such as proliferation and apoptosis in pancreatic ß-cells. Finally, we found direct correlations between DNA methylation of 22,773 (4.9%) CpGs with mRNA expression of 4,876 genes, where 90% of the correlations were negative when CpGs were located in the region surrounding transcription start site. Our study demonstrates for the first time how genome-wide genetic and epigenetic variation interacts to influence gene expression, islet function and potential diabetes risk in humans.

Genome-wide DNA methylation analysis of human pancreatic islets from type 2 diabetic and non-diabetic donors identifies candidate genes that influence insulin secretion.

Impaired insulin secretion is a hallmark of type 2 diabetes (T2D). Epigenetics may affect disease susceptibility. To describe the human methylome in pancreatic islets and determine the epigenetic basis of T2D, we analyzed DNA methylation of 479,927 CpG sites and the transcriptome in pancreatic islets from T2D and non-diabetic donors. We provide a detailed map of the global DNA methylation pattern in human islets, ß- and α-cells. Genomic regions close to the transcription start site showed low degrees of methylation and regions further away from the transcription start site such as the gene body, 3'UTR and intergenic regions showed a higher degree of methylation. While CpG islands were hypomethylated, the surrounding 2 kb shores showed an intermediate degree of methylation, whereas regions further away (shelves and open sea) were hypermethylated in human islets, ß- and α-cells. We identified 1,649 CpG sites and 853 genes, including TCF7L2, FTO and KCNQ1, with differential DNA methylation in T2D islets after correction for multiple testing. The majority of the differentially methylated CpG sites had an intermediate degree of methylation and were underrepresented in CpG islands (∼ 7%) and overrepresented in the open sea (∼ 60%). 102 of the differentially methylated genes, including CDKN1A, PDE7B, SEPT9 and EXOC3L2, were differentially expressed in T2D islets. Methylation of CDKN1A and PDE7B promoters in vitro suppressed their transcriptional activity. Functional analyses demonstrated that identified candidate genes affect pancreatic ß- and α-cells as Exoc3l silencing reduced exocytosis and overexpression of Cdkn1a, Pde7b and Sept9 perturbed insulin and glucagon secretion in clonal ß- and α-cells, respectively. Together, our data can serve as a reference methylome in human islets. We provide new target genes with altered DNA methylation and expression in human T2D islets that contribute to perturbed insulin and glucagon secretion. These results highlight the importance of epigenetics in the pathogenesis of T2D.

A six months exercise intervention influences the genome-wide DNA methylation pattern in human adipose tissue.

Epigenetic mechanisms are implicated in gene regulation and the development of different diseases. The epigenome differs between cell types and has until now only been characterized for a few human tissues. Environmental factors potentially alter the epigenome. Here we describe the genome-wide pattern of DNA methylation in human adipose tissue from 23 healthy men, with a previous low level of physical activity, before and after a six months exercise intervention. We also investigate the differences in adipose tissue DNA methylation between 31 individuals with or without a family history of type 2 diabetes. DNA methylation was analyzed using Infinium HumanMethylation450 BeadChip, an array containing 485,577 probes covering 99% RefSeq genes. Global DNA methylation changed and 17,975 individual CpG sites in 7,663 unique genes showed altered levels of DNA methylation after the exercise intervention (q<0.05). Differential mRNA expression was present in 1/3 of gene regions with altered DNA methylation, including RALBP1, HDAC4 and NCOR2 (q<0.05). Using a luciferase assay, we could show that increased DNA methylation in vitro of the RALBP1 promoter suppressed the transcriptional activity (p = 0.03). Moreover, 18 obesity and 21 type 2 diabetes candidate genes had CpG sites with differences in adipose tissue DNA methylation in response to exercise (q<0.05), including TCF7L2 (6 CpG sites) and KCNQ1 (10 CpG sites). A simultaneous change in mRNA expression was seen for 6 of those genes. To understand if genes that exhibit differential DNA methylation and mRNA expression in human adipose tissue in vivo affect adipocyte metabolism, we silenced Hdac4 and Ncor2 respectively in 3T3-L1 adipocytes, which resulted in increased lipogenesis both in the basal and insulin stimulated state. In conclusion, exercise induces genome-wide changes in DNA methylation in human adipose tissue, potentially affecting adipocyte metabolism.

Reduced insulin secretion correlates with decreased expression of exocytotic genes in pancreatic islets from patients with type 2 diabetes.

Reduced insulin release has been linked to defect exocytosis in ß-cells. However, whether expression of genes suggested to be involved in the exocytotic process (exocytotic genes) is altered in pancreatic islets from patients with type 2 diabetes (T2D), and correlate to insulin secretion, needs to be further investigated. Analysing expression levels of 23 exocytotic genes using microarray revealed reduced expression of five genes in human T2D islets (χ(2)=13.25; p<0.001). Gene expression of STX1A, SYT4, SYT7, SYT11, SYT13, SNAP25 and STXBP1 correlated negatively to in vivo measurements of HbA1c levels and positively to glucose stimulated insulin secretion (GSIS) in vitro in human islets. STX1A, SYT4 and SYT11 protein levels correspondingly decreased in human T2D islets. Moreover, silencing of SYT4 and SYT13 reduced GSIS in INS1-832/13 cells. Our data support that reduced expression of exocytotic genes contributes to impaired insulin secretion, and suggest decreased expression of these genes as part of T2D pathogenesis.

Regulation of core clock genes in human islets.

Nearly all mammalian cells express a set of genes known as clock genes. These regulate the circadian rhythm of cellular processes by means of negative and positive autoregulatory feedback loops of transcription and translation. Recent genomewide association studies have demonstrated an association between a polymorphism near the circadian clock gene CRY2 and elevated fasting glucose. To determine whether clock genes could play a pathogenetic role in the disease, we examined messenger RNA (mRNA) expression of core clock genes in human islets from donors with or without type 2 diabetes mellitus. Microarray and quantitative real-time polymerase chain reaction analyses were used to assess expression of the core clock genes CLOCK, BMAL-1, PER1 to 3, and CRY1 and 2 in human islets. Insulin secretion and insulin content in human islets were measured by radioimmunoassay. The mRNA levels of PER2, PER3, and CRY2 were significantly lower in islets from donors with type 2 diabetes mellitus. To investigate the functional relevance of these clock genes, we correlated their expression to insulin content and glycated hemoglobin levels: mRNA levels of PER2 (ρ = 0.33, P = .012), PER3 (ρ = 0.30, P = .023), and CRY2 (ρ = 0.37, P = .0047) correlated positively with insulin content. Of these genes, expression of PER3 and CRY2 correlated negatively with glycated hemoglobin levels (ρ = -0.44, P = .0012; ρ = -0.28, P = .042). Furthermore, in an in vitro model mimicking pathogenetic conditions, the PER3 mRNA level was reduced in human islets exposed to 16.7 mmol/L glucose per 1 mmol/L palmitate for 48 hours (P = .003). Core clock genes are regulated in human islets. The data suggest that perturbations of circadian clock components may contribute to islet pathophysiology in human type 2 diabetes mellitus.

Decreased expression of genes involved in oxidative phosphorylation in human pancreatic islets from patients with type 2 diabetes.

OBJECTIVE: Gene expression alterations, especially in target tissues of insulin, have been associated with type 2 diabetes (T2D). In this study, we examined if genes involved in oxidative phosphorylation (OXPHOS) show differential gene expression and DNA methylation in pancreatic islets from patients with T2D compared with non-diabetic donors. DESIGN AND METHODS: Gene expression was analyzed in human pancreatic islets from 55 non-diabetic donors and nine T2D donors using microarray. RESULTS: While the expected number of OXPHOS genes with reduced gene expression is 7.21, we identified 21 downregulated OXPHOS genes in pancreatic islets from patients with T2D using microarray analysis. This gives a ratio of observed over expected OXPHOS genes of 26.37 by a χ(2)-test with P=2.81 × 10(-7). The microarray data was validated by qRT-PCR for four selected OXPHOS genes: NDUFA5, NDUFA10, COX11, and ATP6V1H. All four OXPHOS genes were significantly downregulated in islets from patients with T2D compared with non-diabetic donors using qRT-PCR (P ≤ 0.01). Furthermore, HbAlc levels correlated negatively with gene expression of NDUFA5, COX11, and ATP6V1H (P<0.05). Gene expression of NDUFA5, NDUFA10, COX11, and ATP6V1H correlated positively with glucose-stimulated insulin secretion (P<0.03). Finally, DNA methylation was analyzed upstream of the transcription start for NDUFA5, COX11, and ATP6V1H. However, none of the analyzed CpG sites in the three genes showed differences in DNA methylation in islets from donors with T2D compared with non-diabetic donors. CONCLUSION: Pancreatic islets from patients with T2D show decreased expression of a set of OXPHOS genes, which may lead to impaired insulin secretion.

The expression of myosin heavy chain (MHC) genes in human skeletal muscle is related to metabolic characteristics involved in the pathogenesis of type 2 diabetes.

Type 2 diabetes patients exhibit a reduction in oxidative muscle fibres and an increase in glycolytic muscle fibres. In this study, we investigated whether both genetic and non-genetic factors influence the mRNA expression levels of three myosin heavy chain (MHC) genes represented in different fibre types. Specifically, we examined the MHC7 (slow-twitch oxidative fibre), MHCIIa (fast-twitch oxidative fibre) and MHCIIx/d (fast-twitch glycolytic fibre) genes in human skeletal muscle. We further investigated the use of MHC mRNA expression as a proxy to determine fibre-type composition, as measured by traditional ATP staining. Two cohorts of age-matched Swedish men were studied to determine the relationship of muscle mRNA expression of MHC7, MHCIIa, and MHCIIx/d with muscle fibre composition. A classical twin approach, including young and elderly Danish twin pairs, was utilised to examine if differences in expression levels were due to genetic or environmental factors. Although MHCIIx/d mRNA expression correlated positively with the level of type IIx/d muscle fibres in the two cohorts (P<0.05), a relatively low magnitude of correlation suggests that mRNA does not fully correlate with fibre-type composition. Heritability estimates and genetic analysis suggest that the levels of MHC7, MHCIIa and MHCIIx/d expression are primarily under non-genetic influence, and MHCIIa indicated an age-related decline. PGC-1α exhibited a positive relationship with the expression of all three MHC genes (P<0.05); meanwhile, PGC-1ß related positively with MHCIIa expression and negatively with MHCIIx/d expression (P<0.05). While MHCIIa expression related positively with insulin-stimulated glucose uptake (P<0.01), MHCIIx/d expression related negatively with insulin-stimulated glucose uptake (P<0.05). Our findings suggest that the expression levels of the MHC genes are associated with age and both PGC-1α and PGC-1ß and indicate that the MHC genes may to some extent be used to determine fibre-type composition in human skeletal muscle.

Two common genetic variants near nuclear-encoded OXPHOS genes are associated with insulin secretion in vivo.

CONTEXT: Mitochondrial ATP production is important in the regulation of glucose-stimulated insulin secretion. Genetic factors may modulate the capacity of the ß-cells to secrete insulin and thereby contribute to the risk of type 2 diabetes. OBJECTIVE: The aim of this study was to identify genetic loci in or adjacent to nuclear-encoded genes of the oxidative phosphorylation (OXPHOS) pathway that are associated with insulin secretion in vivo. DESIGN AND METHODS: To find polymorphisms associated with glucose-stimulated insulin secretion, data from a genome-wide association study (GWAS) of 1467 non-diabetic individuals, including the Diabetes Genetic Initiative (DGI), was examined. A total of 413 single nucleotide polymorphisms with a minor allele frequency ≥0.05 located in or adjacent to 76 OXPHOS genes were included in the DGI GWAS. A more extensive population-based study of 4323 non-diabetics, the PPP-Botnia, was used as a replication cohort. Insulinogenic index during an oral glucose tolerance test was used as a surrogate marker of glucose-stimulated insulin secretion. Multivariate linear regression analyses were used to test genotype-phenotype associations. RESULTS: Two common variants were identified in the DGI, where the major C-allele of rs606164, adjacent to NADH dehydrogenase (ubiquinone) 1 subunit C2 (NDUFC2), and the minor G-allele of rs1323070, adjacent to cytochrome c oxidase subunit VIIa polypeptide 2 (COX7A2), showed nominal associations with decreased glucose-stimulated insulin secretion (P=0.0009, respective P=0.003). These associations were replicated in PPP-Botnia (P=0.002 and P=0.05). CONCLUSION: Our study shows that genetic variation near genes involved in OXPHOS may influence glucose-stimulated insulin secretion in vivo.

A common variant in TFB1M is associated with reduced insulin secretion and increased future risk of type 2 diabetes.

Type 2 diabetes (T2D) evolves when insulin secretion fails. Insulin release from the pancreatic ß cell is controlled by mitochondrial metabolism, which translates fluctuations in blood glucose into metabolic coupling signals. We identified a common variant (rs950994) in the human transcription factor B1 mitochondrial (TFB1M) gene associated with reduced insulin secretion, elevated postprandial glucose levels, and future risk of T2D. Because islet TFB1M mRNA levels were lower in carriers of the risk allele and correlated with insulin secretion, we examined mice heterozygous for Tfb1m deficiency. These mice displayed lower expression of TFB1M in islets and impaired mitochondrial function and released less insulin in response to glucose in vivo and in vitro. Reducing TFB1M mRNA and protein in clonal ß cells by RNA interference impaired complexes of the mitochondrial oxidative phosphorylation system. Consequently, nutrient-stimulated ATP generation was reduced, leading to perturbed insulin secretion. We conclude that a deficiency in TFB1M and impaired mitochondrial function contribute to the pathogenesis of T2D.

A common variant near the PRL gene is associated with increased adiposity in males.

A common variant (rs4712652) adjacent to the prolactin gene was recently associated with obesity using a genome-wide association study. The aim of this study was to replicate the association between rs4712652 and obesity and further examine if rs4712652 is associated with fat percentage and adiponectin levels in a population based Scandinavian cohort. rs4712652 was genotyped in 4879 participants (mean BMI 26.5±4.5 kg/m(2)) from the population-based PPP-Botnia Study and related to BMI, fat percentage and adiponectin levels. We found that the risk A allele of rs4712652 is associated with increased BMI and fat percentage in males (P=0.0047 and P=0.025, respectively), but not in females (P=0.98, P=0.45). Male A allele carriers have a higher risk of being overweight with an OR of 1.16 (P=0.025). While there was a significant negative correlation between adiponectin levels and fat percentage (r=-0.36; P=0.039) in male carriers of the protective GG genotype, this correlation was lost in male carriers of the risk rs4712652 A allele (P=0.33). Thus, the common SNP rs4712652 near the PRL gene seems to affect body fat and adiposity in a sex-specific fashion. It remains to be shown whether this is mediated by different prolactin concentrations or differences in tissue sensitivity to prolactin.