Genetic Variants Associated With Quantitative Glucose Homeostasis Traits Translate to Type 2 Diabetes in Mexican Americans: The GUARDIAN (Genetics Underlying Diabetes in Hispanics) Consortium

A genome-wide linkage and association analysis of imputed insertions and deletions with cardiometabolic phenotypes in Mexican Americans: The Insulin Resistance Atherosclerosis Family Study

Insertions and deletions (INDELs) represent a significant fraction of interindividual variation in the human genome yet their contribution to phenotypes is poorly understood. To confirm the quality of imputed INDELs and investigate their roles in mediating cardiometabolic phenotypes, genome-wide association and linkage analyses were performed for 15 phenotypes with 1,273,952 imputed INDELs in 1,024 Mexican-origin Americans. Imputation quality was validated using whole exome sequencing with an average kappa of 0.93 in common INDELs (minor allele frequencies [MAFs] ≥ 5%). Association analysis revealed one genome-wide significant association signal for the cholesterylester transfer protein gene (CETP) with high-density lipoprotein levels (rs36229491, P = 3.06 × 10^-12 ); linkage analysis identified two peaks with logarithm of the odds (LOD) > 5 (rs60560566, LOD = 5.36 with insulin sensitivity (S_I ) and rs5825825, LOD = 5.11 with adiponectin levels). Suggestive overlapping signals between linkage and association were observed: rs59849892 in the WSC domain containing 2 gene (WSCD2) was associated and nominally linked with S_I (P = 1.17 × 10^-7 , LOD = 1.99). This gene has been implicated in glucose metabolism in human islet cell expression studies. In addition, rs201606363 was linked and nominally associated with low-density lipoprotein (P = 4.73 × 10^-4 , LOD = 3.67), apolipoprotein B (P = 1.39 × 10^-3 , LOD = 4.64), and total cholesterol (P = 1.35 × 10^-2 , LOD = 3.80) levels. rs201606363 is an intronic variant of the UBE2F-SCLY (where UBE2F is ubiquitin-conjugating enzyme E2F and SCLY is selenocysteine lyase) fusion gene that may regulate cholesterol through selenium metabolism. In conclusion, these results confirm the feasibility of imputing INDELs from array-based single nucleotide polymorphism (SNP) genotypes. Analysis of these variants using association and linkage replicated previously identified SNP signals and identified multiple novel INDEL signals. These results support the inclusion of INDELs into genetic studies to more fully interrogate the spectrum of genetic variation.

The autoimmune risk gene ZMIZ1 is a vitamin D responsive marker of a molecular phenotype of multiple sclerosis

Multiple Sclerosis (MS) is a neurological condition driven in part by immune cells from the peripheral circulation, the targets for current successful therapies. The autoimmune and MS risk gene ZMIZ1 is underexpressed in blood in people with MS. We show that, from three independent sets of transcriptomic data, expression of ZMIZ1 is tightly correlated with that of hundreds of other genes. Further we show expression is partially heritable (heritability 0.26), relatively stable over time, predominantly in plasmacytoid dendritic cells and non-classical monocytes, and that levels of ZMIZ1 protein expression are reduced in MS. ZMIZ1 gene expression is increased in response to calcipotriol (1,25 Vitamin D3) (p < 0.0003) and associated with Epstein Barr Virus (EBV) EBNA-1 antibody titre (p < 0.004). MS therapies fingolimod and dimethyl fumarate altered blood ZMIZ1 gene expression compared to untreated MS. The phenotype indicates susceptibility to MS, and may correspond with clinical response and represent a novel clinical target.

Insulin Clearance Is Associated with Hepatic Lipase Activity and Lipid and Adiposity Traits in Mexican Americans

Reduction in insulin clearance plays an important role in the compensatory response to insulin resistance. Given the importance of this trait to the pathogenesis of diabetes, a deeper understanding of its regulation is warranted. Our goal was to identify metabolic and cardiovascular traits that are independently associated with metabolic clearance rate of insulin (MCRI). We conducted a cross-sectional analysis of metabolic and cardiovascular traits in 765 participants from the Mexican-American Coronary Artery Disease (MACAD) project who had undergone blood sampling, oral glucose tolerance test, euglycemic-hyperinsulinemic clamp, dual-energy X-ray absorptiometry, and carotid ultrasound. We assessed correlations of MCRI with traits from seven domains, including anthropometry, biomarkers, cardiovascular, glucose homeostasis, lipase activity, lipid profile, and liver function tests. We found inverse independent correlations between MCRI and hepatic lipase (P = 0.0004), insulin secretion (P = 0.0002), alanine aminotransferase (P = 0.0045), total fat mass (P = 0.014), and diabetes (P = 0.03). MCRI and apolipoprotein A-I exhibited a positive independent correlation (P = 0.035). These results generate a hypothesis that lipid and adiposity associated traits related to liver function may play a role in insulin clearance.

