Publisher Correction: Whole-genome sequencing of a sporadic primary immunodeficiency cohort.

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

Whole-genome sequencing of a sporadic primary immunodeficiency cohort.

Primary immunodeficiency (PID) is characterized by recurrent and often life-threatening infections, autoimmunity and cancer, and it poses major diagnostic and therapeutic challenges. Although the most severe forms of PID are identified in early childhood, most patients present in adulthood, typically with no apparent family history and a variable clinical phenotype of widespread immune dysregulation: about 25% of patients have autoimmune disease, allergy is prevalent and up to 10% develop lymphoid malignancies1-3. Consequently, in sporadic (or non-familial) PID genetic diagnosis is difficult and the role of genetics is not well defined. Here we address these challenges by performing whole-genome sequencing in a large PID cohort of 1,318 participants. An analysis of the coding regions of the genome in 886 index cases of PID found that disease-causing mutations in known genes that are implicated in monogenic PID occurred in 10.3% of these patients, and a Bayesian approach (BeviMed4) identified multiple new candidate PID-associated genes, including IVNS1ABP. We also examined the noncoding genome, and found deletions in regulatory regions that contribute to disease causation. In addition, we used a genome-wide association study to identify loci that are associated with PID, and found evidence for the colocalization of-and interplay between-novel high-penetrance monogenic variants and common variants (at the PTPN2 and SOCS1 loci). This begins to explain the contribution of common variants to the variable penetrance and phenotypic complexity that are observed in PID. Thus, using a cohort-based whole-genome-sequencing approach in the diagnosis of PID can increase diagnostic yield and further our understanding of the key pathways that influence immune responsiveness in humans.

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.

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.

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.

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.

Pathogenic NFKB2 variant in the ankyrin repeat domain (R635X) causes a variable antibody deficiency.

Genetic studies are identifying an increasing number of monogenic causes of Common Variable Immunodeficiency (CVID). Pathogenic variants in the C-terminus of NFKB2 have been identified in the subset of CVID patients whose immunodeficiency is associated with ectodermal dysplasia and central adrenal insufficiency. We describe 2 unrelated CVID pedigrees with 4 cases of pathogenic stop gain variants (c.1903C > T) in the ankyrin repeat domain (ARD) of NF-κB2, leading to a premature truncation of the protein at p.Arg635Term (R635X). By immunophenotyping and functional ex vivo B- and T-cell experiments we characterized the variant by reduced class-switched memory B-cell counts and immature plasmablasts, unable to produce IgG and IgA. Features of a poor proliferative T-cell response and reduced expansion of CD4+CXCR5+ T cells was only observed in the two clinically affected index cases without any clear clinical correlate. In conclusion, pathogenic stop variants in the ARD of NFKB2 can cause 'infection-only' CVID with an abnormal B-cell phenotype and a variable clinical penetrance.

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.

Loss-of-function nuclear factor κB subunit 1 (NFKB1) variants are the most common monogenic cause of common variable immunodeficiency in Europeans.

BACKGROUND: The genetic cause of primary immunodeficiency disease (PID) carries prognostic information. OBJECTIVE: We conducted a whole-genome sequencing study assessing a large proportion of the NIHR BioResource-Rare Diseases cohort. METHODS: In the predominantly European study population of principally sporadic unrelated PID cases (n = 846), a novel Bayesian method identified nuclear factor κB subunit 1 (NFKB1) as one of the genes most strongly associated with PID, and the association was explained by 16 novel heterozygous truncating, missense, and gene deletion variants. This accounted for 4% of common variable immunodeficiency (CVID) cases (n = 390) in the cohort. Amino acid substitutions predicted to be pathogenic were assessed by means of analysis of structural protein data. Immunophenotyping, immunoblotting, and ex vivo stimulation of lymphocytes determined the functional effects of these variants. Detailed clinical and pedigree information was collected for genotype-phenotype cosegregation analyses. RESULTS: Both sporadic and familial cases demonstrated evidence of the noninfective complications of CVID, including massive lymphadenopathy (24%), unexplained splenomegaly (48%), and autoimmune disease (48%), features prior studies correlated with worse clinical prognosis. Although partial penetrance of clinical symptoms was noted in certain pedigrees, all carriers have a deficiency in B-lymphocyte differentiation. Detailed assessment of B-lymphocyte numbers, phenotype, and function identifies the presence of an increased CD21low B-cell population. Combined with identification of the disease-causing variant, this distinguishes between healthy subjects, asymptomatic carriers, and clinically affected cases. CONCLUSION: We show that heterozygous loss-of-function variants in NFKB1 are the most common known monogenic cause of CVID, which results in a temporally progressive defect in the formation of immunoglobulin-producing B cells.

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.

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.

