Prospective, international, multisite comparison of platelet isolation techniques for genome-wide transcriptomics: communication from the SSC of the ISTH.

Genome-wide platelet transcriptomics is increasingly used to uncover new aspects of platelet biology and as a diagnostic and prognostic tool. Nevertheless, platelet isolation methods for transcriptomic studies are not standardized, introducing challenges for cross-study comparisons, data integration, and replication. In this prospective multicenter study, called "Standardizing Platelet Transcriptomics for Discovery, Diagnostics, and Therapeutics in the Thrombosis and Hemostasis Community (STRIDE)" by the International Society on Thrombosis and Haemostasis Scientific and Standardization Committees, we assessed how 3 of the most commonly used platelet isolation protocols influence metrics from next-generation bulk RNA sequencing and functional assays. Compared with washing alone, more stringent removal of leukocytes by anti-CD45 beads or PALL filters resulted in a sufficient quantity of RNA for next-generation sequencing and similar quality of RNA sequencing metrics. Importantly, stringent removal of leukocytes resulted in the lower relative expression of known leukocyte-specific genes and the higher relative expression of known platelet-specific genes. The results were consistent across enrolling sites, suggesting that the techniques are transferrable and reproducible. Moreover, all 3 isolation techniques did not influence basal platelet reactivity, but agonist-induced integrin αIIbß3 activation is reduced by anti-CD45 bead isolation compared with washing alone. In conclusion, the isolation technique chosen influences genome-wide transcriptional and functional assays in platelets. These results should help the research community make informed choices about platelet isolation techniques in their own platelet studies.

Mutations in the U4 snRNA gene RNU4-2 cause one of the most prevalent monogenic neurodevelopmental disorders.

Most people with intellectual disability (ID) do not receive a molecular diagnosis following genetic testing. To identify new etiologies of ID, we performed a genetic association analysis comparing the burden of rare variants in 41,132 noncoding genes between 5,529 unrelated cases and 46,401 unrelated controls. RNU4-2, which encodes U4 small nuclear RNA, a critical component of the spliceosome, was the most strongly associated gene. We implicated de novo variants among 47 cases in two regions of RNU4-2 in the etiology of a syndrome characterized by ID, microcephaly, short stature, hypotonia, seizures and motor delay. We replicated this finding in three collections, bringing the number of unrelated cases to 73. Analysis of national genomic diagnostic data showed RNU4-2 to be a more common etiological gene for neurodevelopmental abnormality than any previously reported autosomal gene. Our findings add to the growing evidence of spliceosome dysfunction in the etiologies of neurological disorders.

Biallelic loss-of-function variants of SLC12A9 cause lysosome dysfunction and a syndromic neurodevelopmental disorder.

PURPOSE: Pathogenic variants of FIG4 generate enlarged lysosomes and neurological and developmental disorders. To identify additional genes regulating lysosomal volume, we carried out a genome-wide activation screen to detect suppression of enlarged lysosomes in FIG4-/- cells. METHODS: The CRISPR-a gene activation screen utilized sgRNAs from the promoters of protein-coding genes. Fluorescence-activated cell sorting separated cells with correction of the enlarged lysosomes from uncorrected cells. Patient variants of SLC12A9 were identified by exome or genome sequencing and studied by segregation analysis and clinical characterization. RESULTS: Overexpression of SLC12A9, a solute co-transporter, corrected lysosomal swelling in FIG4-/- cells. SLC12A9 (NP_064631.2) colocalized with LAMP2 at the lysosome membrane. Biallelic variants of SLC12A9 were identified in 3 unrelated probands with neurodevelopmental disorders. Common features included intellectual disability, skeletal and brain structural abnormalities, congenital heart defects, and hypopigmented hair. Patient 1 was homozygous for nonsense variant p.(Arg615∗), patient 2 was compound heterozygous for p.(Ser109Lysfs∗20) and a large deletion, and proband 3 was compound heterozygous for p.(Glu290Glyfs∗36) and p.(Asn552Lys). Fibroblasts from proband 1 contained enlarged lysosomes that were corrected by wild-type SLC12A9 cDNA. Patient variant p.(Asn552Lys) failed to correct the lysosomal defect. CONCLUSION: Impaired function of SLC12A9 results in enlarged lysosomes and a recessive disorder with a recognizable neurodevelopmental phenotype.

