The common, VTE-protective, G haplotype of F5 increases factor V-short, TFPI function, and risk of bleeding.

The G haplotype is a group of co-inherited single nucleotide variants in the F5 gene that reduce venous thromboembolism (VTE) risk. Even though seven percent of the population is homozygous for the G haplotype (F5-G/G), the underlying mechanism of VTE protection is poorly understood. Using RNA-seq data from 4,651 blood donors in the INTERVAL study we detected a rare excision event at the FV-short splice sites in 5% of F5-G/Gs as compared with 2.16% of homozygotes for the F5 reference sequence (F5-ref) (p=0.003). Highly elevated (~10-fold) FV-short, an FV isoform lacking most of the B-domain, has been linked with increased tissue factor inhibitor alpha (TFPIα) levels in rare hemorrhagic diathesis including East Texas Bleeding Disorder. To ascertain whether the enhanced FV-short splicing seen in F5-G/G INTERVAL participants translated to increased plasma FV-short levels we analyzed plasma samples from 7 F5-G/G and 13 F5-ref individuals in a recall-by-genotype study. A ~2.2-fold higher amount of FV-short was found in a plasma pool from F5-G/G participants as compared with F5-refs (p=0.029), but no difference in total FV levels. Whilst no significant difference in TFPI levels were found, F5-G/Gs showed a ~1.4-fold TFPI-dependent increase in lag time to thrombin generation compared to F5-refs (p=0.0085). Finally, in an analysis of 117,699 UK Biobank participants we discovered that, while being protective against VTE, the G haplotype also confers an increase in bleeding episodes (p=0.011). Our study provides evidence that the effect of the common G haplotype is mediated by the FV-short/TFPI pathway.

Evaluating the clinical validity of genes related to hemostasis and thrombosis using the Clinical Genome Resource gene curation framework.

BACKGROUND: Inherited bleeding, thrombotic, and platelet disorders (BTPDs) are a heterogeneous set of diseases, many of which are very rare globally. Over the past 5 decades, the genetic basis of some of these disorders has been identified, and recently, high-throughput sequencing has become the primary means of identifying disease-causing genetic variants. OBJECTIVES: Knowledge of the clinical validity of a gene-disease relationship is essential to provide an accurate diagnosis based on results of diagnostic gene panel tests and inform the construction of such panels. The Scientific and Standardization Committee for Genetics in Thrombosis and Hemostasis undertook a curation process for selecting 96 TIER1 genes for BTPDs. The purpose of the process was to evaluate the evidence supporting each gene-disease relationship and provide an expert-reviewed classification for the clinical validity of genes associated with BTPDs. METHODS: The Clinical Genome Resource (ClinGen) Hemostasis/Thrombosis Gene Curation Expert Panel assessed the strength of evidence for TIER1 genes using the semiquantitative ClinGen gene-disease clinical validity framework. ClinGen Lumping and Splitting guidelines were used to determine the appropriate disease entity or entities for each gene, and 101 gene-disease relationships were identified for curation. RESULTS: The final outcome included 68 Definitive (67%), 26 Moderate (26%), and 7 Limited (7%) classifications. The summary of each curation is available on the ClinGen website. CONCLUSION: Expert-reviewed assignment of gene-disease relationships by the ClinGen Hemostasis/Thrombosis Gene Curation Expert Panel facilitates accurate molecular diagnoses of BTPDs by clinicians and diagnostic laboratories. These curation efforts can allow genetic testing to focus on genes with a validated role in disease.

Platelet Receptor Glycoprotein VI-Dimer Is Overexpressed in Patients with Atrial Fibrillation at High Risk of Ischemic Stroke.

Introduction Atrial fibrillation (AF) increases the risk of ischemic stroke (IS). We hypothesized that the functional form of platelet receptor glycoprotein (GP) VI, GPVI-dimer, which binds to collagen and fibrin causing platelet activation, is overexpressed in patients with AF who have not had a stroke. Methods A total of 75 inpatients with AF were recruited. None were admitted with or had previously had thrombotic events, including IS or myocardial infarction. Platelet surface expression of total GPVI, GPVI-dimer, and the platelet activation marker P-selectin were quantitated by whole blood flow cytometry. Serum biomarkers were collected in AF patients. Results were compared against patients contemporaneously admitted to hospital with similar age and vascular risk-factor profiles without AF (noAF, n = 30). Results Patients with AF have similar total GPVI surface expression ( p = 0.58) and P-selectin exposure ( p = 0.73) on their platelets compared with noAF patients but demonstrate significantly higher GPVI-dimer expression ( p = 0.02 ). Patients with paroxysmal AF express similar GPVI-dimer levels compared with permanent AF and GPVI-dimer levels were not different between anticoagulated groups. Serum N-terminal pro b-type natriuretic peptide ( p < 0.0001 ) and high sensitivity C-reactive protein ( p < 0.0001 ) were significantly correlated with GPVI-dimer expression in AF platelets. AF was the only vascular risk factor that was independently associated with higher GPVI-dimer expression in the whole population ( p = 0.02 ) . Conclusion GPVI inhibition is being explored in clinical trials as a novel target for IS treatment. As GPVI-dimer is elevated in AF patients' platelets, the exploration of targeted GPVI-dimer inhibition for stroke prevention in patients at high risk of IS due to AF is supported.

