A genome-wide association study for venous thromboembolism: the extended cohorts for heart and aging research in genomic epidemiology (CHARGE) consortium

Thrombosis-related characteristics of pregnant women with antithrombin deficiency, protein C deficiency and protein S deficiency in Japan

BACKGROUND

 We previously conducted a primary survey of pregnant women with hereditary thrombophilia based on national surveillance in Japan, but did not examine their thrombosis-related characteristics. Antithrombin (AT) deficiency, protein C (PC) deficiency and protein S (PS) deficiency are the major types of hereditary thrombophilia in Japan.

METHODS

We examined their detailed information related to thrombosis, and evaluated peripartum outcomes in comparison with control data obtained from the Japan Society of Obstetrics and Gynecology.

RESULTS

Definite or probable AT deficiency, PC deficiency and PS deficiency were observed in 80, 50, and 317 pregnancies, respectively, from 2014 to 2018 in Japan, with prevalence rates among total deliveries of 0.011%, 0.007%, 0.044%. The number of pregnancies with AT, PC and PS deficiency might have been as many as 27, 17 and 108 every year if complete answers had been provided. In the peripartum period of current pregnancies, 27.5% of women with AT deficiency, 28.0% with PC deficiency and 13.2% with PS deficiency developed thrombosis (p < 0.001 vs. control). Pregnant women with AT and PC deficiency were more susceptible to thrombosis than those with PS deficiency (P < 0.01). Of the thromboses, 92.3% occurred during pregnancy, 62.8% at less than 15 gestational weeks. The earliest onset of thrombosis was 5 gestational weeks. Prophylactic anticoagulation significantly prevented the onset of both antepartum and postpartum thrombosis (p < 0.0001). The rate of recurrent pregnancy loss in women with low PC or PS activities was significantly higher than in controls (p < 0.0001); however, it is unknown whether recurrent pregnancy loss is related to hereditary PS deficiency. There seem to have been few serious maternal or fetal/neonatal complications due to placental insufficiency related to a hypercoagulable state other than growth restriction.

CONCLUSIONS

This survey revealed the thrombosis-related characteristics of pregnant women with hereditary thrombophilia in Japan. We suggest prophylactic anticoagulation to prevent maternal or fetal/neonatal complications.

Association between venous thromboembolism-associated genetic variants, coagulation factor levels, and thrombin generation potential

Recently three large meta-analyses of genome-wide association studies for venous thromboembolism (VTE) identified over 130 genetic variants. However, mechanisms by which newly identified and therefore underexplored VTE-associated genetic variants influence VTE remain unclear. To elucidate the mechanism, we investigated the association between 61 newly identified VTE-associated genetic variants and the levels of coagulation factor (F) VIII, FIX, FXI, and fibrinogen as well as thrombin generation parameters (lag time, peak, endogenous thrombin potential, time-to-peak, and velocity), which are well-known biological traits associated with VTE. This study was conducted on 5341 participants of the Netherlands Epidemiology of Obesity study. The associations between VTE-associated genetic variants and coagulation factor levels and thrombin generation parameters were examined using linear regression analyses, adjusted for age, sex, body mass index, oral contraceptive use, hormone replacement therapy, and menopausal status. Of 61 genetic variants, 33 were associated with one or more of the coagulation factor levels and thrombin generation parameters. Following multiple testing corrections, five genetic variants remained significant, of which MAP1A rs55707100 exhibited the most robust association with thrombin generation parameters and FXI levels (β = -5.33%, 95% confidence interval: -8.44, -2.22). Our findings shed light on the underlying mechanisms by which these genetic variants influence the risk of VTE.

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 GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). 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 may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

Genomic science of risk prediction for venous thromboembolic disease: convenient clarification or compounding complexity

Venous thromboembolism (VTE) refers to abnormal blood clots in veins occurring in 1 to 2 per 1000 individuals every year. While anticoagulant treatment can prevent VTE, it increases the risk of bleeding. This emphasizes the importance of identifying individuals with a high risk of VTE and providing prophylactic interventions to these individuals to reduce both VTE and bleeding risks. Current risk assessment of VTE is based on the combination of mainly clinical risk factors. With the identification of an increasing number of genetic variants associated with the risk of VTE, the addition of genetic findings to clinical prediction models can improve risk prediction for VTE. Especially for individuals in high-risk situations, the added value of genetic findings to clinical prediction models may have benefits such as better prophylaxis of VTE and the reduced side effects of bleeding from unnecessary treatment. Nevertheless, the question of whether these models will eventually have clinical utility remains to be proven. Here, we review the current state of knowledge on genetic risk factors for VTE, explore genetic prediction models for VTE, and discuss their clinical implications and challenges.

Thrombogenesis-associated genetic determinants as predictors of thromboembolism and prognosis in cervical cancer

Venous thromboembolism (VTE) is a leading cause of death among cancer patients. Khorana score (KS) is the most studied tool to predict cancer-related VTE, however, it exerts poor sensitivity. Several single-nucleotide polymorphisms (SNPs) have been associated with VTE risk in the general population, but whether they are predictors of cancer-related VTE is a matter of discussion. Compared to other solid tumours, little is known about VTE in the setting of cervical cancer (CC) and whether thrombogenesis-related polymorphisms could be valuable biomarkers in patients with this neoplasia. This study aims to analyse the effect of VTE occurrence on the prognosis of CC patients, explore the predictive capability of KS and the impact of thrombogenesis-related polymorphisms on CC-related VTE incidence and patients' prognosis regardless of VTE. A profile of eight SNPs was evaluated. A retrospective hospital-based cohort study was conducted with 400 CC patients under chemoradiotherapy. SNP genotyping was carried on by using TaqMan® Allelic Discrimination methodology. Time to VTE occurrence and overall survival were the two measures of clinical outcome evaluated. The results indicated that VTE occurrence (8.5%) had a significant impact on the patient's survival (log-rank test, P < 0.001). KS showed poor performance (KS ≥ 3, χ², P = 0.191). PROCR rs10747514 and RGS7 rs2502448 were significantly associated with the risk of CC-related VTE development (P = 0.021 and P = 0.006, respectively) and represented valuable prognostic biomarkers regardless of VTE (P = 0.004 and P = 0.010, respectively). Thus, thrombogenesis-related genetic polymorphisms may constitute valuable biomarkers among CC patients allowing a more personalized clinical intervention.

A Clinical-Genetic Risk Score for Predicting Cancer-Associated Venous Thromboembolism: A Development and Validation Study Involving Two Independent Prospective Cohorts

PURPOSE

Venous thromboembolism (VTE) is a leading cause of death among patients with cancer. The Khorana score was developed for assessing the risk of VTE in outpatients with cancer receiving chemotherapy, but its accuracy in identifying patients at high risk has been questioned. The aim of this study was to develop and validate a clinical-genetic score that improves the assessment of VTE risk in oncology outpatients within 6 months of diagnosis.

METHODS

