Addressing Ancestry and Sex Bias in Pharmacogenomics.

The association of an individual's genetic makeup with their response to drugs is referred to as pharmacogenomics. By understanding the relationship between genetic variants and drug efficacy or toxicity, we are able to optimize pharmacological therapy according to an individual's genotype. Pharmacogenomics research has historically suffered from bias and underrepresentation of people from certain ancestry groups and of the female sex. These biases can arise from factors such as drugs and indications studied, selection of study participants, and methods used to collect and analyze data. To examine the representation of biogeographical populations in pharmacogenomic data sets, we describe individuals involved in gene-drug response studies from PharmGKB, a leading repository of drug-gene annotations, and showcaseCYP2D6, a gene that metabolizes approximately 25% of all prescribed drugs. We also show how the historical underrepresentation of females in clinical trials has led to significantly more adverse drug reactions in females than in males.

Major cardiovascular events and subsequent risk of kidney failure with replacement therapy: a CKD Prognosis Consortium study.

AIMS: Chronic kidney disease (CKD) increases risk of cardiovascular disease (CVD). Less is known about how CVD associates with future risk of kidney failure with replacement therapy (KFRT). METHODS AND RESULTS: The study included 25 903 761 individuals from the CKD Prognosis Consortium with known baseline estimated glomerular filtration rate (eGFR) and evaluated the impact of prevalent and incident coronary heart disease (CHD), stroke, heart failure (HF), and atrial fibrillation (AF) events as time-varying exposures on KFRT outcomes. Mean age was 53 (standard deviation 17) years and mean eGFR was 89 mL/min/1.73 m2, 15% had diabetes and 8.4% had urinary albumin-to-creatinine ratio (ACR) available (median 13 mg/g); 9.5% had prevalent CHD, 3.2% prior stroke, 3.3% HF, and 4.4% prior AF. During follow-up, there were 269 142 CHD, 311 021 stroke, 712 556 HF, and 605 596 AF incident events and 101 044 (0.4%) patients experienced KFRT. Both prevalent and incident CVD were associated with subsequent KFRT with adjusted hazard ratios (HRs) of 3.1 [95% confidence interval (CI): 2.9-3.3], 2.0 (1.9-2.1), 4.5 (4.2-4.9), 2.8 (2.7-3.1) after incident CHD, stroke, HF and AF, respectively. HRs were highest in first 3 months post-CVD incidence declining to baseline after 3 years. Incident HF hospitalizations showed the strongest association with KFRT [HR 46 (95% CI: 43-50) within 3 months] after adjustment for other CVD subtype incidence. CONCLUSION: Incident CVD events strongly and independently associate with future KFRT risk, most notably after HF, then CHD, stroke, and AF. Optimal strategies for addressing the dramatic risk of KFRT following CVD events are needed.

Competing risks analysis for neutrophil to lymphocyte ratio as a predictor of diabetic retinopathy incidence in the Scottish population.

BACKGROUND: Diabetic retinopathy (DR) is a major sight-threatening microvascular complication in individuals with diabetes. Systemic inflammation combined with oxidative stress is thought to capture most of the complexities involved in the pathology of diabetic retinopathy. A high level of neutrophil-lymphocyte ratio (NLR) is an indicator of abnormal immune system activity. Current estimates of the association of NLR with diabetes and its complications are almost entirely derived from cross-sectional studies, suggesting that the nature of the reported association may be more diagnostic than prognostic. Therefore, in the present study, we examined the utility of NLR as a biomarker to predict the incidence of DR in the Scottish population. METHODS: The incidence of DR was defined as the time to the first diagnosis of R1 or above grade in the Scottish retinopathy grading scheme from type 2 diabetes diagnosis. The effect of NLR and its interactions were explored using a competing risks survival model adjusting for other risk factors and accounting for deaths. The Fine and Gray subdistribution hazard model (FGR) was used to predict the effect of NLR on the incidence of DR. RESULTS: We analysed data from 23,531 individuals with complete covariate information. At 10 years, 8416 (35.8%) had developed DR and 2989 (12.7%) were lost to competing events (death) without developing DR and 12,126 individuals did not have DR. The median (interquartile range) level of NLR was 2.04 (1.5 to 2.7). The optimal NLR cut-off value to predict retinopathy incidence was 3.04. After accounting for competing risks at 10 years, the cumulative incidence of DR and deaths without DR were 50.7% and 21.9%, respectively. NLR was associated with incident DR in both Cause-specific hazard (CSH = 1.63; 95% CI: 1.28-2.07) and FGR models the subdistribution hazard (sHR = 2.24; 95% CI: 1.70-2.94). Both age and HbA1c were found to modulate the association between NLR and the risk of DR. CONCLUSIONS: The current study suggests that NLR has a promising potential to predict DR incidence in the Scottish population, especially in individuals less than 65 years and in those with well-controlled glycaemic status.

