Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Increased genetic contribution to wellbeing during the COVID-19 pandemic.

Physical and mental health are determined by an interplay between nature, for example genetics, and nurture, which encompasses experiences and exposures that can be short or long-lasting. The COVID-19 pandemic represents a unique situation in which whole communities were suddenly and simultaneously exposed to both the virus and the societal changes required to combat the virus. We studied 27,537 population-based biobank participants for whom we have genetic data and extensive longitudinal data collected via 19 questionnaires over 10 months, starting in March 2020. This allowed us to explore the interaction between genetics and the impact of the COVID-19 pandemic on individuals' wellbeing over time. We observe that genetics affected many aspects of wellbeing, but also that its impact on several phenotypes changed over time. Over the course of the pandemic, we observed that the genetic predisposition to life satisfaction had an increasing influence on perceived quality of life. We also estimated heritability and the proportion of variance explained by shared environment using variance components methods based on pedigree information and household composition. The results suggest that people's genetic constitution manifested more prominently over time, potentially due to social isolation driven by strict COVID-19 containment measures. Overall, our findings demonstrate that the relative contribution of genetic variation to complex phenotypes is dynamic rather than static.

Association between 6PPD-quinone exposure and BMI, influenza, and diarrhea in children.

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine-quinone (6PPD-quinone) has received extensive attention due to its ubiquitous distribution and potential toxicity. However, the distribution characteristics of 6PPD-quinone in dust from e-waste recycling areas and the consequential health risks to children are unclear. A total of 183 dust samples were collected from roads (n = 40), homes (n = 91), and kindergartens (n = 52) in Guiyu (the e-waste-exposed group) and Haojiang (the reference group) from 2019 to 2021. The results show that the concentrations of 6PPD-quinone in kindergarten and house dust from the exposed group were significantly higher than those from the reference group (P < 0.001). These findings show that e-waste may be another potential source of 6PPD-quinone, in addition to rubber tires. The exposure risk of 6PPD-quinone in children was assessed using their daily intake. The daily intake of 925 kindergarten children was calculated using the concentration of 6PPD-quinone in kindergarten dust. The daily intake of 6PPD-quinone via ingestion was approximately five orders of magnitude higher than via inhalation. Children in the exposed group had a higher exposure risk to 6PPD-quinone than the reference group. A higher daily intake of 6PPD-quinone from kindergarten dust was associated with a lower BMI and a higher frequency of influenza and diarrhea in children. This study reports the distribution of 6PPD-quinone in an e-waste recycling town and explores the associated health risks to children.

Association of CYP2D6 and CYP2C19 metabolizer status with switching and discontinuing antidepressant drugs: an exploratory study.

BACKGROUND: Tailoring antidepressant drugs (AD) to patients' genetic drug-metabolism profile is promising. However, literature regarding associations of ADs' treatment effect and/or side effects with drug metabolizing genes CYP2D6 and CYP2C19 has yielded inconsistent results. Therefore, our aim was to longitudinally investigate associations between CYP2D6 (poor, intermediate, and normal) and CYP2C19 (poor, intermediate, normal, and ultrarapid) metabolizer-status, and switching/discontinuing of ADs. Next, we investigated whether the number of perceived side effects differed between metabolizer statuses. METHODS: Data came from the multi-site naturalistic longitudinal cohort Netherlands Study of Depression and Anxiety (NESDA). We selected depression- and/or anxiety patients, who used AD at some point in the course of the 9 years follow-up period (n = 928). Medication use was followed to assess patterns of AD switching/discontinuation over time. CYP2D6 and CYP2C19 alleles were derived using genome-wide data of the NESDA samples and haplotype data from the PharmGKB database. Logistic regression analyses were conducted to investigate the association of metabolizer status with switching/discontinuing ADs. Mann-Whitney U-tests were conducted to compare the number of patient-perceived side effects between metabolizer statuses. RESULTS: No significant associations were observed of CYP metabolizer status with switching/discontinuing ADs, nor with the number of perceived side effects. CONCLUSIONS: We found no evidence for associations between CYP metabolizer statuses and switching/discontinuing AD, nor with side effects of ADs, suggesting that metabolizer status only plays a limited role in switching/discontinuing ADs. Additional studies with larger numbers of PM and UM patients are needed to further determine the potential added value of pharmacogenetics to guide pharmacotherapy.

A genome-wide association study of hand eczema identifies locus 20q13.33 and reveals genetic overlap with atopic dermatitis.