The new score was developed using the data of 364 outpatients belonging to the Spanish ONCOTHROMB 12-01 population. In this cohort, clinical data associated with the risk of VTE were collected at the time of diagnosis, including the Khorana score. These patients were also genotyped for the 51 genetic variants known to be associated with VTE. Multivariate logistic regression was performed to determine the weight of each genetic and clinical variable in relation to VTE risk, allowing a clinical-genetic risk score (the ONCOTHROMB score) to be developed. The Khorana and the ONCOTHROMB scores were then compared via the area under the receiver operating characteristic curve (AUC), calibration, and the number of patients needed to treat. The new score was then validated in a study of 263 patients in the Vienna Cancer and Thrombosis Study population.

RESULTS

Nine genetic variants, tumor site, TNM stage, and a body mass index of > 25 kg/m2 were found to be associated with VTE and were used to build the ONCOTHROMB score, which better predicted the overall risk of VTE than did the Khorana score (AUC, 0.781 v 0.580; P < .001). Similar AUC results were recorded in the validation study the Vienna Cancer and Thrombosis Study cohort involving patients with the same type of tumor (AUC for the ONCOTHROMB score v the Khorana score: 0.686 v 0.577; P < .001) and with all type of tumors (AUC for the ONCOTHROMB score v the Khorana score: 0.720 v 0.561; P < .0001).

CONCLUSION

The ONCOTHROMB score for VTE risk in outpatients with cancer, which takes into account both clinical and genetic variables, better identifies patients who might benefit from primary thromboprophylaxis than does the Khorana score.

Genome-wide association analyses identified novel susceptibility loci for pulmonary embolism among Han Chinese population

BACKGROUND

A large proportion of pulmonary embolism (PE) heritability remains unexplained, particularly among the East Asian (EAS) population. Our study aims to expand the genetic architecture of PE and reveal more genetic determinants in Han Chinese.

METHODS

We conducted the first genome-wide association study (GWAS) of PE in Han Chinese, then performed the GWAS meta-analysis based on the discovery and replication stages. To validate the effect of the risk allele, qPCR and Western blotting experiments were used to investigate possible changes in gene expression. Mendelian randomization (MR) analysis was employed to implicate pathogenic mechanisms, and a polygenic risk score (PRS) for PE risk prediction was generated.

RESULTS

After meta-analysis of the discovery dataset (622 cases, 8853 controls) and replication dataset (646 cases, 8810 controls), GWAS identified 3 independent loci associated with PE, including the reported loci FGG rs2066865 (p-value = 3.81 × 10^-14), ABO rs582094 (p-value = 1.16 × 10^-10) and newly reported locus FABP2 rs1799883 (p-value = 7.59 × 10^-17). Previously reported 10 variants were successfully replicated in our cohort. Functional experiments confirmed that FABP2-A163G(rs1799883) promoted the transcription and protein expression of FABP2. Meanwhile, MR analysis revealed that high LDL-C and TC levels were associated with an increased risk of PE. Individuals with the top 10% of PRS had over a fivefold increased risk for PE compared to the general population.

CONCLUSIONS

We identified FABP2, related to the transport of long-chain fatty acids, contributing to the risk of PE and provided more evidence for the essential role of metabolic pathways in PE development.

Whole Genome Analysis of Venous Thromboembolism: the Trans-Omics for Precision Medicine Program

BACKGROUND

Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies.

METHODS

The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants).