Incidence and Outcomes of Patients Receiving Chronic Kidney Replacement Therapy Admitted to Scottish ICUs Between 2009 and 2019-A National Observational Cohort Study.

OBJECTIVES: To determine the incidence and characteristics of ICU admissions in the Scottish population of patients treated with chronic kidney replacement therapy (KRT) over an 11-year period and determine factors associated with post-ICU admission mortality. DESIGN: Retrospective observational cohort study. SETTING: We analyzed admissions to Scottish intensive care environments between January 1, 2009, and December 31, 2019. PATIENTS: All patients receiving chronic KRT-including maintenance dialysis and kidney transplant-in Scotland. INTERVENTION: None. MEASUREMENTS AND MAIN RESULTS: Descriptive statistics and factors associated with mortality using logistic regression and Cox proportional hazard models. From 10,657 unique individuals registered in the Scottish Renal Registry over the 11-year study period and alive as of January 1, 2009, 1,402 adult patients were identified as being admitted to a Scottish critical care setting. Between 2009 and 2019, admissions to ICU increased in a nonlinear manner driven by increases in admissions for renal causes and elective cardiac surgery. The ICU admission rate was higher among patients on chronic dialysis than in kidney transplant recipients (59.1 vs 19.9 per 1,000 person-years), but post-ICU mortality was similar (about 24% at 30 d and 40% at 1 year). Admissions for renal reasons were most common (20.9%) in patients undergoing chronic dialysis, whereas kidney transplant recipients were most frequently admitted for pneumonia (19.3%) or sepsis (12.8%). Adjusted Cox PH models showed that receiving invasive ventilation and vasoactive drugs was associated with an increased risk of death at 30 days post-ICU admission (HR, 1.75; 95% CI, 1.28-2.39 and 1.72; 95% CI, 1.28-2.31, respectively). CONCLUSIONS: With a growing population of kidney transplant recipients and the improved survival of patients on chronic dialysis, the number of ICU admissions is rising in the chronic KRT population. Mortality post-ICU admission is high for these patients.

Weight variability and cardiovascular outcomes: a systematic review and meta-analysis.