BACKGROUND: Twin studies revealed that genetic effects play a role in hand eczema (HE), but the responsible genetic factors are unknown. OBJECTIVES: To identify and characterise genetic loci associated with HE and to provide insight into the genetic overlap between HE and atopic dermatitis (AD). METHODS: We used questionnaire-derived and genotype data from the European population-based Lifelines cohort and biobank. We performed a discovery genome-wide association study (GWAS) of HE (2879 cases and 16 249 controls) and of AD (1706 cases and 17 190 controls). We replicated our findings in an independent Lifelines sample for HE (1188 cases and 6431 controls) and AD (757 cases and 6747 controls). We conducted several post-GWAS analyses and performed genetic correlation analyses between our HE results and independent AD data. RESULTS: The two-step GWAS of HE, regardless of adjusting for AD, identified one independent locus 20q13.33, likely driven by a number of causal single-nucleotide polymorphisms. For the AD GWAS, we replicated a known stop-gained rs61816761 at locus 1q21.3 (FLG, FLGAS1). We found a strong genetic correlation (p < 0.01) between HE and AD (rg = 0.65), regardless of adjusting for AD (rg = 0.63). CONCLUSIONS: Locus 20q13.33 is associated with HE, and there is a large genetic overlap between HE and AD.

External Validation of Mortality Prediction Models for Critical Illness Reveals Preserved Discrimination but Poor Calibration.

OBJECTIVES: In a recent scoping review, we identified 43 mortality prediction models for critically ill patients. We aimed to assess the performances of these models through external validation. DESIGN: Multicenter study. SETTING: External validation of models was performed in the Simple Intensive Care Studies-I (SICS-I) and the Finnish Acute Kidney Injury (FINNAKI) study. PATIENTS: The SICS-I study consisted of 1,075 patients, and the FINNAKI study consisted of 2,901 critically ill patients. MEASUREMENTS AND MAIN RESULTS: For each model, we assessed: 1) the original publications for the data needed for model reconstruction, 2) availability of the variables, 3) model performance in two independent cohorts, and 4) the effects of recalibration on model performance. The models were recalibrated using data of the SICS-I and subsequently validated using data of the FINNAKI study. We evaluated overall model performance using various indexes, including the (scaled) Brier score, discrimination (area under the curve of the receiver operating characteristics), calibration (intercepts and slopes), and decision curves. Eleven models (26%) could be externally validated. The Acute Physiology And Chronic Health Evaluation (APACHE) II, APACHE IV, Simplified Acute Physiology Score (SAPS)-Reduced (SAPS-R)' and Simplified Mortality Score for the ICU models showed the best scaled Brier scores of 0.11' 0.10' 0.10' and 0.06' respectively. SAPS II, APACHE II, and APACHE IV discriminated best; overall discrimination of models ranged from area under the curve of the receiver operating characteristics of 0.63 (0.61-0.66) to 0.83 (0.81-0.85). We observed poor calibration in most models, which improved to at least moderate after recalibration of intercepts and slopes. The decision curve showed a positive net benefit in the 0-60% threshold probability range for APACHE IV and SAPS-R. CONCLUSIONS: In only 11 out of 43 available mortality prediction models, the performance could be studied using two cohorts of critically ill patients. External validation showed that the discriminative ability of APACHE II, APACHE IV, and SAPS II was acceptable to excellent, whereas calibration was poor.

Familial co-aggregation and shared genetics of cardiometabolic disorders and traits: data from the multi-generational Lifelines Cohort Study.