RESULTS

Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P=7.4×10-14) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P=1.6×10-14), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P=1.8×10-6 with the secondary filter), while SERPINC1 did not (minimum P=4.4×10-5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1, became significant (P=4.4×10-7 using all missense variants with minor allele frequency <0.0005).

CONCLUSIONS

Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.

Genomewide Association Study Identifies Copy Number Variants Associated With Warfarin Dose Response and Risk of Venous Thromboembolism in African Americans

The anticoagulant warfarin is commonly used to control and prevent thrombotic disorders, such as venous thromboembolism (VTE), which disproportionately afflicts African Americans. Despite the importance of copy number variants (CNVs), few studies have focused on characterizing and understanding their role in drug response and disease risk among African Americans. In this study, we conduct the first genome-wide analysis of CNVs to more comprehensively account for the contribution of genetic variation in warfarin dose requirement and VTE risk among African Americans. We used hidden Markov models to detect CNVs from high-throughput single-nucleotide polymorphism arrays for 340 African American participants in the International Warfarin Pharmacogenetics Consortium. We identified 11,570 CNVs resulting in 2,038 copy number variable regions (CNVRs) and found 3 CNVRs associated with warfarin dose requirement and 3 CNVRs associated with VTE risk in African Americans. CNVRs 1q31.2del and 6q14.1del were associated with increased warfarin dose requirement (β = 11.18 and 4.94, respectively; P_emp  = < 0.002); CNVR 19p13.31del was associated with decreased warfarin dose requirement (β = -1.41, P_emp  = 0.0004); CNVRs (2p22.1del and 5q35.1-q35.2del) were found to be associated with increased risk of VTE (odds ratios (ORs) = 1.88 and 14.9, respectively; P_emp  ≤0.02); and CNVR 10q26.12del was associated with a decreased risk of VTE (OR = 0.6; P_emp  = 0.05). Modeling of the 10q26.12del in HepG2 cells revealed that this deletion results in decreased fibrinogen gene expression, decreased fibrinogen and WDR11 protein levels, and decreased secretion of fibrinogen into the extracellular matrix. We found robust evidence that CNVRs could contribute to warfarin dose requirement and risk of VTE in African Americans and for 10q26.3del describe a possible pathogenic mechanism.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors

BACKGROUND

Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources.

METHODS

We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations.

RESULTS

In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis.

CONCLUSIONS

Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Peripartum management of hereditary thrombophilia: results of primary surveillance in Japan

This study investigated patients with thrombophilia and current peripartum management practices based on national surveillance in Japan. Between 2014 and 2018, antithrombin (AT), protein C (PC) and protein S (PS) deficiency were observed in 84, 67, and 443 pregnancies, respectively, with incidence rates among total deliveries at 0.012%, 0.009%, and 0.061%. The percentage of institutions that measured both antigens and AT, PC, and PS activity for the diagnosis of thrombophilia was 50.2%, and 46.9% of institutions did not perform gene analysis. Prophylactic anticoagulation therapy was used in the ante- and postpartum management of patients with AT deficiency at 67.1% and 66.3% of institutions, most commonly with 10,000 units of unfractionated heparin. Ante- and postpartum management of PC and PS deficiency was performed at 75.3% and 67.1% of institutions. Approximately half of the institutions performed peripartum prophylactic AT supplementation for AT deficiency. Low trough AT activity before supplementation was most commonly 50 ≤  < 70%, and the highest AT supplementation was 1500 ≤  < 3000 units. The number of pregnancies with AT, PC and PS deficiency might be as many as 29, 23 and 151 every year in Japan if complete answers were provided.

Oral postmenopausal hormone therapy and genetic risk on venous thromboembolism: gene-hormone interaction results from a large prospective cohort study

OBJECTIVE

Oral postmenopausal hormone therapy (HT) has been shown to be associated with venous thromboembolism (VTE), but whether this association is modified by VTE-associated genetic susceptibility is unknown. We examined interactions between oral HT use and a genetic risk score (GRS) of VTE.

METHOD

Eligible women were postmenopausal women who had data on oral HT use, VTE incidence between 1990 and 2012, and genetic data in the Nurses' Health Study. We built a GRS aggregating 16 VTE-related genetic variants. We used Cox regression to estimate associations of HT use with incident VTE and assessed interactions between HT use and VTE GRS. We also estimated incidence of VTE between age 50 and 79 years for groups of women defined by HT use and VTE GRS.

RESULTS

We identified 432 incident VTE cases. Current HT users were at higher risk of VTE than never users (HR: 1.9, 95% CI: 1.5-2.6), with slightly higher risk for estrogen plus progestin HT than estrogen only (HR: 2.4 vs 1.9). The GRS was associated with VTE risk (HR comparing 4th quartile to 1st: 2.0, 95% CI: 1.2-3.4). We did not observe significant multiplicative interactions between HT use and GRS. The estimated VTE risk difference (per 10,000 person-years) comparing 50-year-old current HT users to never users was 22.5 for women in the highest GRS quartile and 9.8 for women in the lowest GRS quartile.

CONCLUSION

The VTE GRS might inform clinical guidance regarding the balance of risks and benefits of HT use, especially among younger women.

Multivariate GWAS of Structural Dental Anomalies and Dental Caries in a Multi-Ethnic Cohort

Odontogenesis is a complex process, where disruption can result in dental anomalies and/or increase the risk of developing dental caries. Based on previous studies, certain dental anomalies tend to co-occur in patients, suggesting that these traits may share common genetic and etiological components. The main goal of this study was to implement a multivariate genome-wide association study approach to identify genetic variants shared between correlated structural dental anomalies and dental caries. Our cohort (N = 3,579) was derived from the Pittsburgh Orofacial Clefts Study, where multiple dental traits were assessed in both the unaffected relatives of orofacial cleft (OFC) cases (n = 2,187) and unaffected controls (n = 1,392). We identified four multivariate patterns of correlated traits in this data: tooth agenesis, impaction, and rotation (AIR); enamel hypoplasia, displacement, and rotation (HDR); displacement, rotation, and mamelon (DRM); and dental caries, tooth agenesis and enamel hypoplasia (CAH). We analyzed each of these four models using genome-wide multivariate tests of association. No genome-wide statistically significant results were found, but we identified multiple suggestive association signals (P < 10-5) near genes with known biological roles during tooth development, including ADAMTS9 and PRICKLE2 associated with AIR; GLIS3, WDR72, and ROR2 associated with HDR and DRM; ROBO2 associated with DRM; BMP7 associated with HDR; and ROBO1, SMAD2, and MSX2 associated with CAH. This is the first study to investigate genetic associations for multivariate patterns of correlated dental anomalies and dental caries. Further studies are needed to replicate these results in independent cohorts.

Today's view of hereditary thrombophilia

Venous thromboembolism (VTE) is still a serious medical problem with the non-decreasing incidence of new cases despite prophylaxis in risky situations. It is a multifactorial disease, in which the hereditary component is also significantly involved. The aim of the current research is to search for new polymorphisms that are involved in thrombogenesis in addition to classical thrombophilia (deficiency of natural coagulation inhibitors and FVL and FII prothrombin mutations). The article provides an overview of the results of already performed genome-wide association studies of VTE and their use for the calculation of the so-called polygenic risk score, which could be used for individualized prevention of VTE after standardization of the method.

Polygenic risk scores: An overview from bench to bedside for personalised medicine

Since the first polygenic risk score (PRS) in 2007, research in this area has progressed significantly. The increasing number of SNPs that have been identified by large scale GWAS analyses has fuelled the development of a myriad of PRSs for a wide variety of diseases and, more recently, to PRSs that potentially identify differential response to specific drugs. PRSs constitute a composite genomic biomarker and potential applications for PRSs in clinical practice encompass risk prediction and disease screening, early diagnosis, prognostication, and drug stratification to improve efficacy or reduce adverse drug reactions. Nevertheless, to our knowledge, no PRSs have yet been adopted into routine clinical practice. Beyond the technical considerations of PRS development, the major challenges that face PRSs include demonstrating clinical utility and circumnavigating the implementation of novel genomic technologies at scale into stretched healthcare systems. In this review, we discuss progress in developing disease susceptibility PRSs across multiple medical specialties, development of pharmacogenomic PRSs, and future directions for the field.

An assessment of the effect of the genotype on postoperative venous thromboembolism risk in 140,831 surgical patients

Background

Postoperative Venous Thromboembolism Events (VTE) constitute a major source of morbidity and mortality after surgery. The aim of this study was to investigate whether commonly occurring Single Nucleotide Polymorphisms (SNPs) are associated with VTE in the surgical setting.

Methods

Retrospective study using data from the United Kingdom (UK) biobank, a genome biobank containing healthcare and genotyping data from more than 500.000 individuals. A cohort of 140,831patients with a registered surgical procedure was identified and used for a discovery genome wide association study (GWAS), with the remainder of the cohort (305,349 non-surgical patients) used as a replication cohort. Primary outcome was associations between SNPs and VTE within 30 days after a surgical procedure. Genome wide significance was set at p = 5 × 10⁻⁸.

Results

In the surgical (discovery) cohort, no SNPs reached genome wide significance. The VTE association of the top candidate SNP in the ABO gene rs505922 (p = 3.33 × 10⁻⁷), was replicated in the general (replication) cohort (p = 2.42 × 10⁻⁵⁹).

Conclusions

and Relevance: This study did not identify associations between SNPs and postoperative VTE events reaching genome-wide significance, although the VTE relevance of top candidates were demonstrated.

•Venous thromboembolisms are common after surgery.

•Large biobanks help identify common genetic risk factors.

•Common variants in ABO may increase the risk of postoperative venous thromboembolisms.

•Genotyping surgical patients may benefit preoperative risk assessment.

Prediction of primary venous thromboembolism based on clinical and genetic factors within the U.K. Biobank

Both clinical and genetic factors drive the risk of venous thromboembolism. However, whether clinically recorded risk factors and genetic variants can be combined into a clinically applicable predictive score remains unknown. Using Cox proportional-hazard models, we analyzed the association of risk factors with the likelihood of venous thromboembolism in U.K. Biobank, a large prospective cohort. We then created a polygenic risk score of 36 single nucleotide polymorphisms and a clinical score determined by age, sex, body mass index, previous cancer diagnosis, smoking status, and fracture in the last 5 years. Participants were at significantly increased risk of venous thromboembolism if they were at high clinical risk (subhazard ratio, 4.37 [95% CI, 3.85-4.97]) or high genetic risk (subhazard ratio, 3.02 [95% CI, 2.63-3.47]) relative to participants at low clinical or genetic risk, respectively. The combined model, consisting of clinical and genetic components, was significantly better than either the clinical or the genetic model alone (P < 0.001). Participants at high risk in the combined score had nearly an eightfold increased risk of venous thromboembolism relative to participants at low risk (subhazard ratio, 7.51 [95% CI, 6.28-8.98]). This risk score can be used to guide decisions regarding venous thromboembolism prophylaxis, although external validation is needed.

Identification of transcriptome alterations in the prefrontal cortex, hippocampus, amygdala and hippocampus of suicide victims

Suicide is one of the leading causes of death globally for all ages, and as such presents a very serious problem for clinicians worldwide. However, the underlying neurobiological pathology remains to a large extent unknown. In order to address this gap, we have carried out a genome-wide investigation of the gene expression in the amygdala, hippocampus, prefrontal cortex and thalamus in post-mortem brain samples obtained from 20 suicide completers and 7 control subjects. By KEGG enrichment analysis indicated we identified novel clusters of downregulated pathways involved in antigen neutralization and autoimmune thyroid disease (amygdala, thalamus), decreased axonal plasticity in the hippocampus. Two upregulated pathways were involved in neuronal death in the hippocampus and olfactory transduction in the thalamus and the prefrontal cortex. Autoimmune thyroid disease pathway was downregulated only in females. Metabolic pathways involved in Notch signaling amino acid metabolism and unsaturated lipid synthesis were thalamus-specific. Suicide-associated changes in the expression of several genes and pseudogenes that point to various functional mechanisms possibly implicated in the pathology of suicide. Two genes (SNORA13 and RNU4-2) involved in RNA processing were common to all brain regions analyzed. Most of the identified gene expression changes were related to region-specific dysregulated manifestation of genetic and epigenetic mechanisms underlying neurodevelopmental disorders (SNORD114-10, SUSd1), motivation, addiction and motor disorders (CHRNA6), long-term depression (RAB3B), stress response, major depression and schizophrenia (GFAP), signal transduction at the neurovascular unit (NEXN) and inhibitory neurotransmission in spatial learning, neural plasticity (CALB2; CLIC6, ENPP1). Some of the differentially expressed genes were brain specific non-coding RNAs involved in the regulation of translation (SNORA13). One, (PARM1) is a potential oncogene and prognostic biomarker for colorectal cancer with no known function in the brain. Disturbed gene expression involved in antigen neutralization, autoimmunity, neural plasticity, stress response, signal transduction at the neurovascular unit, dysregulated nuclear RNA processing and translation and epigenetic imprinting signatures is associated with suicide and point to regulatory non-coding RNAs as potential targets of new drugs development.

Donor genetic variants as risk factors for thrombosis after liver transplantation: A genome-wide association study

Thrombosis after liver transplantation substantially impairs graft- and patient survival. Inevitably, heritable disorders of coagulation originating in the donor liver are transmitted by transplantation. We hypothesized that genetic variants in donor thrombophilia genes are associated with increased risk of posttransplant thrombosis. We genotyped 775 donors for adult recipients and 310 donors for pediatric recipients transplanted between 1993 and 2018. We determined the association between known donor thrombophilia gene variants and recipient posttransplant thrombosis. In addition, we performed a genome-wide association study (GWAS) and meta-analyzed 1085 liver transplantations. In our donor cohort, known thrombosis risk loci were not associated with posttransplant thrombosis, suggesting that it is unnecessary to exclude liver donors based on thrombosis-susceptible polymorphisms. By performing a meta-GWAS from children and adults, we identified 280 variants in 55 loci at suggestive genetic significance threshold. Downstream prioritization strategies identified biologically plausible candidate genes, among which were AK4 (rs11208611-T, p = 4.22 × 10^-05 ) which encodes a protein that regulates cellular ATP levels and concurrent activation of AMPK and mTOR, and RGS5 (rs10917696-C, p = 2.62 × 10^-05 ) which is involved in vascular development. We provide evidence that common genetic variants in the donor, but not previously known thrombophilia-related variants, are associated with increased risk of thrombosis after liver transplantation.

Clinical guidance for peripartum management of patients with hereditary thrombophilia

Hereditary thrombophilia is a condition in which individuals are susceptible to the formation of thrombi due to a hereditary deficiency in anticoagulant factors, antithrombin (AT), protein C (PC), or protein S (PS). Many Japanese thrombophilia patients have PS deficiency, especially PS p.K196E (also called as PS Tokushima), which is exclusive to the Japanese population, and thrombosis sometimes occurs during pregnancy. At present, no management guidelines for pregnancy and delivery in thrombophilia patients have been developed. The Study Group for Hereditary Thrombophilia, one of the research groups of blood coagulation abnormalities in the Research Program on Rare and Intractable Diseases supported with the Research Grants of the Ministry of Health, Labour and Welfare Science, has therefore developed this clinical guidance to provide healthcare workers with necessary information on safe pregnancy, parturition and neonatal management, adopting a format of responses to seven clinical questions (CQ). At the end of each answer, the corresponding Recommendation Level (A, B, C) is indicated.

Cis-Segregation of c.1171C>T Stop Codon (p.R391*) in SERPINC1 Gene and c.1691G>A Transition (p.R506Q) in F5 Gene and Selected GWAS Multilocus Approach in Inherited Thrombophilia

Inherited thrombophilia (e.g., venous thromboembolism, VTE) is due to rare loss-of-function mutations in anticoagulant factors genes (i.e., SERPINC1, PROC, PROS1), common gain-of-function mutations in procoagulant factors genes (i.e., F5, F2), and acquired risk conditions. Genome Wide Association Studies (GWAS) recently recognized several genes associated with VTE though gene defects may unpredictably remain asymptomatic, so calculating the individual genetic predisposition is a challenging task. We investigated a large family with severe, recurrent, early-onset VTE in which two sisters experienced VTE during pregnancies characterized by a perinatal in-utero thrombosis in the newborn and a life-saving pregnancy-interruption because of massive VTE, respectively. A nonsense mutation (CGA > TGA) generating a premature stop-codon (c.1171C>T; p.R391*) in the exon 6 of SERPINC1 gene (1q25.1) causing Antithrombin (AT) deficiency and the common missense mutation (c.1691G>A; p.R506Q) in the exon 10 of F5 gene (1q24.2) (i.e., FV Leiden; rs6025) were coinherited in all the symptomatic members investigated suspecting a cis-segregation further confirmed by STR-linkage-analyses [i.e., SERPINC1 IVS5 (ATT)_5–18, F5 IVS2 (AT)_6–33 and F5 IVS11 (GT)_12–16] and SERPINC1 intragenic variants (i.e., rs5878 and rs677). A multilocus investigation of blood-coagulation balance genes detected the coexistence of FV Leiden (rs6025) in trans with FV HR2-haplotype (p.H1299R; rs1800595) in the aborted fetus, and F11 rs2289252, F12 rs1801020, F13A1 rs5985, and KNG1 rs710446 in the newborn and other members. Common selected gene variants may strongly synergize with less common mutations tuning potential life-threatening conditions when combined with rare severest mutations. Merging classic and newly GWAS-identified gene markers in at risk families is mandatory for VTE risk estimation in the clinical practice, avoiding partial risk score evaluation in unrecognized at risk patients.

Epidemiology and Genetics of Venous Thromboembolism and Chronic Venous Disease

Venous disease is a term that broadly covers both venous thromboembolic disease and chronic venous disease. The basic pathophysiology of venous thromboembolism and chronic venous disease differ as venous thromboembolism results from an imbalance of hemostasis and thrombosis while chronic venous disease occurs in the setting of tissue damage because of prolonged venous hypertension. Both diseases are common and account for significant mortality and morbidity, respectively, and collectively make up a large health care burden. Despite both diseases having well-characterized environmental components, it has been known for decades that family history is an important risk factor, implicating a genetic element to a patient's risk. Our understanding of the pathogenesis of these diseases has greatly benefited from an expansion of population genetic studies from pioneering familial studies to large genome-wide association studies; we now have multiple risk loci for each venous disease. In this review, we will highlight the current state of knowledge on the epidemiology and genetics of venous thromboembolism and chronic venous disease and directions for future research.

Robust Performance of Potentially Functional SNPs in Machine Learning Models for the Prediction of Atorvastatin-Induced Myalgia

Statins can cause muscle symptoms resulting in poor adherence to therapy and increased cardiovascular risk. We hypothesize that combinations of potentially functional SNPs (pfSNPs), rather than individual SNPs, better predict myalgia in patients on atorvastatin. This study assesses the value of potentially functional single nucleotide polymorphisms (pfSNPs) and employs six machine learning algorithms to identify the combination of SNPs that best predict myalgia.

Methods: Whole genome sequencing of 183 Chinese, Malay and Indian patients from Singapore was conducted to identify genetic variants associated with atorvastatin induced myalgia. To adjust for confounding factors, demographic and clinical characteristics were also examined for their association with myalgia. The top factor, sex, was then used as a covariate in the whole genome association analyses. Variants that were highly associated with myalgia from this and previous studies were extracted, assessed for potential functionality (pfSNPs) and incorporated into six machine learning models. Predictive performance of a combination of different models and inputs were compared using the average cross validation area under ROC curve (AUC). The minimum combination of SNPs to achieve maximum sensitivity and specificity as determined by AUC, that predict atorvastatin-induced myalgia in most, if not all the six machine learning models was determined.

Results: Through whole genome association analyses using sex as a covariate, a larger proportion of pfSNPs compared to non-pf SNPs were found to be highly associated with myalgia. Although none of the individual SNPs achieved genome wide significance in univariate analyses, machine learning models identified a combination of 15 SNPs that predict myalgia with good predictive performance (AUC >0.9). SNPs within genes identified in this study significantly outperformed SNPs within genes previously reported to be associated with myalgia. pfSNPs were found to be more robust in predicting myalgia, outperforming non-pf SNPs in the majority of machine learning models tested.

Conclusion: Combinations of pfSNPs that were consistently identified by different machine learning models to have high predictive performance have good potential to be clinically useful for predicting atorvastatin-induced myalgia once validated against an independent cohort of patients.

Venous thromboembolism and prostate cancer: what about genetic markers?

Aim & methods: To review the existing literature concerning the relationship between venous thromboembolism (VTE) and prostate cancer (PC) and explore the putative biological and clinical implications of VTE genetic markers on PC patients by screening the PubMed database. Results: Considering the roles of VTE genome-wide association studies-identified genetic determinants in disease development in the general population, these variants might also underlie the susceptibility for PC-related VTE. Therefore, they could help to identify those with a positive benefit-to-harm ratio for thromboprophylaxis approaches during cancer therapy management, thereby improving patient's prognosis. Conclusion: Future studies are mandatory to explore the relationship between VTE and PC and dissect the predictive value of VTE genome-wide association studies-identified genetic determinants in PC patients, given their clinical implications.

Development and Validation of Epigenetic Signature Predict Survival for Patients with Laryngeal Squamous Cell Carcinoma

Establishing epigenetic signature to improve the accuracy of survival prediction and optimize therapeutic strategies for laryngeal squamous cell carcinoma (LSCC) by a genome-wide integrated analysis of methylation and the transcriptome. LSCC DNA methylation datasets and RNA sequencing datasets were acquired from the Cancer Genome Atlas (TCGA). MethylMix was applied to detect DNA methylation-driven genes (MDGs), which developed an epigenetic signature. The predictive accuracy and clinical value of the epigenetic signature were evaluated by receiver operating characteristic and decision curve analysis, and compared with tumor-node-metastasis (TNM) stage system. In addition, prognostic value of the epigenetic signature was validated by external Gene Expression Omnibus (GEO) database. According to five MDGs of epigenetic signature, the candidate small molecules for LSCC were screen out by the CMap database. A total of 88 DNA MDGs were identified, five of which (MAGEB2, SUSD1, ZNF382, ZNF418, and ZNF732) were chosen to construct an epigenetic signature. The epigenetic signature can effectively divide patients into high-risk and low-risk group, with the area under curve (AUC) of 0.8 (5-year overall survival [OS]) and AUC of 0.745 (3-year OS). Stratification analysis affirmed that the epigenetic signature was still a significant statistical prognostic model in subsets of patients with different clinical variables. Multivariate Cox regression analysis indicated that the efficacy of epigenetic signature appears independent of other clinicopathological characteristics. In terms of predictive capacity and clinical usefulness, the epigenetic signature was superior to traditional TNM stage. In addition, the epigenetic signature was confirmed in external LSCC cohorts from GEO. Finally, CMap matched the 10 most significant small molecules as promising therapeutic drugs to reverse the LSCC gene expression. An epigenetic signature, with five DNA MDGs, was identified and validated in LSCC patients by integrating multidimensional genomic data, which may offer novel research directions and prospects for individualized treatment of patients with LSCC.

Associations between blood type and COVID-19 infection, intubation, and death

The rapid global spread of the novel coronavirus SARS-CoV-2 has strained healthcare and testing resources, making the identification and prioritization of individuals most at-risk a critical challenge. Recent evidence suggests blood type may affect risk of severe COVID-19. Here, we use observational healthcare data on 14,112 individuals tested for SARS-CoV-2 with known blood type in the New York Presbyterian (NYP) hospital system to assess the association between ABO and Rh blood types and infection, intubation, and death. We find slightly increased infection prevalence among non-O types. Risk of intubation was decreased among A and increased among AB and B types, compared with type O, while risk of death was increased for type AB and decreased for types A and B. We estimate Rh-negative blood type to have a protective effect for all three outcomes. Our results add to the growing body of evidence suggesting blood type may play a role in COVID-19.

Recent evidence has suggested that blood type may be associated with severe COVID-19. Here, the authors use data from ~14,000 individuals tested for SARS-CoV-2 at a New York City hospital, and find that certain ABO and Rh blood types are associated with infection, intubation, and death.

Combined effects of five prothrombotic genotypes and cancer on the risk of a first venous thromboembolic event

BACKGROUND

The role of combined prothrombotic genotypes in cancer-related venous thromboembolism (VTE) is scarcely studied. We aimed to investigate the impact of a 5-single nucleotide polymorphism (SNP) score on the risk of VTE in patients with and without cancer using a population-based case-cohort.

METHODS

Cases with a first VTE (n = 1493) and a subcohort (n = 13 072) were derived from the Tromsø Study (1994-2012) and the Nord-Trøndelag Health Study (1995-2008). Five SNPs previously reported as a risk score were genotyped: ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865), and F11 (rs2036914). Hazard ratios (HRs) for VTE were estimated according to cancer status and the number of risk alleles in the 5-SNP score (0-1, 2-3, and ≥4 alleles).