A genome-wide association study of blood cell morphology identifies cellular proteins implicated in disease aetiology.

Blood cells contain functionally important intracellular structures, such as granules, critical to immunity and thrombosis. Quantitative variation in these structures has not been subjected previously to large-scale genetic analysis. We perform genome-wide association studies of 63 flow-cytometry derived cellular phenotypes-including cell-type specific measures of granularity, nucleic acid content and reactivity-in 41,515 participants in the INTERVAL study. We identify 2172 distinct variant-trait associations, including associations near genes coding for proteins in organelles implicated in inflammatory and thrombotic diseases. By integrating with epigenetic data we show that many intracellular structures are likely to be determined in immature precursor cells. By integrating with proteomic data we identify the transcription factor FOG2 as an early regulator of platelet formation and α-granularity. Finally, we show that colocalisation of our associations with disease risk signals can suggest aetiological cell-types-variants in IL2RA and ITGA4 respectively mirror the known effects of daclizumab in multiple sclerosis and vedolizumab in inflammatory bowel disease.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants.

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

A signature of platelet reactivity in CBC scattergrams reveals genetic predictors of thrombotic disease risk.

Genetic studies of platelet reactivity (PR) phenotypes may identify novel antiplatelet drug targets. However, such studies have been limited by small sample sizes (n < 5000) because of the complexity of measuring PR. We trained a model to predict PR from complete blood count (CBC) scattergrams. A genome-wide association study of this phenotype in 29 806 blood donors identified 21 distinct associations implicating 20 genes, of which 6 have been identified previously. The effect size estimates were significantly correlated with estimates from a study of flow cytometry-measured PR and a study of a phenotype of in vitro thrombus formation. A genetic score of PR built from the 21 variants was associated with the incidence rates of myocardial infarction and pulmonary embolism. Mendelian randomization analyses showed that PR was causally associated with the risks of coronary artery disease, stroke, and venous thromboembolism. Our approach provides a blueprint for using phenotype imputation to study the determinants of hard-to-measure but biologically important hematological traits.

Incorporating Alternative Polygenic Risk Scores into the BOADICEA Breast Cancer Risk Prediction Model.

BACKGROUND: The multifactorial risk prediction model BOADICEA enables identification of women at higher or lower risk of developing breast cancer. BOADICEA models genetic susceptibility in terms of the effects of rare variants in breast cancer susceptibility genes and a polygenic component, decomposed into an unmeasured and a measured component - the polygenic risk score (PRS). The current version was developed using a 313 SNP PRS. Here, we evaluated approaches to incorporating this PRS and alternative PRS in BOADICEA. METHODS: The mean, SD, and proportion of the overall polygenic component explained by the PRS (α2) need to be estimated. $\alpha $ was estimated using logistic regression, where the age-specific log-OR is constrained to be a function of the age-dependent polygenic relative risk in BOADICEA; and using a retrospective likelihood (RL) approach that models, in addition, the unmeasured polygenic component. RESULTS: Parameters were computed for 11 PRS, including 6 variations of the 313 SNP PRS used in clinical trials and implementation studies. The logistic regression approach underestimates $\alpha $, as compared with the RL estimates. The RL $\alpha $ estimates were very close to those obtained by assuming proportionality to the OR per 1 SD, with the constant of proportionality estimated using the 313 SNP PRS. Small variations in the SNPs included in the PRS can lead to large differences in the mean. CONCLUSIONS: BOADICEA can be readily adapted to different PRS in a manner that maintains consistency of the model. IMPACT: : The methods described facilitate comprehensive breast cancer risk assessment.

Effectiveness and costs of a stepwise versus an all-in-one approach to diagnose mild bleeding disorders.

