Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis

BACKGROUND

Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes.

METHODS & FINDINGS

Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants.

CONCLUSIONS

As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Genetic analysis of emerging risk factors in coronary artery disease

BACKGROUND AND AIMS

Type 2 diabetes (T2D), low-density lipoprotein-cholesterol (LDL-c), body mass index (BMI), blood pressure and smoking are established risk factors that play a causal role in coronary artery disease (CAD). Numerous common genetic variants associating with these and other risk factors have been identified, but their association with CAD has not been comprehensively examined in a single study. Our goal was to comprehensively evaluate the associations of established and emerging risk factors with CAD using genetic variants identified from Genome-wide Association Studies (GWAS).

METHODS

We tested the effect of 60 traditional and putative risk factors with CAD, using summary statistics obtained in GWAS. We approximated the regression of a response variable onto an additive multi-SNP genetic risk score in the Coronary Artery DIsease Genomewide Replication And Meta-analysis (CARDIoGRAM) consortium dataset weighted by the effect of the SNP on the risk factors.

RESULTS

The strongest association with risk of CAD was for LDL-c SNPs (p = 3.96E-34). For non-established CAD risk factors, we found significant CAD associations for coronary artery calcification (CAC), Lp(a), LP-PLA2 activity, plaque, vWF and FVIII. In an attempt to identify independent associations between risk factors and CAD, only SNPs with an effect on the target trait were included. This identified CAD associations for Lp(a)(p = 1.77E-21), LDL-c (p = 4.16E-06), triglycerides (TG) (p = 1.94E-05), height (p = 2.06E-05), CAC (p = 3.13E-23) and carotid plaque (p = 2.08E-05).

CONCLUSIONS

We identified SNPs associated with the emerging risk factors Lp(a), TG, plaque, height and CAC to be independently associated with risk of CAD. This provides further support for-ongoing clinical trials of Lp(a) and TG, and suggests that CAC and plaque could be used as surrogate markers for CAD in clinical trials.

HBsAg loss in patients treated with peginterferon alfa-2a and adefovir is associated with SLC16A9 gene variation and lower plasma carnitine levels

BACKGROUND & AIMS

Achievement of HBsAg loss remains the hallmark of chronic hepatitis B treatment. In order to identify host factors contributing to treatment-induced HBsAg loss, we performed a genome-wide screen of single nucleotide polymorphisms (SNPs) and studied its immunological consequence.

METHODS

Chronic hepatitis B patients (40 HBeAg-positive and 44 HBeAg-negative) treated with peginterferon alfa-2a and adefovir were genotyped for 999,091 SNPs, which were associated with HBsAg loss at week 96 (n = 9). Plasma carnitine levels were measured by tandem-mass spectrometry, and the effect of carnitine on the proliferative capacity of hepatitis B virus (HBV)-specific and non-specific CD8 T cells was studied in vitro.

RESULTS

One polymorphism, rs12356193 located in the SLC16A9 gene, was genome-wide significantly associated with HBsAg loss at week 96 (p = 1.84 × 10(-8)). The previously reported association of rs12356193 with lower carnitine levels was confirmed in our cohort, and baseline carnitine levels were lower in patients with HBsAg loss compared to patients with HBsAg persistence (p = 0.02). Furthermore, we demonstrated that carnitine suppressed HBV-specific CD8 T cell proliferation.

CONCLUSIONS

In chronic hepatitis B patients treated with peginterferon and adefovir, we identified strong associations of SLC16A9 gene variation and carnitine levels with HBsAg loss. Our results further suggest that a lower baseline plasma carnitine level increases the proliferative capacity of CD8 T cells, making patients more susceptible to the immunological effect of this treatment. These novel findings may provide new insight into factors involved in treatment-induced HBsAg loss, and play a role in the prediction of treatment outcome.

Secretory phospholipase A(2)-IIA and cardiovascular disease: a mendelian randomization study

OBJECTIVES

This study sought to investigate the role of secretory phospholipase A2 (sPLA2)-IIA in cardiovascular disease.

BACKGROUND

