Research Paper Volume 5, Issue 5 pp 373—385

Centenarians as super-controls to assess the biological relevance of genetic risk factors for common age-related diseases: A proof of principle on type 2 diabetes

Paolo Garagnani^1,2,3, , Cristina Giuliani^4, , Chiara Pirazzini^1,2, , Fabiola Olivieri^5,6, , Maria Giulia Bacalini^1,2, , Rita Ostan^1,2, , Daniela Mari^7, , Giuseppe Passarino^8, , Daniela Monti^9, , Anna Rita Bonfigli^10, , Massimo Boemi^10, , Antonio Ceriello^11,12, , Stefano Genovese^13, , Federica Sevini^1,2, , Donata Luiselli^4, , Paolo Tieri^14, , Miriam Capri^1,2, , Stefano Salvioli^1,2, , Jan Vijg^15,17, , Yousin Suh^15,16,17,18, , Massimo Delledonne^19,20, , Roberto Testa^21, , Claudio Franceschi^1, ,

¹ DIMES - Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, Bologna, 40126 Italy
² C.I.G. Interdepartmental Center “L. Galvani”, University of Bologna, Bologna, Italy
³ CRBA - Applied Biomedical Research Center, S. Orsola-Malpighi Polyclinic, Bologna, 40138 Italy
⁴ Department of Biological, Geological and Environmental Sciences, Laboratory of Molecular Anthropology & Centre for Genome Biology, University of Bologna, Bologna 40126, Italy
⁵ Department of Clinical and Molecular Sciences, Università Politecnica delle Marche, Ancona, Italy
⁶ Center of Clinical Pathology and Innovative Therapy, Italian National Research Center on Aging INRCA-IRCCS, Ancona, Italy
⁷ Geriatric Unit IRCCS Ca' Granda Foundation Maggiore Policlinico Hospital Hospital and Department of Clinical Sciences and Community Health, University of Milano, Italy
⁸ Department of Cell Biology, University of Calabria, Rende, Italy
⁹ Department of Experimental and Clinical Biomedical Sciences, University of Florence, Florence, Italy
¹⁰ Metabolic Diseases and Diabetology Unit, IRCCS-INRCA, Ancona, Italy
¹¹ Institut d'Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain
¹² Centro de Investigación Biomédica en Red de Diabetes y Enfermedades Metabólicas Asociadas (CIBERDEM), Barcelona, Spain
¹³ Department of Cardiovascular and Metabolic Diseases, IRCCS Gruppo Multimedica Sesto San Giovanni (MI), Italy
¹⁴ IAC-CNR Istituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, Rome, Italy
¹⁵ Department of Genetics, Albert Einstein College of Medicine, Bronx, NY 10461, USA
¹⁶ Department of Medicine, Albert Einstein College of Medicine, Bronx, NY 10461, USA
¹⁷ Institute for Aging Research, Diabetes Research and Training Center, Albert Einstein College of Medicine, Bronx, NY 10461, USA
¹⁸ Institute of Aging Research, Guangdong Medical College, Dongguan 523808, China
¹⁹ Personal Genomics SRL, Strada le Grazie 15, 37133 Verona - Italy
²⁰ Functional Genomics Center, Dept. of Biotechnologies, University of Verona, Strada le Grazie 15, 37133 Verona - Italy
²¹ Experimental Models in Clinical Pathology, IRCCS-INRCA, Ancona, Italy

* Equal contribution

Received: May 17, 2013 Accepted: May 31, 2013 Published: May 31, 2013

Copyright: © 2013 Garagnani et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Genetic association studies of age-related, chronic human diseases often suffer from a lack of power to detect modest effects. Here we propose an alternative approach of including healthy centenarians as a more homogeneous and extreme control group. As a proof of principle we focused on type 2 diabetes (T2D) and assessed allelic/genotypic associations of 31 SNPs associated with T2D, diabetes complications and metabolic diseases and SNPs of genes relevant for telomere stability and age-related diseases. We hypothesized that the frequencies of risk variants are inversely correlated with decreasing health and longevity. We performed association analyses comparing diabetic patients and non-diabetic controls followed by association analyses with extreme phenotypic groups (T2D patients with complications and centenarians). Results drew attention to rs7903146 (TCF7L2 gene) that showed a constant increase in the frequencies of risk genotype (TT) from centenarians to diabetic patients who developed macro-complications and the strongest genotypic association was detected when diabetic patients were compared to centenarians (p_value = 9.066*10⁻⁷). We conclude that robust and biologically relevant associations can be obtained when extreme phenotypes, even with a small sample size, are compared.

Introduction

In the last ten years the scientific community has devoted a consistent effort to identify the genetic basis of the most common age-related diseases, as they represent one of the most important public health and socio-economical burden all over the world and particularly in Western Countries. This challenge was mainly faced up by genome wide association studies (GWASs) based on microarray technology that allows the simultaneous analyses of hundred thousands of single nucleotide polymorphisms (SNPs), within the framework of the “common variant common disease” theory [1]. So far, more than 1,000 published GWASs reported significant associations of ~4,000 SNPs for more than 200 traits/diseases [2]. GWASs of age-related, chronic human diseases often suffer from a lack of power to detect modest effects, which can to some extent explain why the identified genetic effects comprise only a small fraction of the estimated trait heritability. These limitations can be overcome simply by ever increasing sample size in order to achieve the necessary statistical power to detect variants with small effects, which is not always feasible. Moreover, in most cases the biological role/relevance of the genetic variants emerged from GWASs is still unclear. In past, several strategies were proposed to increase detection power without affecting the cohort size, such as the use of genetically isolated populations characterized by a reduced genetic diversity that facilitate the discovery of relevant loci for complex diseases [3,4]. To date this approach has not been successful to boost the discovery of relevant associations with complex traits [5]. Thus, there is an urgent need for an effective strategy to obtain biological insights from the genetic knowledge derived from GWASs that can be translated into clinical benefits.