BACKGROUND: It is unclear to what extent genetics explain the familial clustering and the co-occurrence of distinct cardiometabolic disorders in the general population. We therefore aimed to quantify the familial (co-)aggregation of various cardiometabolic disorders and to estimate the heritability of cardiometabolic traits and their genetic correlations using the large, multi-generational Lifelines Cohort Study. METHODS: We used baseline data of 162,416 participants from Lifelines. Cardiometabolic disorders including type 2 diabetes (T2D), cardiovascular diseases, hypertension, obesity, hypercholesterolemia, and metabolic syndrome (MetS), were defined in adult participants. Fifteen additional cardiometabolic traits indexing obesity, blood pressure, inflammation, glucose regulation, and lipid levels were measured in all included participants. Recurrence risk ratios (λR) for first-degree relatives (FDR) indexed familial (co-)aggregation of cardiometabolic disorders using modified conditional Cox proportional hazards models and were compared to those of spouses. Heritability (h2), shared environment, and genetic correlation (rg) were estimated using restricted maximum likelihood variance decomposition methods, adjusted for age, age2, and sex. RESULTS: Individuals with a first-degree relative with a cardiometabolic disorder had a higher risk of the same disorder, ranging from λFDR of 1.23 (95% CI 1.20-1.25) for hypertension to λFDR of 2.48 (95% CI 2.15-2.86) for T2D. Most of these were higher than in spouses (λSpouses < λFDR), except for obesity which was slightly higher in spouses. We found moderate heritability for cardiometabolic traits (from h2CRP: 0.26 to h2HDL: 0.50). Cardiometabolic disorders showed positive familial co-aggregation, particularly between T2D, MetS, and obesity (from λFDR obesity-MetS: 1.28 (95% CI 1.24-1.32) to λFDR MetS-T2D: 1.61 (95% CI 1.52-1.70)), consistent with the genetic correlations between continuous intermediate traits (ranging from rg HDL-Triglycerides: - 0.53 to rg LDL-Apolipoprotein B: 0.94). CONCLUSIONS: There is positive familial (co-)aggregation of cardiometabolic disorder, moderate heritability of intermediate traits, and moderate genetic correlations between traits. These results indicate that shared genetics and common genetic architecture contribute to cardiometabolic disease.

Twenty-Five Novel Loci for Carotid Intima-Media Thickness: A Genome-Wide Association Study in >45 000 Individuals and Meta-Analysis of >100 000 Individuals.

OBJECTIVE: Carotid artery intima-media thickness (cIMT) is a widely accepted marker of subclinical atherosclerosis. Twenty susceptibility loci for cIMT were previously identified and the identification of additional susceptibility loci furthers our knowledge on the genetic architecture underlying atherosclerosis. APPROACH AND RESULTS: We performed 3 genome-wide association studies in 45 185 participants from the UK Biobank study who underwent cIMT measurements and had data on minimum, mean, and maximum thickness. We replicated 15 known loci and identified 20 novel loci associated with cIMT at P<5×10-8. Seven novel loci (ZNF385D, ADAMTS9, EDNRA, HAND2, MYOCD, ITCH/EDEM2/MMP24, and MRTFA) were identified in all 3 phenotypes. An additional new locus (LOXL1) was identified in the meta-analysis of the 3 phenotypes. Sex interaction analysis revealed sex differences in 7 loci including a novel locus (SYNE3) in males. Meta-analysis of UK Biobank data with a previous meta-analysis led to identification of three novel loci (APOB, FIP1L1, and LOXL4). Transcriptome-wide association analyses implicated additional genes ARHGAP42, NDRG4, and KANK2. Gene set analysis showed an enrichment in extracellular organization and the PDGF (platelet-derived growth factor) signaling pathway. We found positive genetic correlations of cIMT with coronary artery disease rg=0.21 (P=1.4×10-7), peripheral artery disease rg=0.45 (P=5.3×10-5), and systolic blood pressure rg=0.30 (P=4.0×10-18). A negative genetic correlation between average of maximum cIMT and high-density lipoprotein was found rg=-0.12 (P=7.0×10-4). CONCLUSIONS: Genome-wide association meta-analyses in >100 000 individuals identified 25 novel loci associated with cIMT providing insights into genes and tissue-specific regulatory mechanisms of proatherosclerotic processes. We found evidence for shared biological mechanisms with cardiovascular diseases.

Variants in the GPR146 Gene Are Associated With a Favorable Cardiometabolic Risk Profile.