The association between body weight variability and the risk of cardiovascular disease (CVD) has been investigated previously with mixed findings. However, there has been no extensive study which systematically evaluates the current evidence. Furthermore, the impact of ethnicity and type 2 diabetes on this phenomena has not yet been investigated. Therefore, the aim of this study was to comprehensively evaluate the effect of weight variability on risk of CVD (any cardiovascular (CV) event, composite CV outcome, CV death, Stroke, Myocardial Infarction) and the influence of ethnicity and type 2 diabetes status on the observed association. A systematic review and meta-analysis was performed according to the meta-analyses of observational studies in epidemiology (MOOSE) guidelines. The electronic databases PubMed, Web of Science, and the Cochrane Library were searched for studies that investigated the relationship between body weight or BMI variability and CV diseases using Medical Subject Headings (MeSH) terms and keywords. The relative risks (RRs) for the outcomes were collected from studies, pooled, and analysed using a random-effects model to estimate the overall relative risk. Of 5645 articles screened, 23 studies with a total population of 15,382,537 fulfilled the prespecified criteria and were included. Individuals in the highest strata of body weight variability were found to have significantly increased risk of any CV event (RR = 1.27; 95% Confidence Interval (CI) 1.17-1.38; P < 0.0001; I2 = 97.28%), cardiovascular death (RR = 1.29; 95% CI 1.03-1.60; P < 0.0001; I2 = 55.16%), myocardial infarction (RR = 1.32; 95% CI 1.09-1.59; P = 0.0037; I2 = 97.14%), stroke (RR = 1.21; 95% CI 1.19-1.24; P < 0.0001; I2 = 0.06%), and compound CVD outcomes (RR = 1.36; 95% CI 1.08-1.73; P = 0.01; I2 = 92.41%). Similar RRs were observed regarding BMI variability and per unit standard deviation (SD) increase in body weight variability. Comparable effects were seen in people with and without diabetes, in White Europeans and Asians. In conclusion, body weight variability is associated with increased risk of CV diseases regardless of ethnicity or diabetes status. Future research is needed to prove a causative link between weight variability and CVD risk, as appropriate interventions to maintain stable weight could positively influence CVD.

Diabetes status modifies the long-term effect of lipoprotein-associated phospholipase A2 on major coronary events.

AIMS/HYPOTHESIS: Lipoprotein-associated phospholipase A2 (Lp-PLA2) activity has an independent prognostic association with major coronary events (MCE). However, no study has investigated whether type 2 diabetes status modifies the effect of Lp-PLA2 activity or inhibition on the risk of MCE. We investigate the interaction between diabetes status and Lp-PLA2 activity with risk of MCE. Subsequently, we test the resulting hypothesis that diabetes status will play a role in modifying the efficacy of an Lp-PLA2 inhibitor. METHODS: A retrospective cohort study design was utilised in two study populations. Discovery analyses were performed in the Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) cohort based in Scotland, UK. Participants were categorised by type 2 diabetes control status: poorly controlled (HbA1c ≥ 48 mmol/mol or ≥6.5%) and well-controlled (HbA1c < 48 mmol/mol or <6.5%) diabetes (n = 7420). In a secondary analysis of the Stabilization of Atherosclerotic Plaque by Initiation of Darapladib Therapy (STABILITY) trial of Lp-PLA2 inhibitor (darapladib) efficacy, 15,828 participants were stratified post hoc by type 2 diabetes diagnosis status (diabetes or no diabetes) at time of recruitment. Lp-PLA2 activity was then divided into population-specific quartiles. MCE were determined from linked medical records in GoDARTS and trial records in STABILITY. First, the interaction between diabetes control status and Lp-PLA2 activity on the outcome of MCE was explored in GoDARTS. The effect was replicated in the placebo arm of STABILITY. The effect of Lp-PLA2 on MCE was then examined in models stratified by diabetes status. This helped determine participants at higher risk. Finally, the effect of Lp-PLA2 inhibition was assessed in STABILITY in the higher risk group. Cox proportional hazards models adjusted for confounders were used to assess associations. RESULTS: In GoDARTS, a significant interaction between increased Lp-PLA2 activity (continuous and quartile divided) and diabetes control status was observed in the prediction of MCE (p < 0.0001). These effects were replicated in the placebo arm of STABILITY (p < 0.0001). In GoDARTS, stratified analyses showed that, among individuals with poorly controlled diabetes, the hazards of MCE for those with high (Q4) Lp-PLA2 activity was 1.19 compared with individuals with lower (Q1-3) Lp-PLA2 activity (95% CI 1.11, 1.38; p < 0.0001) and 1.35 (95% CI 1.16, 1.57; p < 0.0001) when compared with those with the lowest activity (Q1). Those in the higher risk group were identified as individuals with the highest Lp-PLA2 activity (Q4) and poorly controlled diabetes or diabetes. Based on these observations in untreated populations, we hypothesised that the Lp-PLA2 inhibitor would have more benefit in this higher risk group. In this risk group, Lp-PLA2 inhibitor use was associated with a 33% reduction in MCE compared with placebo (HR 0.67 [95% CI 0.50, 0.90]; p = 0.008). In contrast, Lp-PLA2 inhibitor showed no efficacy in individuals with low activity, regardless of diabetes status, or among those with no baseline diabetes and high Lp-PLA2 activity. CONCLUSIONS/INTERPRETATION: These results support the hypothesis that diabetes status modifies the association between Lp-PLA2 activity and MCE. These results suggest that cardiovascular morbidity and mortality associated with Lp-PLA2 activity is especially important in patients with type 2 diabetes, particularly those with worse glycaemic control. Further investigation of the effects of Lp-PLA2 inhibition in diabetes appears warranted. DATA AVAILABILITY: STABILITY trial data are available from clinicaltrials.gov repository through the GlaxoSmithKline clinical study register https://clinicaltrials.gov/ct2/show/NCT00799903 . GoDARTS datasets generated during and/or analysed during the current study are available following request to the GoDARTS Access Managements Group https://godarts.org/scientific-community/ .