The aim of this study was to test the hypothesis that GWASs sensitivity can be boosted by including extreme phenotypic groups. As a proof of principle study, we considered one of the major age-related diseases, i.e. type 2 diabetes (T2D). We applied a candidate gene approach to assess associations of a limited set of SNPs (31) in or nearby genes relevant for T2D, diabetes complications, metabolic diseases, telomere stability and age-related diseases. The extreme phenotypic groups we focused on are as follows: i) CENTENARIANS who reached the extreme limit of human life, escaping or largely postponing the major age-related diseases, including T2D, and who can be considered a paradigm of healthy aging [6]. In particular, a major characteristic of centenarians is their extraordinarily well preserved glucose metabolism. Aging is frequently associated with impaired glucose metabolism related to a raise in insulin resistance (IR), not fully compensated by a sufficient β-cell function. Such age-related metabolic changes are key risk factors for T2D and are associated with a variety of intermediate phenotypes (hypertension, atherosclerosis, obesity) strongly affecting morbidity, disability and mortality among elderly. In a pioneering study, Paolisso et al. [7], using an oral glucose tolerance test and euglycemic glucose clamp, showed that centenarians had a 2-h plasma glucose concentration that was lower than that of aged subjects but not different from adults, and an insulin-mediated glucose uptake higher than that reported for aged subjects but not different from that found in adults. In a subsequent study, involving a large cohort of individuals (age range: 28-111 years) carefully selected for health status, we showed that the age-related trajectories of IR and β-cell function increase with age (the increase of β-cell function is necessary to compensate for the raise in IR). However, beyond 85-90 years of age, people with a lower degree of IR and with a lower β-cell function emerge, indicating that in very old people lower IR does not require a compensatory increment in β-cell function, thus allowing to preserve endocrine pancreas secretion and to prevent the development of T2D [8]. On the whole, these results point towards IR and β-cell function as phenotypes under strong selective forces during aging, supporting the hypothesis that effective peripheral glucose disposal is pivotal in determining longevity. These observations are confirmed by studies from our and other groups showing that centenarians offspring and nonagenarian siblings have a better health status than subjects of the same cohort, without the parental extreme longevity [9], and a reduced risk of diabetes and an enhanced insulin sensitivity [10-14]. ii) PATIENTS WITH T2D PLUS ONE OR MORE COMPLICATIONS. The T2D patients included in the present study have been fully characterized from a phenotypical (clinical, biochemical, pharmacological, among others) point of view, and have been extensively studied regarding a variety of genetic [15-18] and non-genetic factors. In particular, we showed that T2D patients undergo accelerated aging [19] and that IL-10 polymorphisms [20], telomere length in peripheral blood cells [21] and mitochondrial DNA variants [22] are able to distinguish between patients with and without complications. Thus, the presence of one or more micro- and macro-vascular complications not only represents the most severe and extreme phenotype of T2D, of relevance from a clinical and therapeutic point of view, but also has a biological counterpart which is still largely unknown. In this study, T2D patients have been compared with age-, gender- and geography-matched control group.

Here we considered a limited set of SNPs belonging to four classes of genetic variants: i) SNPs previously found associated with T2D by GWASs; ii) variants of genes involved in vascular pathology potentially involved in the development of T2D complications; iii) SNPs previously associated with metabolic diseases, including T2D; and iv) SNPs of genes relevant for telomere stability and involved in age-related diseases, including T2D complications.

Results

In this study we analyzed the following cohorts described in depth in Methods: the whole cohort of diabetic patients (D), diabetic patients with complications (D+Co), diabetic patients who developed micro-vascular complications (D+microCo), diabetic patients who developed macro-vascular complications (D+macroCo), non diabetic controls matched for age, gender and geographical origin with diabetic patients (CTR) and centenarians (100+). These samples were tested for 31 SNPs, 22 of which passed the quality check. Results from association analyses are reported in Supplemental Materials (S1, S2, S3, S4). The most significant results involve TCF7L2 and ADIPOQ gene variants, which are reported in Table 1 and summarized in Figure 1. The comparisons that have been performed to test the effectiveness of our approach are reported in the following five sections.

Table 1.

Allelic and genotypic association analyses. Samples considered were diabetic patients (N=562), non diabetic controls (N=558), diabetic individuals with complications (N=241), centenarians (N=229), diabetic patients with macro-vascular complications (N=64) and diabetic patients with micro-vascular complications (N=110).

	Cohorts	Gene	SNP (ID)	Allele/Genotype	Frequency (%)		p_value	OR (95% CI)	Model
	Cohorts	Gene	SNP (ID)	Allele/Genotype	1	2	p_value	OR (95% CI)	Model
Allelic Association	D (1) vs 100+ (2)	ADIPOQ	rs266729	G	24.6	31.6	5.89*10⁻³	0.71 (0.55-0.91)	/
	D (1) vs CTR (2)	TCF7L2	rs7903146	T	45	37	3.1*10⁻⁴	1.37 (1.16-1.63)	/
	D+Co (1) vs CTR (2)	TCF7L2	rs7903146	T	46.1	37	1.66*10⁻³	1.43 (1.14-1.79)	/
	D (1) vs 100+ (2)	TCF7L2	rs7903146	T	44.8	31.4	1.34*10⁻⁶	1.77 (1.40-2.24)	/
	D+Co (1) vs 100+ (2)	TCF7L2	rs7903146	T	45	31.4	1.06*10⁻⁶	1.78 (1.41-2.25)	/
	D+macroCo (1) vs 100+ (2)	TCF7L2	rs7903146	T	50	31.4	1.13*10⁻⁴	2.177 (1.46-3.24)	/
	D+microCo (1) vs 100+ (2)	TCF7L2	rs7903146	T	46	31.4	2.64*10⁻⁴	1.85 (1.33-2.57)	/
Genotypic association	D (1) vs 100+ (2)	ADIPOQ	rs266729	G/C	36.6	50.7	4.19*10⁻⁴	0.56 (0.41-0.77)	Over dominant
	D (1) vs CTR (2)	TCF7L2	rs7903146	TC/TT	58.9	70.0	1.41*10⁻⁴	1.63 (1.26-2.09)	Dominant
	D+Co (1) vs CTR (2)	TCF7L2	rs7903146	TC/TT	58.8	71.9	5.96*10⁻⁴	1.79 (1.28-2.52)	Dominant
	D (1) vs 100+ (2)	TCF7L2	rs7903146	CC	30	47.5	9.07*10⁻⁷	1.78 (1.41-2.26)	Log-additive
				TC	50	42.1
				TT	20	10.4
	D+Co (1) vs 100+ (2)	TCF7L2	rs7903146	CC	28.1	47.5	5.47*10⁻⁶	1.89 (1.43-2.50)	Log-additive
				TC	51.1	42.1
				TT	20.8	10.4
	D+macroCo (1) vs 100+ (2)	TCF7L2	rs7903146	CC	23.4	47.5	1.54*10⁻⁴	2.18 (1.45 −3.29)	Log-additive
				TC	53.1	42.1
				TT	23.4	10.4
	D+microCo (1) vs 100+ (2)	TCF7L2	rs7903146	CC	29.1	47.5	3.30*10⁻⁴	1.84 (1.31 −2.58)	Log-additive
				TC	50.0	42.1
				TT	20.9	10.4

Figure 1. TCF7L2 rs7903146 genotypic and allelic frequencies among the different groups analyzed. TCF7L2 rs7903146 genotypic and allelic frequencies were reported according to the decreasing of health/longevity and the increasing of T2D severity (from left to right).

All diabetic patients (D) vs non diabetic controls (CTR) (classical approach of genetic association studies)

Allelic and genotypic association analyses were performed on all diabetic patients vs the classical control group. Only TCF7L2 rs7903146 T-allele was associated with T2D with an odds ratio (OR) of 1.372 per risk allele (95% CI 1.155-1.629) (Table 1). Genotypic association of the TCF7L2 rs7903146 was significantly detected (p_value = 1.415*10⁻⁴, OR = 1.63 (95% CI 1.26-2.09)) under a dominant model (Table 1).

All diabetic patients (D) vs centenarians (100+)

