Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters.

Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases.

The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease.

Many common variants have been associated with hematological traits, but identification of causal genes and pathways has proven challenging. We performed a genome-wide association analysis in the UK Biobank and INTERVAL studies, testing 29.5 million genetic variants for association with 36 red cell, white cell, and platelet properties in 173,480 European-ancestry participants. This effort yielded hundreds of low frequency (<5%) and rare (<1%) variants with a strong impact on blood cell phenotypes. Our data highlight general properties of the allelic architecture of complex traits, including the proportion of the heritable component of each blood trait explained by the polygenic signal across different genome regulatory domains. Finally, through Mendelian randomization, we provide evidence of shared genetic pathways linking blood cell indices with complex pathologies, including autoimmune diseases, schizophrenia, and coronary heart disease and evidence suggesting previously reported population associations between blood cell indices and cardiovascular disease may be non-causal.

Whole-genome sequencing of patients with rare diseases in a national health system.

Most patients with rare diseases do not receive a molecular diagnosis and the aetiological variants and causative genes for more than half such disorders remain to be discovered1. Here we used whole-genome sequencing (WGS) in a national health system to streamline diagnosis and to discover unknown aetiological variants in the coding and non-coding regions of the genome. We generated WGS data for 13,037 participants, of whom 9,802 had a rare disease, and provided a genetic diagnosis to 1,138 of the 7,065 extensively phenotyped participants. We identified 95 Mendelian associations between genes and rare diseases, of which 11 have been discovered since 2015 and at least 79 are confirmed to be aetiological. By generating WGS data of UK Biobank participants2, we found that rare alleles can explain the presence of some individuals in the tails of a quantitative trait for red blood cells. Finally, we identified four novel non-coding variants that cause disease through the disruption of transcription of ARPC1B, GATA1, LRBA and MPL. Our study demonstrates a synergy by using WGS for diagnosis and aetiological discovery in routine healthcare.

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.

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.

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.

Unravelling the disease mechanism for TSPYL1 deficiency.

We describe a lethal combined nervous and reproductive systems disease in three affected siblings of a consanguineous family. The phenotype was characterized by visceroautonomic dysfunction (neonatal bradycardia/apnea, feeding problems, hyperactive startle reflex), severe postnatal progressive neurological abnormalities (including abnormal neonatal cry, hypotonia, epilepsy, polyneuropathy, cerebral gray matter atrophy), visual impairment, testicular dysgenesis in males and sudden death at infant age by brainstem-mediated cardiorespiratory arrest. Whole-exome sequencing revealed a novel homozygous frameshift variant p.Val242GlufsTer52 in the TSPY-like 1 gene (TSPYL1). The truncated TSPYL1 protein that lacks the nucleosome assembly protein domain was retained in the Golgi of fibroblasts from the three patients, whereas control fibroblasts express full-length TSPYL1 in the nucleus. Proteomic analysis of nuclear extracts from fibroblasts identified 24 upregulated and 20 downregulated proteins in the patients compared with 5 controls with 'regulation of cell cycle' as the highest scored biological pathway affected. TSPYL1-deficient cells had prolonged S and G2 phases with reduced cellular proliferation rates. Tspyl1 depletion in zebrafish mimicked the patients' phenotype with early lethality, defects in neurogenesis and cardiac dilation. In conclusion, this study reports the third pedigree with recessive TSPYL1 variants, confirming that TSPYL1 deficiency leads to a combined nervous and reproductive systems disease, and provides for the first time insights into the disease mechanism.

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.

Mutational and phenotypic characterization of hereditary hemorrhagic telangiectasia.

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant vascular dysplasia. Care delivery for HHT patients is impeded by the need for laborious, repeated phenotyping and gaps in knowledge regarding the relationships between causal DNA variants in ENG, ACVRL1, SMAD4 and GDF2, and clinical manifestations. To address this, we analyzed DNA samples from 183 previously uncharacterized, unrelated HHT and suspected HHT cases using the ThromboGenomics high-throughput sequencing platform. We identified 127 rare variants across 168 heterozygous genotypes. Applying modified American College of Medical Genetics and Genomics Guidelines, 106 variants were classified as pathogenic/likely pathogenic and 21 as nonpathogenic (variant of uncertain significance/benign). Unlike the protein products of ACVRL1 and SMAD4, the extracellular ENG amino acids are not strongly conserved. Our inferences of the functional consequences of causal variants in ENG were therefore informed by the crystal structure of endoglin. We then compared the accuracy of predictions of the causal gene blinded to the genetic data using 2 approaches: subjective clinical predictions and statistical predictions based on 8 Human Phenotype Ontology terms. Both approaches had some predictive power, but they were insufficiently accurate to be used clinically, without genetic testing. The distributions of red cell indices differed by causal gene but not sufficiently for clinical use in isolation from genetic data. We conclude that parallel sequencing of the 4 known HHT genes, multidisciplinary team review of variant calls in the context of detailed clinical information, and statistical and structural modeling improve the prognostication and treatment of HHT.

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.

