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.

100,000 Genomes Pilot on Rare-Disease Diagnosis in Health Care - Preliminary Report.

BACKGROUND: The U.K. 100,000 Genomes Project is in the process of investigating the role of genome sequencing in patients with undiagnosed rare diseases after usual care and the alignment of this research with health care implementation in the U.K. National Health Service. Other parts of this project focus on patients with cancer and infection. METHODS: We conducted a pilot study involving 4660 participants from 2183 families, among whom 161 disorders covering a broad spectrum of rare diseases were present. We collected data on clinical features with the use of Human Phenotype Ontology terms, undertook genome sequencing, applied automated variant prioritization on the basis of applied virtual gene panels and phenotypes, and identified novel pathogenic variants through research analysis. RESULTS: Diagnostic yields varied among family structures and were highest in family trios (both parents and a proband) and families with larger pedigrees. Diagnostic yields were much higher for disorders likely to have a monogenic cause (35%) than for disorders likely to have a complex cause (11%). Diagnostic yields for intellectual disability, hearing disorders, and vision disorders ranged from 40 to 55%. We made genetic diagnoses in 25% of the probands. A total of 14% of the diagnoses were made by means of the combination of research and automated approaches, which was critical for cases in which we found etiologic noncoding, structural, and mitochondrial genome variants and coding variants poorly covered by exome sequencing. Cohortwide burden testing across 57,000 genomes enabled the discovery of three new disease genes and 19 new associations. Of the genetic diagnoses that we made, 25% had immediate ramifications for clinical decision making for the patients or their relatives. CONCLUSIONS: Our pilot study of genome sequencing in a national health care system showed an increase in diagnostic yield across a range of rare diseases. (Funded by the National Institute for Health Research and others.).

First Genotype-Phenotype Study in TBX4 Syndrome: Gain-of-Function Mutations Causative for Lung Disease.

Rationale: Despite the increased recognition of TBX4 (T-BOX transcription factor 4)-associated pulmonary arterial hypertension (PAH), genotype-phenotype associations are lacking and may provide important insights. Objectives: To compile and functionally characterize all TBX4 variants reported to date and undertake a comprehensive genotype-phenotype analysis. Methods: We assembled a multicenter cohort of 137 patients harboring monoallelic TBX4 variants and assessed the pathogenicity of missense variation (n = 42) using a novel luciferase reporter assay containing T-BOX binding motifs. We sought genotype-phenotype correlations and undertook a comparative analysis with patients with PAH with BMPR2 (Bone Morphogenetic Protein Receptor type 2) causal variants (n = 162) or no identified variants in PAH-associated genes (n = 741) genotyped via the National Institute for Health Research BioResource-Rare Diseases. Measurements and Main Results: Functional assessment of TBX4 missense variants led to the novel finding of gain-of-function effects associated with older age at diagnosis of lung disease compared with loss-of-function effects (P = 0.038). Variants located in the T-BOX and nuclear localization domains were associated with earlier presentation (P = 0.005) and increased incidence of interstitial lung disease (P = 0.003). Event-free survival (death or transplantation) was shorter in the T-BOX group (P = 0.022), although age had a significant effect in the hazard model (P = 0.0461). Carriers of TBX4 variants were diagnosed at a younger age (P < 0.001) and had worse baseline lung function (FEV1, FVC) (P = 0.009) than the BMPR2 and no identified causal variant groups. Conclusions: We demonstrated that TBX4 syndrome is not strictly the result of haploinsufficiency but can also be caused by gain of function. The pleiotropic effects of TBX4 in lung disease may be in part explained by the differential effect of pathogenic mutations located in critical protein domains.

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.

Exploring the relevance of NUP93 variants in steroid-resistant nephrotic syndrome using next generation sequencing and a fly kidney model.

BACKGROUND: Variants in genes encoding nuclear pore complex (NPC) proteins are a newly identified cause of paediatric steroid-resistant nephrotic syndrome (SRNS). Recent reports describing NUP93 variants suggest these could be a significant cause of paediatric onset SRNS. We report NUP93 cases in the UK and demonstrate in vivo functional effects of Nup93 depletion in a fly (Drosophila melanogaster) nephrocyte model. METHODS: Three hundred thirty-seven paediatric SRNS patients from the National cohort of patients with Nephrotic Syndrome (NephroS) were whole exome and/or whole genome sequenced. Patients were screened for over 70 genes known to be associated with Nephrotic Syndrome (NS). D. melanogaster Nup93 knockdown was achieved by RNA interference using nephrocyte-restricted drivers. RESULTS: Six novel homozygous and compound heterozygous NUP93 variants were detected in 3 sporadic and 2 familial paediatric onset SRNS characterised histologically by focal segmental glomerulosclerosis (FSGS) and progressing to kidney failure by 12 months from clinical diagnosis. Silencing of the two orthologs of human NUP93 expressed in D. melanogaster, Nup93-1, and Nup93-2 resulted in significant signal reduction of up to 82% in adult pericardial nephrocytes with concomitant disruption of NPC protein expression. Additionally, nephrocyte morphology was highly abnormal in Nup93-1 and Nup93-2 silenced flies surviving to adulthood. CONCLUSION: We expand the spectrum of NUP93 variants detected in paediatric onset SRNS and demonstrate its incidence within a national cohort. Silencing of either D. melanogaster Nup93 ortholog caused a severe nephrocyte phenotype, signaling an important role for the nucleoporin complex in podocyte biology. A higher resolution version of the Graphical abstract is available as Supplementary information.

Comprehensive Cancer-Predisposition Gene Testing in an Adult Multiple Primary Tumor Series Shows a Broad Range of Deleterious Variants and Atypical Tumor Phenotypes.