Allelic and genotypic association analyses were performed on all diabetic patients vs centenarians. Two nominal significant differences in the allelic frequencies of diabetic patients and centenarians were observed in the ADIPOQ rs266729 (OR = 0.706 (95% CI 0.551 - 0.905), p_value = 5.89*10⁻³) and in TCF7L2 rs7903146. In this comparison the TCF7L2 rs7903146 was highly associated with T2D with an OR = 1.775 per risk allele (95% CI 1.404-2.243, p_value = 1.35*10⁻⁶) (Table 1). TCF7L2 association was supported by genotypic association analysis. Indeed, significant differences in genotypic frequencies for rs7903146 were reported (OR = 1.78 (95% CI 1.41-2.26), p_value = 9.066*10⁻⁷) (Table 1). Genotypic association was found also for rs266729 in ADIPOQ gene (p_value = 4.19*10⁻⁴).

Diabetic patients with complications (D+Co) vs non diabetic controls (CTR) Allelic and genotypic association analyses were performed on diabetic patients with at least one complication vs non diabetic controls. Only TCF7L2 rs7903146 showed an allelic association with T2D susceptibility with an OR = 1.434 (95% CI 1.145-1.795, p_value = 1.66*10⁻³). No allelic associations were detected for the other SNPs genotyped (Supplementary Table 2). Genotypic association was observed for the same SNP with an OR = 1.79 (95% CI 1.28-2.52, p_value = 5.96*10⁻⁴) (Table 1).

Table 2.

Diabetic patients and non diabetic controls collected by the Diabetology Unit in Ancona: samples description

		All samples		Males		Females
	Variables	Diabetic Patients (N=562)	Controls (N=558)	Diabetic Patients (N=305)	Controls (N=214)	Diabetic Patients (N=257)	Controls (N=344)
Traits	Age ± STD	65.76± 8.11	58.11 ± 12.40	64.93 ± 8.42	58.35 ± 11.88	66.75 ± 7.64	57.96 ±12.73
	Sex (M/F)	306/258	213/345	…	…	…	…
	BMI ± STD	28.76 ±4.62	27.14 ± 4.46	28.24 ± 4.15	27.29±3.51	29.40 ± 5.07	27.04 ±5.00
Complications	Retinopathy (%)	27.22 %	…	27.21 %	…	27.24 %	…
	Somatic Neuropathy (%)	17.79 %	…	21.64 %	…	13.23 %	…
	Renal failure (%)	3.56 %	…	4.92 %	…	1.94 %	…
	Nephropathy (%)	12.81 %	…	16.06 %	…	8.95 %	…
	Ischemic heart disease (%)	17.08 %	…	20.33 %	…	13.23 %	…
	Acute myocardial infarction (%)	9.07 %	…	13.11 %	…	4.28 %	…
	Ictus (%)	6.76 %	…	6.88 %	…	6.61%	…

D+Covs 100+, D+macroCo vs 100+, and D+microCo vs 100+ (comparison between extreme phenotypes)

Genotypic association analyses were performed on diabetic patients with at least one complication vs centenarians as a non-diabetes control group. When we compared D+Co vs 100+, D+macroCo vs 100+ and D+microCo vs 100+, the TCF7L2 rs7903146 resulted associated with T2D complicated patients with an OR = 1.89 (95% CI 1.43-2.50, p_value = 5.473*10⁻⁶), 2.18 (95% CI 1.45 -3.29, p value = 1.538*10⁻⁴) and 1.84 (95% CI 1.31 - 2.58, p_value = 2.68*10⁻⁴), respectively (Table 1). In Figure 1 allele and genotype frequencies for rs7903146 were reported for each group considered.

TCF7L2 interactome analysis

Since our results confirmed a core role of TCF7L2 in T2D pathogenesis, consistent with its implicated role in T2D [23,24], an interactome analysis was performed. Tcf7l2 is a 619 AA-long transcription factor that participates in the Wnt signaling pathway and modulates MYC expression by binding to its promoter in a sequence-specific manner. The most central elements of the Tcf7l2 interactome are Ctbn1 (ranking 1^st) and Ep300, and, on a second order of importance, Csk21, Myc and Cebpa (Figure 2A). Among the whole Tcf7l2 interactome (including co-localization and physical association) with 30 elements (Figure 2), there is evidence of direct physical interactions (experimentally validated protein-protein interaction) of Tcf7l2 (Figure 2B, red line) with 12 other proteins (Figure 2B, blue lines). Listing OMIM enrichment analysis of the Tcf7l2 interactome was performed and the first ten OMIM diseases significantly connected to the interactome are reported in Figure 3.

Figure 2. Tcf7l2 interactome. Graph of the Tcf7l2 interactome (A) and list of its elements (B). The most central elements in terms of betweenness centrality and node degree [85,89] are Ctnb1, Tfc7l2 (as expected) and Ep300. Non-connected nodes: TCF7L2 gene-listed gene interactions only, and no PPIs, are reported. (A) Node color code: from red to green, from highest to lowest values of betweeness centrality; interaction link color code: from red to green, from higher to lower number of experimental evidences. (B) Direct physical interactors are in blue, other types of interaction, such as colocalization and physical association are in black.