RESULTS

During a median follow-up of 12.3 years, 1496 individuals were diagnosed with cancer, of whom 232 experienced VTE. The VTE risk increased with the number of risk alleles in the 5-SNP score among subjects without and with cancer. In cancer-free subjects, the HR was 2.17 (95% confidence interval [CI] 1.79-2.62) for ≥4 versus 0-1 risk alleles. In cancer patients, the corresponding HR was 1.93 (95% CI 1.28-2.91). The combination of cancer and ≥4 risk alleles yielded a 17-fold (HR 17.1, 95% CI 12.5-23.4) higher risk of VTE compared with cancer-free subjects with 0-1 risk alleles.

CONCLUSION

The risk of VTE increases with the number of prothrombotic risk alleles in subjects with and without cancer, and the combination of prothrombotic risk alleles and cancer leads to a highly elevated risk of VTE.

Genetic association study of fatal pulmonary embolism

Pulmonary embolism (PE) is a complex multi-factorial disease and represents one manifestation of venous thromboembolism (VTE). Most commonly PE constitutes a complication of VTE’s other clinical presentation deep vein thrombosis (DVT). The majority of studies concerning risk factors do not distinguish between PE and DVT. The risk factors are often estimated to be alike, but the prevalence and the risk associated with the major genetic factor Factor V Leiden differ between the two disease states. We have investigated the association of 22 SNPs with PE in 185 PE case and 375 healthy control subjects. At p = 0.05, eight SNPs presented with nominally significant evidence of association (EOA), although no significantly different genotype distributions remained between cases and controls after Bonferroni correction. Three of these variants (rs1800790, rs3813948, rs6025) showed EOA in the main analysis, and five variants (rs169713, rs1801131, rs4524, rs5985 and rs8176592) demonstrated EOAs in subgroups. Genomic variation modulating Factor V, Factor XIII, Beta fibrinogen (FGB), TFPI or HIVEP1 should be worth to be followed in subsequent studies. The findings of this study support the view that PE represents a complex disease with many factors contributing relatively small effects. Larger sample sizes will be required to reliably detect these small effects.