BACKGROUND: In mice, GPR146 (G-protein-coupled receptor 146) deficiency reduces plasma lipids and protects against atherosclerosis. Whether these findings translate to humans is unknown. METHODS: Common and rare genetic variants in the GPR146 gene locus were used as research instruments in the UK Biobank. The Lifelines, The Copenhagen-City Heart Study, and a cohort of individuals with familial hypobetalipoproteinemia were used to find and study rare GPR146 variants. RESULTS: In the UK Biobank, carriers of the common rs2362529-C allele present with lower low-density lipoprotein cholesterol, apo (apolipoprotein) B, high-density lipoprotein cholesterol, apoAI, CRP (C-reactive protein), and plasma liver enzymes compared with noncarriers. Carriers of the common rs1997243-G allele, associated with higher GPR146 expression, present with the exact opposite phenotype. The associations with plasma lipids of the above alleles are allele dose-dependent. Heterozygote carriers of a rare coding variant (p.Pro62Leu; n=2615), predicted to be damaging, show a stronger reductions in the above parameters compared with carriers of the common rs2362529-C allele. The p.Pro62Leu variant is furthermore shown to segregate with low low-density lipoprotein cholesterol in a family with familial hypobetalipoproteinemia. Compared with controls, carriers of the common rs2362529-C allele show a marginally reduced risk of coronary artery disease (P=0.03) concomitant with a small effect size on low-density lipoprotein cholesterol (average decrease of 2.24 mg/dL in homozygotes) of this variant. Finally, mendelian randomization analyses suggest a causal relationship between GPR146 gene expression and plasma lipid and liver enzyme levels. CONCLUSIONS: This study shows that carriers of new genetic GPR146 variants have a beneficial cardiometabolic risk profile, but it remains to be shown whether genetic or pharmaceutical inhibition of GPR146 protects against atherosclerosis in humans.

Mediators of the association between educational attainment and type 2 diabetes mellitus: a two-step multivariable Mendelian randomisation study.

AIMS/HYPOTHESIS: Type 2 diabetes mellitus is a major health burden disproportionately affecting those with lower educational attainment (EA). We aimed to obtain causal estimates of the association between EA and type 2 diabetes and to quantify mediating effects of known modifiable risk factors. METHODS: We applied two-step, two-sample multivariable Mendelian randomisation (MR) techniques using SNPs as genetic instruments for exposure and mediators, thereby minimising bias due to confounding and reverse causation. We leveraged summary data on genome-wide association studies for EA, proposed mediators (i.e. BMI, blood pressure, smoking, television watching) and type 2 diabetes. The total effect of EA on type 2 diabetes was decomposed into a direct effect and indirect effects through multiple mediators. Additionally, traditional mediation analysis was performed in a subset of the National Health and Nutrition Examination Survey 2013-2014. RESULTS: EA was inversely associated with type 2 diabetes (OR 0.53 for each 4.2 years of schooling; 95% CI 0.49, 0.56). Individually, the largest contributors were BMI (51.18% mediation; 95% CI 46.39%, 55.98%) and television watching (50.79% mediation; 95% CI 19.42%, 82.15%). Combined, the mediators explained 83.93% (95% CI 70.51%, 96.78%) of the EA-type 2 diabetes association. Traditional analysis yielded smaller effects but showed consistent direction and priority ranking of mediators. CONCLUSIONS/INTERPRETATION: These results support a potentially causal protective effect of EA against type 2 diabetes, with considerable mediation by a number of modifiable risk factors. Interventions on these factors thus have the potential of substantially reducing the burden of type 2 diabetes attributable to low EA.

Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts.

AIMS/HYPOTHESIS: Type 2 diabetes is a complex metabolic disease with increasing prevalence worldwide. Improving the prediction of incident type 2 diabetes using epigenetic markers could help tailor prevention efforts to those at the highest risk. The aim of this study was to identify predictive methylation markers for incident type 2 diabetes by combining epigenome-wide association study (EWAS) results from five prospective European cohorts. METHODS: We conducted a meta-analysis of EWASs in blood collected 7-10 years prior to type 2 diabetes diagnosis. DNA methylation was measured with Illumina Infinium Methylation arrays. A total of 1250 cases and 1950 controls from five longitudinal cohorts were included: Doetinchem, ESTHER, KORA1, KORA2 and EPIC-Norfolk. Associations between DNA methylation and incident type 2 diabetes were examined using robust linear regression with adjustment for potential confounders. Inverse-variance fixed-effects meta-analysis of cohort-level individual CpG EWAS estimates was performed using METAL. The methylGSA R package was used for gene set enrichment analysis. Confirmation of genome-wide significant CpG sites was performed in a cohort of Indian Asians (LOLIPOP, UK). RESULTS: The meta-analysis identified 76 CpG sites that were differentially methylated in individuals with incident type 2 diabetes compared with control individuals (p values <1.1 × 10-7). Sixty-four out of 76 (84.2%) CpG sites were confirmed by directionally consistent effects and p values <0.05 in an independent cohort of Indian Asians. However, on adjustment for baseline BMI only four CpG sites remained genome-wide significant, and addition of the 76 CpG methylation risk score to a prediction model including established predictors of type 2 diabetes (age, sex, BMI and HbA1c) showed no improvement (AUC 0.757 vs 0.753). Gene set enrichment analysis of the full epigenome-wide results clearly showed enrichment of processes linked to insulin signalling, lipid homeostasis and inflammation. CONCLUSIONS/INTERPRETATION: By combining results from five European cohorts, and thus significantly increasing study sample size, we identified 76 CpG sites associated with incident type 2 diabetes. Replication of 64 CpGs in an independent cohort of Indian Asians suggests that the association between DNA methylation levels and incident type 2 diabetes is robust and independent of ethnicity. Our data also indicate that BMI partly explains the association between DNA methylation and incident type 2 diabetes. Further studies are required to elucidate the underlying biological mechanisms and to determine potential causal roles of the differentially methylated CpG sites in type 2 diabetes development.

