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.

Long-Read Sequencing Identifies the First Retrotransposon Insertion and Resolves Structural Variants Causing Antithrombin Deficiency.

The identification of inherited antithrombin deficiency (ATD) is critical to prevent potentially life-threatening thrombotic events. Causal variants in SERPINC1 are identified for up to 70% of cases, the majority being single-nucleotide variants and indels. The detection and characterization of structural variants (SVs) in ATD remain challenging due to the high number of repetitive elements in SERPINC1. Here, we performed long-read whole-genome sequencing on 10 familial and 9 singleton cases with type I ATD proven by functional and antigen assays, who were selected from a cohort of 340 patients with this rare disorder because genetic analyses were either negative, ambiguous, or not fully characterized. We developed an analysis workflow to identify disease-associated SVs. This approach resolved, independently of its size or type, all eight SVs detected by multiple ligation-dependent probe amplification, and identified for the first time a complex rearrangement previously misclassified as a deletion. Remarkably, we identified the mechanism explaining ATD in 2 out of 11 cases with previous unknown defect: the insertion of a novel 2.4 kb SINE-VNTR-Alu retroelement, which was characterized by de novo assembly and verified by specific polymerase chain reaction amplification and sequencing in the probands and affected relatives. The nucleotide-level resolution achieved for all SVs allowed breakpoint analysis, which revealed repetitive elements and microhomologies supporting a common replication-based mechanism for all the SVs. Our study underscores the utility of long-read sequencing technology as a complementary method to identify, characterize, and unveil the molecular mechanism of disease-causing SVs involved in ATD, and enlarges the catalogue of genetic disorders caused by retrotransposon insertions.

Non-coding genetic variation in regulatory elements determines thrombosis and hemostasis phenotypes.

Since the early inception of genome-wide association studies (GWAS), it became clear that, in all diseases or traits studied, most genetic variants are likely to exert their effect on gene expression mainly by altering the function of regulatory elements. At the same time, the regulation of the gene expression field broadened its boundaries, from the univocal relationship between regulatory elements and genes to include genome organization, long-range DNA interactions, and epigenetics. Next-generation sequencing has introduced genome-wide approaches that have greatly improved our understanding of the general principles of gene expression. However, elucidating how these apply in every single genomic locus still requires painstaking experimental work, in which several independent lines of evidence are required, and often this is helped by rare genetic variants in individuals with rare diseases. This review will focus on the non-coding features of the genome involved in transcriptional regulation, that when altered, leads to known cases of inherited (familial) thrombotic and hemostatic phenotypes, emphasizing the role of enhancers and super-enhancers.

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.

Familial pseudohyperkalemia induces significantly higher levels of extracellular potassium in early storage of red cell concentrates without affecting other standard measures of quality: A case control and allele frequency study.

BACKGROUND: Familial pseudohyperkalemia (FP) is characterized by an increased rate of potassium leakage in refrigerated red cells and is associated with the minor allele of the single nucleotide polymorphism rs148211042 (R723Q) in the ABCB6 gene. The study aims were to obtain the minor allele frequencies of ABCB6 variants and to measure supernatant potassium accumulation, and other red cell storage parameters, in red cell concentrates (RCC) from carriers of variant rs148211042 under standard blood bank conditions. STUDY DESIGN: Whole blood units were collected from 6 FP individuals and 11 controls and processed into RCC in additive solution. RCC were sampled and tested over cold storage for full blood count, extracellular potassium, glucose, lactate, microvesicle release, deformability, hemolysis, pH, adenosine triphosphate, and 2,3-diphosphoglycerate. RESULTS: Screening of genotyped cohorts identified that variant rs148211042 is present in 1 in 394 British citizens of European ancestry. FP RCC had significantly higher supernatant potassium at all time points from day 3 onwards (p < .001) and higher mean cell volume (p = .032) than controls. The initial rate of potassium release was higher in FP RCC; supernatant potassium reached 46.0 (23.8-57.6) mmol/L (mean [range]) by day 5, increasing to 68.9 (58.8-73.7) mmol/L by day 35. Other quality parameters were not significantly different between FP RCC and controls. CONCLUSION: These data suggest that if a blood donor has FP, reducing the RCC shelf-life to 5 days may be insufficient to reduce the risk of hyperkalemia in clinical scenarios such as neonatal large volume transfusion.

A rare coding mutation in the MAST2 gene causes venous thrombosis in a French family with unexplained thrombophilia: The Breizh MAST2 Arg89Gln variant.