A Deep Intronic PKHD1 Variant Identified by SpliceAI in a Deceased Neonate With Autosomal Recessive Polycystic Kidney Disease.

The etiologies of newborn deaths in neonatal intensive care units usually remain unknown, even after genetic testing. Whole-genome sequencing, combined with artificial intelligence-based methods for predicting the effects of non-coding variants, provide an avenue for resolving these deaths. Using one such method, SpliceAI, we identified a maternally inherited deep intronic PKHD1 splice variant (chr6:52030169T>C), in trans with a pathogenic missense variant (p.Thr36Met), in a newborn who died of autosomal recessive polycystic kidney disease at age 2 days. We validated the deep intronic variant's impact in maternal urine-derived cells expressing PKHD1. Reverse transcription polymerase chain reaction followed by Sanger sequencing showed that the variant causes inclusion of 147bp of the canonical intron between exons 29 and 30 of PKHD1 into the mRNA, including a premature stop codon. Allele-specific expression analysis at a heterozygous site in the mother showed that the mutant allele completely suppresses canonical splicing. In an unrelated healthy control, there was no evidence of transcripts including the novel splice junction. We returned a diagnostic report to the parents, who underwent in vitro embryo selection.

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.

Genetic association analysis of 77,539 genomes reveals rare disease etiologies.

The genetic etiologies of more than half of rare diseases remain unknown. Standardized genome sequencing and phenotyping of large patient cohorts provide an opportunity for discovering the unknown etiologies, but this depends on efficient and powerful analytical methods. We built a compact database, the 'Rareservoir', containing the rare variant genotypes and phenotypes of 77,539 participants sequenced by the 100,000 Genomes Project. We then used the Bayesian genetic association method BeviMed to infer associations between genes and each of 269 rare disease classes assigned by clinicians to the participants. We identified 241 known and 19 previously unidentified associations. We validated associations with ERG, PMEPA1 and GPR156 by searching for pedigrees in other cohorts and using bioinformatic and experimental approaches. We provide evidence that (1) loss-of-function variants in the Erythroblast Transformation Specific (ETS)-family transcription factor encoding gene ERG lead to primary lymphoedema, (2) truncating variants in the last exon of transforming growth factor-ß regulator PMEPA1 result in Loeys-Dietz syndrome and (3) loss-of-function variants in GPR156 give rise to recessive congenital hearing impairment. The Rareservoir provides a lightweight, flexible and portable system for synthesizing the genetic and phenotypic data required to study rare disease cohorts with tens of thousands of participants.

MitoPhen database: a human phenotype ontology-based approach to identify mitochondrial DNA diseases.

Diagnosing mitochondrial disorders remains challenging. This is partly because the clinical phenotypes of patients overlap with those of other sporadic and inherited disorders. Although the widespread availability of genetic testing has increased the rate of diagnosis, the combination of phenotypic and genetic heterogeneity still makes it difficult to reach a timely molecular diagnosis with confidence. An objective, systematic method for describing the phenotypic spectra for each variant provides a potential solution to this problem. We curated the clinical phenotypes of 6688 published individuals with 89 pathogenic mitochondrial DNA (mtDNA) mutations, collating 26 348 human phenotype ontology (HPO) terms to establish the MitoPhen database. This enabled a hypothesis-free definition of mtDNA clinical syndromes, an overview of heteroplasmy-phenotype relationships, the identification of under-recognized phenotypes, and provides a publicly available reference dataset for objective clinical comparison with new patients using the HPO. Studying 77 patients with independently confirmed positive mtDNA diagnoses and 1083 confirmed rare disease cases with a non-mitochondrial nuclear genetic diagnosis, we show that HPO-based phenotype similarity scores can distinguish these two classes of rare disease patients with a false discovery rate <10% at a sensitivity of 80%. Enriching the MitoPhen database with more patients will improve predictions for increasingly rare variants.