Diagnostic high-throughput sequencing of 2396 patients with bleeding, thrombotic, and platelet disorders.

A targeted high-throughput sequencing (HTS) panel test for clinical diagnostics requires careful consideration of the inclusion of appropriate diagnostic-grade genes, the ability to detect multiple types of genomic variation with high levels of analytic sensitivity and reproducibility, and variant interpretation by a multidisciplinary team (MDT) in the context of the clinical phenotype. We have sequenced 2396 index patients using the ThromboGenomics HTS panel test of diagnostic-grade genes known to harbor variants associated with rare bleeding, thrombotic, or platelet disorders (BTPDs). The molecular diagnostic rate was determined by the clinical phenotype, with an overall rate of 49.2% for all thrombotic, coagulation, platelet count, and function disorder patients and a rate of 3.2% for patients with unexplained bleeding disorders characterized by normal hemostasis test results. The MDT classified 745 unique variants, including copy number variants (CNVs) and intronic variants, as pathogenic, likely pathogenic, or variants of uncertain significance. Half of these variants (50.9%) are novel and 41 unique variants were identified in 7 genes recently found to be implicated in BTPDs. Inspection of canonical hemostasis pathways identified 29 patients with evidence of oligogenic inheritance. A molecular diagnosis has been reported for 894 index patients providing evidence that introducing an HTS genetic test is a valuable addition to laboratory diagnostics in patients with a high likelihood of having an inherited BTPD.

Identification of a homozygous recessive variant in PTGS1 resulting in a congenital aspirin-like defect in platelet function.

We have identified a rare missense variant on chromosome 9, position 125145990 (GRCh37), in exon 8 in PTGS1 (the gene encoding cyclo-oxygenase 1, COX-1, the target of anti-thrombotic aspirin therapy). We report that in the homozygous state within a large consanguineous family this variant is associated with a bleeding phenotype and alterations in platelet reactivity and eicosanoid production. Western blotting and confocal imaging demonstrated that COX-1 was absent in the platelets of three family members homozygous for the PTGS1 variant but present in their leukocytes. Platelet reactivity, as assessed by aggregometry, lumi-aggregometry and flow cytometry, was impaired in homozygous family members, as were platelet adhesion and spreading. The productions of COX-derived eicosanoids by stimulated platelets were greatly reduced but there were no changes in the levels of urinary metabolites of COX-derived eicosanoids. The proband exhibited additional defects in platelet aggregation and spreading which may explain why her bleeding phenotype was slightly more severe than those of other homozygous affected relatives. This is the first demonstration in humans of the specific loss of platelet COX-1 activity and provides insight into its consequences for platelet function and eicosanoid metabolism. Notably despite the absence of thromboxane A2 (TXA2) formation by platelets, urinary TXA2 metabolites were in the normal range indicating these cannot be assumed as markers of in vivo platelet function. Results from this study are important benchmarks for the effects of aspirin upon platelet COX-1, platelet function and eicosanoid production as they define selective platelet COX-1 ablation within humans.

Cell type-specific novel long non-coding RNA and circular RNA in the BLUEPRINT hematopoietic transcriptomes atlas.

Transcriptional profiling of hematopoietic cell subpopulations has helped to characterize the developmental stages of the hematopoietic system and the molecular bases of malignant and non-malignant blood diseases. Previously, only the genes targeted by expression microarrays could be profiled genome-wide. High-throughput RNA sequencing, however, encompasses a broader repertoire of RNA molecules, without restriction to previously annotated genes. We analyzed the BLUEPRINT consortium RNA-sequencing data for mature hematopoietic cell types. The data comprised 90 total RNA-sequencing samples, each composed of one of 27 cell types, and 32 small RNA-sequencing samples, each composed of one of 11 cell types. We estimated gene and isoform expression levels for each cell type using existing annotations from Ensembl. We then used guided transcriptome assembly to discover unannotated transcripts. We identified hundreds of novel non-coding RNA genes and showed that the majority have cell type-dependent expression. We also characterized the expression of circular RNA and found that these are also cell type-specific. These analyses refine the active transcriptional landscape of mature hematopoietic cells, highlight abundant genes and transcriptional isoforms for each blood cell type, and provide a valuable resource for researchers of hematologic development and diseases. Finally, we made the data accessible via a web-based interface: https://blueprint.haem.cam.ac.uk/bloodatlas/.