Higher circulating levels of sPLA2-IIA mass or sPLA2 enzyme activity have been associated with increased risk of cardiovascular events. However, it is not clear if this association is causal. A recent phase III clinical trial of an sPLA2 inhibitor (varespladib) was stopped prematurely for lack of efficacy.

METHODS

We conducted a Mendelian randomization meta-analysis of 19 general population studies (8,021 incident, 7,513 prevalent major vascular events [MVE] in 74,683 individuals) and 10 acute coronary syndrome (ACS) cohorts (2,520 recurrent MVE in 18,355 individuals) using rs11573156, a variant in PLA2G2A encoding the sPLA2-IIA isoenzyme, as an instrumental variable.

RESULTS

PLA2G2A rs11573156 C allele associated with lower circulating sPLA2-IIA mass (38% to 44%) and sPLA2 enzyme activity (3% to 23%) per C allele. The odds ratio (OR) for MVE per rs11573156 C allele was 1.02 (95% confidence interval [CI]: 0.98 to 1.06) in general populations and 0.96 (95% CI: 0.90 to 1.03) in ACS cohorts. In the general population studies, the OR derived from the genetic instrumental variable analysis for MVE for a 1-log unit lower sPLA2-IIA mass was 1.04 (95% CI: 0.96 to 1.13), and differed from the non-genetic observational estimate (OR: 0.69; 95% CI: 0.61 to 0.79). In the ACS cohorts, both the genetic instrumental variable and observational ORs showed a null association with MVE. Instrumental variable analysis failed to show associations between sPLA2 enzyme activity and MVE.

CONCLUSIONS

Reducing sPLA2-IIA mass is unlikely to be a useful therapeutic goal for preventing cardiovascular events.

Common genetic variants do not associate with CAD in familial hypercholesterolemia

In recent years, multiple loci dispersed on the genome have been shown to be associated with coronary artery disease (CAD). We investigated whether these common genetic variants also hold value for CAD prediction in a large cohort of patients with familial hypercholesterolemia (FH). We genotyped a total of 41 single-nucleotide polymorphisms (SNPs) in 1701 FH patients, of whom 482 patients (28.3%) had at least one coronary event during an average follow up of 66 years. The association of each SNP with event-free survival time was calculated with a Cox proportional hazard model. In the cardiovascular disease risk factor adjusted analysis, the most significant SNP was rs1122608:G>T in the SMARCA4 gene near the LDL-receptor (LDLR) gene, with a hazard ratio for CAD risk of 0.74 (95% CI 0.49-0.99; P-value 0.021). However, none of the SNPs reached the Bonferroni threshold. Of all the known CAD loci analyzed, the SMARCA4 locus near the LDLR had the strongest negative association with CAD in this high-risk FH cohort. The effect is contrary to what was expected. None of the other loci showed association with CAD.

Gene-centric meta-analysis of lipid traits in African, East Asian and Hispanic populations

Meta-analyses of European populations has successfully identified genetic variants in over 100 loci associated with lipid levels, but our knowledge in other ethnicities remains limited. To address this, we performed dense genotyping of ∼2,000 candidate genes in 7,657 African Americans, 1,315 Hispanics and 841 East Asians, using the IBC array, a custom ∼50,000 SNP genotyping array. Meta-analyses confirmed 16 lipid loci previously established in European populations at genome-wide significance level, and found multiple independent association signals within these lipid loci. Initial discovery and in silico follow-up in 7,000 additional African American samples, confirmed two novel loci: rs5030359 within ICAM1 is associated with total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) (p = 8.8×10(-7) and p = 1.5×10(-6) respectively) and a nonsense mutation rs3211938 within CD36 is associated with high-density lipoprotein cholesterol (HDL-C) levels (p = 13.5×10(-12)). The rs3211938-G allele, which is nearly absent in European and Asian populations, has been previously found to be associated with CD36 deficiency and shows a signature of selection in Africans and African Americans. Finally, we have evaluated the effect of SNPs established in European populations on lipid levels in multi-ethnic populations and show that most known lipid association signals span across ethnicities. However, differences between populations, especially differences in allele frequency, can be leveraged to identify novel signals, as shown by the discovery of ICAM1 and CD36 in the current report.