Polygenic Risk for Type 2 Diabetes in African Americans.

African Americans (AAs) have been underrepresented in polygenic risk score (PRS) studies. Here, we integrated genome-wide data from multiple observational studies on type 2 diabetes (T2D), encompassing a total of 101,987 AAs, to train and optimize an AA-focused T2D PRS (PRSAA), using a Bayesian polygenic modeling method. We further tested the score in three independent studies with a total of 7,275 AAs and compared the PRSAA with other published scores. Results show that a 1-SD increase in the PRSAA was associated with 40-60% increase in the odds of T2D (odds ratio [OR] 1.60, 95% CI 1.37-1.88; OR 1.40, 95% CI 1.16-1.70; and OR 1.45, 95% CI 1.30-1.62) across three testing cohorts. These models captured 1.0-2.6% of the variance (R2) in T2D on the liability scale. The positive predictive values for three calculated score thresholds (the top 2%, 5%, and 10%) ranged from 14 to 35%. The PRSAA, in general, performed similarly to existing T2D PRS. The need remains for larger data sets to continue to evaluate the utility of within-ancestry scores in the AA population.

Methods for estimating insulin resistance from untargeted metabolomics data.

CONTEXT: Insulin resistance is associated with multiple complex diseases; however, precise measures of insulin resistance are invasive, expensive, and time-consuming. OBJECTIVE: Develop estimation models for measures of insulin resistance, including insulin sensitivity index (SI) and homeostatic model assessment of insulin resistance (HOMA-IR) from metabolomics data. DESIGN: Insulin Resistance Atherosclerosis Family Study (IRASFS). SETTING: Community based. PARTICIPANTS: Mexican Americans (MA) and African Americans (AA). MAIN OUTCOME: Estimation models for measures of insulin resistance, i.e. SI and HOMA-IR. RESULTS: Least Absolute Shrinkage and Selection Operator (LASSO) and Elastic Net regression were used to build insulin resistance estimation models from 1274 metabolites combined with clinical data, e.g. age, sex, body mass index (BMI). Metabolite data were transformed using three approaches, i.e. inverse normal transformation, standardization, and Box Cox transformation. The analysis was performed in one MA recruitment site (San Luis Valley, Colorado (SLV); N = 450) and tested in another MA recruitment site (San Antonio, Texas (SA); N = 473). In addition, the two MA recruitment sites were combined and estimation models tested in the AA recruitment sample (Los Angeles, California; N = 495). Estimated and empiric SI were correlated in the SA (r2 = 0.77) and AA (r2 = 0.74) testing datasets. Further, estimated and empiric SI were consistently associated with BMI, low-density lipoprotein cholesterol (LDL), and triglycerides. We applied similar approaches to estimate HOMA-IR with similar results. CONCLUSIONS: We have developed a method for estimating insulin resistance with metabolomics data that has the potential for application to a wide range of biomedical studies and conditions.

