Hemoglobin A1c Genetics and Disparities in Risk of Diabetic Retinopathy in Individuals of Genetically Inferred African American/African British and European Ancestries.

OBJECTIVE: Individuals with diabetes who carry genetic variants that lower hemoglobin A1c (HbA1c) independently of glycemia may have higher real, but undetected, hyperglycemia compared with those without these variants despite achieving similar HbA1c targets, potentially placing them at greater risk for diabetes-related complications. We sought to determine whether these genetic variants, aggregated in a polygenic score, and the large-effect African ancestry-specific missense variant in G6PD (rs1050828) that lower HbA1c were associated with higher retinopathy risk. RESEARCH DESIGN AND METHODS: Using data from 29,828 type 2 diabetes cases of genetically inferred African American/African British and European ancestries, we calculated ancestry-specific nonglycemic HbA1c polygenic scores (ngA1cPS) composed of 122 variants associated with HbA1c at genome-wide significance, but not with glucose. We tested the association of the ngA1cPS and the G6PD variant with retinopathy, adjusting for measured HbA1c and retinopathy risk factors. RESULTS: Participants in the bottom quintile of the ngA1cPS showed between 20% and 50% higher retinopathy prevalence, compared with those above this quintile, despite similar levels of measured HbA1c. The adjusted meta-analytic odds ratio for the bottom quintile was 1.31 (95% CI 1.0, 1.73; P = 0.05) in African ancestry and 1.31 (95% CI 1.15, 1.50; P = 6.5 × 10-5) in European ancestry. Among individuals of African ancestry with HbA1c below 7% units, retinopathy prevalence was higher in individuals below, compared with above, the 50th percentile of the ngA1cPS regardless of sex or G6PD carrier status. CONCLUSIONS: Genetic effects need to be considered to personalize HbA1c targets and improve outcomes of people with diabetes from diverse ancestries.

Multi-omics characterization of type 2 diabetes associated genetic variation.

Discerning the mechanisms driving type 2 diabetes (T2D) pathophysiology from genome-wide association studies (GWAS) remains a challenge. To this end, we integrated omics information from 16 multi-tissue and multi-ancestry expression, protein, and metabolite quantitative trait loci (QTL) studies and 46 multi-ancestry GWAS for T2D-related traits with the largest, most ancestrally diverse T2D GWAS to date. Of the 1,289 T2D GWAS index variants, 716 (56%) demonstrated strong evidence of colocalization with a molecular or T2D-related trait, implicating 657 cis-effector genes, 1,691 distal-effector genes, 731 metabolites, and 43 T2D-related traits. We identified 773 of these cis- and distal-effector genes using either expression QTL data from understudied ancestry groups or inclusion of T2D index variants enriched in underrepresented populations, emphasizing the value of increasing population diversity in functional mapping. Linking these variants, genes, metabolites, and traits into a network, we elucidated mechanisms through which T2D-associated variation may impact disease risk. Finally, we showed that drugs targeting effector proteins were enriched in those approved to treat T2D, highlighting the potential of these results to prioritize drug targets for T2D. These results represent a leap in the molecular characterization of T2D-associated genetic variation and will aid in translating genetic findings into novel therapeutic strategies.

Adaptive selection at G6PD and disparities in diabetes complications.

Diabetes complications occur at higher rates in individuals of African ancestry. Glucose-6-phosphate dehydrogenase deficiency (G6PDdef), common in some African populations, confers malaria resistance, and reduces hemoglobin A1c (HbA1c) levels by shortening erythrocyte lifespan. In a combined-ancestry genome-wide association study of diabetic retinopathy, we identified nine loci including a G6PDdef causal variant, rs1050828 -T (Val98Met), which was also associated with increased risk of other diabetes complications. The effect of rs1050828 -T on retinopathy was fully mediated by glucose levels. In the years preceding diabetes diagnosis and insulin prescription, glucose levels were significantly higher and HbA1c significantly lower in those with versus without G6PDdef. In the Action to Control Cardiovascular Risk in Diabetes (ACCORD) trial, participants with G6PDdef had significantly higher hazards of incident retinopathy and neuropathy. At the same HbA1c levels, G6PDdef participants in both ACCORD and the Million Veteran Program had significantly increased risk of retinopathy. We estimate that 12% and 9% of diabetic retinopathy and neuropathy cases, respectively, in participants of African ancestry are due to this exposure. Across continentally defined ancestral populations, the differences in frequency of rs1050828 -T and other G6PDdef alleles contribute to disparities in diabetes complications. Diabetes management guided by glucose or potentially genotype-adjusted HbA1c levels could lead to more timely diagnoses and appropriate intensification of therapy, decreasing the risk of diabetes complications in patients with G6PDdef alleles.