The diagnostic work-up of patients referred to the haematologist for bleeding evaluation is performed in a stepwise way: bleeding history and results of screening laboratory tests guide further diagnostic evaluation. This can be ineffective, time-consuming and burdensome for patients. To improve this strategy, the initial laboratory investigation can be extended. In a model-based approach, effectiveness and costs of a conventional stepwise versus a newly proposed all-in-one diagnostic approach for bleeding evaluation were evaluated and compared, using data from an observational patient cohort study, including adult patients referred for bleeding evaluation. In the all-in-one approach, specialized platelet function tests, coagulation factors, and fibrinolysis tests were included in the initial investigation. Final diagnosis, hospital resource use and costs and patient burden were compared. A total of 150 patients were included. Compared to the stepwise approach, in the all-in-one approach, 19 additional patients reached a diagnosis and patient burden was lower, but total costs per patient were higher [€359, 95% bootstrapped confidence interval (BCI) 283-518, p = 0.001]. For bleeding evaluation of patients referred to the haematologist, an all-in-one diagnostic approach has a higher diagnostic yield and reduces patient burden, at a higher cost. This raises the question what costs justify the diagnosis of a bleeding disorder and a less burdensome diagnostic strategy.

Detection of enterovirus RNA in peripheral blood mononuclear cells correlates with the presence of the predisposing allele of the type 1 diabetes risk gene IFIH1 and with disease stage.

AIMS/HYPOTHESIS: Enteroviral infection has been implicated consistently as a key environmental factor correlating with the appearance of autoimmunity and/or the presence of overt type 1 diabetes, in which pancreatic insulin-producing beta cells are destroyed by an autoimmune response. Genetic predisposition through variation in the type 1 diabetes risk gene IFIH1 (interferon induced with helicase C domain 1), which encodes the viral pattern-recognition receptor melanoma differentiation-associated protein 5 (MDA5), supports a potential link between enterovirus infection and type 1 diabetes. METHODS: We used molecular techniques to detect enterovirus RNA in peripheral blood samples (in separated cellular compartments or plasma) from two cohorts comprising 79 children or 72 adults that include individuals with and without type 1 diabetes who had multiple autoantibodies. We also used immunohistochemistry to detect the enteroviral protein VP1 in the pancreatic islets of post-mortem donors (n=43) with type 1 diabetes. RESULTS: We observed enhanced detection sensitivity when sampling the cellular compartment compared with the non-cellular compartment of peripheral blood (OR 21.69; 95% CI 3.64, 229.20; p<0.0001). In addition, we show that children with autoimmunity are more likely to test positive for enterovirus RNA than those without autoimmunity (OR 11.60; 95% CI 1.89, 126.90; p=0.0065). Furthermore, we found that individuals carrying the predisposing allele (946Thr) of the common variant in IFIH1 (rs1990760, Thr946Ala) are more likely to test positive for enterovirus in peripheral blood (OR 3.07; 95% CI 1.02, 8.58; p=0.045). In contrast, using immunohistochemistry, there was no correlation between the common variant in IFIH1 and detection of enteroviral VP1 protein in the pancreatic islets of donors with type 1 diabetes. CONCLUSIONS/INTERPRETATION: Our data indicate that, in peripheral blood, antigen-presenting cells are the predominant source of enterovirus infection, and that infection is correlated with disease stage and genetic predisposition, thereby supporting a role for enterovirus infection prior to disease onset.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome.

BACKGROUND: The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts. METHODS: We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups. RESULTS: We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants. CONCLUSIONS: These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

Personalised Risk Prediction in Hereditary Breast and Ovarian Cancer: A Protocol for a Multi-Centre Randomised Controlled Trial.

Women who test positive for an inherited pathogenic/likely pathogenic gene variant in BRCA1, BRCA2, PALB2, CHEK2 and ATM are at an increased risk of developing certain types of cancer-specifically breast (all) and epithelial ovarian cancer (only BRCA1, BRCA2, PALB2). Women receive broad cancer risk figures that are not personalised (e.g., 44-63% lifetime risk of breast cancer for those with PALB2). Broad, non-personalised risk estimates may be problematic for women when they are considering how to manage their risk. Multifactorial-risk-prediction tools have the potential to deliver personalised risk estimates. These may be useful in the patient's decision-making process and impact uptake of risk-management options. This randomised control trial (registration number to follow), based in genetic centres in the UK and US, will randomise participants on a 1:1 basis to either receive conventional cancer risk estimates, as per routine clinical practice, or to receive a personalised risk estimate. This personalised risk estimate will be calculated using the CanRisk risk prediction tool, which combines the patient's genetic result, family history and polygenic risk score (PRS), along with hormonal and lifestyle factors. Women's decision-making around risk management will be monitored using questionnaires, completed at baseline (pre-appointment) and follow-up (one, three and twelve months after receiving their risk assessment). The primary outcome for this study is the type and timing of risk management options (surveillance, chemoprevention, surgery) taken up over the course of the study (i.e., 12 months). The type of risk-management options planned to be taken up in the future (i.e., beyond the end of the study) and the potential impact of personalised risk estimates on women's psychosocial health will be collected as secondary-outcome measures. This study will also assess the acceptability, feasibility and cost-effectiveness of using personalised risk estimates in clinical care.