Testing the association between blood type and COVID-19 infection, intubation, and death

The rapid global spread of the novel coronavirus SARS-CoV-2 has strained healthcare and testing resources, making the identification and prioritization of individuals most at-risk a critical challenge. Recent evidence suggests blood type may affect risk of severe COVID-19. We used observational healthcare data on 14,112 individuals tested for SARS-CoV-2 with known blood type in the New York Presbyterian (NYP) hospital system to assess the association between ABO and Rh blood types and infection, intubation, and death. We found slightly increased infection prevalence among non-O types. Risk of intubation was decreased among A and increased among AB and B types, compared with type O, while risk of death was increased for type AB and decreased for types A and B. We estimated Rh-negative blood type to have a protective effect for all three outcomes. Our results add to the growing body of evidence suggesting blood type may play a role in COVID-19.

Whole-exome sequencing identifies rare variants in STAB2 associated with venous thromboembolic disease

Deep vein thrombosis and pulmonary embolism, collectively defined as venous thromboembolism (VTE), are the third leading cause of cardiovascular death in the United States. Common genetic variants conferring increased varying degrees of VTE risk have been identified by genome-wide association studies (GWAS). Rare mutations in the anticoagulant genes PROC, PROS1 and SERPINC1 result in perinatal lethal thrombosis in homozygotes and markedly increased VTE risk in heterozygotes. However, currently described VTE variants account for an insufficient portion of risk to be routinely used for clinical decision making. To identify new rare VTE risk variants, we performed a whole-exome study of 393 individuals with unprovoked VTE and 6114 controls. This study identified 4 genes harboring an excess number of rare damaging variants in patients with VTE: PROS1, STAB2, PROC, and SERPINC1. At STAB2, 7.8% of VTE cases and 2.4% of controls had a qualifying rare variant. In cell culture, VTE-associated variants of STAB2 had a reduced surface expression compared with reference STAB2. Common variants in STAB2 have been previously associated with plasma von Willebrand factor and coagulation factor VIII levels in GWAS, suggesting that haploinsufficiency of stabilin-2 may increase VTE risk through elevated levels of these procoagulants. In an independent cohort, we found higher von Willebrand factor levels and equivalent propeptide levels in individuals with rare STAB2 variants compared with controls. Taken together, this study demonstrates the utility of gene-based collapsing analyses to identify loci harboring an excess of rare variants with functional connections to a complex thrombotic disease.

