Bayesian Inference Associates Rare KDR Variants with Specific Phenotypes in Pulmonary Arterial Hypertension.

Background - Approximately 25% of patients with pulmonary arterial hypertension (PAH) have been found to harbor rare mutations in disease-causing genes. To identify missing heritability in PAH we integrated deep phenotyping with whole-genome sequencing data using Bayesian statistics. Methods - We analyzed 13,037 participants enrolled in the NIHR BioResource - Rare Diseases (NBR) study, of which 1,148 were recruited to the PAH domain. To test for genetic associations between genes and selected phenotypes of pulmonary hypertension (PH), we used the Bayesian rare-variant association method BeviMed. Results - Heterozygous, high impact, likely loss-of-function variants in the Kinase Insert Domain Receptor (KDR) gene were strongly associated with significantly reduced transfer coefficient for carbon monoxide (KCO, posterior probability (PP)=0.989) and older age at diagnosis (PP=0.912). We also provide evidence for familial segregation of a rare nonsense KDR variant with these phenotypes. On computed tomographic imaging of the lungs, a range of parenchymal abnormalities were observed in the five patients harboring these predicted deleterious variants in KDR. Four additional PAH cases with rare likely loss-of-function variants in KDR were independently identified in the US PAH Biobank cohort with similar phenotypic characteristics. Conclusions - The Bayesian inference approach allowed us to independently validate KDR, which encodes for the Vascular Endothelial Growth Factor Receptor 2 (VEGFR2), as a novel PAH candidate gene. Furthermore, this approach specifically associated high impact likely loss-of-function variants in the genetically constrained gene with distinct phenotypes. These findings provide evidence for KDR being a clinically actionable PAH gene and further support the central role of the vascular endothelium in the pathobiology of PAH.

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.

ADA2 deficiency complicated by EBV-driven lymphoproliferative disease.

A 29-year old male with recurrent respiratory and skin infections, anaemia and neutropaenia during childhood required immunoglobulin replacement for antibody deficiency from age 16. He remained relatively well until age 28 when he presented with a two-week history of fatigue, sore throat, fever and productive cough. He was found to have EBV viraemia and splenomegaly and a diagnosis of EBV-driven lymphoproliferative disease was made following bone marrow trephine. Family history was notable with three siblings: a healthy sister and two brothers with anaemia and neutropaenia; one who succumbed to septicaemia secondary to neutropaenic enterocolitis age 5 and another who developed intestinal vasculitis and antibody deficiency and had a successful haemopoetic stem cell transplant. The proband's DNA underwent targeted sequencing of 279 genes associated with immunodeficiency (GRID panel). The best candidates were two ADA2 variants, p.Arg169Gln (R169Q) and p.Asn370Lys (N370K). Sanger sequencing and co-segregation of variants in the parents, unaffected sister and all three affected brothers was fully consistent with compound heterozygous inheritance. Subsequent whole genome sequencing of the proband identified no other potential causal variants. ADA2 activity was consistent with a diagnosis of ADA2 deficiency in affected family members. This is the first description of EBV-driven lymphoproliferative disease in ADA2 deficiency. ADA2 deficiency may cause susceptibility to severe EBV-induced disease and we would recommend that EBV status and viral load is monitored in patients with this diagnosis and allogeneic SCT is considered at an early stage for patients whose ADA2 deficiency is associated with significant complications.

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.

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.

A high-throughput sequencing test for diagnosing inherited bleeding, thrombotic, and platelet disorders.

Inherited bleeding, thrombotic, and platelet disorders (BPDs) are diseases that affect ∼300 individuals per million births. With the exception of hemophilia and von Willebrand disease patients, a molecular analysis for patients with a BPD is often unavailable. Many specialized tests are usually required to reach a putative diagnosis and they are typically performed in a step-wise manner to control costs. This approach causes delays and a conclusive molecular diagnosis is often never reached, which can compromise treatment and impede rapid identification of affected relatives. To address this unmet diagnostic need, we designed a high-throughput sequencing platform targeting 63 genes relevant for BPDs. The platform can call single nucleotide variants, short insertions/deletions, and large copy number variants (though not inversions) which are subjected to automated filtering for diagnostic prioritization, resulting in an average of 5.34 candidate variants per individual. We sequenced 159 and 137 samples, respectively, from cases with and without previously known causal variants. Among the latter group, 61 cases had clinical and laboratory phenotypes indicative of a particular molecular etiology, whereas the remainder had an a priori highly uncertain etiology. All previously detected variants were recapitulated and, when the etiology was suspected but unknown or uncertain, a molecular diagnosis was reached in 56 of 61 and only 8 of 76 cases, respectively. The latter category highlights the need for further research into novel causes of BPDs. The ThromboGenomics platform thus provides an affordable DNA-based test to diagnose patients suspected of having a known inherited BPD.