Thrombomodulin in patients with mild to moderate bleeding tendency.

INTRODUCTION: A massive increase of soluble thrombomodulin (sTM) due to variants in the thrombomodulin gene (THBD) has recently been identified as a novel bleeding disorder. AIM: To investigate sTM levels and underlying genetic variants as a cause for haemostatic impairment and bleeding in a large number of patients with a mild to moderate bleeding disorder (MBD), including patients with bleeding of unknown cause (BUC). PATIENTS AND METHODS: In 507 MBD patients, sTM levels, thrombin generation and plasma clot formation were measured and compared to 90 age- and sex-matched healthy controls. In patients, genetic analysis of the THBD gene was performed. RESULTS: No difference in sTM levels between patients and controls was found overall (median ([IQR] 5.0 [3.8-6.3] vs. 5.1 [3.7-6.4] ng/ml, p = .762), and according to specific diagnoses of MBD or BUC, and high sTM levels (≥95th percentile of healthy controls) were not overrepresented in patients. Soluble TM levels had no impact on bleeding severity or global tests of haemostasis, including thrombin generation or plasma clot formation. In the THBD gene, no known pathogenic or novel disease-causing variants affecting sTM plasma levels were identified in our patient cohort. CONCLUSION: TM-associated coagulopathy appears to be rare, as it was not identified in our large cohort of patients with MBD. Soluble TM did not arise as a risk factor for bleeding or altered haemostasis in these patients.

Next-generation sequencing for the diagnosis of MYH9-RD: Predicting pathogenic variants.

The heterogeneous manifestations of MYH9-related disorder (MYH9-RD), characterized by macrothrombocytopenia, Döhle-like inclusion bodies in leukocytes, bleeding of variable severity with, in some cases, ear, eye, kidney, and liver involvement, make the diagnosis for these patients still challenging in clinical practice. We collected phenotypic data and analyzed the genetic variants in more than 3,000 patients with a bleeding or platelet disorder. Patients were enrolled in the BRIDGE-BPD and ThromboGenomics Projects and their samples processed by high throughput sequencing (HTS). We identified 50 patients with a rare variant in MYH9. All patients had macrothrombocytes and all except two had thrombocytopenia. Some degree of bleeding diathesis was reported in 41 of the 50 patients. Eleven patients presented hearing impairment, three renal failure and two elevated liver enzymes. Among the 28 rare variants identified in MYH9, 12 were novel. HTS was instrumental in diagnosing 23 patients (46%). Our results confirm the clinical heterogeneity of MYH9-RD and show that, in the presence of an unclassified platelet disorder with macrothrombocytes, MYH9-RD should always be considered. A HTS-based strategy is a reliable method to reach a conclusive diagnosis of MYH9-RD in clinical practice.

Bleeding of unknown cause and unclassified bleeding disorders; diagnosis, pathophysiology and management.

Bleeding of unknown cause (BUC), also known as unclassified bleeding disorders (UBD), has been defined as a clear bleeding tendency in the presence of normal haemostatic tests. There are challenges in the diagnosis and management of these patients. BUC/UBD encompasses a heterogenous group of disorders which may include undiagnosed rare monogenic diseases, polygenic reasons for bleeding; and patients without a clear bleeding disorder but with a previous bleeding event. Nevertheless, these patients may have heavy menstrual bleeding or be at risk of bleeding when undergoing surgical procedures, or childbirth; optimizing haemostasis and establishing a mode of inheritance is important to minimize morbidity. The bleeding score has been used to clinically assess and describe these patients, but its value remains uncertain. In addition, accurate distinction between normal and pathological bleeding remains difficult. Several studies have investigated cohorts of these patients using research haemostasis tests, including thrombin generation and fibrinolytic assays, yet no clear characteristics have consistently emerged. Thus far, detailed genetic analysis of these patients has not been fruitful in unravelling the cause of bleeding. There is a need for standardization of diagnosis and management guidelines for these patients. This review gives an overview of this field with some suggestions for future research.

Expanded repertoire of RASGRP2 variants responsible for platelet dysfunction and severe bleeding.

Heritable platelet function disorders (PFDs) are genetically heterogeneous and poorly characterized. Pathogenic variants in RASGRP2, which encodes calcium and diacylglycerol-regulated guanine exchange factor I (CalDAG-GEFI), have been reported previously in 3 pedigrees with bleeding and reduced platelet aggregation responses. To better define the phenotype associated with pathogenic RASGRP2 variants, we compared high-throughput sequencing and phenotype data from 2042 cases in pedigrees with unexplained bleeding or platelet disorders to data from 5422 controls. Eleven cases harbored 11 different, previously unreported RASGRP2 variants that were biallelic and likely pathogenic. The variants included 5 high-impact variants predicted to prevent CalDAG-GEFI expression and 6 missense variants affecting the CalDAG-GEFI CDC25 domain, which mediates Rap1 activation during platelet inside-out αIIbß3 signaling. Cases with biallelic RASGRP2 variants had abnormal mucocutaneous, surgical, and dental bleeding from childhood, requiring ≥1 blood or platelet transfusion in 78% of cases. Platelets displayed reduced aggregation in response to adenosine 5'-diphosphate and epinephrine, but variable aggregation defects with other agonists. There were no other consistent clinical or laboratory features. These data enable definition of human CalDAG-GEFI deficiency as a nonsyndromic, recessive PFD associated with a moderate or severe bleeding phenotype and complex defects in platelet aggregation.

Rare variants in GP1BB are responsible for autosomal dominant macrothrombocytopenia.

The von Willebrand receptor complex, which is composed of the glycoproteins Ibα, Ibß, GPV, and GPIX, plays an essential role in the earliest steps in hemostasis. During the last 4 decades, it has become apparent that loss of function of any 1 of 3 of the genes encoding these glycoproteins (namely, GP1BA, GP1BB, and GP9) leads to autosomal recessive macrothrombocytopenia complicated by bleeding. A small number of variants in GP1BA have been reported to cause a milder and dominant form of macrothrombocytopenia, but only 2 tentative reports exist of such a variant in GP1BB By analyzing data from a collection of more than 1000 genome-sequenced patients with a rare bleeding and/or platelet disorder, we have identified a significant association between rare monoallelic variants in GP1BB and macrothrombocytopenia. To strengthen our findings, we sought further cases in 2 additional collections in the United Kingdom and Japan. Across 18 families exhibiting phenotypes consistent with autosomal dominant inheritance of macrothrombocytopenia, we report on 27 affected cases carrying 1 of 9 rare variants in GP1BB.

High-throughput elucidation of thrombus formation reveals sources of platelet function variability.

In combination with microspotting, whole-blood microfluidics can provide high-throughput information on multiple platelet functions in thrombus formation. Based on assessment of the inter- and intra-subject variability in parameters of microspot-based thrombus formation, we aimed to determine the platelet factors contributing to this variation. Blood samples from 94 genotyped healthy subjects were analyzed for conventional platelet phenotyping: i.e. hematologic parameters, platelet glycoprotein (GP) expression levels and activation markers (24 parameters). Furthermore, platelets were activated by ADP, CRP-XL or TRAP. Parallel samples were investigated for whole-blood thrombus formation (6 microspots, providing 48 parameters of adhesion, aggregation and activation). Microspots triggered platelet activation through GP Ib-V-IX, GPVI, CLEC-2 and integrins. For most thrombus parameters, inter-subject variation was 2-4 times higher than the intra-subject variation. Principal component analyses indicated coherence between the majority of parameters for the GPVI-dependent microspots, partly linked to hematologic parameters, and glycoprotein expression levels. Prediction models identified parameters per microspot that were linked to variation in agonist-induced αIIbß3 activation and secretion. Common sequence variation of GP6 and FCER1G, associated with GPVI-induced αIIbß3 activation and secretion, affected parameters of GPVI-and CLEC-2-dependent thrombus formation. Subsequent analysis of blood samples from patients with Glanzmann thrombasthenia or storage pool disease revealed thrombus signatures of aggregation-dependent parameters that were subject-dependent, but not linked to GPVI activity. Taken together, this high-throughput elucidation of thrombus formation revealed patterns of inter-subject differences in platelet function, which were partly related to GPVI-induced activation and common genetic variance linked to GPVI, but also included a distinct platelet aggregation component.