Epigenome-wide association study of triglyceride postprandial responses to a high-fat dietary challenge

Postprandial lipemia (PPL), the increased plasma TG concentration after consuming a high-fat meal, is an independent risk factor for CVD. Individual responses to a meal high in fat vary greatly, depending on genetic and lifestyle factors. However, only a few loci have been associated with TG-PPL response. Heritable epigenomic changes may be significant contributors to the unexplained inter-individual PPL variability. We conducted an epigenome-wide association study on 979 subjects with DNA methylation measured from CD4⁺ T cells, who were challenged with a high-fat meal as a part of the Genetics of Lipid Lowering Drugs and Diet Network study. Eight methylation sites encompassing five genes, LPP, CPT1A, APOA5, SREBF1, and ABCG1, were significantly associated with PPL response at an epigenome-wide level (P < 1.1 × 10⁻⁷), but no methylation site reached epigenome-wide significance after adjusting for baseline TG levels. Higher methylation at LPP, APOA5, SREBF1, and ABCG1, and lower methylation at CPT1A methylation were correlated with an increased TG-PPL response. These PPL-associated methylation sites, also correlated with fasting TG, account for a substantially greater amount of phenotypic variance (14.9%) in PPL and fasting TG (16.3%) when compared with the genetic contribution of loci identified by our previous genome-wide association study (4.5%). In summary, the epigenome is a large contributor to the variation in PPL, and this has the potential to be used to modulate PPL and reduce CVD.

Genome-wide linkage and association analysis of cardiometabolic phenotypes in Hispanic Americans

Linkage studies of complex genetic diseases have been largely replaced by genome-wide association studies, due in part to limited success in complex trait discovery. However, recent interest in rare and low-frequency variants motivates re-examination of family-based methods. In this study, we investigated the performance of two-point linkage analysis for over 1.6 million single-nucleotide polymorphisms (SNPs) combined with single variant association analysis to identify high impact variants, which are both strongly linked and associated with cardiometabolic traits in up to 1414 Hispanics from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of all 50 phenotypes yielded 83 557 000 LOD (logarithm of the odds) scores, with 9214 LOD scores ⩾3.0, 845 ⩾4.0 and 89 ⩾5.0, with a maximal LOD score of 6.49 (rs12956744 in the LAMA1 gene for tumor necrosis factor-α (TNFα) receptor 2). Twenty-seven variants were associated with P<0.005 as well as having an LOD score >4, including variants in the NFIB gene under a linkage peak with TNFα receptor 2 levels on chromosome 9. Linkage regions of interest included a broad peak (31 Mb) on chromosome 1q with acute insulin response (max LOD=5.37). This region was previously documented with type 2 diabetes in family-based studies, providing support for the validity of these results. Overall, we have demonstrated the utility of two-point linkage and association in comprehensive genome-wide array-based SNP genotypes.

Improved Performance of Dynamic Measures of Insulin Response Over Surrogate Indices to Identify Genetic Contributors of Type 2 Diabetes: The GUARDIAN Consortium

Type 2 diabetes (T2D) is a heterogeneous disorder with contributions from peripheral insulin resistance and β-cell dysfunction. For minimization of phenotypic heterogeneity, quantitative intermediate phenotypes characterizing basal glucose homeostasis (insulin resistance and HOMA of insulin resistance [HOMAIR] and of β-cell function [HOMAB]) have shown promise in relatively large samples. We investigated the utility of dynamic measures of glucose homeostasis (insulin sensitivity [SI] and acute insulin response [AIRg]) evaluating T2D-susceptibility variants (n = 57) in Hispanic Americans from the GUARDIAN Consortium (n = 2,560). Basal and dynamic measures were genetically correlated (HOMAB-AIRg: ρG = 0.28-0.73; HOMAIR-SI: ρG = -0.73 to -0.83) with increased heritability for the dynamic measure AIRg Significant association of variants with dynamic measures (P < 8.77 × 10(-4)) was observed. A pattern of superior performance of AIRg was observed for well-established loci including MTNR1B (P = 9.46 × 10(-12)), KCNQ1 (P = 1.35 × 10(-4)), and TCF7L2 (P = 5.10 × 10(-4)) with study-wise statistical significance. Notably, significant association of MTNR1B with AIRg (P < 1.38 × 10(-9)) was observed in a population one-fourteenth the size of the initial discovery cohort. These observations suggest that basal and dynamic measures provide different views and levels of sensitivity to discrete elements of glucose homeostasis. Although more costly to obtain, dynamic measures yield significant results that could be considered physiologically "closer" to causal pathways and provide insight into the discrete mechanisms of action.

Genome-Wide Studies of Type 2 Diabetes and Lipid Traits in Hispanics