Predicting the Risk of Recurrent Venous Thromboembolism: Current Challenges and Future Opportunities

Acute venous thromboembolism (VTE) is a commonly diagnosed condition and requires treatment with anticoagulation to reduce the risk of embolisation as well as recurrent venous thrombotic events. In many cases, cessation of anticoagulation is associated with an unacceptably high risk of recurrent VTE, precipitating the use of indefinite anticoagulation. In contrast, however, continuing anticoagulation is associated with increased major bleeding events. As a consequence, it is essential to accurately predict the subgroup of patients who have the highest probability of experiencing recurrent VTE, so that treatment can be appropriately tailored to each individual. To this end, the development of clinical prediction models has aided in calculating the risk of recurrent thrombotic events; however, there are several limitations with regards to routine use for all patients with acute VTE. More recently, focus has shifted towards the utility of novel biomarkers in the understanding of disease pathogenesis as well as their application in predicting recurrent VTE. Below, we review the current strategies used to predict the development of recurrent VTE, with emphasis on the application of several promising novel biomarkers in this field.

Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality. To advance our understanding of the biology contributing to VTE, we conducted a genome-wide association study (GWAS) of VTE and a transcriptome-wide association study (TWAS) based on imputed gene expression from whole blood and liver. We meta-analyzed GWAS data from 18 studies for 30 234 VTE cases and 172 122 controls and assessed the association between 12 923 718 genetic variants and VTE. We generated variant prediction scores of gene expression from whole blood and liver tissue and assessed them for association with VTE. Mendelian randomization analyses were conducted for traits genetically associated with novel VTE loci. We identified 34 independent genetic signals for VTE risk from GWAS meta-analysis, of which 14 are newly reported associations. This included 11 newly associated genetic loci (C1orf198, PLEK, OSMR-AS1, NUGGC/SCARA5, GRK5, MPHOSPH9, ARID4A, PLCG2, SMG6, EIF5A, and STX10) of which 6 replicated, and 3 new independent signals in 3 known genes. Further, TWAS identified 5 additional genetic loci with imputed gene expression levels differing between cases and controls in whole blood (SH2B3, SPSB1, RP11-747H7.3, RP4-737E23.2) and in liver (ERAP1). At some GWAS loci, we found suggestive evidence that the VTE association signal for novel and previously known regions colocalized with expression quantitative trait locus signals. Mendelian randomization analyses suggested that blood traits may contribute to the underlying risk of VTE. To conclude, we identified 16 novel susceptibility loci for VTE; for some loci, the association signals are likely mediated through gene expression of nearby genes.