Increased Genetic Risk for ß-Cell Failure Is Associated With ß-Cell Function Decline in People With Prediabetes.

Partitioned polygenic scores (pPS) have been developed to capture pathophysiologic processes underlying type 2 diabetes (T2D). We investigated the association of T2D pPS with diabetes-related traits and T2D incidence in the Diabetes Prevention Program. We generated five T2D pPS (ß-cell, proinsulin, liver/lipid, obesity, lipodystrophy) in 2,647 participants randomized to intensive lifestyle, metformin, or placebo arms. Associations were tested with general linear models and Cox regression with adjustment for age, sex, and principal components. Sensitivity analyses included adjustment for BMI. Higher ß-cell pPS was associated with lower insulinogenic index and corrected insulin response at 1-year follow-up with adjustment for baseline measures (effect per pPS SD -0.04, P = 9.6 × 10-7, and -8.45 µU/mg, P = 5.6 × 10-6, respectively) and with increased diabetes incidence with adjustment for BMI at nominal significance (hazard ratio 1.10 per SD, P = 0.035). The liver/lipid pPS was associated with reduced 1-year baseline-adjusted triglyceride levels (effect per SD -4.37, P = 0.001). There was no significant interaction between T2D pPS and randomized groups. The remaining pPS were associated with baseline measures only. We conclude that despite interventions for diabetes prevention, participants with a high genetic burden of the ß-cell cluster pPS had worsening in measures of ß-cell function.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes.

OBJECTIVE: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD. RESULTS: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Polygenic scores for longitudinal prediction of incident type 2 diabetes in an ancestrally and medically diverse primary care physician network: a patient cohort study.