Figure 3. Major OMIM diseases linked to the Tcf7l2 interactome. OMIM enrichment analysis of the Tcf7l2 interactome was performed using the Enrichr web service (http://amp.pharm.mssm.edu/Enrichr/) and the first ten OMIM diseases significantly linked to the interactome (“Enrichr combined score”) are reported.

Discussion

In the last few years GWASs provided a great amount of data regarding genetic risk factors for major age-related diseases. Most variants identified so far confer relatively small increments in risk, and explain only a small proportion of the estimated trait heritability. Increasing sample size has been suggested as an answer to the missing heritability to increase sensitivity to detect variants with smaller effects. As an alternative for dramatically increasing cohort size, we tested a strategy that involves extreme phenotypes replicating a limited set of SNPs emerged and validated in GWASs. Individuals that are at the extreme ends of a trait distribution have been proposed for identifying variants that are rare but not private and that have modest to high effect sizes.

We hypothesized that the frequencies of GWASs risk variants are inversely correlated with increasing health and longevity, and extreme phenotypic differences between cases and controls would allow detection of the genetic variants with the most relevant biological and clinical value at relatively small sample size. For an extreme-trait approach, accurate phenotyping will be of vital importance. Central to our approach is the “control” group. The proper “normal” controls, especially in relation to common age-related diseases, are difficult to define, especially when large numbers of “healthy” subjects are collected from different regions, characterized by peculiar lifestyle habits, population genetics, among others. To circumvent these difficulties, we propose to use centenarians as an additional, more informative and defined control group for studies on the genetic determinants of age-related diseases, and particularly to establish the biological and clinical relevance of the genetic variants emerged in GWASs. Indeed, this group can be considered as a “super-control”, as most centenarians achieved their remarkable age avoiding or largely postponing major age-related diseases, and thus, at variance with much younger controls, we can be sure that they never developed such diseases. We reasoned that if a given genetic variant emerged from GWASs and considered a statistically significant but “weak” risk factor for the disease of interest, is present at similar frequency in centenarians and in patients, its biological relevance is reasonably negligible. On the contrary, if a given genetic variant considered a risk factor for an age-associated disease has a frequency much higher in patients than in centenarians, we can reasonably argue that this variant is a strong candidate to play a consistent biological role in the pathogenesis of the disease. We also assumed that centenarians, as the best model of healthy aging in humans [6], are a control group less blurred and undefined than the usually employed “age-matched controls”. Moreover, in order to properly take into account the above-mentioned largely neglected clinical heterogeneity/severity of age-related diseases in studies on their genetic determinants, we focused on T2D and more specifically on two main groups of patients, i.e. those who did not develop complications and can enjoy a relatively long and acceptable quality of life, and those who developed one or more (micro- and macro-vascular) complications with a poorer quality of life. Accordingly, this approach allowed us to compare extreme phenotypes, i.e. centenarians on one side and T2D patients with severe complications on the other, thus maximizing possible genetic differences among patients and controls.

Using this methodological approach on a total of 1,349 individuals we tested 31 SNPs in or nearby 16 different genes that were previously strongly associated with T2D and metabolic diseases and that could be considered good candidate to play a role in T2D complications, on the basis of other large genetic studies.

Our results indicate that most of the SNPs analyzed were not significantly associated with T2D, even when extreme groups of cases and controls were used. This somewhat unexpected result could be due to a lack of biological relevance in spite of their previously strong association signal in GWASs. Indeed, a major drawback of GWASs is the observation of many significant association signals from loci of weak effect. Therefore, we have to consider the possibility that such SNPs, despite having reached statistical significance in large GWASs, could have a weak biological relevance. Indeed, the presence of risk alleles at the same frequency in T2D patients with or without complications and in centenarians suggests that per se they do not represent a strong biological risk, being compatible with exceptional longevity, and that they likely necessitate to interact (epistatic effects) with (many) other risk alleles and/or specific environmental conditions to give a certain phenotype. Further studies, possibly in more than one population, are needed to test this possibility.

The two SNPs that reached statistical significance in our study were ADIPOQ rs266729 and TCF7L2 rs7903146, the last one being the strongest. ADIPOQ encodes for adiponectin, the most abundant adipokine in human plasma. Low adiponectin levels are associated to IR and T2D [26-28] and the chromosome locus 3q27, where ADIPOQ is located, is linked to metabolic syndrome and T2D [29]. Several studies have investigated the association between ADIPOQ polymorphisms and T2D, but their results are controversial, probably due to differences in the size and in the composition of the analyzed cohorts [30]. To this regard it is interesting to note that the statistical significance of ADIPOQ rs266729, a polymorphism associated to decreased adiponectin levels, emerged only when T2D patients were compared with centenarians. Regarding TCF7L2 rs7903146, localizedin an intronicregion of the gene, the most interesting findings were that the risk allele T and the genotype TT showed an increasing statistical significance, from modest genotypic association with the disease when all T2D patients (D) were compared with the age-matched controls (CTR) (OR=1.63; 95% CI 1.26-2.09) to intermediate values (OR=1.79; 95% CI 1.28-2.52) when T2D patients with complications (D+Co) were compared with CTR. These OR values perfectly overlap OR previously described in literature [31] in GWASs that analyzed thousands of samples. Finally, the OR and p_values showed the most significant results when centenarians were compared with T2D patients, and a further increase in OR was observed when centenarians and patients with complications, the extreme phenotypes, were compared (D vs 100+: OR=1.78 (95% CI 1.41-2.26), p_value=9.066*10⁻⁷; D+Co vs 100+: OR=1.89 (95% CI 1.43-2.50), p_value=5.473*10⁻⁶). It is important to note that the reverse is true for the C allele and the CC genotype. Notably, centenarians show not only the lowest frequency of T allele and TT genotype, but also the highest frequency of C allele and CC genotype, which emerges as a robust protective longevity variant. These data also suggest that TCF7L2 rs7903146 is a robust genetic risk variant and indicate that further analyses based on the use of centenarians as super-controls, possibly at genome wide level, are worthwhile to pursue in T2D study, and likely in other age-related diseases. TT genotype frequency was found to be higher not only in diabetic patients but also in patients with macro-vascular complications (Figure 1), reaching a frequency of 23.4%. These results demonstrate an important role of rs7903146 not only in T2D susceptibility but also in susceptibility to develop macro-vascular complications. Different studies investigated the relationship between macro-vascular complications, cardiovascular diseases and diabetes but the molecular mechanisms are still unknown [32-35].

The TCF7L2 polymorphism rs 7903146 was one of the first tobe related to diabetes [36]and subsequently its association has been confirmed in several studies, placing it among the most reproducible markers of T2D [37-40]. Recent studies have shown that this SNP leads to conformational changes inchromatin structure according to the allele type. In particular the risk allele (T) inducesan open conformation of the chromatin, resulting in an increased transcriptional activityand TCF7L2 mRNA levels [39,41,42]. It has been demonstrated that the increase of expression of TCF7L2 leads to a progressive reduction of glucose to lerance in vivo [43].

In order to gain information about the biological relevance of TCF7L2 gene, we performed an interactome analysis, which identified a complex network of interactions with Ctnb1 and Ep300 as the central elements. Ctnb1 is a Beta-catenin (adherens junction protein), critical for the establishment and maintenance of epithelial layers, and a key downstream component of the canonical Wnt signaling pathway. Ep300 is a histone acetyltransferase and regulates transcription via chromatin remodeling. Chromatin remodeling and epigenetic mechanisms seem to play a fundamental role in susceptibility to diabetes and related complications [44-47] even if mechanisms are not completely understood and further studies should be designed to clarify them. The in silico analysis suggests that Tcf7l2 interactome is involved in a variety of age-related pathologies, such as Alzheimer's disease, cardiovascular diseases and colorectal cancer [48-52], all sharing vascular alterations, as a further support to the biological relevance of this gene and to the emerging idea that the same gene can play a major role in apparently different diseases [53].

On the whole, our data on TCF7L2 rs7903146 indicate that the frequency of T2D risk genotypes decreases according to the severity of the disease phenotype, supporting the assumption that centenarians represent a powerful and informative control group in association studies on T2D, a pathology characterized by significant clinical heterogeneity. Different authors have emphasized the problem of phenotypical heterogeneity of cohorts, which can seriously affect the ability todetect genetic associations [54]. Our data suggest that the use of relatively small numbers of cases and controls with extreme phenotypes (thus reducing sample heterogeneity) in genetic association studies of age-related diseases can be successful in detecting significant associations as suggested by other authors [55].

Methods

Samples

A total of 1,349 individuals from Northern/Central Italy, including 562 T2D patients (D) (mean age: 65.76 ± 8.11) and 558 unrelated age-, gender- and geographically-matched controls (CTR) (mean age: 58.11 ± 12.40), and 229 centenarians (100+) (mean age: 105.0 ± 2.9) have been considered in this study. D and CTR subjects were collected by the Diabetology Unit, INRCA (National Institute on Health and Science on Aging) in Ancona (Italy). T2D diagnosis was made according to the American Diabetes Association Criteria (www.diabetes.org) All patients were fully characterized from a clinical point of view and a large number of biochemical/endocronological parameters were measured. To avoid population stratification effects, only patients and controls with at least two generations of maternal ancestry from the Marche region (Central Italy) were included in this study. CTR subjects were carefully assessed and fully characterized according to the same protocol used for T2D patients, and a detailed clinical history was recorded in order to exclude the presence of T2D and of any other overt illness. Within the group of diabetic patients, N = 241 were affected by complications (D+Co), and in particular N = 110 had only micro-vascular (D+microCo) and N = 64 had only macro-vascular complications (D+macroCo). The main phenotypical and clinical characteristics of D and CTR are summarized in Table 2.

The presence/absence of microvascular and macrovascular diabetic complications was assessed according to the following criteria: i) microvascular: diabetic retinopathy by fundoscopy through dilated pupils and/or fluorescence angiography; incipient nephropathy, defined by an excessive urinary albumin excretion (>30 mg/24 h) and a normal creatinine clearance; renal failure, detected as an estimated glomerular filtration rate >60 mL/min per 1.73 m2; neuropathy established by electromyography; ii) macrovascular: ictus by clinical history; ischemic heart disease and acute myocardial infarction by clinical history and by resting electrocardiogram; peripheral vascular disease by clinical history and, for lower limbs, by ankle-brachial index. Hypertension was defined as a systolic blood pressure >140 mmHg and/or a diastolic blood pressure >90 mmHg. The values were measured while the subjects were sitting and confirmed at least three times. Overnight fasting venous blood samples from all subjects were collected from 8:00 to 9:00 a.m. Blood concentrations for HDL cholesterol, triglycerides, HbA1c, fasting insulin, fibrinogen, high-sensitivity C reactive protein (hsCRP), creatinine, urea nitrogen, and white blood cells count were measured by standard procedures.