Sex Hormones and Risk of Aneurysmal Subarachnoid Hemorrhage: A Mendelian Randomization Study.

BACKGROUND: The risk of aneurysmal subarachnoid hemorrhage (aSAH) is increased in postmenopausal women compared with men of similar age, suggesting a role for sex hormones. We aimed to explore whether sex hormones, and age at menarche/menopause have a causal effect on aSAH risk by conducting a 2-sample MR study (Mendelian randomization). METHODS: We obtained sex-specific genetic instruments for serum estradiol, bioavailable testosterone (BioT), SHBG (sex hormone-binding globulin), and age at menarche/menopause from genome-wide association studies. The associated sex-specific aSAH risk was estimated with inverse-variance weighted MR analyses with various statistical sensitivity analyses. Multivariable and cluster MR analyses were performed for BioT and SHBG to account for a genetic and phenotypic correlation between the 2 exposures. The clusters represented (1) single-nucleotide polymorphisms primarily increasing SHBG, with secondary decreasing effects on BioT, and (2) single-nucleotide polymorphisms affecting BioT without affecting SHBG. RESULTS: Univariable MR analyses showed an 18% increased aSAH risk among women per 1-SD increase in genetically determined SHBG levels (odds ratio, 1.18 [95% CI, 1.05-1.34]; P=0.007). Suggestive evidence was identified for a 27% decreased risk of aSAH among women per 1-SD increase in BioT (odds ratio, 0.73 [95% CI, 0.55-0.95]; P=0.02). The latter association disappeared in cluster analysis when only using SHBG-independent variants. MR analyses with variants from the cluster with primary SHBG effects and secondary (opposite) BioT-effects yielded a statistically significant association (odds ratio, 1.21 [95% CI, 1.05-1.40]; P=0.008). No other causal associations were identified. CONCLUSIONS: Genetic predisposition to elevated serum levels of SHBG, with secondary lower serum BioT levels, is associated with an increased aSAH risk among women, suggesting that SHBG and BioT causally elevate aSAH risk. Further studies are required to elucidate the underlying mechanisms and their potential as an interventional target to lower aSAH incidence.

Genomic Mismatch at LIMS1 Locus and Kidney Allograft Rejection.

BACKGROUND: In the context of kidney transplantation, genomic incompatibilities between donor and recipient may lead to allosensitization against new antigens. We hypothesized that recessive inheritance of gene-disrupting variants may represent a risk factor for allograft rejection. METHODS: We performed a two-stage genetic association study of kidney allograft rejection. In the first stage, we performed a recessive association screen of 50 common gene-intersecting deletion polymorphisms in a cohort of kidney transplant recipients. In the second stage, we replicated our findings in three independent cohorts of donor-recipient pairs. We defined genomic collision as a specific donor-recipient genotype combination in which a recipient who was homozygous for a gene-intersecting deletion received a transplant from a nonhomozygous donor. Identification of alloantibodies was performed with the use of protein arrays, enzyme-linked immunosorbent assays, and Western blot analyses. RESULTS: In the discovery cohort, which included 705 recipients, we found a significant association with allograft rejection at the LIMS1 locus represented by rs893403 (hazard ratio with the risk genotype vs. nonrisk genotypes, 1.84; 95% confidence interval [CI], 1.35 to 2.50; P = 9.8×10-5). This effect was replicated under the genomic-collision model in three independent cohorts involving a total of 2004 donor-recipient pairs (hazard ratio, 1.55; 95% CI, 1.25 to 1.93; P = 6.5×10-5). In the combined analysis (discovery cohort plus replication cohorts), the risk genotype was associated with a higher risk of rejection than the nonrisk genotype (hazard ratio, 1.63; 95% CI, 1.37 to 1.95; P = 4.7×10-8). We identified a specific antibody response against LIMS1, a kidney-expressed protein encoded within the collision locus. The response involved predominantly IgG2 and IgG3 antibody subclasses. CONCLUSIONS: We found that the LIMS1 locus appeared to encode a minor histocompatibility antigen. Genomic collision at this locus was associated with rejection of the kidney allograft and with production of anti-LIMS1 IgG2 and IgG3. (Funded by the Columbia University Transplant Center and others.).

