Introduction

African Americans are at increased risk for type 2 diabetes.

Objectives

This work aimed to examine metabolomic signature of glucose homeostasis in African Americans.

Methods

We used an untargeted liquid chromatography-mass spectrometry metabolomic approach to comprehensively profile 727 plasma metabolites among 571 African Americans from the Insulin Resistance Atherosclerosis Family Study (IRAS-FS) and investigate the associations between these metabolites and both the dynamic (S_I, insulin sensitivity; AIR, acute insulin response; DI, disposition index; and S_G, glucose effectiveness) and basal (HOMA-IR and HOMA-B) measures of glucose homeostasis using univariate and regularized regression models. We also compared the results with our previous findings in the IRAS-FS Mexican Americans.

Results

We confirmed increased plasma metabolite levels of branched-chain amino acids and their metabolic derivatives, 2-aminoadipate, 2-hydroxybutyrate, glutamate, arginine and its metabolic derivatives, carbohydrate metabolites, and medium- and long-chain fatty acids were associated with insulin resistance, while increased plasma metabolite levels in the glycine, serine and threonine metabolic pathway were associated with insulin sensitivity. We also observed a differential ancestral effect of glutamate on glucose homeostasis with significantly stronger effects observed in African Americans than those previously observed in Mexican Americans.

Conclusion

We extended the observations that metabolites are useful biomarkers in the identification of prediabetes in individuals at risk of type 2 diabetes in African Americans. We revealed, for the first time, differential ancestral effect of certain metabolites (i.e., glutamate) on glucose homeostasis traits. Our study highlights the need for additional comprehensive metabolomic studies in well-characterized multiethnic cohorts.

Linkage analysis of complex traits has had limited success in identifying trait-influencing loci. Recently, coding variants have been implicated as the basis for some biomedical associations. We tested whether coding variants are the basis for linkage peaks of complex traits in 42 African-American (n = 596) and 90 Hispanic (n = 1,414) families in the Insulin Resistance Atherosclerosis Family Study (IRASFS) using Illumina HumanExome Beadchips. A total of 92,157 variants in African Americans (34%) and 81,559 (31%) in Hispanics were polymorphic and tested using two-point linkage and association analyses with 37 cardiometabolic phenotypes. In African Americans 77 LOD scores greater than 3 were observed. The highest LOD score was 4.91 with the APOE SNP rs7412 (MAF = 0.13) with plasma apolipoprotein B (ApoB). This SNP was associated with ApoB (P-value = 4 × 10(-19)) and accounted for 16.2% of the variance in African Americans. In Hispanic families, 104 LOD scores were greater than 3. The strongest evidence of linkage (LOD = 4.29) was with rs5882 (MAF = 0.46) in CETP with HDL. CETP variants were strongly associated with HDL (0.00049 < P-value <4.6 × 10(-12)), accounting for up to 4.5% of the variance. These loci have previously been shown to have effects on the biomedical traits evaluated here. Thus, evidence of strong linkage in this genome wide survey of primarily coding variants was uncommon. Loci with strong evidence of linkage was characterized by large contributions to the variance, and, in these cases, are common variants. Less compelling evidence of linkage and association was observed with additional loci that may require larger family sets to confirm.

Objective

Adiponectin is found in human serum in three groups of multimers (high molecular weight [HMW], medium molecular weight [MMW], and low molecular weight [LMW]). Two ethnic-specific variants in ADIPOQ, G45R (Hispanic-Americans) and R55C (African-Americans), were previously reported. Although carriers of both variants had mean adiponectin levels ≤ 20% of those of noncarriers, they were not clinically different from noncarriers. To compare carriers of both variants and noncarriers, relative quantification of adiponectin isoforms to total adiponectin was performed on serum samples.

Methods

The multimeric patterns of serum adiponectin in G45R carriers (n = 23), R55C carriers (n = 3), and Hispanic- and African-American noncarriers (n = 84 and 44, respectively) from the Insulin Resistance Atherosclerosis Family Study were explored using native Western blotting and densitometry.

Results

Serum samples from carriers showed an absence of the HMW isoform and a marked reduction in the MMW isoform but an approximate twofold increase in the amount of the LMW isoform. Thus, individuals making only LMW adiponectin are metabolically normal.