Although disproportionately affected by increasing rates of type 2 diabetes and dyslipidemias, Hispanic populations are underrepresented in efforts to understand genetic susceptibility to these disorders. Where research has been undertaken, these populations have provided substantial insight into identification of novel risk-associated genes and have aided in the ability to fine map previously described risk loci. Genome-wide analyses in Hispanic and trans-ethnic populations have resulted in identification of more than 40 replicated or novel genes with significant effects for type 2 diabetes or lipid traits. Initial investigations into rare variant effects have identified new risk-associated variants private to Hispanic populations, and preliminary results suggest metagenomic approaches in Hispanic populations, such as characterizing the gut microbiome, will enable the development of new predictive tools and therapeutic targets for type 2 diabetes. Future genome-wide studies in expanded cohorts of Hispanics are likely to result in new insights into the genetic etiology of metabolic health.

A novel type 2 diabetes risk allele increases the promoter activity of the muscle-specific small ankyrin 1 gene

Genome-wide association studies have identified Ankyrin-1 (ANK1) as a common type 2 diabetes (T2D) susceptibility locus. However, the underlying causal variants and functional mechanisms remain unknown. We screened for 8 tag single nucleotide polymorphisms (SNPs) in ANK1 between 2 case-control studies. Genotype analysis revealed significant associations of 3 SNPs, rs508419 (first identified here), rs515071, and rs516946 with T2D (P < 0.001). These SNPs were in linkage disequilibrium (r² > 0.80); subsequent analysis indicated that the CCC haplotype associated with increased T2D susceptibility (OR 1.447, P < 0.001). Further mapping showed that rs508419 resides in the muscle-specific ANK1 gene promoter. Allele-specific mRNA and protein level measurements confirmed association of the C allele with increased small ANK1 (sAnk1) expression in human skeletal muscle (P = 0.018 and P < 0.001, respectively). Luciferase assays showed increased rs508419-C allele transcriptional activity in murine skeletal muscle C2C12 myoblasts, and electrophoretic mobility-shift assays demonstrated altered rs508419 DNA-protein complex formation. Glucose uptake was decreased with excess sAnk1 expression upon insulin stimulation. Thus, the ANK1 rs508419-C T2D-risk allele alters DNA-protein complex binding leading to increased promoter activity and sAnk1 expression; thus, increased sAnk1 expression in skeletal muscle might contribute to T2D susceptibility.

Genetic Diversity and Association Studies in US Hispanic/Latino Populations: Applications in the Hispanic Community Health Study/Study of Latinos

US Hispanic/Latino individuals are diverse in genetic ancestry, culture, and environmental exposures. Here, we characterized and controlled for this diversity in genome-wide association studies (GWASs) for the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). We simultaneously estimated population-structure principal components (PCs) robust to familial relatedness and pairwise kinship coefficients (KCs) robust to population structure, admixture, and Hardy-Weinberg departures. The PCs revealed substantial genetic differentiation within and among six self-identified background groups (Cuban, Dominican, Puerto Rican, Mexican, and Central and South American). To control for variation among groups, we developed a multi-dimensional clustering method to define a "genetic-analysis group" variable that retains many properties of self-identified background while achieving substantially greater genetic homogeneity within groups and including participants with non-specific self-identification. In GWASs of 22 biomedical traits, we used a linear mixed model (LMM) including pairwise empirical KCs to account for familial relatedness, PCs for ancestry, and genetic-analysis groups for additional group-associated effects. Including the genetic-analysis group as a covariate accounted for significant trait variation in 8 of 22 traits, even after we fit 20 PCs. Additionally, genetic-analysis groups had significant heterogeneity of residual variance for 20 of 22 traits, and modeling this heteroscedasticity within the LMM reduced genomic inflation for 19 traits. Furthermore, fitting an LMM that utilized a genetic-analysis group rather than a self-identified background group achieved higher power to detect previously reported associations. We expect that the methods applied here will be useful in other studies with multiple ethnic groups, admixture, and relatedness.

The Relationship between Native American Ancestry, Body Mass Index and Diabetes Risk among Mexican-Americans