Metabolomic profiling of glucose homeostasis in African Americans: the Insulin Resistance Atherosclerosis Family Study (IRAS-FS).

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 (SI, insulin sensitivity; AIR, acute insulin response; DI, disposition index; and SG, 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.

Genomics and phenomics of body mass index reveals a complex disease network.

Elevated body mass index (BMI) is heritable and associated with many health conditions that impact morbidity and mortality. The study of the genetic association of BMI across a broad range of common disease conditions offers the opportunity to extend current knowledge regarding the breadth and depth of adiposity-related diseases. We identify 906 (364 novel) and 41 (6 novel) genome-wide significant loci for BMI among participants of European (N~1.1 million) and African (N~100,000) ancestry, respectively. Using a BMI genetic risk score including 2446 variants, 316 diagnoses are associated in the Million Veteran Program, with 96.5% showing increased risk. A co-morbidity network analysis reveals seven disease communities containing multiple interconnected diseases associated with BMI as well as extensive connections across communities. Mendelian randomization analysis confirms numerous phenotypes across a breadth of organ systems, including conditions of the circulatory (heart failure, ischemic heart disease, atrial fibrillation), genitourinary (chronic renal failure), respiratory (respiratory failure, asthma), musculoskeletal and dermatologic systems that are deeply interconnected within and across the disease communities. This work shows that the complex genetic architecture of BMI associates with a broad range of major health conditions, supporting the need for comprehensive approaches to prevent and treat obesity.

Development and validation of a trans-ancestry polygenic risk score for type 2 diabetes in diverse populations.

BACKGROUND: Type 2 diabetes (T2D) is a worldwide scourge caused by both genetic and environmental risk factors that disproportionately afflicts communities of color. Leveraging existing large-scale genome-wide association studies (GWAS), polygenic risk scores (PRS) have shown promise to complement established clinical risk factors and intervention paradigms, and improve early diagnosis and prevention of T2D. However, to date, T2D PRS have been most widely developed and validated in individuals of European descent. Comprehensive assessment of T2D PRS in non-European populations is critical for equitable deployment of PRS to clinical practice that benefits global populations. METHODS: We integrated T2D GWAS in European, African, and East Asian populations to construct a trans-ancestry T2D PRS using a newly developed Bayesian polygenic modeling method, and assessed the prediction accuracy of the PRS in the multi-ethnic Electronic Medical Records and Genomics (eMERGE) study (11,945 cases; 57,694 controls), four Black cohorts (5137 cases; 9657 controls), and the Taiwan Biobank (4570 cases; 84,996 controls). We additionally evaluated a post hoc ancestry adjustment method that can express the polygenic risk on the same scale across ancestrally diverse individuals and facilitate the clinical implementation of the PRS in prospective cohorts. RESULTS: The trans-ancestry PRS was significantly associated with T2D status across the ancestral groups examined. The top 2% of the PRS distribution can identify individuals with an approximately 2.5-4.5-fold of increase in T2D risk, which corresponds to the increased risk of T2D for first-degree relatives. The post hoc ancestry adjustment method eliminated major distributional differences in the PRS across ancestries without compromising its predictive performance. CONCLUSIONS: By integrating T2D GWAS from multiple populations, we developed and validated a trans-ancestry PRS, and demonstrated its potential as a meaningful index of risk among diverse patients in clinical settings. Our efforts represent the first step towards the implementation of the T2D PRS into routine healthcare.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation.

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

Identification of genetic effects underlying type 2 diabetes in South Asian and European populations.

South Asians are at high risk of developing type 2 diabetes (T2D). We carried out a genome-wide association meta-analysis with South Asian T2D cases (n = 16,677) and controls (n = 33,856), followed by combined analyses with Europeans (neff = 231,420). We identify 21 novel genetic loci for significant association with T2D (P = 4.7 × 10-8 to 5.2 × 10-12), to the best of our knowledge at the point of analysis. The loci are enriched for regulatory features, including DNA methylation and gene expression in relevant tissues, and highlight CHMP4B, PDHB, LRIG1 and other genes linked to adiposity and glucose metabolism. A polygenic risk score based on South Asian-derived summary statistics shows ~4-fold higher risk for T2D between the top and bottom quartile. Our results provide further insights into the genetic mechanisms underlying T2D, and highlight the opportunities for discovery from joint analysis of data from across ancestral populations.

Gene Set Enrichment Analsyes Identify Pathways Involved in Genetic Risk for Diabetic Retinopathy.

To identify functionally related genes associated with diabetic retinopathy (DR) risk using gene set enrichment analyses applied to genome-wide association study meta-analyses. METHODS: We analyzed DR GWAS meta-analyses performed on 3246 Europeans and 2611 African Americans with type 2 diabetes. Gene sets relevant to 5 key DR pathophysiology processes were investigated: tissue injury, vascular events, metabolic events and glial dysregulation, neuronal dysfunction, and inflammation. Keywords relevant to these processes were queried in 4 pathway and ontology databases. Two GSEA methods, Meta-Analysis Gene set Enrichment of variaNT Associations (MAGENTA) and Multi-marker Analysis of GenoMic Annotation (MAGMA), were used. Gene sets were defined to be enriched for gene associations with DR if the P value corrected for multiple testing (Pcorr) was <.05. RESULTS: Five gene sets were significantly enriched for numerous modest genetic associations with DR in one method (MAGENTA or MAGMA) and also at least nominally significant (uncorrected P < .05) in the other method. These pathways were regulation of the lipid catabolic process (2-fold enrichment, Pcorr = .014); nitric oxide biosynthesis (1.92-fold enrichment, Pcorr = .022); lipid digestion, mobilization, and transport (1.6-fold enrichment, P = .032); apoptosis (1.53-fold enrichment, P = .041); and retinal ganglion cell degeneration (2-fold enrichment, Pcorr = .049). The interferon gamma (IFNG) gene, previously implicated in DR by protein-protein interactions in our GWAS, was among the top ranked genes in the nitric oxide pathway (best variant P = .0001). CONCLUSIONS: These GSEA indicate that variants in genes involved in oxidative stress, lipid transport and catabolism, and cell degeneration are enriched for genes associated with DR risk. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Genetic discovery and risk characterization in type 2 diabetes across diverse populations.

Genomic discovery and characterization of risk loci for type 2 diabetes (T2D) have been conducted primarily in individuals of European ancestry. We conducted a multiethnic genome-wide association study of T2D among 53,102 cases and 193,679 control subjects from African, Hispanic, Asian, Native Hawaiian, and European population groups in the Population Architecture Genomics and Epidemiology (PAGE) and Diabetes Genetics Replication and Meta-analysis (DIAGRAM) Consortia. In individuals of African ancestry, we discovered a risk variant in the TGFB1 gene (rs11466334, risk allele frequency (RAF) = 6.8%, odds ratio [OR] = 1.27, p = 2.06 × 10-8), which replicated in independent studies of African ancestry (p = 6.26 × 10-23). We identified a multiethnic risk variant in the BACE2 gene (rs13052926, RAF = 14.1%, OR = 1.08, p = 5.75 × 10-9), which also replicated in independent studies (p = 3.45 × 10-4). We also observed a significant difference in the performance of a multiethnic genetic risk score (GRS) across population groups (pheterogeneity = 3.85 × 10-20). Comparing individuals in the top GRS risk category (40%-60%), the OR was highest in Asians (OR = 3.08) and European (OR = 2.94) ancestry populations, followed by Hispanic (OR = 2.39), Native Hawaiian (OR = 2.02), and African ancestry (OR = 1.57) populations. These findings underscore the importance of genetic discovery and risk characterization in diverse populations and the urgent need to further increase representation of non-European ancestry individuals in genetics research to improve genetic-based risk prediction across populations.

Metabolomic architecture of obesity implicates metabolonic lactone sulfate in cardiometabolic disease.

OBJECTIVE: Identify and characterize circulating metabolite profiles associated with adiposity to inform precision medicine. METHODS: Untargeted plasma metabolomic profiles in the Insulin Resistance Atherosclerosis Family Study (IRASFS) Mexican American cohort (n = 1108) were analyzed for association with anthropometric (body mass index, BMI; waist circumference, WC; waist-to-hip ratio, WHR) and computed tomography measures (visceral adipose tissue, VAT; subcutaneous adipose tissue, SAT; visceral-to-subcutaneous ratio, VSR) of adiposity. Genetic data, inclusive of genome-wide array-based genotyping, whole exome sequencing (WES) and whole genome sequencing (WGS), were evaluated to identify the genetic contributors. Phenotypic and genetic association signals were replicated across ancestries. Transcriptomic data were analyzed to explore the relationship between genetic and metabolomic data. RESULTS: A partially characterized metabolite, tentatively named metabolonic lactone sulfate (X-12063), was consistently associated with BMI, WC, WHR, VAT, and SAT in IRASFS Mexican Americans (PMA <2.02 × 10-27). Trait associations were replicated in IRASFS African Americans (PAA < 1.12 × 10-07). Expanded analyses revealed associations with multiple phenotypic measures of cardiometabolic health, e.g. insulin sensitivity (SI), triglycerides (TG), diastolic blood pressure (DBP) and plasminogen activator inhibitor-1 (PAI-1) in both ancestries. Metabolonic lactone sulfate levels were heritable (h2 > 0.47), and a significant genetic signal at the ZSCAN25/CYP3A5 locus (PMA = 9.00 × 10-41, PAA = 2.31 × 10-10) was observed, highlighting a putative functional variant (rs776746, CYP3A5∗3). Transcriptomic analysis in the African American Genetics of Metabolism and Expression (AAGMEx) cohort supported the association of CYP3A5 with metabolonic lactone sulfate levels (PFDR = 6.64 × 10-07). CONCLUSIONS: Variant rs776746 is associated with a decrease in the transcript levels of CYP3A5, which in turn is associated with increased metabolonic lactone sulfate levels and poor cardiometabolic health.

Multiethnic Genome-Wide Association Study of Subclinical Atherosclerosis in Individuals With Type 2 Diabetes.

BACKGROUND: Coronary artery calcification (CAC) and carotid artery intima-media thickness (cIMT) are measures of subclinical atherosclerosis in asymptomatic individuals and strong risk factors for cardiovascular disease. Type 2 diabetes (T2D) is an independent cardiovascular disease risk factor that accelerates atherosclerosis. METHODS: We performed meta-analyses of genome-wide association studies in up to 2500 T2D individuals of European ancestry (EA) and 1590 T2D individuals of African ancestry with or without exclusion of prevalent cardiovascular disease, for CAC measured by cardiac computed tomography, and 3608 individuals of EA and 838 individuals of African ancestry with T2D for cIMT measured by ultrasonography within the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium. RESULTS: We replicated 2 loci (rs9369640 and rs9349379 near PHACTR1 and rs10757278 near CDKN2B) for CAC and one locus for cIMT (rs7412 and rs445925 near APOE-APOC1) that were previously reported in the general EA populations. We identified one novel CAC locus (rs8000449 near CSNK1A1L/LINC00547/POSTN at 13q13.3) at P=2.0×10-8 in EA. No additional loci were identified with the meta-analyses of EA and African ancestry. The expression quantitative trait loci analysis with nearby expressed genes derived from arterial wall and metabolic tissues from the Genotype-Tissue Expression project pinpoints POSTN, encoding a matricellular protein involved in bone formation and bone matrix organization, as the potential candidate gene at this locus. In addition, we found significant associations (P<3.1×10-4) for 3 previously reported coronary artery disease loci for these subclinical atherosclerotic phenotypes (rs2891168 near CDKN2B-AS1 and rs11170820 near FLJ12825 for CAC, and rs7412 near APOE for cIMT). CONCLUSIONS: Our results provide potential biological mechanisms that could link CAC and cIMT to increased cardiovascular disease risk in individuals with T2D.

Genome-wide association study of vitamin D concentrations and bone mineral density in the African American-Diabetes Heart Study.

Relative to European Americans, African Americans have lower 25-hydroxyvitamin D (25OHD) and vitamin D binding protein (VDBP) concentrations, higher 1,25-dihydroxyvitamin D (1,25(OH)2D3) concentrations and bone mineral density (BMD), and paradoxically reduced burdens of calcified atherosclerotic plaque (subclinical atherosclerosis). To identify genetic factors contributing to vitamin D and BMD measures, association analysis of >14M variants was conducted in a maximum of 697 African American-Diabetes Heart Study participants with type 2 diabetes (T2D). The most significant association signals were detected for VDBP on chromosome 4; variants rs7041 (ß = 0.44, SE = 0.019, P = 9.4x10-86) and rs4588 (ß = 0.17, SE = 0.021, P = 3.5x10-08) in the group-specific component (vitamin D binding protein) gene (GC). These variants were found to be independently associated. In addition, rs7041 was also associated with bioavailable vitamin D (BAVD; ß = 0.16, SE = 0.02, P = 3.3x10-19). Six rare variants were significantly associated with 25OHD, including a non-synonymous variant in HSPG2 (rs116788687; ß = -1.07, SE = 0.17, P = 2.2x10-10) and an intronic variant in TNIK (rs143555701; ß = -1.01, SE = 0.18, P = 9.0x10-10), both biologically related to bone development. Variants associated with 25OHD failed to replicate in African Americans from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of vitamin D metabolism and bone mineral density phenotypes in an African American population enriched for T2D could provide insight into ethnic specific differences in vitamin D metabolism and bone mineral density.

Development of genome-wide polygenic risk scores for lipid traits and clinical applications for dyslipidemia, subclinical atherosclerosis, and diabetes cardiovascular complications among East Asians.

BACKGROUND: The clinical utility of personal genomic information in identifying individuals at increased risks for dyslipidemia and cardiovascular diseases remains unclear. METHODS: We used data from Biobank Japan (n = 70,657-128,305) and developed novel East Asian-specific genome-wide polygenic risk scores (PRSs) for four lipid traits. We validated (n = 4271) and subsequently tested associations of these scores with 3-year lipid changes in adolescents (n = 620), carotid intima-media thickness (cIMT) in adult women (n = 781), dyslipidemia (n = 7723), and coronary heart disease (CHD) (n = 2374 cases and 6246 controls) in type 2 diabetes (T2D) patients. RESULTS: Our PRSs aggregating 84-549 genetic variants (0.251 < correlation coefficients (r) < 0.272) had comparably stronger association with lipid variations than the typical PRSs derived based on the genome-wide significant variants (0.089 < r < 0.240). Our PRSs were robustly associated with their corresponding lipid levels (7.5 × 10- 103 < P < 1.3 × 10- 75) and 3-year lipid changes (1.4 × 10- 6 < P < 0.0130) which started to emerge in childhood and adolescence. With the adjustments for principal components (PCs), sex, age, and body mass index, there was an elevation of 5.3% in TC (ß ± SE = 0.052 ± 0.002), 11.7% in TG (ß ± SE = 0.111 ± 0.006), 5.8% in HDL-C (ß ± SE = 0.057 ± 0.003), and 8.4% in LDL-C (ß ± SE = 0.081 ± 0.004) per one standard deviation increase in the corresponding PRS. However, their predictive power was attenuated in T2D patients (0.183 < r < 0.231). When we included each PRS (for TC, TG, and LDL-C) in addition to the clinical factors and PCs, the AUC for dyslipidemia was significantly increased by 0.032-0.057 in the general population (7.5 × 10- 3 < P < 0.0400) and 0.029-0.069 in T2D patients (2.1 × 10- 10 < P < 0.0428). Moreover, the quintile of TC-related PRS was moderately associated with cIMT in adult women (ß ± SE = 0.011 ± 0.005, Ptrend = 0.0182). Independent of conventional risk factors, the quintile of PRSs for TC [OR (95% CI) = 1.07 (1.03-1.11)], TG [OR (95% CI) = 1.05 (1.01-1.09)], and LDL-C [OR (95% CI) = 1.05 (1.01-1.09)] were significantly associated with increased risk of CHD in T2D patients (4.8 × 10- 4 < P < 0.0197). Further adjustment for baseline lipid drug use notably attenuated the CHD association. CONCLUSIONS: The PRSs derived and validated here highlight the potential for early genomic screening and personalized risk assessment for cardiovascular disease.

Identification of type 2 diabetes loci in 433,540 East Asian individuals.

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genome-wide association study identifies novel loci for type 2 diabetes-attributed end-stage kidney disease in African Americans.

BACKGROUND: End-stage kidney disease (ESKD) is a significant public health concern disproportionately affecting African Americans (AAs). Type 2 diabetes (T2D) is the leading cause of ESKD in the USA, and efforts to uncover genetic susceptibility to diabetic kidney disease (DKD) have had limited success. A prior genome-wide association study (GWAS) in AAs with T2D-ESKD was expanded with additional AA cases and controls and genotypes imputed to the higher density 1000 Genomes reference panel. The discovery analysis included 3432 T2D-ESKD cases and 6977 non-diabetic non-nephropathy controls (N = 10,409), followed by a discrimination analysis in 2756 T2D non-nephropathy controls to exclude T2D-associated variants. RESULTS: Six independent variants located in or near RND3/RBM43, SLITRK3, ENPP7, GNG7, and APOL1 achieved genome-wide significant association (P < 5 × 10-8) with T2D-ESKD. Following extension analyses in 1910 non-diabetic ESKD cases and 908 non-diabetic non-nephropathy controls, a meta-analysis of 5342 AA all-cause ESKD cases and 6977 AA non-diabetic non-nephropathy controls revealed an additional novel all-cause ESKD locus at EFNB2 (rs77113398; P = 9.84 × 10-9; OR = 1.94). Exclusion of APOL1 renal-risk genotype carriers identified two additional genome-wide significant T2D-ESKD-associated loci at GRAMD3 and MGAT4C. A second variant at GNG7 (rs373971520; P = 2.17 × 10-8, OR = 1.46) remained associated with all-cause ESKD in the APOL1-negative analysis. CONCLUSIONS: Findings provide further evidence for genetic factors associated with advanced kidney disease in AAs with T2D.