Conclusions

The results contrast with the proposed biological importance of the HMW multimer. This suggests that the LMW isoform may functionally compensate for some of the loss or reduction of the higher-order multimers in carriers of the G45R and R55C mutations.

Objective

This study aimed to explore the genetic mechanisms of regional fat deposition, which is a strong risk factor for metabolic diseases beyond total adiposity.

Methods

A genome-wide association study of 7,757,139 single-nucleotide polymorphisms (SNPs) in 983 Mexican Americans (n_male  = 403; n_female  = 580) from the Insulin Resistance Atherosclerosis Family Study was performed. Association analyses were performed with and without sex stratification for subcutaneous adipose tissue, visceral adipose tissue (VAT), and visceral-subcutaneous ratio (VSR) obtained from computed tomography.

Results

The strongest signal identified was SNP rs2185405 (minor allele frequencies [MAF] = 40%; P_VAT = 1.98 × 10^-8 ) with VAT. It is an intronic variant of the GLIS family zinc finger 3 gene (GLIS3). In addition, SNP rs12657394 (MAF = 19%) was associated with VAT in males (P_male = 2.39×10^-8 ; P_female = 2.5 × 10^-3 ). It is located intronically in the serum response factor binding protein 1 gene (SRFBP1). On average, male carriers of the variant had 24.6 cm² increased VAT compared with noncarriers. Subsequently, genome-wide SNP-sex interaction analysis was performed. SNP rs10913233 (MAF = 14%; P_int = 3.07 × 10^-8 ) in PAPPA2 and rs10923724 (MAF = 38%; P_int = 2.89 × 10^-8 ) upstream of TBX15 were strongly associated with the interaction effect for VSR.

Conclusions

Six loci were identified with genome-wide significant associations with fat deposition and interactive effects. These results provided genetic evidence for a differential basis of fat deposition between genders.

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.

Objective

To identify and investigate the susceptibility genes of Kashin-Beck disease (KBD) in Chinese population.

Methods

Whole-exome capturing and sequencing technology was used for the detection of genetic variations in 19 individuals from six families with high incidence of KBD. A total of 44 polymorphisms from 41 genes were genotyped from a total of 144 cases and 144 controls by using MassARRAY under the standard protocol from Sequenom. Association was applied on the data by using PLINK1.07.

Results

In the sequencing stage, each sample showed approximately 70-fold coverage, thus covering more than 99% of the target regions. Among the single nucleotide polymorphisms (SNPs) used in the transmission disequilibrium test, 108 had a p-value of <0.01, whereas 1056 had a p-value of <0.05. Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway analysis indicates that these SNPs focus on three major pathways: regulation of actin cytoskeleton, focal adhesion, and metabolic pathways. In the validation stage, single locus effects revealed that two of these polymorphisms (rs7745040 and rs9275295) in the human leukocyte antigen (HLA)-DRB1 gene and one polymorphism (rs9473132) in CD2-associated protein (CD2AP) gene have a significant statistical association with KBD.

Conclusions

HLA-DRB1 and CD2AP gene were identified to be among the susceptibility genes of KBD, thus supporting the role of the autoimmune response in KBD and the possibility of shared etiology between osteoarthritis, rheumatoid arthritis, and KBD.

