Epigenome-wide association of DNA methylation markers in peripheral blood from Indian Asians and Europeans with incident type 2 diabetes: a nested case-control study

BACKGROUND

Indian Asians, who make up a quarter of the world's population, are at high risk of developing type 2 diabetes. We investigated whether DNA methylation is associated with future type 2 diabetes incidence in Indian Asians and whether differences in methylation patterns between Indian Asians and Europeans are associated with, and could be used to predict, differences in the magnitude of risk of developing type 2 diabetes.

METHODS

We did a nested case-control study of DNA methylation in Indian Asians and Europeans with incident type 2 diabetes who were identified from the 8-year follow-up of 25 372 participants in the London Life Sciences Prospective Population (LOLIPOP) study. Patients were recruited between May 1, 2002, and Sept 12, 2008. We did epigenome-wide association analysis using samples from Indian Asians with incident type 2 diabetes and age-matched and sex-matched Indian Asian controls, followed by replication testing of top-ranking signals in Europeans. For both discovery and replication, DNA methylation was measured in the baseline blood sample, which was collected before the onset of type 2 diabetes. Epigenome-wide significance was set at p<1 × 10(-7). We compared methylation levels between Indian Asian and European controls without type 2 diabetes at baseline to estimate the potential contribution of DNA methylation to increased risk of future type 2 diabetes incidence among Indian Asians.

FINDINGS

1608 (11·9%) of 13 535 Indian Asians and 306 (4·3%) of 7066 Europeans developed type 2 diabetes over a mean of 8·5 years (SD 1·8) of follow-up. The age-adjusted and sex-adjusted incidence of type 2 diabetes was 3·1 times (95% CI 2·8-3·6; p<0·0001) higher among Indian Asians than among Europeans, and remained 2·5 times (2·1-2·9; p<0·0001) higher after adjustment for adiposity, physical activity, family history of type 2 diabetes, and baseline glycaemic measures. The mean absolute difference in methylation level between type 2 diabetes cases and controls ranged from 0·5% (SD 0·1) to 1·1% (0·2). Methylation markers at five loci were associated with future type 2 diabetes incidence; the relative risk per 1% increase in methylation was 1·09 (95% CI 1·07-1·11; p=1·3 × 10(-17)) for ABCG1, 0·94 (0·92-0·95; p=4·2 × 10(-11)) for PHOSPHO1, 0·94 (0·92-0·96; p=1·4 × 10(-9)) for SOCS3, 1·07 (1·04-1·09; p=2·1 × 10(-10)) for SREBF1, and 0·92 (0·90-0·94; p=1·2 × 10(-17)) for TXNIP. A methylation score combining results for the five loci was associated with future type 2 diabetes incidence (relative risk quartile 4 vs quartile 1 3·51, 95% CI 2·79-4·42; p=1·3 × 10(-26)), and was independent of established risk factors. Methylation score was higher among Indian Asians than Europeans (p=1 × 10(-34)).

INTERPRETATION

DNA methylation might provide new insights into the pathways underlying type 2 diabetes and offer new opportunities for risk stratification and prevention of type 2 diabetes among Indian Asians.

FUNDING

The European Union, the UK National Institute for Health Research, the Wellcome Trust, the UK Medical Research Council, Action on Hearing Loss, the UK Biotechnology and Biological Sciences Research Council, the Oak Foundation, the Economic and Social Research Council, Helmholtz Zentrum Munchen, the German Research Center for Environmental Health, the German Federal Ministry of Education and Research, the German Center for Diabetes Research, the Munich Center for Health Sciences, the Ministry of Science and Research of the State of North Rhine-Westphalia, and the German Federal Ministry of Health.

Genetic influences on type 2 diabetes and metabolic syndrome related quantitative traits in Mauritius

Epidemiological studies report a high prevalence of type 2 diabetes and metabolic syndrome in the island nation of Mauritius. The Mauritius Family Study was initiated to examine heritable factors that contribute to these high rates of prevalence and consists of 400 individuals in 24 large extended multigenerational pedigrees. Anthropometric and biochemical measurements relating to the metabolic syndrome were undertaken in addition to family and lifestyle based information for each individual. Variance components methods were used to determine the heritability of the type 2 diabetes and metabolic syndrome related quantitative traits. The cohort was made up of 218 females (55%) and 182 males with 22% diagnosed with type 2 diabetes and a further 30% having impaired glucose tolerance or impaired fasting glucose. Notably BMI was not significantly increased in those with type 2 diabetes (P= .12), however a significant increase in waist circumference was observed in these groups (P= .02). The heritable proportion of trait variance was substantial and greater than values previously published for hip circumference, LDL and total cholesterol, diastolic and systolic blood pressure and serum creatinine. Height, weight and BMI heritabilities were all in the upper range of those previously reported. The phenotypic characteristics of the Mauritius family cohort are similar to those previously reported in the Mauritian population with a high observed prevalence rate of type 2 diabetes. A high heritability for key type 2 diabetes and metabolic syndrome related phenotypes (range 0.23 to 0.68), suggest the cohort will have utility in identifying genes that influence these quantitative traits.

Comparison of TAP2 frequencies in type 1 diabetes patients and healthy controls from three ethnic groups indicates an African origin for the TAP2 G allele