Transcriptional, epigenetic and metabolic signatures in cardiometabolic syndrome defined by extreme phenotypes.

BACKGROUND: This work is aimed at improving the understanding of cardiometabolic syndrome pathophysiology and its relationship with thrombosis by generating a multi-omic disease signature. METHODS/RESULTS: We combined classic plasma biochemistry and plasma biomarkers with the transcriptional and epigenetic characterisation of cell types involved in thrombosis, obtained from two extreme phenotype groups (morbidly obese and lipodystrophy) and lean individuals to identify the molecular mechanisms at play, highlighting patterns of abnormal activation in innate immune phagocytic cells. Our analyses showed that extreme phenotype groups could be distinguished from lean individuals, and from each other, across all data layers. The characterisation of the same obese group, 6 months after bariatric surgery, revealed the loss of the abnormal activation of innate immune cells previously observed. However, rather than reverting to the gene expression landscape of lean individuals, this occurred via the establishment of novel gene expression landscapes. NETosis and its control mechanisms emerge amongst the pathways that show an improvement after surgical intervention. CONCLUSIONS: We showed that the morbidly obese and lipodystrophy groups, despite some differences, shared a common cardiometabolic syndrome signature. We also showed that this could be used to discriminate, amongst the normal population, those individuals with a higher likelihood of presenting with the disease, even when not displaying the classic features.

Elucidating mechanisms of genetic cross-disease associations at the PROCR vascular disease locus.

Many individual genetic risk loci have been associated with multiple common human diseases. However, the molecular basis of this pleiotropy often remains unclear. We present an integrative approach to reveal the molecular mechanism underlying the PROCR locus, associated with lower coronary artery disease (CAD) risk but higher venous thromboembolism (VTE) risk. We identify PROCR-p.Ser219Gly as the likely causal variant at the locus and protein C as a causal factor. Using genetic analyses, human recall-by-genotype and in vitro experimentation, we demonstrate that PROCR-219Gly increases plasma levels of (activated) protein C through endothelial protein C receptor (EPCR) ectodomain shedding in endothelial cells, attenuating leukocyte-endothelial cell adhesion and vascular inflammation. We also associate PROCR-219Gly with an increased pro-thrombotic state via coagulation factor VII, a ligand of EPCR. Our study, which links PROCR-219Gly to CAD through anti-inflammatory mechanisms and to VTE through pro-thrombotic mechanisms, provides a framework to reveal the mechanisms underlying similar cross-phenotype associations.

Platelet surface receptor glycoprotein VI-dimer is overexpressed in stroke: The Glycoprotein VI in Stroke (GYPSIE) study results.

OBJECTIVES: Platelet activation underpins thrombus formation in ischemic stroke. The active, dimeric form of platelet receptor glycoprotein (GP) VI plays key roles by binding platelet ligands collagen and fibrin, leading to platelet activation. We investigated whether patients presenting with stroke expressed more GPVI on their platelet surface and had more active circulating platelets as measured by platelet P-selectin exposure. METHODS: 129 ischemic or hemorrhagic stroke patients were recruited within 8h of symptom onset. Whole blood was analyzed for platelet-surface expression of total GPVI, GPVI-dimer, and P-selectin by flow cytometry at admission and day-90 post-stroke. Results were compared against a healthy control population (n = 301). RESULTS: The platelets of stroke patients expressed significantly higher total GPVI and GPVI-dimer (P<0.0001) as well as demonstrating higher resting P-selectin exposure (P<0.0001), a measure of platelet activity, compared to the control group, suggesting increased circulating platelet activation. GPVI-dimer expression was strongly correlated circulating platelet activation [r2 = 0.88, P<0.0001] in stroke patients. Furthermore, higher platelet surface GPVI expression was associated with increased stroke severity at admission. At day-90 post-stroke, GPVI-dimer expression and was further raised compared to the level at admission (P<0.0001) despite anti-thrombotic therapy. All ischemic stroke subtypes and hemorrhagic strokes expressed significantly higher GPVI-dimer compared to controls (P<0.0001). CONCLUSIONS: Stroke patients express more GPVI-dimer on their platelet surface at presentation, lasting at least until day-90 post-stroke. Small molecule GPVI-dimer inhibitors are currently in development and the results of this study validate that GPVI-dimer as an anti-thrombotic target in ischemic stroke.