A dominant gain-of-function mutation in universal tyrosine kinase SRC causes thrombocytopenia, myelofibrosis, bleeding, and bone pathologies.

The Src family kinase (SFK) member SRC is a major target in drug development because it is activated in many human cancers, yet deleterious SRC germline mutations have not been reported. We used genome sequencing and Human Phenotype Ontology patient coding to identify a gain-of-function mutation in SRC causing thrombocytopenia, myelofibrosis, bleeding, and bone pathologies in nine cases. Modeling of the E527K substitution predicts loss of SRC's self-inhibitory capacity, which we confirmed with in vitro studies showing increased SRC kinase activity and enhanced Tyr(419) phosphorylation in COS-7 cells overexpressing E527K SRC. The active form of SRC predominates in patients' platelets, resulting in enhanced overall tyrosine phosphorylation. Patients with myelofibrosis have hypercellular bone marrow with trilineage dysplasia, and their stem cells grown in vitro form more myeloid and megakaryocyte (MK) colonies than control cells. These MKs generate platelets that are dysmorphic, low in number, highly variable in size, and have a paucity of α-granules. Overactive SRC in patient-derived MKs causes a reduction in proplatelet formation, which can be rescued by SRC kinase inhibition. Stem cells transduced with lentiviral E527K SRC form MKs with a similar defect and enhanced tyrosine phosphorylation levels. Patient-derived and E527K-transduced MKs show Y419 SRC-positive stained podosomes that induce altered actin organization. Expression of mutated src in zebrafish recapitulates patients' blood and bone phenotypes. Similar studies of platelets and MKs may reveal the mechanism underlying the severe bleeding frequently observed in cancer patients treated with next-generation SFK inhibitors.

A gain-of-function variant in DIAPH1 causes dominant macrothrombocytopenia and hearing loss.

Macrothrombocytopenia (MTP) is a heterogeneous group of disorders characterized by enlarged and reduced numbers of circulating platelets, sometimes resulting in abnormal bleeding. In most MTP, this phenotype arises because of altered regulation of platelet formation from megakaryocytes (MKs). We report the identification of DIAPH1, which encodes the Rho-effector diaphanous-related formin 1 (DIAPH1), as a candidate gene for MTP using exome sequencing, ontological phenotyping, and similarity regression. We describe 2 unrelated pedigrees with MTP and sensorineural hearing loss that segregate with a DIAPH1 R1213* variant predicting partial truncation of the DIAPH1 diaphanous autoregulatory domain. The R1213* variant was linked to reduced proplatelet formation from cultured MKs, cell clustering, and abnormal cortical filamentous actin. Similarly, in platelets, there was increased filamentous actin and stable microtubules, indicating constitutive activation of DIAPH1. Overexpression of DIAPH1 R1213* in cells reproduced the cytoskeletal alterations found in platelets. Our description of a novel disorder of platelet formation and hearing loss extends the repertoire of DIAPH1-related disease and provides new insight into the autoregulation of DIAPH1 activity.

Human phenotype ontology annotation and cluster analysis to unravel genetic defects in 707 cases with unexplained bleeding and platelet disorders.