Genome-wide association analysis of venous thromboembolism identifies new risk loci and genetic overlap with arterial vascular disease

Venous thromboembolism is a significant cause of mortality1, yet its genetic determinants are incompletely defined. We performed a discovery genome-wide association study in the Million Veteran Program and UK Biobank, with testing of approximately 13 million DNA sequence variants for association with venous thromboembolism (26,066 cases and 624,053 controls) and meta-analyzed both studies, followed by independent replication with up to 17,672 venous thromboembolism cases and 167,295 controls. We identified 22 previously unknown loci, bringing the total number of venous thromboembolism-associated loci to 33, and subsequently fine-mapped these associations. We developed a genome-wide polygenic risk score for venous thromboembolism that identifies 5% of the population at an equivalent incident venous thromboembolism risk to carriers of the established factor V Leiden p.R506Q and prothrombin G20210A mutations. Our data provide mechanistic insights into the genetic epidemiology of venous thromboembolism and suggest a greater overlap among venous and arterial cardiovascular disease than previously thought.

Identifying mutation-driven changes in gene functionality that lead to venous thromboembolism

Venous thromboembolism (VTE) is a common hematological disorder. VTE affects millions of people around the world each year and can be fatal. Earlier studies have revealed the possible VTE genetic risk factors in Europeans. The 2018 Critical Assessment of Genome Interpretation (CAGI) challenge had asked participants to distinguish between 66 VTE and 37 non-VTE African American (AA) individuals based on their exome sequencing data. We used variants from AA VTE association studies and VTE genes from DisGeNET database to evaluate VTE risk via four different approaches; two of these methods were most successful at the task. Our best performing method represented each exome as a vector of predicted functional effect scores of variants within the known genes. These exome vectors were then clustered with k-means. This approach achieved 70.8% precision and 69.7% recall in identifying VTE patients. Our second-best ranked method had collapsed the variant effect scores into gene-level function changes, using the same vector clustering approach for patient/control identification. These results show predictability of VTE risk in AA population and highlight the importance of variant-driven gene functional changes in judging disease status. Of course, more in-depth understanding of AA VTE pathogenicity is still needed for more precise predictions.

CAGI5: Objective performance assessments of predictions based on the Evolutionary Action equation

Many computational approaches estimate the effect of coding variants, but their predictions often disagree with each other. These contradictions confound users and raise questions regarding reliability. Performance assessments can indicate the expected accuracy for each method and highlight advantages and limitations. The Critical Assessment of Genome Interpretation (CAGI) community aims to organize objective and systematic assessments: They challenge predictors on unpublished experimental and clinical data and assign independent assessors to evaluate the submissions. We participated in CAGI experiments as predictors, using the Evolutionary Action (EA) method to estimate the fitness effect of coding mutations. EA is untrained, uses homology information, and relies on a formal equation: The fitness effect equals the functional sensitivity to residue changes multiplied by the magnitude of the substitution. In previous CAGI experiments (between 2011 and 2016), our submissions aimed to predict the protein activity of single mutants. In 2018 (CAGI5), we also submitted predictions regarding clinical associations, folding stability, and matching genomic data with phenotype. For all these diverse challenges, we used EA to predict the fitness effect of variants, adjusted to specifically address each question. Our submissions had consistently good performance, suggesting that EA predicts reliably the effects of genetic variants.

Influence of ABO Locus on PFA-100 Collagen-ADP Closure Time Is Not Totally Dependent on the Von Willebrand Factor. Results of a GWAS on GAIT-2 Project Phenotypes