αIIbß3 variants defined by next-generation sequencing: predicting variants likely to cause Glanzmann thrombasthenia.

Next-generation sequencing is transforming our understanding of human genetic variation but assessing the functional impact of novel variants presents challenges. We analyzed missense variants in the integrin αIIbß3 receptor subunit genes ITGA2B and ITGB3 identified by whole-exome or -genome sequencing in the ThromboGenomics project, comprising ∼32,000 alleles from 16,108 individuals. We analyzed the results in comparison with 111 missense variants in these genes previously reported as being associated with Glanzmann thrombasthenia (GT), 20 associated with alloimmune thrombocytopenia, and 5 associated with aniso/macrothrombocytopenia. We identified 114 novel missense variants in ITGA2B (affecting ∼11% of the amino acids) and 68 novel missense variants in ITGB3 (affecting ∼9% of the amino acids). Of the variants, 96% had minor allele frequencies (MAF) < 0.1%, indicating their rarity. Based on sequence conservation, MAF, and location on a complete model of αIIbß3, we selected three novel variants that affect amino acids previously associated with GT for expression in HEK293 cells. αIIb P176H and ß3 C547G severely reduced αIIbß3 expression, whereas αIIb P943A partially reduced αIIbß3 expression and had no effect on fibrinogen binding. We used receiver operating characteristic curves of combined annotation-dependent depletion, Polyphen 2-HDIV, and sorting intolerant from tolerant to estimate the percentage of novel variants likely to be deleterious. At optimal cut-off values, which had 69-98% sensitivity in detecting GT mutations, between 27% and 71% of the novel αIIb or ß3 missense variants were predicted to be deleterious. Our data have implications for understanding the evolutionary pressure on αIIbß3 and highlight the challenges in predicting the clinical significance of novel missense variants.

Distinct epigenomic features in end-stage failing human hearts.

BACKGROUND: The epigenome refers to marks on the genome, including DNA methylation and histone modifications, that regulate the expression of underlying genes. A consistent profile of gene expression changes in end-stage cardiomyopathy led us to hypothesize that distinct global patterns of the epigenome may also exist. METHODS AND RESULTS: We constructed genome-wide maps of DNA methylation and histone-3 lysine-36 trimethylation (H3K36me3) enrichment for cardiomyopathic and normal human hearts. More than 506 Mb sequences per library were generated by high-throughput sequencing, allowing us to assign methylation scores to ≈28 million CG dinucleotides in the human genome. DNA methylation was significantly different in promoter CpG islands, intragenic CpG islands, gene bodies, and H3K36me3-enriched regions of the genome. DNA methylation differences were present in promoters of upregulated genes but not downregulated genes. H3K36me3 enrichment itself was also significantly different in coding regions of the genome. Specifically, abundance of RNA transcripts encoded by the DUX4 locus correlated to differential DNA methylation and H3K36me3 enrichment. In vitro, Dux gene expression was responsive to a specific inhibitor of DNA methyltransferase, and Dux siRNA knockdown led to reduced cell viability. CONCLUSIONS: Distinct epigenomic patterns exist in important DNA elements of the cardiac genome in human end-stage cardiomyopathy. The epigenome may control the expression of local or distal genes with critical functions in myocardial stress response. If epigenomic patterns track with disease progression, assays for the epigenome may be useful for assessing prognosis in heart failure. Further studies are needed to determine whether and how the epigenome contributes to the development of cardiomyopathy.

Xenopus oocytes reactivate muscle gene transcription in transplanted somatic nuclei independently of myogenic factors.

Transplantation into eggs or oocytes is an effective means of achieving the reprogramming of somatic cell nuclei. We ask here whether the provision of gene-specific transcription factors forms part of the mechanism by which a gene that is repressed in somatic cells is transcribed in oocytes. We find that M1 oocytes have an extremely strong transcription-inducing activity. They cause muscle genes of nuclei from non-muscle somatic cells, after injection into oocytes, to be transcribed to nearly the same extent as muscle genes in muscle cells. We show, surprisingly, that the myogenic factor MyoD and other known myogenic factors are not required to induce the transcription of muscle genes in a range of non-muscle somatic cell nuclei after transplantation to Xenopus oocytes. The overexpression of Id, a dominant-negative repressor of MyoD, prevents maternal MyoD from binding to its consensus sequences; nevertheless, muscle genes are activated in somatic nuclei to the same extent as without Id. We conclude that M1 oocytes can reprogram somatic nuclei in a different way to other experimental procedures: oocytes do not suppress the transcription of inappropriate genes and they activate a gene without the help of its known transcription factors. We suggest that these characteristics might be a special property of amphibian oocytes, and possibly of oocytes in general.