All centenarians were born in Italy between the 1900 and 1908. Trained physicians and nursing staff collected demographic and lifestyle data, anthropometric measurements, functional, cognitive and health status, clinical anamnesis, and details on drug use, as previously described [9].

SNPs selection

A total of 31 SNPs mapping within and nearby 16 different genes (Table S5) were genotyped in the 1,349 individuals included in the study. The SNPs were selected according to the following criteria: i) SNPs reported to be risk factor for T2D [56-62]such as Transcription factor 7-like 2 (TCF7L2), Insulin-like growth factor 2 mRNA-binding proteins (IGF2BP), Potassium inwardly-rectifying channel, subfamily J, member 11 (KCNJ11), Potassium voltage-gated channel KQT-like sub-family, member 1 (KCNQ1); ii)SNPs of candidate genes for the development of T2D complications, such as Catalase (CAT), erythropoietin (EPO) [61], hypoxia-inducible factor 1 α subunit (HIF1A) [63] and DDAH1. Interestingly, significant association of CAT and SOD polymorphisms with T2D were recently reported [64]. Genetic variations in the DDAH1 and DDAH2 genes are significantly associated with serum Asymmetric DiMethylArginine (ADMA) levels in T2D patients [65]. Moreover, it was reported that DDAH gene polymorphisms play a central role in determining ADMA in diabetic renal impairment [66]. iii)SNPs of genes reported to be involved in metabolic diseases, including T2D, such as ADIPOQ and IRS1 encoding for adiponectin and insulin receptor substrate 1 respectively, and FTO, a major gene for obesity susceptibility [62-65]. Adiponectin is an adipocyte-produced protein involved in regulating glucose, lipid and energy metabolism, and ADIPOQ polymorphisms were previously associated with T2D in Caucasian and non-Caucasian populations [59,71,72]. IRS1 polymorphisms have been associated with risk of T2D and adiposity in GWASs. Recently, it was reported that IRS1 rs2943641 interacts with carbohydrate and fat intakes in incident T2D in a sex-specific manner [73]; iv)SNPs of genes relevant for aging and age-related diseases, such as hTERT and TERC, involved intelomere stability/attrition and in T2D complications [21,74,75].

SNPs genotyping

DNA was extracted from whole blood (QIAmp 96 DNA Blood kit, QIAGEN). Genotyping analysis was performed by using SEQUENOM MassArray iPLEX technology, following the manufacturer's instructions and as previously described [76,77]. Genotype calls were analyzed by using SEQUENOM Typer 4.0 software and the individual spectrograms were checked in order to evaluate the presence of calling errors.

SNP analysis and quality control

Four SNPs (rs16889462, rs4880, rs8047395 and rs10434) that did not satisfy the Hardy-Weinberg equilibrium (HWE) in the control group and five SNPs (rs669173, rs2853669, rs3025021, rs13266634 and rs7901695) that showed 25% missing call rates were excluded from the analysis (Table S5). Diabetic patients were also analyzed after stratifications for the type of complications (micro and macro-vascular).