NRG1 fusions in breast cancer.

BACKGROUND: NRG1 gene fusions may be clinically actionable, since cancers carrying the fusion transcripts can be sensitive to tyrosine kinase inhibitors. The NRG1 gene encodes ligands for the HER2(ERBB2)-ERBB3 heterodimeric receptor tyrosine kinase, and the gene fusions are thought to lead to autocrine stimulation of the receptor. The NRG1 fusion expressed in the breast cancer cell line MDA-MB-175 serves as a model example of such fusions, showing the proposed autocrine loop and exceptional drug sensitivity. However, its structure has not been properly characterised, its oncogenic activity has not been fully explained, and there is limited data on such fusions in breast cancer. METHODS: We analysed genomic rearrangements and transcripts of NRG1 in MDA-MB-175 and a panel of 571 breast cancers. RESULTS: We found that the MDA-MB-175 fusion-originally reported as a DOC4(TENM4)-NRG1 fusion, lacking the cytoplasmic tail of NRG1-is in reality a double fusion, PPP6R3-TENM4-NRG1, producing multiple transcripts, some of which include the cytoplasmic tail. We hypothesise that many NRG1 fusions may be oncogenic not for lacking the cytoplasmic domain but because they do not encode NRG1's nuclear-localised form. The fusion in MDA-MB-175 is the result of a very complex genomic rearrangement, which we partially characterised, that creates additional expressed gene fusions, RSF1-TENM4, TPCN2-RSF1, and MRPL48-GAB2. We searched for NRG1 rearrangements in 571 breast cancers subjected to genome sequencing and transcriptome sequencing and found four cases (0.7%) with fusions, WRN-NRG1, FAM91A1-NRG1, ARHGEF39-NRG1, and ZNF704-NRG1, all splicing into NRG1 at the same exon as in MDA-MB-175. However, the WRN-NRG1 and ARHGEF39-NRG1 fusions were out of frame. We identified rearrangements of NRG1 in many more (8% of) cases that seemed more likely to inactivate than to create activating fusions, or whose outcome could not be predicted because they were complex, or both. This is not surprising because NRG1 can be pro-apoptotic and is inactivated in some breast cancers. CONCLUSIONS: Our results highlight the complexity of rearrangements of NRG1 in breast cancers and confirm that some do not activate but inactivate. Careful interpretation of NRG1 rearrangements will therefore be necessary for appropriate patient management.

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/.

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.

Author Correction: Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Novel manifestations of immune dysregulation and granule defects in gray platelet syndrome.

Gray platelet syndrome (GPS) is a rare recessive disorder caused by biallelic variants in NBEAL2 and characterized by bleeding symptoms, the absence of platelet α-granules, splenomegaly, and bone marrow (BM) fibrosis. Due to the rarity of GPS, it has been difficult to fully understand the pathogenic processes that lead to these clinical sequelae. To discern the spectrum of pathologic features, we performed a detailed clinical genotypic and phenotypic study of 47 patients with GPS and identified 32 new etiologic variants in NBEAL2. The GPS patient cohort exhibited known phenotypes, including macrothrombocytopenia, BM fibrosis, megakaryocyte emperipolesis of neutrophils, splenomegaly, and elevated serum vitamin B12 levels. Novel clinical phenotypes were also observed, including reduced leukocyte counts and increased presence of autoimmune disease and positive autoantibodies. There were widespread differences in the transcriptome and proteome of GPS platelets, neutrophils, monocytes, and CD4 lymphocytes. Proteins less abundant in these cells were enriched for constituents of granules, supporting a role for Nbeal2 in the function of these organelles across a wide range of blood cells. Proteomic analysis of GPS plasma showed increased levels of proteins associated with inflammation and immune response. One-quarter of plasma proteins increased in GPS are known to be synthesized outside of hematopoietic cells, predominantly in the liver. In summary, our data show that, in addition to the well-described platelet defects in GPS, there are immune defects. The abnormal immune cells may be the drivers of systemic abnormalities such as autoimmune disease.

A coagulation defect arising from heterozygous premature termination of tissue factor.