BACKGROUND: Heritable bleeding and platelet disorders (BPD) are heterogeneous and frequently have an unknown genetic basis. The BRIDGE-BPD study aims to discover new causal genes for BPD by high throughput sequencing using cluster analyses based on improved and standardised deep, multi-system phenotyping of cases. METHODS: We report a new approach in which the clinical and laboratory characteristics of BPD cases are annotated with adapted Human Phenotype Ontology (HPO) terms. Cluster analyses are then used to characterise groups of cases with similar HPO terms and variants in the same genes. RESULTS: We show that 60% of index cases with heritable BPD enrolled at 10 European or US centres were annotated with HPO terms indicating abnormalities in organ systems other than blood or blood-forming tissues, particularly the nervous system. Cases within pedigrees clustered closely together on the bases of their HPO-coded phenotypes, as did cases sharing several clinically suspected syndromic disorders. Cases subsequently found to harbour variants in ACTN1 also clustered closely, even though diagnosis of this recently described disorder was not possible using only the clinical and laboratory data available to the enrolling clinician. CONCLUSIONS: These findings validate our novel HPO-based phenotype clustering methodology for known BPD, thus providing a new discovery tool for BPD of unknown genetic basis. This approach will also be relevant for other rare diseases with significant genetic heterogeneity.

Common genetic variants do not associate with CAD in familial hypercholesterolemia.

In recent years, multiple loci dispersed on the genome have been shown to be associated with coronary artery disease (CAD). We investigated whether these common genetic variants also hold value for CAD prediction in a large cohort of patients with familial hypercholesterolemia (FH). We genotyped a total of 41 single-nucleotide polymorphisms (SNPs) in 1701 FH patients, of whom 482 patients (28.3%) had at least one coronary event during an average follow up of 66 years. The association of each SNP with event-free survival time was calculated with a Cox proportional hazard model. In the cardiovascular disease risk factor adjusted analysis, the most significant SNP was rs1122608:G>T in the SMARCA4 gene near the LDL-receptor (LDLR) gene, with a hazard ratio for CAD risk of 0.74 (95% CI 0.49-0.99; P-value 0.021). However, none of the SNPs reached the Bonferroni threshold. Of all the known CAD loci analyzed, the SMARCA4 locus near the LDLR had the strongest negative association with CAD in this high-risk FH cohort. The effect is contrary to what was expected. None of the other loci showed association with CAD.

SMIM1 underlies the Vel blood group and influences red blood cell traits.

The blood group Vel was discovered 60 years ago, but the underlying gene is unknown. Individuals negative for the Vel antigen are rare and are required for the safe transfusion of patients with antibodies to Vel. To identify the responsible gene, we sequenced the exomes of five individuals negative for the Vel antigen and found that four were homozygous and one was heterozygous for a low-frequency 17-nucleotide frameshift deletion in the gene encoding the 78-amino-acid transmembrane protein SMIM1. A follow-up study showing that 59 of 64 Vel-negative individuals were homozygous for the same deletion and expression of the Vel antigen on SMIM1-transfected cells confirm SMIM1 as the gene underlying the Vel blood group. An expression quantitative trait locus (eQTL), the common SNP rs1175550 contributes to variable expression of the Vel antigen (P = 0.003) and influences the mean hemoglobin concentration of red blood cells (RBCs; P = 8.6 × 10(-15)). In vivo, zebrafish with smim1 knockdown showed a mild reduction in the number of RBCs, identifying SMIM1 as a new regulator of RBC formation. Our findings are of immediate relevance, as the homozygous presence of the deletion allows the unequivocal identification of Vel-negative blood donors.

Compound inheritance of a low-frequency regulatory SNP and a rare null mutation in exon-junction complex subunit RBM8A causes TAR syndrome.

The exon-junction complex (EJC) performs essential RNA processing tasks. Here, we describe the first human disorder, thrombocytopenia with absent radii (TAR), caused by deficiency in one of the four EJC subunits. Compound inheritance of a rare null allele and one of two low-frequency SNPs in the regulatory regions of RBM8A, encoding the Y14 subunit of EJC, causes TAR. We found that this inheritance mechanism explained 53 of 55 cases (P < 5 × 10(-228)) of the rare congenital malformation syndrome. Of the 53 cases with this inheritance pattern, 51 carried a submicroscopic deletion of 1q21.1 that has previously been associated with TAR, and two carried a truncation or frameshift null mutation in RBM8A. We show that the two regulatory SNPs result in diminished RBM8A transcription in vitro and that Y14 expression is reduced in platelets from individuals with TAR. Our data implicate Y14 insufficiency and, presumably, an EJC defect as the cause of TAR syndrome.