Only SNPs with call rates higher than 95%, without significant deviation from HWE in controls and with minor allele frequency (MAF) exceeding 5% were retained for the association analysis. Genotypic analysis was performed using R 2.15.2 and SNPassoc package (http://www.r-project.org) [78] that implements binary logistic regression methods under five different genetic models. Allelic association analyses were performed using PLINK 1.07, an open-source whole genome association analysis toolset (http://pngu.mgh.harvard.edu/purcell/plink/) [71].

Interactome analysis

Protein-protein interaction (PPI) data have been retrieved from the Agile Protein Interactions Database APID [80] accessed through the dedicated plugin APID2NET [81] and from the metadatabase InnateDB [82] and then integrated to form the Tcf7l2 interactome (i.e. the set of molecular interactions related to Tcf7l2) and analyzed with the network analysis platform Cytoscape [83,84]. Topological measures such as betweenness centrality and node degree are used to rank the importance of the interactome elements [85,86]. The betweenness centrality of a node reflects the amount of control that this node exerts over the interactions of other nodes in the network. This measure favors nodes that join communities (dense subnetworks), rather than nodes that lie inside a community [87,88]. OMIM diseases overrepresentation analyses and ranking have been performed using Enrichr (http://amp.pharm.mssm.edu/Enrichr).

Supplementary Materials

Supplemental Data

Acknowledgments

We thank Vilma Mantovani and Elena Marasco for their technical support during experimental procedure of genotyping at CRBA (Applied Biomedical Research Center, S. Orsola-Malpighi Polyclinic, Bologna, Italy). We also thank and mention the European Union's Seventh Framework Programme (grant agreement no. 259679, “IDEAL” and grand agreement no. 600803, “MISSION T2D”) and the Italian Ministry of University and Research (Project PRIN 2006 to DM and Project PRIN 2009 to CF).

Conflicts of Interest

The authors declare no conflicts of interest.

References

1. Bodmer W and Bonilla C. Common and rare variants in multifactorial susceptibility to common diseases. Nature Genetics. 2008; 40: 695 -701. [PubMed] .
2. Tazearslan C, Cho M, Suh Y. Discovery of functional gene variants associated with human longevity: opportunities and challenges. J. Gerontol. A Biol. Sci. Med. Sci. 2012; 67: 376 -383. [PubMed] .
3. Heutink P and Oostra BA. Gene finding in genetically isolated populations. Hum. Mol. Genet. 2002; 11: 2507 -2515. [PubMed] .
4. Garagnani P, Laayouni H, González-Neira A, Sikora M, Luiselli D, Bertranpetit J, et al. Isolated populations as treasure troves in genetic epidemiology: the case of the Basques. Eur. J. Hum. Genet. 2009; 17: 1490 -1494. [PubMed] .
5. Lowe JK, Maller JB, Pe'er I, Neale BM, Salit J, Kenny EE, et al. Genome-wide association studies in an isolated founder population from the Pacific Island of Kosrae. PLoS Genet. 2009; 5: e1000365 [PubMed] .
6. Franceschi C and Bonafè M. Centenarians as a model for healthy aging. Biochem. Soc. Trans. 2003; 31: 457 -461. [PubMed] .
7. Paolisso G, Gambardella A, Ammendola S, D'Amore A, Balbi V, Varricchio M, et al. Glucose tolerance and insulin action in healty centenarians. Am. J. Physiol. 1996; 270: E890 -894. [PubMed] .
8. Paolisso G, Barbieri M, Rizzo MR, Carella C, Rotondi M, Bonafè M, et al. Low insulin resistance and preserved beta-cell function contribute to human longevity but are not associated with TH-INS genes. Exp. Gerontol. 2001; 37: 149 -156. [PubMed] .
9. Vitale G, Brugts MP, Ogliari G, Castaldi D, Fatti LM, Varewijck AJ, et al. Low circulating IGF-I bioactivity is associated with human longevity: findings in centenarians' offspring. Aging (Albany NY). 2012; 4: 580 -589. [PubMed] .
10. Evert J, Lawler E, Bogan H, Perls T. Morbidity profiles of centenarians: survivors, delayers, and escapers. J. Gerontol. A Biol. Sci. Med. Sci. 2003; 58232 -52837. .
11. Atzmon G, Schechter C, Greiner W, Davidson D, Rennert G, Barzilai N. Clinical phenotype of families with longevity. J Am Geriatr Soc. 2004; 52: 274 -277. [PubMed] .
12. Westendorp RGJ, van Heemst D, Rozing MP, Frölich M, Mooijaart SP, Blauw G-J, et al. Nonagenarian siblings and their offspring display lower risk of mortality and morbidity than sporadic nonagenarians: The Leiden Longevity Study. J Am Geriatr Soc. 2009; 57: 1634 -1637. [PubMed] .
13. Wijsman CA, Rozing MP, Streefland TCM, le Cessie S, Mooijaart SP, Slagboom PE, et al. Familial longevity is marked by enhanced insulin sensitivity. Aging Cell. 2011; 10: 114 -121. [PubMed] .
14. Kowall B, Peters A, Thorand B, Rathmann W, Meisinger C. Associations between longevity of parents and glucose regulation in their offspring: the KORA S4/F4 Study. Diabetologia. 2013; 56: 268 -274. [PubMed] .
15. Testa R, Olivieri F, Bonfigli AR, Sirolla C, Boemi M, Marchegiani F, et al. Interleukin-6-174 G > C polymorphism affects the association between IL-6 plasma levels and insulin resistance in type 2 diabetic patients. Diabetes Res. Clin. Pract. 2006; 71: 299 -305. [PubMed] .
16. Cormio A, Milella F, Marra M, Pala M, Lezza AMS, Bonfigli AR, et al. Variations at the H-strand replication origins of mitochondrial DNA and mitochondrial DNA content in the blood of type 2 diabetes patients. Biochim. Biophys. Acta. 2009; 1787: 547 -552. [PubMed] .
17. Bonfigli AR, Sirolla C, Testa R, Cucchi M, Spazzafumo L, Salvioli S, et al. The p53 codon 72 (Arg72Pro) polymorphism is associated with the degree of insulin resistance in type 2 diabetic subjects: a cross-sectional study. Acta Diabetol. 2012; in press .
18. Marra M, Marchegiani F, Ceriello A, Sirolla C, Boemi M, Franceschi C, et al. Chronic renal impairment and DDAH2-1151 A/C polymorphism determine ADMA levels in type 2 diabetic subjects. Nephrol. Dial. Transplant. 2013; 28: 964 -971. [PubMed] .
19. Spazzafumo L, Olivieri F, Abbatecola AM, Castellani G, Monti D, Lisa R, et al. Remodelling of biological parameters during human ageing: evidence for complex regulation in longevity and in type 2 diabetes. Age (Dordr). 2013; 35: 419 -429. [PubMed] .
20. Forte GI, Pilato G, Vaccarino L, Sanacore M, Candore G, Romano GC, et al. Risk profiles in type 2 diabetes (metabolic syndrome): integration of IL-10 polymorphisms and laboratory parameters to identify vascular damages related complications. Curr. Pharm. 2010; 16: 898 -903. .
21. Testa R, Olivieri F, Sirolla C, Spazzafumo L, Rippo MR, Marra M, et al. Leukocyte telomere length is associated with complications of type 2 diabetes mellitus. Diabet. Med. 2011; 28: 1388 -1394. [PubMed] .
22. Achilli A, Olivieri A, Pala M, Hooshiar Kashani B, Carossa V, Perego UA, et al. Mitochondrial DNA backgrounds might modulate diabetes complications rather than T2DM as a whole. PLoS ONE. 2011; 6: e21029 [PubMed] .
23. Grant SFA, Thorleifsson G, Reynisdottir I, Benediktsson R, Manolescu A, Sainz J, et al. Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes. Nat. Genet. 2006; 38: 320 -323. [PubMed] .
24. Zhao J, Schug J, Li M, Kaestner KH, Grant SFA. Disease-associated loci are significantly over-represented among genes bound by transcription factor 7-like 2 (TCF7L2) in vivo. Diabetologia. 2010; 53: 2340 -2346. [PubMed] .
25. Franceschi C and Bonafè M. Centenarians as a model for healthy aging. Biochem. Soc. Trans. 2003; 31: 457 -461. [PubMed] .
26. Hotta K, Funahashi T, Arita Y, Takahashi M, Matsuda M, Okamoto Y, et al. Plasma concentrations of a novel, adipose-specific protein, adiponectin, in type 2 diabetic patients. Arterioscler. Thromb. Vasc. Biol. 2000; 20: 1595 -1599. [PubMed] .
27. Daimon M, Oizumi T, Saitoh T, Kameda W, Hirata A, Yamaguchi H, et al. Decreased serum levels of adiponectin are a risk factor for the progression to type 2 diabetes in the Japanese Population: the Funagata study. Diabetes Care. 2003; 26: 2015 -2020. [PubMed] .
28. Lindsay RS, Funahashi T, Hanson RL, Matsuzawa Y, Tanaka S, Tataranni PA, et al. Adiponectin and development of type 2 diabetes in the Pima Indian population. Lancet. 2002; 360: 57 -58. [PubMed] .
29. Vionnet N, Hani EH, Dupont S, Gallina S, Francke S, Dotte S, et al. Genomewide search for type 2 diabetes-susceptibility genes in French whites: evidence for a novel susceptibility locus for early-onset diabetes on chromosome 3q27-qter and independent replication of a type 2-diabetes locus on chromosome 1q21-q24. Am. J. Hum. Genet. 2000; 67: 1470 -1480. [PubMed] .
30. Han LY, Wu QH, Jiao ML, Hao YH, Liang LB, Gao LJ, et al. Associations between single-nucleotide polymorphisms (+45T>G, +276G>T, -11377C>G, -11391G>A) of adiponectin gene and type 2 diabetes mellitus: a systematic review and meta-analysis. Diabetologia. 2011; 54: 2303 -2314. [PubMed] .
31. Scott RA, Lagou V, Welch RP, Wheeler E, Montasser ME, Luan J, et al. Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways. Nat. Genet. 2012; 44: 991 -1005. [PubMed] .
32. Mohlke KL, Boehnke M, Abecasis GR. Metabolic and cardiovascular traits: an abundance of recently identified common genetic variants. Hum. Mol. Genet. 2008; 17: R102 -108. [PubMed] .
33. Sousa AG, Selvatici L, Krieger JE, Pereira AC. Association between genetics of diabetes, coronary artery disease, and macrovascular complications: exploring a common ground hypothesis. Rev Diabet Stud. 2011; 8: 230 -244. [PubMed] .
34. Ciccacci C, Di Fusco D, Cacciotti L, Morganti R, D'Amato C, Novelli G, et al. TCF7L2 gene polymorphisms and type 2 diabetes: association with diabetic retinopathy and cardiovascular autonomic neuropathy. Acta Diabetol. 2012; in press .
35. Muendlein A, Saely CH, Geller-Rhomberg S, Sonderegger G, Rein P, Winder T, et al. Single nucleotide polymorphisms of TCF7L2 are linked to diabetic coronary atherosclerosis. PLoS ONE. 2011; 6: e17978 [PubMed] .
36. Reynisdottir I, Thorleifsson G, Benediktsson R, Sigurdsson G, Emilsson V, Einarsdottir AS, et al. Localization of a susceptibility gene for type 2 diabetes to chromosome 5q34-q35.2. Am. J. Hum. Genet. 2003; 73: 323 -335. [PubMed] .
37. Mayans S, Lackovic K, Lindgren P, Ruikka K, Agren A, Eliasson M, et al. TCF7L2 polymorphisms are associated with type 2 diabetes in northern Sweden. Eur. J. Hum. Genet. 2007; 15: 342 -346. [PubMed] .
38. Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, et al. A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007; 445: 881 -885. [PubMed] .
39. Stitzel ML, Sethupathy P, Pearson DS, Chines PS, Song L, Erdos MR, et al. Global epigenomic analysis of primary human pancreatic islets provides insights into type 2 diabetes susceptibility loci. Cell Metab. 2010; 12: 443 -455. [PubMed] .
40. Van Vliet-Ostaptchouk JV, Shiri-Sverdlov R, Zhernakova A, Strengman E, van Haeften TW, Hofker MH, et al. Association of variants of transcription factor 7-like 2 (TCF7L2) with susceptibility to type 2 diabetes in the Dutch Breda cohort. Diabetologia. 2007; 50: 59 -62. [PubMed] .
41. Lyssenko V, Lupi R, Marchetti P, Del Guerra S, Orho-Melander M, Almgren P, et al. Mechanisms by which common variants in the TCF7L2 gene increase risk of type 2 diabetes. J. Clin. Invest. 2007; 117: 2155 -2163. [PubMed] .
42. Gaulton KJ, Nammo T, Pasquali L, Simon JM, Giresi PG, Fogarty MP, et al. A map of open chromatin in human pancreatic islets. Nat. Genet. 2010; 42: 255 -259. [PubMed] .
43. Ip W, Chiang Y-TA, Jin T. The involvement of the wnt signaling pathway and TCF7L2 in diabetes mellitus: The current understanding, dispute, and perspective. Cell Biosci. 2012; 2: 28 [PubMed] .
44. Ling C and Groop L. Epigenetics: a molecular link between environmental factors and type 2 diabetes. Diabetes. 2009; 58: 2718 -2725. [PubMed] .
45. Pinney SE and Simmons RA. Epigenetic mechanisms in the development of type 2 diabetes. Trends Endocrinol. Metab. 2010; 21: 223 -229. [PubMed] .
46. Pons D, de Vries FR, van den Elsen PJ, Heijmans BT, Quax PHA, Jukema JW. Epigenetic histone acetylation modifiers in vascular remodelling: new targets for therapy in cardiovascular disease. Eur. Heart J. 2009; 30: 266 -277. [PubMed] .
47. Villeneuve LM and Natarajan R. The role of epigenetics in the pathology of diabetic complications. Am. J. Physiol. Renal Physiol. 2010; 299: F14 -25. [PubMed] .
48. Pimentel-Coelho PM and Rivest S. The early contribution of cerebrovascular factors to the pathogenesis of Alzheimer's disease. Eur. J. Neurosci. 2012; 35: 1917 -1937. [PubMed] .
49. Camici PG. Coronary microvascular dysfunction in patients with cardiomyopathies. Circ Heart Fail. settembre. 2008; 1: 150 -152. .
50. Sainz J, Rudolph A, Hoffmeister M, Frank B, Brenner H, Chang-Claude J, et al. Effect of type 2 diabetes predisposing genetic variants on colorectal cancer risk. J. Clin. Endocrinol. Metab. 2012; 97: E845 -851. [PubMed] .
51. Sun J, Wang D, Jin T. Insulin alters the expression of components of the Wnt signaling pathway including TCF-4 in the intestinal cells. Biochim. Biophys. Acta. 2010; 1800: 344 -351. [PubMed] .
52. Jin T. Why diabetes patients are more prone to the development of colon cancer? Med. Hypotheses. 2008; 71: 241 -244. [PubMed] .
53. Vidal M, Cusick ME, Barabási A-L. Interactome networks and human disease. Cell. 2011; 144: 986 -998. [PubMed] .
54. Moonesinghe R, Khoury MJ, Liu T, Ioannidis JPA. Required sample size and nonreplicability thresholds for heterogeneous genetic associations. Proc. Natl. Acad. Sci. U.S.A. 2008; 105: 617 -622. [PubMed] .
55. Perry JRB, Voight BF, Yengo L, Amin N, Dupuis J, Ganser M, et al. Stratifying type 2 diabetes cases by BMI identifies genetic risk variants in LAMA1 and enrichment for risk variants in lean compared to obese cases. PLoS Genet. 2012; 8: e1002741 [PubMed] .
56. Qi L, Cornelis MC, Kraft P, Stanya KJ, Linda Kao WH, Pankow JS, et al. Genetic variants at 2q24 are associated with susceptibility to type 2 diabetes. Hum. Mol. Genet. 2010; 19: 2706 -2715. [PubMed] .
57. Tsai F-J, Yang C-F, Chen C-C, Chuang L-M, Lu C-H, Chang C-T, et al. A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese. PLoS Genet. 2010; 6: e1000847 [PubMed] .
58. Yamauchi T, Hara K, Maeda S, Yasuda K, Takahashi A, Horikoshi M, et al. A genome-wide association study in the Japanese population identifies susceptibility loci for type 2 diabetes at UBE2E2 and C2CD4A-C2CD4B. Nat. Genet. 2010; 42: 864 -868. [PubMed] .
59. Mtiraoui N, Ezzidi I, Turki A, Chaieb A, Mahjoub T, Almawi WY. Single-nucleotide polymorphisms and haplotypes in the adiponectin gene contribute to the genetic risk for type 2 diabetes in Tunisian Arabs. Diabetes Res. Clin. Pract. 2012; 97: 290 -297. [PubMed] .
60. Dupuis J, Langenberg C, Prokopenko I, Saxena R, Soranzo N, Jackson AU, et al. New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat. Genet. 2010; 42: 105 -116. [PubMed] .
61. Williams WW, Salem RM, McKnight AJ, Sandholm N, Forsblom C, Taylor A, et al. Association testing of previously reported variants in a large case-control meta-analysis of diabetic nephropathy. Diabetes. 2012; 61: 2187 -2194. [PubMed] .
62. Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, et al. Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat. Genet. 2010; 42: 579 -589. [PubMed] .
63. Gu HF, Zheng X, Abu Seman N, Gu T, Botusan IR, Sunkari VG, et al. Impact of the Hypoxia-Inducible Factor-1 α (HIF-1α) Pro582Ser Polymorphism on Diabetes Nephropathy. Diabetes Care. 2013; 36: 415 -421. [PubMed] .
64. Ghattas MH and Abo-Elmatty DM. Association of polymorphic markers of the catalase and superoxide dismutase genes with type 2 diabetes mellitus. DNA Cell Biol. 2012; 31: 1598 -1603. [PubMed] .
65. Abhary S, Burdon KP, Kuot A, Javadiyan S, Whiting MJ, Kasmeridis N, et al. Sequence variation in DDAH1 and DDAH2 genes is strongly and additively associated with serum ADMA concentrations in individuals with type 2 diabetes. PLoS ONE. 2010; 5: e9462 [PubMed] .
66. Marra M, Marchegiani F, Ceriello A, Sirolla C, Boemi M, Franceschi C, et al. Chronic renal impairment and DDAH2-1151 A/C polymorphism determine ADMA levels in type 2 diabetic subjects. Nephrol. Dial. Transplant. 2013; 28: 964 -971. [PubMed] .
67. Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, Lindgren CM, et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science. 2007; 316: 889 -894. [PubMed] .
68. Sällman Almén M, Rask-Andersen M, Jacobsson JA, Ameur A, Kalnina I, Moschonis G, et al. Determination of the obesity-associated gene variants within the entire FTO gene by ultra-deep targeted sequencing in obese and lean children. Int J Obes (Lond). 2013; 37: 424 -314. [PubMed] .
69. Jacobsson JA, Schiöth HB, Fredriksson R. The impact of intronic single nucleotide polymorphisms and ethnic diversity for studies on the obesity gene FTO. Obes Rev. 2012; 13: 1096 -1109. [PubMed] .
70. Dina C, Meyre D, Gallina S, Durand E, Körner A, Jacobson P, et al. Variation in FTO contributes to childhood obesity and severe adult obesity. Nat. Genet. 2007; 39: 724 -726. [PubMed] .
71. Han LY, Wu QH, Jiao ML, Hao YH, Liang LB, Gao LJ, et al. Associations between single-nucleotide polymorphisms (+45T>G, +276G>T, -11377C>G, -11391G>A) of adiponectin gene and type 2 diabetes mellitus: a systematic review and meta-analysis. Diabetologia. 2011; 54: 2303 -2314. [PubMed] .
72. Hivert M-F, Manning AK, McAteer JB, Florez JC, Dupuis J, Fox CS, et al. Common variants in the adiponectin gene (ADIPOQ) associated with plasma adiponectin levels, type 2 diabetes, and diabetes-related quantitative traits: the Framingham Offspring Study. Diabetes. 2008; 57: 3353 -3359. [PubMed] .
73. Ericson U, Rukh G, Stojkovic I, Sonestedt E, Gullberg B, Wirfält E, et al. Sex-specific interactions between the IRS1 polymorphism and intakes of carbohydrates and fat on incident type 2 diabetes. Am. J. Clin. Nutr. 2013; 97: 208 -216. [PubMed] .
74. Olivieri F, Lazzarini R, Recchioni R, Marcheselli F, Rippo MR, Di Nuzzo S, et al. MiR-146a as marker of senescence-associated pro-inflammatory status in cells involved in vascular remodelling. Age (Dordr). 2012; in press .
75. Olivieri F, Lorenzi M, Antonicelli R, Testa R, Sirolla C, Cardelli M, et al. Leukocyte telomere shortening in elderly Type2DM patients with previous myocardial infarction. Atherosclerosis. 2009; 206: 588 -593. [PubMed] .
76. Hasan SK, Buttari F, Ottone T, Voso MT, Hohaus S, Marasco E, et al. Risk of acute promyelocytic leukemia in multiple sclerosis: coding variants of DNA repair genes. Neurology. 22 marzo. 2011; 76: 1059 -1065. .
77. Mantovani V, Garagnani P, Selva P, Rossi C, Ferrari S, Cenci M, et al. Simple method for haplotyping the poly(TG) repeat in individuals carrying the IVS8 5T allele in the CFTR gene. Clin. Chem. 2007; 53: 531 -533. [PubMed] .
78. González JR, Armengol L, Solé X, Guinó E, Mercader JM, Estivill X, et al. SNPassoc: an R package to perform whole genome association studies. Bioinformatics. 2007; 23: 644 -645. [PubMed] .
79. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007; 81: 559 -575. [PubMed] .
80. Prieto C and De Las Rivas J. APID: Agile Protein Interaction DataAnalyzer. Nucleic Acids Res. 2006; 34: W298 -302. [PubMed] .
81. Hernandez-Toro J, Prieto C, De las Rivas J. APID2NET: unified interactome graphic analyzer. Bioinformatics. 2007; 23: 2495 -2497. [PubMed] .
82. Lynn DJ, Winsor GL, Chan C, Richard N, Laird MR, Barsky A, et al. InnateDB: facilitating systems-level analyses of the mammalian innate immune response. Mol. Syst. Biol. 2008; 4: 218 [PubMed] .
83. Killcoyne S, Carter GW, Smith J, Boyle J. Cytoscape: a community-based framework for network modeling. Methods Mol. Biol. 2009; 563: 219 -239. [PubMed] .
84. Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011; 27: 431 -432. [PubMed] .
85. Platzer A, Perco P, Lukas A, Mayer B. Characterization of protein-interaction networks in tumors. BMC Bioinformatics. 2007; 8: 224 [PubMed] .
86. Barabási A-L and Oltvai ZN. Network biology: understanding the cell's functional organization. Nat. Rev. Genet. 2004; 5: 101 -113. [PubMed] .
87. Yoon J, Blumer A, Lee K. An algorithm for modularity analysis of directed and weighted biological networks based on edge-betweenness centrality. Bioinformatics. 2006; 22: 3106 -3108. [PubMed] .
88. Assenov Y, Ramírez F, Schelhorn S-E, Lengauer T, Albrecht M. Computing topological parameters of biological networks. Bioinformatics. 2008; 24: 282 -284. [PubMed] .
89. Potapov AP, Voss N, Sasse N, Wingender E. Topology of mammalian transcription networks. Genome Inform. 2005; 16: 270 -278. [PubMed] .