GWASinspector: comprehensive quality control of genome-wide association study results.

SUMMARY: Quality control (QC) of genome wide association study (GWAS) result files has become increasingly difficult due to advances in genomic technology. The main challenges include continuous increases in the number of polymorphic genetic variants contained in recent GWASs and reference panels, the rising number of cohorts participating in a GWAS consortium, and inclusion of new variant types. Here, we present GWASinspector, a flexible R package for comprehensive QC of GWAS results. This package is compatible with recent imputation reference panels, handles insertion/deletion and multi-allelic variants, provides extensive QC reports and efficiently processes big data files. Reference panels covering three human genome builds (NCBI36, GRCh37 and GRCh38) are available. GWASinspector has a user friendly design and allows easy set-up of the QC pipeline through a configuration file. In addition to checking and reporting on individual files, it can be used in preparation of a meta-analysis by testing for systemic differences between studies and generating cleaned, harmonized GWAS files. Comparison with existing GWAS QC tools shows that the main advantages of GWASinspector are its ability to more effectively deal with insertion/deletion and multi-allelic variants and its relatively low memory use. AVAILABILITY AND IMPLEMENTATION: Our package is available at The Comprehensive R Archive Network (CRAN): https://CRAN.R-project.org/package=GWASinspector. Reference datasets and a detailed tutorial can be found at the package website at http://gwasinspector.com/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genomic analysis of diet composition finds novel loci and associations with health and lifestyle.

We conducted genome-wide association studies (GWAS) of relative intake from the macronutrients fat, protein, carbohydrates, and sugar in over 235,000 individuals of European ancestries. We identified 21 unique, approximately independent lead SNPs. Fourteen lead SNPs are uniquely associated with one macronutrient at genome-wide significance (P < 5 × 10-8), while five of the 21 lead SNPs reach suggestive significance (P < 1 × 10-5) for at least one other macronutrient. While the phenotypes are genetically correlated, each phenotype carries a partially unique genetic architecture. Relative protein intake exhibits the strongest relationships with poor health, including positive genetic associations with obesity, type 2 diabetes, and heart disease (rg ≈ 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|rg| ≈ 0.1-0.3) and positive genetic correlations with physical activity (rg ≈ 0.1 and 0.2). Relative fat intake has no consistent pattern of genetic correlations with poor health but has a negative genetic correlation with educational attainment (rg ≈-0.1). Although our analyses do not allow us to draw causal conclusions, we find no evidence of negative health consequences associated with relative carbohydrate, sugar, or fat intake. However, our results are consistent with the hypothesis that relative protein intake plays a role in the etiology of metabolic dysfunction.

Large-Scale Multi-Omics Studies Provide New Insights into Blood Pressure Regulation.

Recent genome-wide association studies uncovered part of blood pressure's heritability. However, there is still a vast gap between genetics and biology that needs to be bridged. Here, we followed up blood pressure genome-wide summary statistics of over 750,000 individuals, leveraging comprehensive epigenomic and transcriptomic data from blood with a follow-up in cardiovascular tissues to prioritise likely causal genes and underlying blood pressure mechanisms. We first prioritised genes based on coding consequences, multilayer molecular associations, blood pressure-associated expression levels, and coregulation evidence. Next, we followed up the prioritised genes in multilayer studies of genomics, epigenomics, and transcriptomics, functional enrichment, and their potential suitability as drug targets. Our analyses yielded 1880 likely causal genes for blood pressure, tens of which are targets of the available licensed drugs. We identified 34 novel genes for blood pressure, supported by more than one source of biological evidence. Twenty-eight (82%) of these new genes were successfully replicated by transcriptome-wide association analyses in a large independent cohort (n = ~220,000). We also found a substantial mediating role for epigenetic regulation of the prioritised genes. Our results provide new insights into genetic regulation of blood pressure in terms of likely causal genes and involved biological pathways offering opportunities for future translation into clinical practice.