Immunodeficiency, autoimmunity, and increased risk of B cell malignancy in humans with TRAF3 mutations.

Tumor necrosis factor receptor-associated factor 3 (TRAF3) is a central regulator of immunity. TRAF3 is often somatically mutated in B cell malignancies, but its role in human immunity is not defined. Here, in five unrelated families, we describe an immune dysregulation syndrome of recurrent bacterial infections, autoimmunity, systemic inflammation, B cell lymphoproliferation, and hypergammaglobulinemia. Affected individuals each had monoallelic mutations in TRAF3 that reduced TRAF3 expression. Immunophenotyping showed that patients' B cells were dysregulated, exhibiting increased nuclear factor-κB 2 activation, elevated mitochondrial respiration, and heightened inflammatory responses. Patients had mild CD4+ T cell lymphopenia, with a reduced proportion of naïve T cells but increased regulatory T cells and circulating T follicular helper cells. Guided by this clinical phenotype, targeted analyses demonstrated that common genetic variants, which also reduce TRAF3 expression, are associated with an increased risk of B cell malignancies, systemic lupus erythematosus, higher immunoglobulin levels, and bacterial infections in the wider population. Reduced TRAF3 conveys disease risks by driving B cell hyperactivity via intrinsic activation of multiple intracellular proinflammatory pathways and increased mitochondrial respiration, with a likely contribution from dysregulated T cell help. Thus, we define monogenic TRAF3 haploinsufficiency syndrome and demonstrate how common TRAF3 variants affect a range of human diseases.

Development and validation of a universal blood donor genotyping platform: a multinational prospective study.

Each year, blood transfusions save millions of lives. However, under current blood-matching practices, sensitization to non-self-antigens is an unavoidable adverse side effect of transfusion. We describe a universal donor typing platform that could be adopted by blood services worldwide to facilitate a universal extended blood-matching policy and reduce sensitization rates. This DNA-based test is capable of simultaneously typing most clinically relevant red blood cell (RBC), human platelet (HPA), and human leukocyte (HLA) antigens. Validation was performed, using samples from 7927 European, 27 South Asian, 21 East Asian, and 9 African blood donors enrolled in 2 national biobanks. We illustrated the usefulness of the platform by analyzing antibody data from patients sensitized with multiple RBC alloantibodies. Genotyping results demonstrated concordance of 99.91%, 99.97%, and 99.03% with RBC, HPA, and HLA clinically validated typing results in 89 371, 3016, and 9289 comparisons, respectively. Genotyping increased the total number of antigen typing results available from 110 980 to >1 200 000. Dense donor typing allowed identification of 2 to 6 times more compatible donors to serve 3146 patients with multiple RBC alloantibodies, providing at least 1 match for 176 individuals for whom previously no blood could be found among the same donors. This genotyping technology is already being used to type thousands of donors taking part in national genotyping studies. Extraction of dense antigen-typing data from these cohorts provides blood supply organizations with the opportunity to implement a policy of genomics-based precision matching of blood.

PIGO deficiency: palmoplantar keratoderma and novel mutations.

BACKGROUND: Several genetic defects have been identified in the glycosylphosphatidylinositol (GPI) anchor synthesis, including mutations in PIGO encoding phosphatidylinositol glycan anchor biosynthesis class O protein. These defects constitute a subgroup of the congenital disorders of glycosylation (CDG). Seven patients from five families have been reported carrying variants in PIGO that cause an autosomal recessive syndrome characterised by dysmorphism, psychomotor disability, epilepsy and hyperphosphatasemia. METHODS: Whole exome sequencing was performed in a boy with dysmorphism, psychomotor disability, epilepsy, palmoplantar keratoderma, hyperphosphatasemia and platelet dysfunction without a clinical bleeding phenotype. RESULTS: Two novel variants in PIGO were detected. The missense variant encoding p. His871Pro was inherited from the boy's father while the frameshift variant encoding p. Arg604ProfsTer40 was maternally inherited. CONCLUSION: A boy with two novel PIGO variants is reported. The skin phenotype and platelet dysfunction in this patient have not been described in previously reported patients with PIGO deficiency but it is of course uncertain whether these are caused by this disorder. The literature on PIGO deficiency is reviewed.