Multiparameter platelet function analysis of bleeding patients with a prolonged platelet function analyser closure time.

Patients referred for evaluation of bleeding symptoms occasionally have a prolonged platelet function analyser (PFA) closure time, without evidence for von Willebrand disease or impaired platelet aggregation. The aim of this study was to establish a shear-dependent platelet function defect in these patients. Patients were included based on high bleeding score and prior PFA prolongation. Common tests of von Willebrand factor (VWF) and platelet function and exome sequencing were performed. Microfluidic analysis of shear-dependent collagen-induced whole-blood thrombus formation was performed. In 14 PFA-only patients, compared to healthy volunteers, microfluidic tests showed significantly lower platelet adhesion and thrombus formation parameters. This was accompanied by lower integrin activation, phosphatidylserine exposure and P-selectin expression. Principal components analysis indicated VWF as primary explaining variable of PFA prolongation, whereas conventional platelet aggregation primarily explained the reduced thrombus parameters under shear. In five patients with severe microfluidic abnormalities, conventional platelet aggregation was in the lowest range of normal. No causal variants in Mendelian genes known to cause bleeding or platelet disorders were identified. Multiparameter assessment of whole-blood thrombus formation under shear indicates single or combined effects of low-normal VWF and low-normal platelet aggregation in these patients, suggesting a shear-dependent platelet function defect, not detected by static conventional haemostatic tests.

A novel RUNX1 exon 3 - 7 deletion causing a familial platelet disorder.

Familial Platelet Disorder with associated Myeloid Malignancy (FPDMM) is a rare inherited disorder confirmed with the presence of a pathogenic germline RUNX1 variant and is thought to be heavily underdiagnosed. RUNX1 has also been found to be mutated in up to 10% of adult AML cases and other cell malignancies. We performed targeted next-generation sequencing and subsequent MLPA analysis in a kindred with multiple affected individuals with low platelet counts and a bleeding history. We detected a novel heterozygous exon 3-7 large deletion in the RUNX1 gene in all affected family members which is predicted to remove all of the Runt-homology DNA-binding domain and a portion of the Activation domain. Our results show that the combination of targeted NGS and MLPA analysis is an effective way to detect copy number variants (CNVs) which would be missed by conventional sequencing methods. This precise diagnosis offers the possibility of accurate counseling and clinical management in such patients who could go onto develop other cell malignancies.

G protein-coupled receptor kinase 5 regulates thrombin signaling in platelets via PAR-1.

The interindividual variation in the functional response of platelets to activation by agonists is heritable. Genome-wide association studies (GWASs) of quantitative measures of platelet function have identified fewer than 20 distinctly associated variants, some with unknown mechanisms. Here, we report GWASs of pathway-specific functional responses to agonism by adenosine 5'-diphosphate, a glycoprotein VI-specific collagen mimetic, and thrombin receptor-agonist peptides, each specific to 1 of the G protein-coupled receptors PAR-1 and PAR-4, in subsets of 1562 individuals. We identified an association (P = 2.75 × 10-40) between a common intronic variant, rs10886430, in the G protein-coupled receptor kinase 5 gene (GRK5) and the sensitivity of platelets to activate through PAR-1. The variant resides in a megakaryocyte-specific enhancer that is bound by the transcription factors GATA1 and MEIS1. The minor allele (G) is associated with fewer GRK5 transcripts in platelets and the greater sensitivity of platelets to activate through PAR-1. We show that thrombin-mediated activation of human platelets causes binding of GRK5 to PAR-1 and that deletion of the mouse homolog Grk5 enhances thrombin-induced platelet activation sensitivity and increases platelet accumulation at the site of vascular injury. This corroborates evidence that the human G allele of rs10886430 is associated with a greater risk for cardiovascular disease. In summary, by combining the results of pathway-specific GWASs and expression quantitative trait locus studies in humans with the results from platelet function studies in Grk5-/- mice, we obtain evidence that GRK5 regulates the human platelet response to thrombin via the PAR-1 pathway.