Multiple primary tumors (MPTs) affect a substantial proportion of cancer survivors and can result from various causes, including inherited predisposition. Currently, germline genetic testing of MPT-affected individuals for variants in cancer-predisposition genes (CPGs) is mostly targeted by tumor type. We ascertained pre-assessed MPT individuals (with at least two primary tumors by age 60 years or at least three by 70 years) from genetics centers and performed whole-genome sequencing (WGS) on 460 individuals from 440 families. Despite previous negative genetic assessment and molecular investigations, pathogenic variants in moderate- and high-risk CPGs were detected in 67/440 (15.2%) probands. WGS detected variants that would not be (or were not) detected by targeted resequencing strategies, including low-frequency structural variants (6/440 [1.4%] probands). In most individuals with a germline variant assessed as pathogenic or likely pathogenic (P/LP), at least one of their tumor types was characteristic of variants in the relevant CPG. However, in 29 probands (42.2% of those with a P/LP variant), the tumor phenotype appeared discordant. The frequency of individuals with truncating or splice-site CPG variants and at least one discordant tumor type was significantly higher than in a control population (χ2 = 43.642; p ≤ 0.0001). 2/67 (3%) probands with P/LP variants had evidence of multiple inherited neoplasia allele syndrome (MINAS) with deleterious variants in two CPGs. Together with variant detection rates from a previous series of similarly ascertained MPT-affected individuals, the present results suggest that first-line comprehensive CPG analysis in an MPT cohort referred to clinical genetics services would detect a deleterious variant in about a third of individuals.

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.

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.

Investigating genotype-phenotype relationship of extreme neuropathic pain disorders in a UK national cohort.

The aims of our study were to use whole genome sequencing in a cross-sectional cohort of patients to identify new variants in genes implicated in neuropathic pain, to determine the prevalence of known pathogenic variants and to understand the relationship between pathogenic variants and clinical presentation. Patients with extreme neuropathic pain phenotypes (both sensory loss and gain) were recruited from secondary care clinics in the UK and underwent whole genome sequencing as part of the National Institute for Health and Care Research Bioresource Rare Diseases project. A multidisciplinary team assessed the pathogenicity of rare variants in genes previously known to cause neuropathic pain disorders and exploratory analysis of research candidate genes was completed. Association testing for genes carrying rare variants was completed using the gene-wise approach of the combined burden and variance-component test SKAT-O. Patch clamp analysis was performed on transfected HEK293T cells for research candidate variants of genes encoding ion channels. The results include the following: (i) Medically actionable variants were found in 12% of study participants (205 recruited), including known pathogenic variants: SCN9A(ENST00000409672.1): c.2544T>C, p.Ile848Thr that causes inherited erythromelalgia, and SPTLC1(ENST00000262554.2):c.340T>G, p.Cys133Tr variant that causes hereditary sensory neuropathy type-1. (ii) Clinically relevant variants were most common in voltage-gated sodium channels (Nav). (iii) SCN9A(ENST00000409672.1):c.554G>A, pArg185His variant was more common in non-freezing cold injury participants than controls and causes a gain of function of NaV1.7 after cooling (the environmental trigger for non-freezing cold injury). (iv) Rare variant association testing showed a significant difference in distribution for genes NGF, KIF1A, SCN8A, TRPM8, KIF1A, TRPA1 and the regulatory regions of genes SCN11A, FLVCR1, KIF1A and SCN9A between European participants with neuropathic pain and controls. (v) The TRPA1(ENST00000262209.4):c.515C>T, p.Ala172Val variant identified in participants with episodic somatic pain disorder demonstrated gain-of-channel function to agonist stimulation. Whole genome sequencing identified clinically relevant variants in over 10% of participants with extreme neuropathic pain phenotypes. The majority of these variants were found in ion channels. Combining genetic analysis with functional validation can lead to a better understanding as to how rare variants in ion channels lead to sensory neuron hyper-excitability, and how cold, as an environmental trigger, interacts with the gain-of-function NaV1.7 p.Arg185His variant. Our findings highlight the role of ion channel variants in the pathogenesis of extreme neuropathic pain disorders, likely mediated through changes in sensory neuron excitability and interaction with environmental triggers.

How common are single gene mutations as a cause for lacunar stroke? A targeted gene panel study.

OBJECTIVES: To determine the frequency of rare and pertinent disease-causing variants in small vessel disease (SVD)-associated genes (such as NOTCH3, HTRA1, COL4A1, COL4A2, FOXC1, TREX1, and GLA) in cerebral SVD, we performed targeted gene sequencing in 950 patients with younger-onset apparently sporadic SVD stroke using a targeted sequencing panel. METHODS: We designed a high-throughput sequencing panel to identify variants in 15 genes (7 known SVD genes, 8 SVD-related disorder genes). The panel was used to screen a population of 950 patients with younger-onset (≤70 years) MRI-confirmed SVD stroke, recruited from stroke centers across the United Kingdom. Variants were filtered according to their frequency in control databases, predicted effect, presence in curated variant lists, and combined annotation dependent depletion scores. Whole genome sequencing and genotyping were performed on a subset of patients to provide a direct comparison of techniques. The frequency of known disease-causing and pertinent variants of uncertain significance was calculated. RESULTS: We identified previously reported variants in 14 patients (8 cysteine-changing NOTCH3 variants in 11 patients, 2 HTRA1 variants in 2 patients, and 1 missense COL4A1 variant in 1 patient). In addition, we identified 29 variants of uncertain significance in 32 patients. CONCLUSION: Rare monogenic variants account for about 1.5% of younger onset lacunar stroke. Most are cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy variants, but the second most common gene affected is HTRA1. A high-throughput sequencing technology platform is an efficient, reliable method to screen for such mutations.