Rare variants outside the classical coagulation cascade might cause inherited thrombosis. We aimed to identify the variant(s) causing venous thromboembolism (VTE) in a family with multiple relatives affected with unprovoked VTE and no thrombophilia defects. We identified by whole exome sequencing an extremely rare Arg to Gln variant (Arg89Gln) in the Microtubule Associated Serine/Threonine Kinase 2 (MAST2) gene that segregates with VTE in the family. Free-tissue factor pathway inhibitor (f-TFPI) plasma levels were significantly decreased in affected family members compared to healthy relatives. Conversely, plasminogen activator inhibitor-1 (PAI-1) levels were significantly higher in affected members than in healthy relatives. RNA sequencing analysis of RNA interference experimental data conducted in endothelial cells revealed that, of the 13,387 detected expressed genes, 2,354 have their level of expression modified by MAST2 knockdown, including SERPINE1 coding for PAI-1 and TFPI. In HEK293 cells overexpressing the MAST2 Gln89 variant, TFPI and SERPINE1 promoter activities were respectively lower and higher than in cells overexpressing the MAST2 wild type. This study identifies a novel thrombophilia-causing Arg89Gln variant in the MAST2 gene that is here proposed as a new molecular player in the etiology of VTE by interfering with hemostatic balance of endothelial cells.

Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13 037 individuals enrolled in the National Institute for Health Research (NIHR) BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor-encoding gene IKZF5 and thrombocytopenia. We report 5 causal missense variants in or near IKZF5 zinc fingers, of which 2 occurred de novo and 3 co-segregated in 3 pedigrees. A canonical DNA-zinc finger binding model predicts that 3 of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared with wild-type IKZF5, and electron microscopy revealed a reduced quantity of α granules in normally sized platelets. Proplatelet formation was reduced in megakaryocytes from 7 cases relative to 6 controls. Comparison of RNA-sequencing data from platelets, monocytes, neutrophils, and CD4+ T cells from 3 cases and 14 healthy controls showed 1194 differentially expressed genes in platelets but only 4 differentially expressed genes in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans.

Complex structural variants in Mendelian disorders: identification and breakpoint resolution using short- and long-read genome sequencing.

BACKGROUND: Studies have shown that complex structural variants (cxSVs) contribute to human genomic variation and can cause Mendelian disease. We aimed to identify cxSVs relevant to Mendelian disease using short-read whole-genome sequencing (WGS), resolve the precise variant configuration and investigate possible mechanisms of cxSV formation. METHODS: We performed short-read WGS and analysis of breakpoint junctions to identify cxSVs in a cohort of 1324 undiagnosed rare disease patients. Long-read WGS and gene expression analysis were used to resolve one case. RESULTS: We identified three pathogenic cxSVs: a de novo duplication-inversion-inversion-deletion affecting ARID1B, a de novo deletion-inversion-duplication affecting HNRNPU and a homozygous deletion-inversion-deletion affecting CEP78. Additionally, a de novo duplication-inversion-duplication overlapping CDKL5 was resolved by long-read WGS demonstrating the presence of both a disrupted and an intact copy of CDKL5 on the same allele, and gene expression analysis showed both parental alleles of CDKL5 were expressed. Breakpoint analysis in all the cxSVs revealed both microhomology and longer repetitive elements. CONCLUSIONS: Our results corroborate that cxSVs cause Mendelian disease, and we recommend their consideration during clinical investigations. We show that resolution of breakpoints can be critical to interpret pathogenicity and present evidence of replication-based mechanisms in cxSV formation.

Selection-free gene repair after adenoviral vector transduction of designer nucleases: rescue of dystrophin synthesis in DMD muscle cell populations.

Duchenne muscular dystrophy (DMD) is a fatal X-linked muscle-wasting disorder caused by mutations in the 2.4 Mb dystrophin-encoding DMD gene. The integration of gene delivery and gene editing technologies based on viral vectors and sequence-specific designer nucleases, respectively, constitutes a potential therapeutic modality for permanently repairing defective DMD alleles in patient-derived myogenic cells. Therefore, we sought to investigate the feasibility of combining adenoviral vectors (AdVs) with CRISPR/Cas9 RNA-guided nucleases (RGNs) alone or together with transcriptional activator-like effector nucleases (TALENs), for endogenous DMD repair through non-homologous end-joining (NHEJ). The strategies tested involved; incorporating small insertions or deletions at out-of-frame sequences for reading frame resetting, splice acceptor knockout for DNA-level exon skipping, and RGN-RGN or RGN-TALEN multiplexing for targeted exon(s) removal. We demonstrate that genome editing based on the activation and recruitment of the NHEJ DNA repair pathway after AdV delivery of designer nuclease genes, is a versatile and robust approach for repairing DMD mutations in bulk populations of patient-derived muscle progenitor cells (up to 37% of corrected DMD templates). These results open up a DNA-level genetic medicine strategy in which viral vector-mediated transient designer nuclease expression leads to permanent and regulated dystrophin synthesis from corrected native DMD alleles.