BACKGROUND: The clinical utility of genetic information for type 2 diabetes (T2D) prediction with polygenic scores (PGS) in ancestrally diverse, real-world US healthcare systems is unclear, especially for those at low clinical phenotypic risk for T2D. METHODS: We tested the association of PGS with T2D incidence in patients followed within a primary care practice network over 16 years in four hypothetical scenarios that varied by clinical data availability (N = 14,712): (1) age and sex; (2) age, sex, body mass index (BMI), systolic blood pressure, and family history of T2D; (3) all variables in (2) and random glucose; and (4) all variables in (3), HDL, total cholesterol, and triglycerides, combined in a clinical risk score (CRS). To determine whether genetic effects differed by baseline clinical risk, we tested for interaction with the CRS. RESULTS: PGS was associated with incident T2D in all models. Adjusting for age and sex only, the Hazard Ratio (HR) per PGS standard deviation (SD) was 1.76 (95% CI 1.68, 1.84) and the HR of top 5% of PGS vs interquartile range (IQR) was 2.80 (2.39, 3.28). Adjusting for the CRS, the HR per SD was 1.48 (1.40, 1.57) and HR of the top 5% of PGS vs IQR was 2.09 (1.72, 2.55). Genetic effects differed by baseline clinical risk ((PGS-CRS interaction p = 0.05; CRS below the median: HR 1.60 (1.43, 1.79); CRS above the median: HR 1.45 (1.35, 1.55)). CONCLUSIONS: Genetic information can help identify high-risk patients even among those perceived to be low risk in a clinical evaluation.

Multi-ancestry polygenic mechanisms of type 2 diabetes.

Type 2 diabetes (T2D) is a multifactorial disease with substantial genetic risk, for which the underlying biological mechanisms are not fully understood. In this study, we identified multi-ancestry T2D genetic clusters by analyzing genetic data from diverse populations in 37 published T2D genome-wide association studies representing more than 1.4 million individuals. We implemented soft clustering with 650 T2D-associated genetic variants and 110 T2D-related traits, capturing known and novel T2D clusters with distinct cardiometabolic trait associations across two independent biobanks representing diverse genetic ancestral populations (African, n = 21,906; Admixed American, n = 14,410; East Asian, n =2,422; European, n = 90,093; and South Asian, n = 1,262). The 12 genetic clusters were enriched for specific single-cell regulatory regions. Several of the polygenic scores derived from the clusters differed in distribution among ancestry groups, including a significantly higher proportion of lipodystrophy-related polygenic risk in East Asian ancestry. T2D risk was equivalent at a body mass index (BMI) of 30 kg m-2 in the European subpopulation and 24.2 (22.9-25.5) kg m-2 in the East Asian subpopulation; after adjusting for cluster-specific genetic risk, the equivalent BMI threshold increased to 28.5 (27.1-30.0) kg m-2 in the East Asian group. Thus, these multi-ancestry T2D genetic clusters encompass a broader range of biological mechanisms and provide preliminary insights to explain ancestry-associated differences in T2D risk profiles.

Assessing the genetic contribution of cumulative behavioral factors associated with longitudinal type 2 diabetes risk highlights adiposity and the brain-metabolic axis.

While genetic factors, behavior, and environmental exposures form a complex web of interrelated associations in type 2 diabetes (T2D), their interaction is poorly understood. Here, using data from ~500K participants of the UK Biobank, we identify the genetic determinants of a "polyexposure risk score" (PXS) a new risk factor that consists of an accumulation of 25 associated individual-level behaviors and environmental risk factors that predict longitudinal T2D incidence. PXS-T2D had a non-zero heritability (h2 = 0.18) extensive shared genetic architecture with established clinical and biological determinants of T2D, most prominently with body mass index (genetic correlation [rg] = 0.57) and Homeostatic Model Assessment for Insulin Resistance (rg = 0.51). Genetic loci associated with PXS-T2D were enriched for expression in the brain. Biobank scale data with genetic information illuminates how complex and cumulative exposures and behaviors as a whole impact T2D risk but whose biology have been elusive in genome-wide studies of T2D.

Genetic architecture and biology of youth-onset type 2 diabetes.

The prevalence of youth-onset type 2 diabetes (T2D) and childhood obesity has been rising steadily1, producing a growing public health concern1 that disproportionately affects minority groups2. The genetic basis of youth-onset T2D and its relationship to other forms of diabetes are unclear3. Here we report a detailed genetic characterization of youth-onset T2D by analysing exome sequences and common variant associations for 3,005 individuals with youth-onset T2D and 9,777 adult control participants matched for ancestry, including both males and females. We identify monogenic diabetes variants in 2.4% of individuals and three exome-wide significant (P < 2.6 × 10-6) gene-level associations (HNF1A, MC4R, ATXN2L). Furthermore, we report rare variant association enrichments within 25 gene sets related to obesity, monogenic diabetes and ß-cell function. Many youth-onset T2D associations are shared with adult-onset T2D, but genetic risk factors of all frequencies-and rare variants in particular-are enriched within youth-onset T2D cases (5.0-fold increase in the rare variant and 3.4-fold increase in common variant genetic liability relative to adult-onset cases). The clinical presentation of participants with youth-onset T2D is influenced in part by the frequency of genetic risk factors within each individual. These findings portray youth-onset T2D as a heterogeneous disease situated on a spectrum between monogenic diabetes and adult-onset T2D.

Rare variant association analysis in 51,256 type 2 diabetes cases and 370,487 controls informs the spectrum of pathogenicity of monogenic diabetes genes.

We meta-analyzed array data imputed with the TOPMed reference panel and whole-genome sequence (WGS) datasets and performed the largest, rare variant (minor allele frequency as low as 5×10-5) GWAS meta-analysis of type 2 diabetes (T2D) comprising 51,256 cases and 370,487 controls. We identified 52 novel variants at genome-wide significance (p<5 × 10-8), including 8 novel variants that were either rare or ancestry-specific. Among them, we identified a rare missense variant in HNF4A p.Arg114Trp (OR=8.2, 95% confidence interval [CI]=4.6-14.0, p = 1.08×10-13), previously reported as a variant implicated in Maturity Onset Diabetes of the Young (MODY) with incomplete penetrance. We demonstrated that the diabetes risk in carriers of this variant was modulated by a T2D common variant polygenic risk score (cvPRS) (carriers in the top PRS tertile [OR=18.3, 95%CI=7.2-46.9, p=1.2×10-9] vs carriers in the bottom PRS tertile [OR=2.6, 95% CI=0.97-7.09, p = 0.06]. Association results identified eight variants of intermediate penetrance (OR>5) in monogenic diabetes (MD), which in aggregate as a rare variant PRS were associated with T2D in an independent WGS dataset (OR=4.7, 95% CI=1.86-11.77], p = 0.001). Our data also provided support evidence for 21% of the variants reported in ClinVar in these MD genes as benign based on lack of association with T2D. Our work provides a framework for using rare variant imputation and WGS analyses in large-scale population-based association studies to identify large-effect rare variants and provide evidence for informing variant pathogenicity.

Multi-ancestry Polygenic Mechanisms of Type 2 Diabetes Elucidate Disease Processes and Clinical Heterogeneity.

We identified genetic subtypes of type 2 diabetes (T2D) by analyzing genetic data from diverse groups, including non-European populations. We implemented soft clustering with 650 T2D-associated genetic variants, capturing known and novel T2D subtypes with distinct cardiometabolic trait associations. The twelve genetic clusters were distinctively enriched for single-cell regulatory regions. Polygenic scores derived from the clusters differed in distribution between ancestry groups, including a significantly higher proportion of lipodystrophy-related polygenic risk in East Asian ancestry. T2D risk was equivalent at a BMI of 30 kg/m2 in the European subpopulation and 24.2 (22.9-25.5) kg/m2 in the East Asian subpopulation; after adjusting for cluster-specific genetic risk, the equivalent BMI threshold increased to 28.5 (27.1-30.0) kg/m2 in the East Asian group, explaining about 75% of the difference in BMI thresholds. Thus, these multi-ancestry T2D genetic subtypes encompass a broader range of biological mechanisms and help explain ancestry-associated differences in T2D risk profiles.

Multi-ancestry Polygenic Mechanisms of Type 2 Diabetes Elucidate Disease Processes and Clinical Heterogeneity.

We identified genetic subtypes of type 2 diabetes (T2D) by analyzing genetic data from diverse groups, including non-European populations. We implemented soft clustering with 650 T2D-associated genetic variants, capturing known and novel T2D subtypes with distinct cardiometabolic trait associations. The twelve genetic clusters were distinctively enriched for single-cell regulatory regions. Polygenic scores derived from the clusters differed in distribution between ancestry groups, including a significantly higher proportion of lipodystrophy-related polygenic risk in East Asian ancestry. T2D risk was equivalent at a BMI of 30 kg/m2 in the European subpopulation and 24.2 (22.9-25.5) kg/m2 in the East Asian subpopulation; after adjusting for cluster-specific genetic risk, the equivalent BMI threshold increased to 28.5 (27.1-30.0) kg/m2 in the East Asian group, explaining about 75% of the difference in BMI thresholds. Thus, these multi-ancestry T2D genetic subtypes encompass a broader range of biological mechanisms and help explain ancestry-associated differences in T2D risk profiles.

A Type 1 Diabetes Polygenic Score Is Not Associated With Prevalent Type 2 Diabetes in Large Population Studies.

Context: Both type 1 diabetes (T1D) and type 2 diabetes (T2D) have significant genetic contributions to risk and understanding their overlap can offer clinical insight. Objective: We examined whether a T1D polygenic score (PS) was associated with a diagnosis of T2D in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium. Methods: We constructed a T1D PS using 79 known single nucleotide polymorphisms associated with T1D risk. We analyzed 13 792 T2D cases and 14 169 controls from CHARGE cohorts to determine the association between the T1D PS and T2D prevalence. We validated findings in an independent sample of 2256 T2D cases and 27 052 controls from the Mass General Brigham Biobank (MGB Biobank). As secondary analyses in 5228 T2D cases from CHARGE, we used multivariable regression models to assess the association of the T1D PS with clinical outcomes associated with T1D. Results: The T1D PS was not associated with T2D both in CHARGE (P = .15) and in the MGB Biobank (P = .87). The partitioned human leukocyte antigens only PS was associated with T2D in CHARGE (OR 1.02 per 1 SD increase in PS, 95% CI 1.01-1.03, P = .006) but not in the MGB Biobank. The T1D PS was weakly associated with insulin use (OR 1.007, 95% CI 1.001-1.012, P = .03) in CHARGE T2D cases but not with other outcomes. Conclusion: In large biobank samples, a common variant PS for T1D was not consistently associated with prevalent T2D. However, possible heterogeneity in T2D cannot be ruled out and future studies are needed do subphenotyping.

Polygenic Scores for Longitudinal Prediction of Incident Type 2 Diabetes in an Ancestrally and Medically Diverse Primary Care Network.

OBJECTIVE: The clinical utility of genetic information for type 2 diabetes (T2D) prediction with polygenic score (PGS) in ancestrally diverse, real-world US healthcare systems is unclear, especially for those at low clinical phenotypic risk for T2D. RESEARCH DESIGN AND METHODS: We tested the association of PGS with T2D incidence in patients followed within a primary care practice network over 16 years in four hypothetical scenarios that varied by clinical data availability (N = 14,712): 1) age and sex, 2) age, sex, BMI, systolic blood pressure, and family history of diabetes; 3) all variables in (2) and random glucose; 4) all variables in (3), HDL, total cholesterol, and triglycerides, combined in a clinical risk score (CRS). To determine whether genetic effects differed by baseline clinical risk, we tested for interaction with the CRS. RESULTS: PGS was associated with incident diabetes in all models. Adjusting for age and sex only, the Hazard Ratio (HR) per PGS standard deviation (SD) was 1.76 (95% CI 1.68, 1.84) and the HR of top 5% of PGS vs interquartile range (IQR) was 2.80 (2.39, 3.28). Adjusting for the CRS, the HR per SD was 1.48 (1.40, 1.57) and HR of top 5% of PGS vs IQR was 2.09 (1.72, 2.55). Genetic effects differed by baseline clinical risk [(PGS-CRS interaction p =0.05; CRS below the median: HR 1.60 (1.43, 1.79); CRS above the median: HR 1.45 (1.35, 1.55)]. CONCLUSIONS: Genetic information can help identify high-risk patients even among those perceived to be low risk in a clinical evaluation.

Algorithms for the identification of prevalent diabetes in the All of Us Research Program validated using polygenic scores - a new resource for diabetes precision medicine.

OBJECTIVE: The study aimed to develop and validate algorithms for identifying people with type 1 and type 2 diabetes in the All of Us Research Program (AoU) cohort, using electronic health record (EHR) and survey data. RESEARCH DESIGN AND METHODS: Two sets of algorithms were developed, one using only EHR data (EHR), and the other using a combination of EHR and survey data (EHR+). Their performance was evaluated by testing their association with polygenic scores for both type 1 and type 2 diabetes. RESULTS: For type 1 diabetes, the EHR-only algorithm showed a stronger association with T1D polygenic score (p=3×10-5) than the EHR+. For type 2 diabetes, the EHR+ algorithm outperformed both the EHR-only and the existing AoU definition, identifying additional cases (25.79% and 22.57% more, respectively) and showing stronger association with T2D polygenic score (DeLong p=0.03 and 1×10-4, respectively). CONCLUSIONS: We provide new validated definitions of type 1 and type 2 diabetes in AoU, and make them available for researchers. These algorithms, by ensuring consistent diabetes definitions, pave the way for high-quality diabetes research and future clinical discoveries.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus.

BACKGROUND: Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD. METHODS: From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D. RESULTS: A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase (P=1.0×10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler.

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.

Type 2 Diabetes Polygenic Score Predicts the Risk of Glucocorticoid-Induced Hyperglycemia in Patients Without Diabetes.

OBJECTIVE: To assess whether increased genetic risk of type 2 diabetes (T2D) is associated with the development of hyperglycemia after glucocorticoid treatment. RESEARCH DESIGN AND METHODS: We performed a retrospective analysis of individuals with no diagnosis of diabetes who received a glucocorticoid dose of ≥10 mg prednisone. We analyzed the association between hyperglycemia and a T2D global extended polygenic score, which was constructed through a meta-analysis of two published genome-wide association studies. RESULTS: Of 546 individuals who received glucocorticoids, 210 developed hyperglycemia and 336 did not. T2D polygenic score was significantly associated with glucocorticoid-induced hyperglycemia (odds ratio 1.4 per SD of polygenic score; P = 0.038). CONCLUSIONS: Individuals with increased genetic risk of T2D have a higher risk of glucocorticoid-induced hyperglycemia. This finding offers a mechanism for risk stratification as part of a precision approach to medical treatment.

Genome-wide association analysis identifies ancestry-specific genetic variation associated with acute response to metformin and glipizide in SUGAR-MGH.

AIMS/HYPOTHESIS: Characterisation of genetic variation that influences the response to glucose-lowering medications is instrumental to precision medicine for treatment of type 2 diabetes. The Study to Understand the Genetics of the Acute Response to Metformin and Glipizide in Humans (SUGAR-MGH) examined the acute response to metformin and glipizide in order to identify new pharmacogenetic associations for the response to common glucose-lowering medications in individuals at risk of type 2 diabetes. METHODS: One thousand participants at risk for type 2 diabetes from diverse ancestries underwent sequential glipizide and metformin challenges. A genome-wide association study was performed using the Illumina Multi-Ethnic Genotyping Array. Imputation was performed with the TOPMed reference panel. Multiple linear regression using an additive model tested for association between genetic variants and primary endpoints of drug response. In a more focused analysis, we evaluated the influence of 804 unique type 2 diabetes- and glycaemic trait-associated variants on SUGAR-MGH outcomes and performed colocalisation analyses to identify shared genetic signals. RESULTS: Five genome-wide significant variants were associated with metformin or glipizide response. The strongest association was between an African ancestry-specific variant (minor allele frequency [MAFAfr]=0.0283) at rs149403252 and lower fasting glucose at Visit 2 following metformin (p=1.9×10-9); carriers were found to have a 0.94 mmol/l larger decrease in fasting glucose. rs111770298, another African ancestry-specific variant (MAFAfr=0.0536), was associated with a reduced response to metformin (p=2.4×10-8), where carriers had a 0.29 mmol/l increase in fasting glucose compared with non-carriers, who experienced a 0.15 mmol/l decrease. This finding was validated in the Diabetes Prevention Program, where rs111770298 was associated with a worse glycaemic response to metformin: heterozygous carriers had an increase in HbA1c of 0.08% and non-carriers had an HbA1c increase of 0.01% after 1 year of treatment (p=3.3×10-3). We also identified associations between type 2 diabetes-associated variants and glycaemic response, including the type 2 diabetes-protective C allele of rs703972 near ZMIZ1 and increased levels of active glucagon-like peptide 1 (GLP-1) (p=1.6×10-5), supporting the role of alterations in incretin levels in type 2 diabetes pathophysiology. CONCLUSIONS/INTERPRETATION: We present a well-phenotyped, densely genotyped, multi-ancestry resource to study gene-drug interactions, uncover novel variation associated with response to common glucose-lowering medications and provide insight into mechanisms of action of type 2 diabetes-related variation. DATA AVAILABILITY: The complete summary statistics from this study are available at the Common Metabolic Diseases Knowledge Portal ( https://hugeamp.org ) and the GWAS Catalog ( www.ebi.ac.uk/gwas/ , accession IDs: GCST90269867 to GCST90269899).