Validating the doubly weighted genetic risk score for the prediction of type 2 diabetes in the Lifelines and Estonian Biobank cohorts.

As many cases of type 2 diabetes (T2D) are likely to remain undiagnosed, better tools for early detection of high-risk individuals are needed to prevent or postpone the disease. We investigated the value of the doubly weighted genetic risk score (dwGRS) for the prediction of incident T2D in the Lifelines and Estonian Biobank (EstBB) cohorts. The dwGRS uses an additional weight for each single nucleotide polymorphism in the risk score, to correct for "Winner's curse" bias in the effect size estimates. The traditional (single-weighted genetic risk score; swGRS) and dwGRS were calculated for participants in Lifelines (n = 12,018) and EstBB (n = 34,129). The dwGRS was found to have stronger association with incident T2D (hazard ratio [HR] = 1.26 [95% confidence interval: 1.10-1.43] and HR = 1.35 [1.28-1.42]) compared to the swGRS (HR = 1.21 [1.07-1.38] and HR = 1.25 [1.19-1.32]) in Lifelines and EstBB, respectively. Comparing the 5-year predicted risks from the models with and without the dwGRS, the continuous net reclassification index was 0.140 (0.034-0.243; p = .009 Lifelines), and 0.257 (0.194-0.319; p < 2 × 10-16 EstBB). The dwGRS provided incremental value to the T2D prediction model with established phenotypic predictors. It clearly distinguished the risk groups for incident T2D in both biobanks thereby showing its clinical relevance.

Genome-wide CNV investigation suggests a role for cadherin, Wnt, and p53 pathways in primary open-angle glaucoma.

BACKGROUND: To investigate whether copy number variations (CNVs) are implicated in molecular mechanisms underlying primary open-angle glaucoma (POAG), we used genotype data of POAG individuals and healthy controls from two case-control studies, AGS (n = 278) and GLGS-UGLI (n = 1292). PennCNV, QuantiSNP, and cnvPartition programs were used to detect CNV. Stringent quality controls at both sample and marker levels were applied. The identified CNVs were intersected in CNV region (CNVR). After, we performed burden analysis, CNV-genome-wide association analysis, gene set overrepresentation and pathway analysis. In addition, in human eye tissues we assessed the expression of the genes lying within significant CNVRs. RESULTS: We reported a statistically significant greater burden of CNVs in POAG cases compared to controls (p-value = 0,007). In common between the two cohorts, CNV-association analysis identified statistically significant CNVRs associated with POAG that span 11 genes (APC, BRCA2, COL3A1, HLA-DRB1, HLA-DRB5, HLA-DRB6, MFSD8, NIPBL, SCN1A, SDHB, and ZDHHC11). Functional annotation and pathway analysis suggested the involvement of cadherin, Wnt signalling, and p53 pathways. CONCLUSIONS: Our data suggest that CNVs may have a role in the susceptibility of POAG and they can reveal more information on the mechanism behind this disease. Additional genetic and functional studies are warranted to ascertain the contribution of CNVs in POAG.

Associations of Genetic Factors, Educational Attainment, and Their Interaction With Kidney Function Outcomes.

Both genetic predisposition and low educational attainment (EA) are associated with higher risk of chronic kidney disease. We examined the interaction of EA and genetic risk in kidney function outcomes. We included 3,597 participants from the Prevention of Renal and Vascular End-Stage Disease Cohort Study, a longitudinal study in a community-based sample from Groningen, the Netherlands (median follow-up, 11 years; 1997-2012). Kidney function was approximated by obtaining estimated glomerular filtration rate (eGFR) from serum creatinine and cystatin C. Individual longitudinal linear eGFR trajectories were derived from linear mixed models. Genotype data on 63 single-nucleotide polymorphisms, with known associations with eGFR, were used to calculate an allele-weighted genetic score (WGS). EA was categorized into high, medium, and low. In ordinary least squares analysis, higher WGS and lower EA showed additive effects on reduced baseline eGFR; the interaction term was nonsignificant. In analysis of eGFR decline, the significant interaction term suggested amplification of genetic risk by low EA. Adjustment for known renal risk factors did not affect our results. This study presents the first evidence of gene-environment interaction between EA and a WGS for eGFR decline and provides population-level insights into the mechanisms underlying socioeconomic disparities in chronic kidney disease.