In order to determine the ethnic origin of the transporter associated with antigen processing 2 (TAP2) G allele, initially discovered by us in a group of type 1 diabetes (insulin-dependent diabetes mellitus) patients living on Reunion Island, HLA TAP2 typing was performed using the polymerase chain reaction-amplification refractory mutation system (PCR-ARMS) method in type 1 diabetes patients and unrelated healthy controls of three different ethnic groups (Caucasians, Indians and black Africans from Senegal and Mauritius). The comparison of TAP2 allele frequencies in controls showed significant racial (ethnic) differences. The TAP2*0101 and TAP2 C alleles were increased, respectively, in the Caucasian (50% in Caucasians vs. 40% in other groups) and Senegalese (27% in Senegalese vs. 10% in other groups) populations. In comparison with Caucasians, the TAP2*0201 variant was significantly increased in the Indian population and decreased in the Senegalese black population. In addition, the TAP2 G allele was observed in the two African populations studied but not in the Caucasian or Indian population. This observation is consistent with the view that this allele is restricted to populations of African origin. In addition, we have determined the large extended haplotype DQA1-DQB1-DRB1 associated with TAP2 G. We found that this allele is preferentially associated with the large conserved haplotype HLA DQA1*0501-DQB1*0201-DRB1*0301.

A multi-center study in order to further define the molecular basis of beta-thalassemia in Thailand, Pakistan, Sri Lanka, Mauritius, Syria, and India, and to develop a simple molecular diagnostic strategy by amplification refractory mutation system-polymerase chain reaction

The spectrum of the beta-thalassemia mutations of Thailand, Pakistan, India, Sri Lanka, Mauritius and Syria has been further characterized by a multi-center study of 1,235 transfusion-dependent patients, and the mutations discovered used to assess the fidelity of a simple diagnostic strategy. A total of 44 beta-thalassemia mutations were identified either by allele-specific oligonucleotide hybridization, amplification with allele-specific primers, or DNA sequencing of amplified product. The results confirm and extend earlier findings for Thailand, Pakistan, India, Mauritius and Syria. This is the first detailed report of the spectrum of mutations for Sri Lanka. Two novel mutations were identified, codon 55 (-A) and IVS-I-129 (A-->C), both found in Sri Lankan patients. Two beta-thalassemia mutations were found to coexist in one beta-globin gene: Sri Lankan patients homozygous for the beta0 codon 16 (-C) frameshift were also homozygous for the beta+ codon 10 (C-->A) mutation. Studies of Sri Lankan, Pakistani, and Indian carriers suggest the codon 10 (C-->A) mutation is just a rare polymorphism on an ancestral allele, on which the beta0 codon 16 (-C) mutation has arisen. Each country was found to have only a few common mutations accounting for 70% or more of the beta-thalassemia alleles. A panel of primers to diagnose the majority of the mutations by the amplification refractory mutation system was developed, enabling a simple molecular diagnostic strategy to be introduced for each country participating in the multi-center study.

Spectrum of beta thalassemia mutations and their linkage to beta-globin gene haplotypes in the Indo-Mauritians

The beta thalassemia alleles in 53 thalassemic Indo-Mauritian patients and their families consisting of 23 homozygous beta-thalassemia, 9 HbE/beta-thalassemia, 18 HbS/beta-thalassemia, 1 HbD/beta-thalassemia, 1 deltabeta/beta-thalassemia and 1 HbH/beta-thalassemia from the island of Mauritius were studied. Characterization by polymerase chain reaction-based reverse dot blot hybridization technique revealed that the IVS1-5 (G-->C) mutation accounted for 74% of the beta thalassemic alleles, while six other mutations occurred at much lower frequencies: HbE codon 26 (G-->A); 10.4%, codon 8/9 (+G); 3.5%, codon 30 (AGG-->ACG) also called IVSI (-1).G-->C; 3.5%, codon 15 (G-->A); 3.5%, codon 41/42 (-CTTT); 2.4% and -28 (A-->G); 2.4%. Association of these mutations to specific beta globin gene sequence framework and haplotype allowed to trace their ancestral link. These data are useful in future molecular screening of the population in view of implementing a thalassemia prevention and control program in Mauritius.

Lack of association of angiotensin I-converting enzyme gene polymorphism and premature myocardial infarction in Mauritian Indians

Eighty-five young Mauritian Indians, male survivors of premature myocardial infarction (MI) and thus belonging to a high risk group, were compared with 108 stringently selected controls for a possible association between premature MI and an insertion/deletion (I/D) polymorphism in the gene encoding angiotensin I-converting enzyme (ACE). The frequency of the D allele was 0.42 in the MI group and 0.43 in the control group, and thus no association between I/D polymorphism of ACE with susceptibility to early-onset MI was found in this population group. Other gene components of the renin-angiotensin system and lipid metabolism need to be explored to understand the genetic factors involved in causing MI at an early age.

Insertion/deletion polymorphism within a polyadenylate stretch at the human atrial natriuretic peptides (hANP) gene locus

A novel 8-bp bi-allelic insertion/deletion polymorphism is described within a polyadenylate stretch in the second intron of the human atrial natriuretic peptide gene locus. This new marker is located in the candidate gene for familial susceptibility to hypertension.

A new polymorphic restriction site at the human atrial natriuretic peptide (hANP) gene locus

A unique two allele polymorphism for both HpaII and SmaI is described in the second intron of the human atrial natriuretic peptide gene. It should be a useful marker of this candidate gene in familial susceptibility to hypertension.

Investigation of the polymorphic ScaI site by a PCR-based assay at the human atrial natriuretic peptides (hANP) gene locus

The ScaI polymorphic site within the stop codon of the human atrial natriuretic peptides (hANP) gene was investigated in Mauritian Indian, black African and French Caucasian populations. A distinct distribution pattern is observed in these three populations.