(1) Background: In a previous study, we found that two phenotypes related to platelet reactivity, measured with the PFA-100 system, were highly heritable. The aim of the present study was to identify genetic determinants that influence the variability of these phenotypes: closure time of collagen-ADP (Col-ADP) and of collagen-epinephrine (Col-Epi). (2) Methods: As part of the GAIT-2 (Genetic Analysis of Idiopathic Thrombophilia (2) Project, 935 individuals from 35 large Spanish families were studied. A genome-wide association study (GWAS) with ≈ 10 M single nucleotide polymorphisms (SNPs) was carried out with Col-ADP and Col-Epi phenotypes. (3) Results: The study yielded significant genetic signals that mapped to the ABO locus. After adjusting both phenotypes for the ABO genotype, these signals disappeared. After adjusting for von Willebrand factor (VWF) or for coagulation factor VIII (FVIII), the significant signals disappeared totally for Col-Epi phenotype but only partially for Col-ADP phenotype. (4) Conclusion: Our results suggest that the ABO locus exerts the main genetic influence on PFA-100 phenotypes. However, while the effect of the ABO locus on Col-Epi phenotype is mediated through VWF and/or FVIII, the effect of the ABO locus on Col-ADP phenotype is partly produced through VWF and/or FVIII, and partly through other mechanisms.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Genetic determinants of activity and antigen levels of contact system factors

Essentials Genetic variation may provide valuable insight into the role of the contact system in thrombosis. Explored associations of genetic variants with activity, antigen, and disease in RATIO study. Two novel loci were identified: KLKB1 rs4253243 for prekallikrein; KNG1 rs5029980 for HMWK levels. Contact system variants and haplotypes were not associated with myocardial infarction or stroke. SUMMARY: Background The complex, interdependent contact activation system has been implicated in thrombotic disease, although few genetic determinants of levels of proteins from this system are known. Objectives Our primary aim was to study the influence of common F11, F12, KLKB1, and KNG1 variants on factor (F) XI activity and FXI, FXII, prekallikrein (PK) and high-molecular-weight kininogen (HMWK) antigen levels, as well as the risk of myocardial infarction and ischemic stroke. Patients/methods We analyzed samples from all 630 healthy participants, 182 ischemic stroke patients and 216 myocardial infarction patients in the RATIO case-control study of women aged < 50 years. Forty-three tagging single nucleotide variants (SNVs) were genotyped to represent common genetic variation in the contact system genes. Antigen and activity levels were measured with sandwich-ELISA-based and one-stage clotting assays. We performed single variant, age-adjusted, linear regression analyses per trait and disease phenotype, assuming additive inheritance and determined conditionally independent associations. Haplotypes based on the lead SNV and all conditionally independent SNVs were tested for association with traits and disease. Results We identified two novel associations of KLKB1 SNV rs4253243 with PK antigen (βconditional = -12.38; 95% CI, -20.07 to -4.69) and KNG1 SNV rs5029980 with HMWK antigen (βconditional = 5.86; 95% CI, 2.40-9.32) and replicated previously reported associations in a single study. Further analyses probed whether the observed associations were indicative of linkage, pleiotropic effects or mediation. No individual SNVs or haplotypes were associated with the disease outcomes. Conclusion This study adds to current knowledge of how genetic variation influences contact system protein levels and clarifies interdependencies.

Metabolites Associated With the Risk of Incident Venous Thromboembolism: A Metabolomic Analysis

Background Venous thromboembolism ( VTE ) is a complex thrombotic disorder that constitutes a major source of mortality and morbidity. To improve understanding of the cause of VTE , we conducted a metabolomic analysis in a case-control study including 240 incident VTE cases and 6963 controls nested within 3 large prospective population-based cohorts, the Nurses' Health Study, the Nurses' Health Study II , and the Health Professionals Follow-Up Study. Methods and Results For each individual, we measured 211 metabolites and collected detailed information on lifestyle factors. We performed logistic regression and enrichment analysis to identify metabolites and biological categories associated with incident VTE risk, accounting for key confounders, such as age, sex, smoking, alcohol consumption, body mass index, and comorbid diseases (eg, cancers). We performed analyses of all VTEs and separate analyses of pulmonary embolism. Using the basic model controlling for age, sex, and primary disease, we identified 60 nominally significant VTE - or pulmonary embolism-associated metabolites ( P<0.05). These metabolites were enriched for diacylglycerols ( P_permutation<0.05). However, after controlling for multiple testing, only 1 metabolite (C5 carnitine; odds ratio, 1.25; 95% confidence interval, 1.10-1.41; P_corrected=0.03) remained significantly associated with VTE . After further adjustment for body mass index, no metabolites were significantly associated with disease after accounting for multiple testing, and no metabolite classes were enriched for nominally significant associations. Conclusions Although our findings suggest that circulating metabolites may influence the risk of incident VTE , the associations we observed were confounded by body mass index. Larger studies involving additional individuals and with broader metabolomics coverage are needed to confirm our findings.

Phenome-wide association studies across large population cohorts support drug target validation

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

Testing the association between genetic variants and a range of phenotypes can assist drug development. Here, in a phenome-wide association study in up to 697,815 individuals, Diogo et al. identify genotype–phenotype associations predicting efficacy, alternative indications or adverse drug effects.

Association of the CYP4V2 polymorphism rs13146272 with venous thromboembolism in a Chinese population

Genome-wide association studies have identified the CYP4V2 polymorphism (rs13146272) as a risk factor associated with venous thromboembolism (VTE). However, due to the small sample size and variance in genetic analysis models, the relationship between VTE and rs13146272 remains unclear. Here, we performed a case-control study to analyse the associations between rs13146272 and VTE in a Chinese population and to compare the differences among various ethnicities. In this study, 226 VTE patients and 205 healthy controls were recruited, and the allele frequency of variant rs13146272 was analysed by a MassARRAY SNP genotyping assay. In addition, 9 case-control cohorts from 5 studies involving 6667 VTE-affected individuals and 8716 control subjects were included in this meta-analysis. Pooled ORs and 95% CIs were calculated to assess the association between rs13146272 and VTE by using different genetic models. Our case-control study results showed that there was no significant association between VTE and rs13146272 under the additive model (OR = 0.92, 95% CIs: 0.70-1.21, p = 0.55) in this Chinese population. However, the results of the meta-analysis performed by merging all cohorts showed that rs13146272 was significantly associated with VTE under the additive model, recessive model and dominant model. In the additive and recessive models, the association reached the threshold for genome-wide significance (p < 5.0e^-08). In conclusion, our pooled systematic study results indicated that individuals with the A allele had a higher risk of developing VTE than those with the C allele of the rs13146272 variant, but the risk was inconsistent among different ethnicities. Further validation of this association with larger sample sizes and multiple ethnicities is warranted.

Case-control studies of gene-environment interactions. When a case might not be the case

Case-control Genome-Wide Association Studies (GWAS) provide a rich resource for studying the genetic architecture of complex diseases. A key is to elucidate how the genetic effects vary by the environment, what is traditionally defined by Gene-Environment interactions (GxE). The overlooked complication is that multiple, distinct pathophysiologic mechanisms may lead to the same clinical diagnosis and often these mechanisms have distinct genetic bases. In this paper, we first show that using the clinically diagnosed status can lead to severely biased estimates of GxE interactions in situations when the frequency of the pathologic diagnosis of interest, as compared to other diagnoses, depends on the environment. We then propose a pseudo-likelihood solution to correct the bias. Finally, we demonstrate our method in extensive simulations and in a GWAS of Alzheimer's disease.

Deep learning sequence-based ab initio prediction of variant effects on expression and disease risk

Key challenges for human genetics, precision medicine and evolutionary biology include deciphering the regulatory code of gene expression and understanding the transcriptional effects of genome variation. However, this is extremely difficult because of the enormous scale of the noncoding mutation space. We developed a deep learning-based framework, ExPecto, that can accurately predict, ab initio from a DNA sequence, the tissue-specific transcriptional effects of mutations, including those that are rare or that have not been observed. We prioritized causal variants within disease- or trait-associated loci from all publicly available genome-wide association studies and experimentally validated predictions for four immune-related diseases. By exploiting the scalability of ExPecto, we characterized the regulatory mutation space for human RNA polymerase II-transcribed genes by in silico saturation mutagenesis and profiled > 140 million promoter-proximal mutations. This enables probing of evolutionary constraints on gene expression and ab initio prediction of mutation disease effects, making ExPecto an end-to-end computational framework for the in silico prediction of expression and disease risk.

The pharmacoepigenomics informatics pipeline defines a pathway of novel and known warfarin pharmacogenomics variants

AIM

'Pharmacoepigenomics' methods informed by omics datasets and pre-existing knowledge have yielded discoveries in neuropsychiatric pharmacogenomics. Now we evaluate the generality of these methods by discovering an extended warfarin pharmacogenomics pathway.

MATERIALS & METHODS

We developed the pharmacoepigenomics informatics pipeline, a scalable multi-omics variant screening pipeline for pharmacogenomics, and conducted an experiment in the genomics of warfarin.

RESULTS

We discovered known and novel pharmacogenomics variants and genes, both coding and regulatory, for warfarin response, including adverse events. Such genes and variants cluster in a warfarin response pathway consolidating known and novel warfarin response variants and genes.

CONCLUSION

These results can inform a new warfarin test. The pharmacoepigenomics informatics pipeline may be able to discover new pharmacogenomics markers in other drug-disease systems.

Genetic variants in mRNA untranslated regions

Genome Wide Association Studies (GWAS) have mapped thousands of genetic variants associated with complex disease risk and regulating quantitative traits, thus exploiting an unprecedented high-resolution genetic characterization of the human genome. A small fraction (3.7%) of the identified associations is located in untranslated regions (UTRs), and the molecular mechanism has been elucidated for few of them. Genetic variations at UTRs may modify regulatory elements affecting the interaction of the UTRs with proteins and microRNAs. The overall functional consequences include modulation of messenger RNA (mRNA) transcription, secondary structure, stability, localization, translation, and access to regulators like microRNAs (miRNAs) and RNA-binding proteins (RBPs). Alterations of these regulatory mechanisms are known to modify molecular pathways and cellular processes, potentially leading to disease processes. Here, we analyze some examples of genetic risk variants mapping in the UTR regulatory elements. We describe a recently identified genetic variant localized in the 3'UTR of the TNFSF13B gene, associated with autoimmunity risk and responsible of an increased stability and translation of TNFSF13B mRNA. We discuss how the correct use and interpretation of public GWAS repositories could lead to a better understanding of etiopathogenetic mechanisms and the generation of robust biological hypothesis as starting point for further functional studies. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Disease.