Young-onset diabetes in Asian Indians is associated with lower measured and genetically determined beta cell function.

AIMS/HYPOTHESIS: South Asians in general, and Asian Indians in particular, have higher risk of type 2 diabetes compared with white Europeans, and a younger age of onset. The reasons for the younger age of onset in relation to obesity, beta cell function and insulin sensitivity are under-explored. METHODS: Two cohorts of Asian Indians, the ICMR-INDIAB cohort (Indian Council of Medical Research-India Diabetes Study) and the DMDSC cohort (Dr Mohan's Diabetes Specialties Centre), and one of white Europeans, the ESDC (East Scotland Diabetes Cohort), were used. Using a cross-sectional design, we examined the comparative prevalence of healthy, overweight and obese participants with young-onset diabetes, classified according to their BMI. We explored the role of clinically measured beta cell function in diabetes onset in Asian Indians. Finally, the comparative distribution of a partitioned polygenic score (pPS) for risk of diabetes due to poor beta cell function was examined. Replication of the genetic findings was sought using data from the UK Biobank. RESULTS: The prevalence of young-onset diabetes with normal BMI was 9.3% amongst white Europeans and 24-39% amongst Asian Indians. In Asian Indians with young-onset diabetes, after adjustment for family history of type 2 diabetes, sex, insulin sensitivity and HDL-cholesterol, stimulated C-peptide was 492 pmol/ml (IQR 353-616, p<0.0001) lower in lean compared with obese individuals. Asian Indians in our study, and South Asians from the UK Biobank, had a higher number of risk alleles than white Europeans. After weighting the pPS for beta cell function, Asian Indians have lower genetically determined beta cell function than white Europeans (p<0.0001). The pPS was associated with age of diagnosis in Asian Indians but not in white Europeans. The pPS explained 2% of the variation in clinically measured beta cell function, and 1.2%, 0.97%, and 0.36% of variance in age of diabetes amongst Asian Indians with normal BMI, or classified as overweight and obese BMI, respectively. CONCLUSIONS/INTERPRETATION: The prevalence of lean BMI in young-onset diabetes is over two times higher in Asian Indians compared with white Europeans. This phenotype of lean, young-onset diabetes appears driven in part by lower beta cell function. We demonstrate that Asian Indians with diabetes also have lower genetically determined beta cell function.

Cardiovascular outcomes in patients with chronic kidney disease and COVID-19: a multi-regional data-linkage study.

BACKGROUND: Data describing cardiovascular outcomes in patients with coronavirus disease 2019 (COVID-19) and chronic kidney disease (CKD) are lacking. We compared cardiovascular outcomes of patients with and without COVID-19, stratified by CKD status. METHODS: This retrospective, multi-regional data-linkage study utilised individual patient-level data from two Scottish cohorts. All patients tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Cohort 1 between 1 February 2020 and 31 March 2021 and in Cohort 2 between 28 February 2020 and 8 February 2021 were included. RESULTS: Overall, 86 964 patients were tested for SARS-CoV-2. There were 36 904 patients (mean±sd age 61±21 years; 58.1% women; 15.9% CKD; 10.1% COVID-19 positive) in Cohort 1 and 50 060 patients (mean±sd age 63±20 years; 62.0% women; 16.4% CKD; 9.1% COVID-19 positive) in Cohort 2. In CKD patients, COVID-19 increased the risk of cardiovascular death by more than two-fold within 30 days (cause-specific hazard ratio (csHR) meta-estimate 2.34, 95% CI 1.83-2.99) and by 57% at the end of study follow-up (csHR meta-estimate 1.57, 95% CI 1.31-1.89). Similarly, the risk of all-cause death in COVID-19 positive versus negative CKD patients was greatest within 30 days (HR 4.53, 95% CI 3.97-5.16). Compared with patients without CKD, those with CKD had a higher risk of testing positive (11.5% versus 9.3%). Following a positive test, CKD patients had higher rates of cardiovascular death (11.1% versus 2.7%), cardiovascular complications and cardiovascular hospitalisations (7.1% versus 3.3%) than those without CKD. CONCLUSIONS: COVID-19 increases the risk of cardiovascular and all-cause death in CKD patients, especially in the short-term. CKD patients with COVID-19 are also at a disproportionate risk of cardiovascular complications than those without CKD.

The Relationship between AKI and CKD in Patients with Type 2 Diabetes: An Observational Cohort Study.

BACKGROUND: There are few observational studies evaluating the risk of AKI in people with type 2 diabetes, and even fewer simultaneously investigating AKI and CKD in this population. This limits understanding of the interplay between AKI and CKD in people with type 2 diabetes compared with the nondiabetic population. METHODS: In this retrospective, cohort study of participants with or without type 2 diabetes, we used electronic healthcare records to evaluate rates of AKI and various statistical methods to determine their relationship to CKD status and further renal function decline. RESULTS: We followed the cohort of 16,700 participants (9417 with type 2 diabetes and 7283 controls without diabetes) for a median of 8.2 years. Those with diabetes were more likely than controls to develop AKI (48.6% versus 17.2%, respectively) and have preexisting CKD or CKD that developed during follow-up (46.3% versus 17.2%, respectively). In the absence of CKD, the AKI rate among people with diabetes was nearly five times that of controls (121.5 versus 24.6 per 1000 person-years). Among participants with CKD, AKI rate in people with diabetes was more than twice that of controls (384.8 versus 180.0 per 1000 person-years after CKD diagnostic date, and 109.3 versus 47.4 per 1000 person-years before CKD onset in those developing CKD after recruitment). Decline in eGFR slope before AKI episodes was steeper in people with diabetes versus controls. After AKI episodes, decline in eGFR slope became steeper in people without diabetes, but not among those with diabetes and preexisting CKD. CONCLUSIONS: Patients with diabetes have significantly higher rates of AKI compared with patients without diabetes, and this remains true for individuals with preexisting CKD.

Genome-wide association study of angioedema induced by angiotensin-converting enzyme inhibitor and angiotensin receptor blocker treatment.

Angioedema in the mouth or upper airways is a feared adverse reaction to angiotensin-converting enzyme inhibitor (ACEi) and angiotensin receptor blocker (ARB) treatment, which is used for hypertension, heart failure and diabetes complications. This candidate gene and genome-wide association study aimed to identify genetic variants predisposing to angioedema induced by these drugs. The discovery cohort consisted of 173 cases and 4890 controls recruited in Sweden. In the candidate gene analysis, ETV6, BDKRB2, MME, and PRKCQ were nominally associated with angioedema (p < 0.05), but did not pass Bonferroni correction for multiple testing (p < 2.89 × 10-5). In the genome-wide analysis, intronic variants in the calcium-activated potassium channel subunit alpha-1 (KCNMA1) gene on chromosome 10 were significantly associated with angioedema (p < 5 × 10-8). Whilst the top KCNMA1 hit was not significant in the replication cohort (413 cases and 599 ACEi-exposed controls from the US and Northern Europe), a meta-analysis of the replication and discovery cohorts (in total 586 cases and 1944 ACEi-exposed controls) revealed that each variant allele increased the odds of experiencing angioedema 1.62 times (95% confidence interval 1.05-2.50, p = 0.030). Associated KCNMA1 variants are not known to be functional, but are in linkage disequilibrium with variants in transcription factor binding sites active in relevant tissues. In summary, our data suggest that common variation in KCNMA1 is associated with risk of angioedema induced by ACEi or ARB treatment. Future whole exome or genome sequencing studies will show whether rare variants in KCNMA1 or other genes contribute to the risk of ACEi- and ARB-induced angioedema.

Lp-PLA2 activity is associated with increased risk of diabetic retinopathy: a longitudinal disease progression study.

AIMS/HYPOTHESIS: The aim of the study was to examine the association between lipoprotein-associated phospholipase A2 (Lp-PLA2) activity levels and incident diabetic retinopathy and change in retinopathy grade. METHODS: This was a cohort study of diabetic participants with serum collected at baseline and routinely collected diabetic retinal screening data. Participants with type 2 diabetes from the GoDARTS (Genetics of Diabetes Audit and Research in Tayside Scotland) cohort were used. This cohort is composed of individuals of white Scottish ancestry from the Tayside region of Scotland. Survival analysis accounting for informative censoring by modelling death as a competing risk was performed for the development of incident diabetic retinopathy from a disease-free state in a 3 year follow-up period (n = 1364) by stratified Lp-PLA2 activity levels (in quartiles). The same analysis was performed for transitions to more severe grades. RESULTS: The hazard of developing incident diabetic retinopathy was 2.08 times higher (95% CI 1.64, 2.63) for the highest quartile of Lp-PLA2 activity compared with the lowest. Higher Lp-PLA2 activity levels were associated with a significantly increased risk for transitions to all grades. The hazards of developing observable (or more severe) and referable (or more severe) retinopathy were 2.82 (95% CI 1.71, 4.65) and 1.87 (95% CI 1.26, 2.77) times higher for the highest quartile of Lp-PLA2 activity compared with the lowest, respectively. CONCLUSIONS/INTERPRETATION: Higher Lp-PLA2 levels are associated with increased risk of death and the development of incident diabetic retinopathy, as well as transitions to more severe grades of diabetic retinopathy. These associations are independent of calculated LDL-cholesterol and other traditional risk factors. Further, this biomarker study shows that the association is temporally sensitive to the proximity of the event to measurement of Lp-PLA2.

A common missense variant of LILRB5 is associated with statin intolerance and myalgia.

Aims: A genetic variant in LILRB5 (leukocyte immunoglobulin-like receptor subfamily-B) (rs12975366: T > C: Asp247Gly) has been reported to be associated with lower creatine phosphokinase (CK) and lactate dehydrogenase (LDH) levels. Both biomarkers are released from injured muscle tissue, making this variant a potential candidate for susceptibility to muscle-related symptoms. We examined the association of this variant with statin intolerance ascertained from electronic medical records in the GoDARTS study. Methods and results: In the GoDARTS cohort, the LILRB5 Asp247 variant was associated with statin intolerance (SI) phenotypes; one defined as having raised CK and being non-adherent to therapy [odds ratio (OR) 1.81; 95% confidence interval (CI): 1.34-2.45] and the other as being intolerant to the lowest approved dose of a statin before being switched to two or more other statins (OR 1.36; 95% CI: 1.07-1.73). Those homozygous for Asp247 had increased odds of developing both definitions of intolerance. Importantly the second definition did not rely on CK elevations. These results were replicated in adjudicated cases of statin-induced myopathy in the PREDICTION-ADR consortium (OR1.48; 95% CI: 1.05-2.10) and for the development of myalgia in the JUPITER randomized clinical trial of rosuvastatin (OR1.35, 95% CI: 1.10-1.68). A meta-analysis across the studies showed a consistent association between Asp247Gly and outcomes associated with SI (OR1.34; 95% CI: 1.16-1.54). Conclusion: This study presents a novel immunogenetic factor associated with statin intolerance, an important risk factor for cardiovascular outcomes. The results suggest that true statin-induced myalgia and non-specific myalgia are distinct, with a potential role for the immune system in their development. We identify a genetic group that is more likely to be intolerant to their statins.

Replication confirms the association of loci in FOXE1, PDE8B, CAPZB and PDE10A with thyroid traits: a Genetics of Diabetes Audit and Research Tayside study (GoDARTS).

OBJECTIVE: Replication of associations in genome-wide association studies is desirable to ensure that such signals are potentially clinically meaningful. This study aimed to replicate associations of selected single-nucleotide polymorphisms (SNPs) with hypothyroidism and serum thyroid-stimulating hormone (TSH) using electronic medical records (EMRs). PATIENTS AND METHODS: A cross-sectional study was carried out among patients of European Caucasian ethnicity from the Genetics of Diabetes Audit and Research Tayside recruited in Tayside (Scotland, UK). EMRs (biochemistry, prescribing, hospital admissions and demographics) were used to ascertain patients with hypothyroidism and their controls as well as average serum TSH concentration, and linked to genetic biobank data. Genetic tests of association were performed using logistic and linear regression models. RESULTS: We analysed 1703 cases of hypothyroidism and 9457 controls. All four SNPs located on chromosome 9 at FOXE1 were associated with hypothyroidism with similar effect estimates (odds ratio=0.75-0.76, P<5e-08). Also, loci on chromosomes 1 (PTPN22), six (HLA-E/HLA-C) and 12 (SH2B3) were replicated. For serum TSH, we confirmed 12 SNPs previously reported at PDE8B, CAPZB, PDE10A, LOC105371356, NR3C2, VEGFA, IGFBP5, INSR, PRDM11, NFIA, ITPK1 and ABO. Overall, these SNPs accounted for 6.8% of the serum TSH variation (P<1e-04). CONCLUSION: EMRs linked to genomic data in large populations enable validation of genome-wide association studies discoveries without additional genotyping costs. Our replication confirmed at genome-wide significance the association of loci at FOXE1 with hypothyroidism, and PDE8B, CAPZB and PDE10A with serum TSH. A total of 12 SNPs seemed to explain nearly 7% of the serum TSH variation.

Meta-analysis of genome-wide association studies on the intolerance of angiotensin-converting enzyme inhibitors.

OBJECTIVES: To identify single nucleotide polymorphisms (SNPs) associated with switching from an angiotensin-converting enzyme (ACE)-inhibitor to an angiotensin receptor blocker. METHODS: Two cohorts of patients starting ACE-inhibitors were identified within the Rotterdam Study in the Netherlands and the Genetics of Diabetes Audit and Research in Tayside Scotland study in Scotland. Cases were intolerant patients who switched from an ACE-inhibitor to an angiotensin receptor blocker and controls were individuals who used ACE-inhibitors continuously for at least 2 years and did not switch. Genome-wide association study (GWAS) using an additive model was run in these sets and the results were meta-analysed using Genome-Wide Association Meta Analysis software. RESULTS: A total of 972 cases out of 5161 ACE-inhibitor starters were identified. Eight SNPs within four genes reached the genome-wide association study significance level (P<5×10) in the meta-analysis [RNA binding protein, Fox-1 homolog (Caenorhabditis elegans), γ-aminobutyric acid receptor subunit γ-2, sarcoma (Src) homology 2 (SH2) B adaptor protein 1 and membrane bound O-acyltransferase domain containing 1]. The strongest associated SNP was located in an intron of RNA binding protein, Fox-1 homolog (Caenorhabditis elegans), which contains an RNA binding protein [rs2061538: minor allele frequency=0.16, odds ratio=1.52 (95% confidence interval: 1.32-1.76), P=6.2×10]. CONCLUSION: These results indicate that genetic variation in the above-mentioned genes may increase the risk of ACE-inhibitor-induced adverse reactions.

Common and Distinct Genetic Architecture of Age at Diagnosis of Diabetes in South Indian and European Populations.

OBJECTIVE: South Asians are diagnosed with type 2 diabetes (T2D) more than a decade earlier in life than seen in European populations. We hypothesized that studying the genomics of age of diagnosis in these populations may give insight into the earlier age diagnosis of T2D among individuals of South Asian descent. RESEARCH DESIGN AND METHODS: We conducted a meta-analysis of genome-wide association studies (GWAS) of age at diagnosis of T2D in 34,001 individuals from four independent cohorts of European and South Asian Indians. RESULTS: We identified two signals near the TCF7L2 and CDKAL1 genes associated with age at the onset of T2D. The strongest genome-wide significant variants at chromosome 10q25.3 in TCF7L2 (rs7903146; P = 2.4 × 10-12, ß = -0.436; SE 0.02) and chromosome 6p22.3 in CDKAL1 (rs9368219; P = 2.29 × 10-8; ß = -0.053; SE 0.01) were directionally consistent across ethnic groups and present at similar frequencies; however, both loci harbored additional independent signals that were only present in the South Indian cohorts. A genome-wide signal was also obtained at chromosome 10q26.12 in WDR11 (rs3011366; P = 3.255 × 10-8; ß = 1.44; SE 0.25), specifically in the South Indian cohorts. Heritability estimates for the age at diagnosis were much stronger in South Indians than Europeans, and a polygenic risk score constructed based on South Indian GWAS explained ∼2% trait variance. CONCLUSIONS: Our findings provide a better understanding of ethnic differences in the age at diagnosis and indicate the potential importance of ethnic differences in the genetic architecture underpinning T2D.

A gene risk score using missense variants in SLCO1B1 is associated with earlier onset statin intolerance.

BACKGROUND AND AIMS: The efficacy of statin therapy is hindered by intolerance to the therapy, leading to discontinuation. Variants in SLCO1B1, which encodes the hepatic transporter OATB1B1, influence statin pharmacokinetics, resulting in altered plasma concentrations of the drug and its metabolites. Current pharmacogenetic guidelines require sequencing of the SLCO1B1 gene, which is more expensive and less accessible than genotyping. In this study, we aimed to develop an easy, clinically implementable functional gene risk score (GRS) of common variants in SLCO1B1 to identify patients at risk of statin intolerance. METHODS AND RESULTS: A GRS was developed from four common variants in SLCO1B1. In statin users from Tayside, Scotland, UK, those with a high-risk GRS had increased odds across three phenotypes of statin intolerance [general statin intolerance (GSI): ORGSI 2.42; 95% confidence interval (CI): 1.29-4.31, P = 0.003; statin-related myopathy: ORSRM 2.51; 95% CI: 1.28-4.53, P = 0.004; statin-related suspected rhabdomyolysis: ORSRSR 2.85; 95% CI: 1.03-6.65, P = 0.02]. In contrast, using the Val174Ala genotype alone or the recommended OATP1B1 functional phenotypes produced weaker and less reliable results. A meta-analysis with results from adjudicated cases of statin-induced myopathy in the PREDICTION-ADR Consortium confirmed these findings (ORVal174Ala 1.99; 95% CI: 1.01-3.95, P = 0.048; ORGRS 1.76; 95% CI: 1.16-2.69, P = 0.008). For those requiring high-dose statin therapy, the high-risk GRS was more consistently associated with the time to onset of statin intolerance amongst the three phenotypes compared with Val174Ala (GSI: HRVal174Ala 2.49; 95% CI: 1.09-5.68, P = 0.03; HRGRS 2.44; 95% CI: 1.46-4.08, P < 0.001). Finally, sequence kernel association testing confirmed that rare variants in SLCO1B1 are associated with the risk of intolerance (P = 0.02). CONCLUSION: We provide evidence that a GRS based on four common SLCO1B1 variants provides an easily implemented genetic tool that is more reliable than the current recommended practice in estimating the risk and predicting early-onset statin intolerance.