Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Advanced variant classification framework reduces the false positive rate of predicted loss-of-function variants in population sequencing data.

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Transcript expression-aware annotation improves rare variant interpretation.

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)1, we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

The mutational constraint spectrum quantified from variation in 141,456 humans.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Treatment of Multisystem Inflammatory Syndrome in Children.

BACKGROUND: Evidence is urgently needed to support treatment decisions for children with multisystem inflammatory syndrome (MIS-C) associated with severe acute respiratory syndrome coronavirus 2. METHODS: We performed an international observational cohort study of clinical and outcome data regarding suspected MIS-C that had been uploaded by physicians onto a Web-based database. We used inverse-probability weighting and generalized linear models to evaluate intravenous immune globulin (IVIG) as a reference, as compared with IVIG plus glucocorticoids and glucocorticoids alone. There were two primary outcomes: the first was a composite of inotropic support or mechanical ventilation by day 2 or later or death; the second was a reduction in disease severity on an ordinal scale by day 2. Secondary outcomes included treatment escalation and the time until a reduction in organ failure and inflammation. RESULTS: Data were available regarding the course of treatment for 614 children from 32 countries from June 2020 through February 2021; 490 met the World Health Organization criteria for MIS-C. Of the 614 children with suspected MIS-C, 246 received primary treatment with IVIG alone, 208 with IVIG plus glucocorticoids, and 99 with glucocorticoids alone; 22 children received other treatment combinations, including biologic agents, and 39 received no immunomodulatory therapy. Receipt of inotropic or ventilatory support or death occurred in 56 patients who received IVIG plus glucocorticoids (adjusted odds ratio for the comparison with IVIG alone, 0.77; 95% confidence interval [CI], 0.33 to 1.82) and in 17 patients who received glucocorticoids alone (adjusted odds ratio, 0.54; 95% CI, 0.22 to 1.33). The adjusted odds ratios for a reduction in disease severity were similar in the two groups, as compared with IVIG alone (0.90 for IVIG plus glucocorticoids and 0.93 for glucocorticoids alone). The time until a reduction in disease severity was similar in the three groups. CONCLUSIONS: We found no evidence that recovery from MIS-C differed after primary treatment with IVIG alone, IVIG plus glucocorticoids, or glucocorticoids alone, although significant differences may emerge as more data accrue. (Funded by the European Union's Horizon 2020 Program and others; BATS ISRCTN number, ISRCTN69546370.).

A gene pathogenicity tool "GenePy" identifies missed biallelic diagnoses in the 100,000 Genomes Project.

PURPOSE: The 100,000 Genomes Project diagnosed a quarter of affected participants, but 26% of diagnoses were not on the applied gene panel(s); with many being de novo variants. Assessing biallelic variants without a gene panel is more challenging. METHODS: We sought to identify missed biallelic diagnoses using GenePy, which incorporates allele frequency, zygosity, and a user-defined deleterious metric, generating an aggregate GenePy score per gene, per participant. We calculated GenePy scores for 2862 recessive disease genes in 78,216 100,000 Genomes Project participants. For each gene, we ranked participant GenePy scores and scrutinized affected participants without a diagnosis, whose scores ranked among the top 5 for each gene. In cases which participant phenotypes overlapped with the disease gene of interest, we extracted rare variants and applied phase, ClinVar, and ACMG classification. RESULTS: 3184 affected individuals without a molecular diagnosis had a top-5-ranked GenePy score and 682 of 3184 (21%) had phenotypes overlapping with a top-ranking gene. In 122 of 669 (18%) phenotype-matched cases (excluding 13 withdrawn participants), we identified a putative missed diagnosis (2.2% of all undiagnosed participants). A further 334 of 669 (50%) cases have a possible missed diagnosis but require functional validation. CONCLUSION: Applying GenePy at scale has identified 456 potential diagnoses, demonstrating the value of novel diagnostic strategies.

Targeting de novo loss-of-function variants in constrained disease genes improves diagnostic rates in the 100,000 Genomes Project.

BACKGROUND: Genome sequencing was first offered clinically in the UK through the 100,000 Genomes Project (100KGP). Analysis was restricted to predefined gene panels associated with the patient's phenotype. However, panels rely on clearly characterised phenotypes and risk missing diagnoses outside of the panel(s) applied. We propose a complementary method to rapidly identify pathogenic variants, including those missed by 100KGP methods. METHODS: The Loss-of-function Observed/Expected Upper-bound Fraction (LOEUF) score quantifies gene constraint, with low scores correlated with haploinsufficiency. We applied DeNovoLOEUF, a filtering strategy to sequencing data from 13,949 rare disease trios in the 100KGP, by filtering for rare, de novo, loss-of-function variants in disease genes with a LOEUF score < 0.2. We compared our findings with the corresponding patient's diagnostic reports. RESULTS: 324/332 (98%) of the variants identified using DeNovoLOEUF were diagnostic or partially diagnostic (whereby the variant was responsible for some of the phenotype). We identified 39 diagnoses that were "missed" by 100KGP standard analyses, which are now being returned to patients. CONCLUSION: We have demonstrated a highly specific and rapid method with a 98% positive predictive value that has good concordance with standard analysis, low false-positive rate, and can identify additional diagnoses. Globally, as more patients are being offered genome sequencing, we anticipate that DeNovoLOEUF will rapidly identify new diagnoses and facilitate iterative analyses when new disease genes are discovered.

Delineation of a Human Mendelian Disorder of the DNA Demethylation Machinery: TET3 Deficiency.

Germline pathogenic variants in chromatin-modifying enzymes are a common cause of pediatric developmental disorders. These enzymes catalyze reactions that regulate epigenetic inheritance via histone post-translational modifications and DNA methylation. Cytosine methylation (5-methylcytosine [5mC]) of DNA is the quintessential epigenetic mark, yet no human Mendelian disorder of DNA demethylation has yet been delineated. Here, we describe in detail a Mendelian disorder caused by the disruption of DNA demethylation. TET3 is a methylcytosine dioxygenase that initiates DNA demethylation during early zygote formation, embryogenesis, and neuronal differentiation and is intolerant to haploinsufficiency in mice and humans. We identify and characterize 11 cases of human TET3 deficiency in eight families with the common phenotypic features of intellectual disability and/or global developmental delay; hypotonia; autistic traits; movement disorders; growth abnormalities; and facial dysmorphism. Mono-allelic frameshift and nonsense variants in TET3 occur throughout the coding region. Mono-allelic and bi-allelic missense variants localize to conserved residues; all but one such variant occur within the catalytic domain, and most display hypomorphic function in an assay of catalytic activity. TET3 deficiency and other Mendelian disorders of the epigenetic machinery show substantial phenotypic overlap, including features of intellectual disability and abnormal growth, underscoring shared disease mechanisms.

A novel variant in GATM causes idiopathic renal Fanconi syndrome and predicts progression to end-stage kidney disease.

Renal Fanconi syndrome (RFS) is a generalised disorder of the proximal convoluted tubule. Many genes have been associated with RFS including those that cause systemic disorders such as cystinosis, as well as isolated RFS. We discuss the case of a 10-year-old female who presented with leg pain and raised creatinine on a screening blood test. Her mother has RFS and required a kidney transplant in her thirties. Further investigations confirmed RFS in the daughter. Exome sequencing was performed on the affected mother, child, and unaffected father. We identified a novel variant in GATM; c.965G>C p.(Arg322Pro) segregating dominantly in the mother and daughter. We validated our finding with molecular dynamics simulations and demonstrated a dynamic signature that differentiates our variant and two previously identified pathogenic variants in GATM from wildtype. Genetic testing has uncovered a novel pathogenic variant that predicts progression to end stage kidney failure and has important implications for family planning and cascade testing. We recommend that GATM is screened for in children presenting with RFS, in addition to adults, particularly with kidney failure, who may have had previous negative gene testing.

The Pediatric Crohn Disease Morbidity Index (PCD-MI): Development of a Tool to Assess Long-Term Disease Burden Using a Data-Driven Approach.

BACKGROUND/OBJECTIVE: Heterogeneity and chronicity of Crohn disease (CD) make prediction of outcomes difficult. To date, no longitudinal measure can quantify burden over a patient's disease course, preventing assessment and integration into predictive modeling. Here, we aimed to demonstrate the feasibility of constructing a data driven, longitudinal disease burden score. METHODS: Literature was reviewed for tools used in assessment of CD activity. Themes were identified to construct a pediatric CD morbidity index (PCD-MI). Scores were assigned to variables. Data were extracted automatically from the electronic patient records at Southampton Children's Hospital, diagnosed from 2012 to 2019 (inclusive). PCD-MI scores were calculated, adjusted for duration of follow up and assessed for variation (ANOVA) and distribution (Kolmogorov-Smirnov). RESULTS: Nineteen clinical/biological features across five themes were included in the PCD-MI including blood/fecal/radiological/endoscopic results, medication usage, surgery, growth parameters, and extraintestinal manifestations. Maximal score was 100 after accounting for follow-up duration. PCD-MI was assessed in 66 patients, mean age 12.5 years. Following quality filtering, 9528 blood/fecal test results and 1309 growth measures were included. Mean PCD-MI score was 14.95 (range 2.2-32.5); data were normally distributed ( P = 0.2) with 25% of patients having a PCD-MI < 10. There was no difference in the mean PCD-MI when split by year of diagnosis, F -statistic 1.625, P = 0.147. CONCLUSIONS: PCD-MI is a calculatable measure for a cohort of patients diagnosed over an 8-year period, integrating a wide-range of data with potential to determine high or low disease burden. Future iterations of the PCD-MI require refinement of included features, optimized scores, and validation on external cohorts.

De novo putative loss-of-function variants in TAF4 are associated with a neuro-developmental disorder.

TATA-binding protein associated factor 4 (TAF4) is a subunit of the Transcription Factor IID (TFIID) complex, a central player in transcription initiation. Other members of this multimeric complex have been implicated previously as monogenic disease genes in human developmental disorders. TAF4 has not been described to date as a monogenic disease gene. We here present a cohort of eight individuals, each carrying de novo putative loss-of-function (pLoF) variants in TAF4 and expressing phenotypes consistent with a neuro-developmental disorder (NDD). Common features include intellectual disability, abnormal behavior, and facial dysmorphisms. We propose TAF4 as a novel dominant disease gene for NDD, and coin this novel disorder "TAF4-related NDD" (T4NDD). We place T4NDD in the context of other disorders related to TFIID subunits, revealing shared features of T4NDD with other TAF-opathies.

A gene-to-patient approach uplifts novel disease gene discovery and identifies 18 putative novel disease genes.

PURPOSE: Exome and genome sequencing have drastically accelerated novel disease gene discoveries. However, discovery is still hindered by myriad variants of uncertain significance found in genes of undetermined biological function. This necessitates intensive functional experiments on genes of equal predicted causality, leading to a major bottleneck. METHODS: We apply the loss-of-function observed/expected upper-bound fraction metric of intolerance to gene inactivation to curate a list of predicted haploinsufficient disease genes. Using data from the 100,000 Genomes Project, we adopt a gene-to-patient approach that matches de novo loss-of-function variants in constrained genes to patients with rare disease. Through large-scale aggregation of data, we reduce excess analytical noise currently hindering novel discoveries. RESULTS: Results from 13,949 trios revealed 643 rare, de novo predicted loss-of-function events filtered from 1044 loss-of-function observed/expected upper-bound fraction-constrained genes. A total of 168 variants occurred within 126 genes without a known disease-gene relationship. Of these, 27 genes had >1 kindred affected, and for 18 of these genes, multiple kindreds had overlapping phenotypes. Two years after initial analysis, 11 of 18 (61%) of these genes have been independently published as novel disease gene discoveries. CONCLUSION: Using large cohorts and adopting gene-based approaches can rapidly and objectively accelerate dominantly inherited novel gene discovery by targeting the most appropriate genes for functional validation.

Rare pathogenic variants in WNK3 cause X-linked intellectual disability.

PURPOSE: WNK3 kinase (PRKWNK3) has been implicated in the development and function of the brain via its regulation of the cation-chloride cotransporters, but the role of WNK3 in human development is unknown. METHOD: We ascertained exome or genome sequences of individuals with rare familial or sporadic forms of intellectual disability (ID). RESULTS: We identified a total of 6 different maternally-inherited, hemizygous, 3 loss-of-function or 3 pathogenic missense variants (p.Pro204Arg, p.Leu300Ser, p.Glu607Val) in WNK3 in 14 male individuals from 6 unrelated families. Affected individuals had ID with variable presence of epilepsy and structural brain defects. WNK3 variants cosegregated with the disease in 3 different families with multiple affected individuals. This included 1 large family previously diagnosed with X-linked Prieto syndrome. WNK3 pathogenic missense variants localize to the catalytic domain and impede the inhibitory phosphorylation of the neuronal-specific chloride cotransporter KCC2 at threonine 1007, a site critically regulated during the development of synaptic inhibition. CONCLUSION: Pathogenic WNK3 variants cause a rare form of human X-linked ID with variable epilepsy and structural brain abnormalities and implicate impaired phospho-regulation of KCC2 as a pathogenic mechanism.

Blood RNA analysis can increase clinical diagnostic rate and resolve variants of uncertain significance.

PURPOSE: Diagnosis of genetic disorders is hampered by large numbers of variants of uncertain significance (VUSs) identified through next-generation sequencing. Many such variants may disrupt normal RNA splicing. We examined effects on splicing of a large cohort of clinically identified variants and compared performance of bioinformatic splicing prediction tools commonly used in diagnostic laboratories. METHODS: Two hundred fifty-seven variants (coding and noncoding) were referred for analysis across three laboratories. Blood RNA samples underwent targeted reverse transcription polymerase chain reaction (RT-PCR) analysis with Sanger sequencing of PCR products and agarose gel electrophoresis. Seventeen samples also underwent transcriptome-wide RNA sequencing with targeted splicing analysis based on Sashimi plot visualization. Bioinformatic splicing predictions were obtained using Alamut, HSF 3.1, and SpliceAI software. RESULTS: Eighty-five variants (33%) were associated with abnormal splicing. The most frequent abnormality was upstream exon skipping (39/85 variants), which was most often associated with splice donor region variants. SpliceAI had greatest accuracy in predicting splicing abnormalities (0.91) and outperformed other tools in sensitivity and specificity. CONCLUSION: Splicing analysis of blood RNA identifies diagnostically important splicing abnormalities and clarifies functional effects of a significant proportion of VUSs. Bioinformatic predictions are improving but still make significant errors. RNA analysis should therefore be routinely considered in genetic disease diagnostics.

Correction: Blood RNA analysis can increase clinical diagnostic rate and resolve variants of uncertain significance.

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

Inactivation of AMMECR1 is associated with growth, bone, and heart alterations.

We report five individuals with loss-of-function of the X-linked AMMECR1: a girl with a balanced X-autosome translocation and inactivation of the normal X-chromosome; two boys with maternally inherited and de novo nonsense variants; and two half-brothers with maternally inherited microdeletion variants. They present with short stature, cardiac and skeletal abnormalities, and hearing loss. Variants of unknown significance in AMMECR1 in four male patients from two families with partially overlapping phenotypes were previously reported. AMMECR1 is coexpressed with genes implicated in cell cycle regulation, five of which were previously associated with growth and bone alterations. Our knockdown of the zebrafish orthologous gene resulted in phenotypes reminiscent of patients' features. The increased transcript and encoded protein levels of AMMECR1L, an AMMECR1 paralog, in the t(X;9) patient's cells indicate a possible partial compensatory mechanism. AMMECR1 and AMMECR1L proteins dimerize and localize to the nucleus as suggested by their nucleic acid-binding RAGNYA folds. Our results suggest that AMMECR1 is potentially involved in cell cycle control and linked to a new syndrome with growth, bone, heart, and kidney alterations with or without elliptocytosis.

Thrombotic microangiopathy following haematopoietic stem cell transplant.

Thrombotic microangiopathy is a potentially lethal complication of haematopoietic stem cell (bone marrow) transplantation. The pathophysiology is incompletely understood, although endothelial damage appears to be central. Platelet activation, neutrophil extracellular traps and complement activation appear to play key roles. Diagnosis may be difficult and universally accepted diagnostic criteria are not available. Treatment remains controversial. In some cases, withdrawal of calcineurin inhibitors is adequate. Rituximab and defibrotide also appear to have been used successfully. In severe cases, complement inhibitors such as eculizumab may play a valuable role. Further research is required to define the pathophysiology and determine both robust diagnostic criteria and the optimal treatment.