Genome-wide association studies have identified numerous variants associated with lipid levels; yet, the majority are located in non-coding regions with unclear mechanisms. In the Insulin Resistance Atherosclerosis Family Study (IRASFS), heritability estimates suggest a strong genetic basis: low-density lipoprotein (LDL, h² = 0.50), high-density lipoprotein (HDL, h² = 0.57), total cholesterol (TC, h² = 0.53), and triglyceride (TG, h² = 0.42) levels. Exome sequencing of 1,205 Mexican Americans (90 pedigrees) from the IRASFS identified 548,889 variants and association and linkage analyses with lipid levels were performed. One genome-wide significant signal was detected in APOA5 with TG (rs651821, P_TG = 3.67 × 10^-10, LOD_TG = 2.36, MAF = 14.2%). In addition, two correlated SNPs (r² = 1.0) rs189547099 (P_TG = 6.31 × 10^-08, LOD_TG = 3.13, MAF = 0.50%) and chr4:157997598 (P_TG = 6.31 × 10^-08, LOD_TG = 3.13, MAF = 0.50%) reached exome-wide significance (P < 9.11 × 10^-08). rs189547099 is an intronic SNP in FNIP2 and SNP chr4:157997598 is intronic in GLRB. Linkage analysis revealed 46 SNPs with a LOD > 3 with the strongest signal at rs1141070 (LOD_LDL = 4.30, P_LDL = 0.33, MAF = 21.6%) in DFFB. A total of 53 nominally associated variants (P < 5.00 × 10^-05, MAF ≥ 1.0%) were selected for replication in six Mexican-American cohorts (N = 3,280). The strongest signal observed was a synonymous variant (rs1160983, P_LDL = 4.44 × 10^-17, MAF = 2.7%) in TOMM40. Beyond primary findings, previously reported lipid loci were fine-mapped using exome sequencing in IRASFS. These results support that exome sequencing complements and extends insights into the genetics of lipid levels.

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.

We previously identified a low-frequency (1.1 %) coding variant (G45R; rs200573126) in the adiponectin gene (ADIPOQ) which was the basis for a multipoint microsatellite linkage signal (LOD = 8.2) for plasma adiponectin levels in Hispanic families. We have empirically evaluated the ability of data from targeted common variants, exome chip genotyping, and genome-wide association study data to detect linkage and association to adiponectin protein levels at this locus. Simple two-point linkage and association analyses were performed in 88 Hispanic families (1,150 individuals) using 10,958 SNPs on chromosome 3. Approaches were compared for their ability to map the functional variant, G45R, which was strongly linked (two-point LOD = 20.98) and powerfully associated (p value = 8.1 × 10(-50)). Over 450 SNPs within a broad 61 Mb interval around rs200573126 showed nominal evidence of linkage (LOD > 3) but only four other SNPs in this region were associated with p values < 1.0 × 10(-4). When G45R was accounted for, the maximum LOD score across the interval dropped to 4.39 and the best p value was 1.1 × 10(-5). Linked and/or associated variants ranged in frequency (0.0018-0.50) and type (coding, non-coding) and had little detectable linkage disequilibrium with rs200573126 (r (2) < 0.20). In addition, the two-point linkage approach empirically outperformed multipoint microsatellite and multipoint SNP analysis. In the absence of data for rs200573126, family-based linkage analysis using a moderately dense SNP dataset, including both common and low-frequency variants, resulted in stronger evidence for an adiponectin locus than association data alone. Thus, linkage analysis can be a useful tool to facilitate identification of high-impact genetic variants.

Although type 2 diabetes (T2D) results from metabolic defects in insulin secretion and insulin sensitivity, most of the genetic risk loci identified to date relates to insulin secretion. We reported that T2D loci influencing insulin sensitivity may be identified through interactions with insulin secretion loci, thereby leading to T2D. Here, we hypothesize that joint testing of variant main effects and interaction effects with an insulin secretion locus increases power to identify genetic interactions leading to T2D. We tested this hypothesis with an intronic MTNR1B SNP, rs10830963, which is associated with acute insulin response to glucose, a dynamic measure of insulin secretion. rs10830963 was tested for interaction and joint (main + interaction) effects with genome-wide data in African Americans (2,452 cases and 3,772 controls) from five cohorts. Genome-wide genotype data (Affymetrix Human Genome 6.0 array) was imputed to a 1000 Genomes Project reference panel. T2D risk was modeled using logistic regression with rs10830963 dosage, age, sex, and principal component as predictors. Joint effects were captured using the Kraft two degrees of freedom test. Genome-wide significant (P < 5 × 10^-8 ) interaction with MTNR1B and joint effects were detected for CMIP intronic SNP rs17197883 (P_interaction  = 1.43 × 10^-8 ; P_joint  = 4.70 × 10^-8 ). CMIP variants have been nominally associated with T2D, fasting glucose, and adiponectin in individuals of East Asian ancestry, with high-density lipoprotein, and with waist-to-hip ratio adjusted for body mass index in Europeans. These data support the hypothesis that additional genetic factors contributing to T2D risk, including insulin sensitivity loci, can be identified through interactions with insulin secretion loci.