Tissue factor (TF) is the primary initiator of blood coagulation in vivo and the only blood coagulation factor for which a human genetic defect has not been described. As there are no routine clinical assays that capture the contribution of endogenous TF to coagulation initiation, the extent to which reduced TF activity contributes to unexplained bleeding is unknown. Using whole genome sequencing, we identified a heterozygous frameshift variant (p.Ser117HisfsTer10) in F3, the gene encoding TF, causing premature termination of TF (TFshort) in a woman with unexplained bleeding. Routine hematological laboratory evaluation of the proposita was normal. CRISPR-edited human induced pluripotent stem cells recapitulating the variant were differentiated into vascular smooth muscle and endothelial cells that demonstrated haploinsufficiency of TF. The variant F3 transcript is eliminated by nonsense-mediated decay. Neither overexpression nor addition of exogenous recombinant TFshort inhibited factor Xa or thrombin generation, excluding a dominant-negative mechanism. F3+/- mice provide an animal model of TF haploinsufficiency and exhibited prolonged bleeding times, impaired thrombus formation, and reduced survival following major injury. Heterozygous TF deficiency is present in at least 1 in 25,000 individuals and could limit coagulation initiation in undiagnosed individuals with abnormal bleeding but a normal routine laboratory evaluation.

A new pedigree with thrombomodulin-associated coagulopathy in which delayed fibrinolysis is partially attenuated by co-inherited TAFI deficiency.

BACKGROUND: Thrombomodulin-associated coagulopathy (TM-AC) is a rare bleeding disorder in which a single reported p.Cys537* variant in the thrombomodulin gene THBD causes high plasma thrombomodulin (TM) levels. High TM levels attenuate thrombin generation and delay fibrinolysis. OBJECTIVES: To report the characteristics of pedigree with a novel THBD variant causing TM-AC, and co-inherited deficiency of thrombin-activatable fibrinolysis inhibitor (TAFI). PATIENTS/METHODS: Identification of pathogenic variants in hemostasis genes by next-generation sequencing and case recall for deep phenotyping. RESULTS: Pedigree members with a previously reported THBD variant predicting p.Pro496Argfs*10 and chain truncation in TM transmembrane domain had abnormal bleeding and greatly increased plasma TM levels. Affected cases had attenuated thrombin generation and delayed fibrinolysis similar to previous reported TM_AC cases with THBD p.Cys537*. Coincidentally, some pedigree members also harbored a stop-gain variant in CPB2 encoding TAFI. This reduced plasma TAFI levels but was asymptomatic. Pedigree members with TM-AC caused by the p.Pro496Argfs*10 THBD variant and also TAFI deficiency had a partially attenuated delay in fibrinolysis, but no change in the defective thrombin generation. CONCLUSIONS: These data extend the reported genetic repertoire of TM-AC and establish a common molecular pathogenesis arising from high plasma levels of TM extra-cellular domain. The data further confirm that the delay in fibrinolysis associated with TM-AC is directly linked to increased TAFI activation. The combination of the rare variants in the pedigree members provides a unique genetic model to develop understanding of the thrombin-TM system and its regulation of TAFI.

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.

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.

Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans.

Several strands of evidence question the dogma that human mitochondrial DNA (mtDNA) is inherited exclusively down the maternal line, most recently in three families where several individuals harbored a 'heteroplasmic haplotype' consistent with biparental transmission. Here we report a similar genetic signature in 7 of 11,035 trios, with allelic fractions of 5-25%, implying biparental inheritance of mtDNA in 0.06% of offspring. However, analysing the nuclear whole genome sequence, we observe likely large rare or unique nuclear-mitochondrial DNA segments (mega-NUMTs) transmitted from the father in all 7 families. Independently detecting mega-NUMTs in 0.13% of fathers, we see autosomal transmission of the haplotype. Finally, we show the haplotype allele fraction can be explained by complex concatenated mtDNA-derived sequences rearranged within the nuclear genome. We conclude that rare cryptic mega-NUMTs can resemble paternally mtDNA heteroplasmy, but find no evidence of paternal transmission of mtDNA in humans.

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.

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.