Generation and characterization of a live attenuated enterotoxigenic Escherichia coli combination vaccine expressing six colonization factors and heat-labile toxin subunit B.

Live attenuated oral enterotoxigenic Escherichia coli (ETEC) vaccines have been demonstrated to be safe and immunogenic in human volunteers and to provide a viable approach to provide protection against this important pathogen. This report describes the construction of new ETEC vaccine candidate strains from recent clinical isolates and their characterization. All known genes for ETEC toxins were removed, and attenuating deletion mutations were made in the aroC, ompC, and ompF chromosomal genes. An isolate expressing coli surface antigen 2 (CS2), CS3, heat-labile toxin (LT), heat-stable toxin (ST), and enteroaggregative Escherichia coli heat-stable toxin 1 (EAST1) was attenuated to generate ACAM2007. The subsequent insertion of the operon encoding CS1 created ACAM2017, and this was further modified by the addition of an expression cassette containing the eltB gene, encoding a pentamer of B subunits of LT (LTB), to generate ACAM2027. Another isolate expressing CS5, CS6, LT, ST, and EAST1 was attenuated to generate ACAM2006, from which a lysogenic prophage was deleted to create ACAM2012 and an LTB gene was introduced to form ACAM2022. Finally, a previously described vaccine strain, ACAM2010, had the eltB gene incorporated to generate ACAM2025. All recombinant genes were incorporated into the chromosomal sites of the attenuating mutations to ensure maximal genetic stability. The expression of the recombinant antigens and the changes in plasmids accompanying the deletion of toxin genes are described. Strains ACAM2025, ACAM2022, and ACAM2027 have been combined to create the ETEC vaccine formulation ACE527, which has recently successfully completed a randomized, double-blind, placebo-controlled phase I trial and is currently undergoing a randomized, double-blind placebo-controlled phase II challenge trial, both in healthy adult volunteers.

Stabilization of a plasmid coding for a heterologous antigen in Salmonella enterica serotype typhi vaccine strain CVD908-htrA by using site-specific recombination.

A gene cassette incorporating the crs-rsd site-specific recombination system from the Salmonella enterica subsp. enterica serovar Dublin virulence plasmid improved the inheritance in S. enterica serotype Typhi strain CVD908-htrA of a multicopy plasmid expression vector. Use of this recombination cassette may improve expression of heterologous antigens from multicopy plasmid expression vectors in attenuated bacterial vaccine strains.

Construction and phase I clinical evaluation of the safety and immunogenicity of a candidate enterotoxigenic Escherichia coli vaccine strain expressing colonization factor antigen CFA/I.

Oral delivery of toxin-negative derivatives of enterotoxigenic Escherichia coli (ETEC) that express colonization factor antigens (CFA) with deletions of the aroC, ompC, ompF, and toxin genes may be an effective approach to vaccination against ETEC-associated diarrhea. We describe the creation and characterization of an attenuated CFA/I-expressing ETEC vaccine candidate, ACAM2010, from a virulent isolate in which the heat-stable enterotoxin (ST) and CFA/I genes were closely linked and on the same virulence plasmid as the enteroaggregative E. coli heat-stable toxin (EAST1) gene. A new suicide vector (pJCB12) was constructed and used to delete the ST and EAST1 genes and to introduce defined deletion mutations into the aroC, ompC, and ompF chromosomal genes. A phase I trial, consisting of an open-label dose escalation phase in 18 adult outpatient volunteers followed by a placebo-controlled double-blind phase in an additional 31 volunteers, was conducted. The vaccine was administered in two formulations, fresh culture and frozen suspension. These were both well tolerated, with no evidence of significant adverse events related to vaccination. Immunoglobulin A (IgA) and IgG antibody-secreting cells specific for CFA/I were assayed by ELISPOT. Positive responses (greater than twofold increase) were seen in 27 of 37 (73%) subjects who received the highest dose level of vaccine (nominally 5 x 10(9) CFU). Twenty-nine of these volunteers were secreting culturable vaccine organisms at day 3 following vaccination; five were still positive on day 7, with a single isolation on day 13. This live attenuated bacterial vaccine is safe and immunogenic in healthy adult volunteers.