Higher body mass index (BMI) is a well-established risk factor for type 2 diabetes, and rates of obesity and type 2 diabetes are substantially higher among Mexican-Americans relative to non-Hispanic European Americans. Mexican-Americans are genetically diverse, with a highly variable distribution of Native American, European, and African ancestries. Here, we evaluate the role of Native American ancestry on BMI and diabetes risk in a well-defined Mexican-American population. Participants were randomly selected among individuals residing in the Houston area who are enrolled in the Mexican-American Cohort study. Using a custom Illumina GoldenGate Panel, we genotyped DNA from 4,662 cohort participants for 87 Ancestry-Informative Markers. On average, the participants were of 50.2% Native American ancestry, 42.7% European ancestry and 7.1% African ancestry. Using multivariate linear regression, we found BMI and Native American ancestry were inversely correlated; individuals with <20% Native American ancestry were 2.5 times more likely to be severely obese compared to those with >80% Native American ancestry. Furthermore, we demonstrated an interaction between BMI and Native American ancestry in diabetes risk among women; Native American ancestry was a strong risk factor for diabetes only among overweight and obese women (OR = 1.190 for each 10% increase in Native American ancestry). This study offers new insight into the complex relationship between obesity, genetic ancestry, and their respective effects on diabetes risk. Findings from this study may improve the diabetes risk prediction among Mexican-American individuals thereby facilitating targeted prevention strategies.

An Islet-Targeted Genome-Wide Association Scan Identifies Novel Genes Implicated in Cytokine-Mediated Islet Stress in Type 2 Diabetes

Genome-wide association studies in human type 2 diabetes (T2D) have renewed interest in the pancreatic islet as a contributor to T2D risk. Chronic low-grade inflammation resulting from obesity is a risk factor for T2D and a possible trigger of β-cell failure. In this study, microarray data were collected from mouse islets after overnight treatment with cytokines at concentrations consistent with the chronic low-grade inflammation in T2D. Genes with a cytokine-induced change of >2-fold were then examined for associations between single nucleotide polymorphisms and the acute insulin response to glucose (AIRg) using data from the Genetics Underlying Diabetes in Hispanics (GUARDIAN) Consortium. Significant evidence of association was found between AIRg and single nucleotide polymorphisms in Arap3 (5q31.3), F13a1 (6p25.3), Klhl6 (3q27.1), Nid1 (1q42.3), Pamr1 (11p13), Ripk2 (8q21.3), and Steap4 (7q21.12). To assess the potential relevance to islet function, mouse islets were exposed to conditions modeling low-grade inflammation, mitochondrial stress, endoplasmic reticulum (ER) stress, glucotoxicity, and lipotoxicity. RT-PCR revealed that one or more forms of stress significantly altered expression levels of all genes except Arap3. Thapsigargin-induced ER stress up-regulated both Pamr1 and Klhl6. Three genes confirmed microarray predictions of significant cytokine sensitivity: F13a1 was down-regulated 3.3-fold by cytokines, Ripk2 was up-regulated 1.5- to 3-fold by all stressors, and Steap4 was profoundly cytokine sensitive (167-fold up-regulation). Three genes were thus closely associated with low-grade inflammation in murine islets and also with a marker for islet function (AIRg) in a diabetes-prone human population. This islet-targeted genome-wide association scan identified several previously unrecognized candidate genes related to islet dysfunction during the development of T2D.

Genetics of diabetic nephropathy: a long road of discovery

The global prevalence of diabetic nephropathy is rising in parallel with the increasing incidence of diabetes in most countries. Unfortunately, up to 40 % of persons diagnosed with diabetes may develop kidney complications. Diabetic nephropathy is associated with substantially increased risks of cardiovascular disease and premature mortality. An inherited susceptibility to diabetic nephropathy exists, and progress is being made unravelling the genetic basis for nephropathy thanks to international research collaborations, shared biological resources and new analytical approaches. Multiple epidemiological studies have highlighted the clinical heterogeneity of nephropathy and the need for better phenotyping to help define important subgroups for analysis and increase the power of genetic studies. Collaborative genome-wide association studies for nephropathy have reported unique genes, highlighted novel biological pathways and suggested new disease mechanisms, but progress towards clinically relevant risk prediction models for diabetic nephropathy has been slow. This review summarises the current status, recent developments and ongoing challenges elucidating the genetics of diabetic nephropathy.

Connecting SNPs in Diabetes: A Spatial Analysis of Meta-GWAS Loci

Meta-analyses of genome-wide association studies (GWAS) have improved our understanding of the genetic foundations of a number of diseases, including diabetes. However, single nucleotide polymorphisms (SNPs) that are identified by GWAS, especially those that fall outside of gene regions, do not always clearly link to the underlying biology. Despite this, these SNPs have often been validated through re-sequencing efforts as not just tag SNPs, but as causative SNPs, and so must play a role in disease development or progression. In this study, we show how the 3D genome (spatial connections) and trans-expression Quantitative Trait Loci connect diabetes loci from different GWAS meta-analyses, informing the backbone of regulatory networks. Our findings include a three-way functional-spatial connection between the TM6SF2, CTRB1-BCAR1, and CELSR2-PSRC1 loci (rs201189528, rs7202844, and rs7202844, respectively) connected through the KCNIP3 and BCAR1/BCAR3 loci, respectively. These spatial hubs serve as an example of how loci in genes with little biological connection to disease come together to contribute to the diabetes phenotype.