Cryptic intronic NBAS variant reveals the genetic basis of recurrent liver failure in a child

RNA sequencing driven diagnosis expands the phenotypic spectrum of NBAS deficiency

One in 10 individuals has a rare disease, with exome and genome sequencing yielding an overall diagnostic rate of approximately 30 %. RNA sequencing can augment genome analysis and improve diagnosis. We present a young woman with global developmental delay, poor growth, distinctive facial features, osteopenia, premature ovarian insufficiency, and ocular abnormalities who had non-diagnostic genome sequencing. RNAseq performed on her skin fibroblasts showed that NBAS gene expression was significantly reduced compared with controls. Manual inspection of the binary alignment map (BAM) files revealed compound heterozygous variants in NBAS: a rare deep intronic variant NM_015909.4:c.2423 + 403G > C which creates a hypomorphic pseudoexon not seen in control samples (gnomad allele frequency (AF) 0.000006572); and a rare premature termination codon (PTC) NM_015909.4:c.4753C > T; p.Arg1585Ter (gnomad AF 0.000006572). Both variants are predicted to cause nonsense mediated decay of transcripts, as the pseudoexon contains a PTC. Biallelic variants in NBAS are associated with two major phenotypes, i.e., infantile liver failure syndrome 2 (MIM # 616483) and short stature, optic nerve atrophy, and Pelger-Huet anomaly (MIM # 614800). Our patient, the first reported with one loss of function and one splice variant resulting in an out of frame transcript in NBAS, manifested a severe phenotype compared with previously reported individuals. This case demonstrates the utility of incorporating RNAseq to generate diagnostic candidates and expands the phenotypic spectrum of NBAS deficiency.

Data-driven insights to inform splice-altering variant assessment

Disease-causing genetic variants often disrupt mRNA splicing, an intricate process that is incompletely understood. Thus, accurate inference of which genetic variants will affect splicing and what their functional consequences will be is challenging, particularly for variants outside of the essential splice sites. Here, we describe a set of data-driven heuristics that inform the interpretation of human splice-altering variants (SAVs) based on the analysis of annotated exons, experimentally validated SAVs, and the currently understood principles of splicing biology. We defined requisite splicing criteria by examining around 202,000 canonical protein-coding exons and 19,000 experimentally validated splicing branchpoints. This analysis defined the sequence, spacing, and motif strength required for splicing, with 95.9% of the exons examined meeting these criteria. By considering over 12,000 experimentally validated variants from the SpliceVarDB, we defined a set of heuristics that inform the evaluation of putative SAVs. To ensure the applicability of each heuristic, only those supported by at least 10 experimentally validated variants were considered. This allowed us to establish a measure of spliceogenicity: the proportion of variants at a location (or motif site) that affected splicing in a given context. This study makes considerable advances toward bridging the gap between computational predictions and the biological process of splicing, offering an evidence-based approach to identifying SAVs and evaluating their impact. Our splicing heuristics enhance the current framework for genetic variant evaluation with a robust, detailed, and comprehensible analysis by adding valuable context over traditional binary prediction tools.

Distinct diagnostic trajectories in NBAS-associated acute liver failure highlights the need for timely functional studies

Variants of uncertain significance (VUS) are commonly found following genomic sequencing, particularly in ethnically diverse populations that are underrepresented in large population databases. Functional characterization of VUS may assist in variant reclassification, however these studies are not readily available and often rely on research funding and good will. We present four individuals from three families at different stages of their diagnostic trajectory with recurrent acute liver failure (RALF) and biallelic NBAS variants, confirmed by either trio analysis or cDNA studies. Functional characterization was undertaken, measuring NBAS and p31 levels by Western blotting, demonstrating reduced NBAS levels in two of three families, and reduced p31 levels in all three families. These results provided functional characterization of the molecular impact of a missense VUS, allowing reclassification of the variant and molecular confirmation of NBAS-associated RALF. Importantly, p31 was decreased in all individuals, including an individual with two missense variants where NBAS protein levels were preserved. These results highlight the importance of access to timely functional studies after identification of putative variants, and the importance of considering a range of assays to validate variants whose pathogenicity is uncertain. We suggest that funding models for genomic sequencing should consider incorporating capabilities for adjunct RNA, protein, biochemical, and other specialized tests to increase the diagnostic yield which will lead to improved medical care, increased equity, and access to molecular diagnoses for all patients.

Splice-disrupt genomic variants in prostate cancer

Background

Splice-disrupt genomic variants are one of the causes of cancer-causing errors in gene expression. Little is known about splice-disrupt genomic variants.

Methods and results

Here, pattern of splice-disrupt variants was investigated using 21,842,764 genomic variants in different types of prostate cancer. A particular attention was paid to genomic locations of splice-disrupt variants on target genes. HLA-A in prostate cancer, MSR1 in familial prostate cancer, and EGFR in both castration-resistant prostate cancer and metastatic castration-resistant had the highest allele frequencies of splice-disrupt variations. Some splice-disrupt variants, located on coding sequences of NCOR2, PTPRC, and CRP, were solely present in the advanced metastatic castration-resistant prostate cancer. High-risk splice-disrupt variants were identified based on computationally calculated Polymorphism Phenotyping (PolyPhen), Sorting Intolerant From Tolerant (SIFT), and Genomic Evolutionary Rate Profiling (GERP) + + scores as well as the recorded clinical significance in dbSNP database of NCBI. Functional annotation of damaging splice-disrupt variants highlighted important cancer-associated functions, including endocrine resistance, lipid metabolic process, steroid metabolic process, regulation of mitotic cell cycle, and regulation of metabolic process. This is the first study that profiles the splice-disrupt genomic variants and their target genes in prostate cancer. Literature mining based variant analysis highlighted the importance of rs1800716 variant, located on the CYP2D6 gene, involved in a range of important functions, such as RNA spicing, drug interaction, death, and urotoxicity.

Conclusions

This is the first study that profiles the splice-disrupt genomic variants and their target genes in different types of prostate cancer. Unravelling alternative splicing opens a new avenue towards the establishment of new diagnostic and prognostic markers for prostate cancer progression and metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-022-07257-9.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Application of Genome Sequencing from Blood to Diagnose Mitochondrial Diseases

Mitochondrial diseases can be caused by pathogenic variants in nuclear or mitochondrial DNA-encoded genes that often lead to multisystemic symptoms and can have any mode of inheritance. Using a single test, Genome Sequencing (GS) can effectively identify variants in both genomes, but it has not yet been universally used as a first-line approach to diagnosing mitochondrial diseases due to related costs and challenges in data analysis. In this article, we report three patients with mitochondrial disease molecularly diagnosed through GS performed on DNA extracted from blood to demonstrate different diagnostic advantages of this technology, including the detection of a low-level heteroplasmic pathogenic variant, an intragenic nuclear DNA deletion, and a large mtDNA deletion. Current technical improvements and cost reductions are likely to lead to an expanded routine diagnostic usage of GS and of the complementary “Omic” technologies in mitochondrial diseases.

Genome-first approach for the characterization of a complex phenotype with combined NBAS and CUL4B deficiency

Biallelic variants in neuroblastoma-amplified sequence (NBAS) cause an extremely broad spectrum of phenotypes. Clinical features range from isolated recurrent episodes of liver failure to multisystemic syndrome including short stature, skeletal osteopenia and dysplasia, optic atrophy, and a variable immunological, cutaneous, muscular, and neurological abnormalities. Hemizygous variants in CUL4B cause syndromic X-linked intellectual disability characterized by limitations in intellectual functions, developmental delays in gait, cognitive, and speech functioning, and other features including short stature, dysmorphism, and cerebral malformations. In this study, we report on a 4.5-month-old preterm infant with a complex phenotype mainly characterized by placental-related severe intrauterine growth restriction, post-natal growth failure with spontaneous bone fractures, which led to a suspicion of osteogenesis imperfecta, and lethal bronchopulmonary dysplasia with pulmonary hypertension. Whole exome sequencing identified compound heterozygosity for a known frameshift and a novel missense variant in NBAS and hemizygosity for a known CUL4B nonsense mutation. In vitro functional studies on the novel NBAS missense substitution demonstrated altered Golgi-to-endoplasmic reticulum retrograde vesicular trafficking and reduced collagen secretion, likely explaining part of the patient's phenotype. We also provided a comprehensive overview of the phenotypic features of NBAS and CUL4B deficiency, thus updating the recently emerging NBAS genotype-phenotype correlations. Our findings highlight the power of a genome-first approach for an early diagnosis of complex phenotypes.

NBAS disease: 14 new patients, a recurrent mutation, and genotype-phenotype correlation among 24 Chinese patients

AIM

Neuroblastoma amplified sequence (NBAS)-associated disease has a wide phenotypic spectrum, including infantile liver failure syndrome type 2 (ILFS2, OMIM #616483), short stature with optic nerve atrophy and Pelger-Huët anomaly (SOPH) syndrome (OMIM #614800), and a combined phenotype overlapping ILFS2 and SOPH syndrome. The mutation spectra of NBAS and its genotype-phenotype correlation among Chinese were not clear.

METHODS

Clinical and genetic data were retrospectively collected from the medical charts of patients with biallelic NBAS mutations, as well as from Chinese patients in previously published reports.

RESULTS

Fourteen new patients were identified, including 10 novel mutations: c.648-1G>A, c.2563_c.2577+5del/p.His855_Gln859del, c.3115C>T/p.Gln1039Ter, c.3284G>A/p.Trp1095Ter, c.2570C>T/p.Ala857Val, c.6859G>T/p.Asp2287Tyr, c.1028G>A/p.Ser343Asn, c.1177_1182delinsAGATAGA/p.Val393ArgfsTer2, c.3432_3435dupCAGT/p.Ala1146GlnfsTer14, and c.680_690dupACTGTTTCAGC/p.Phe231ThrfsTer35. All 14 patients presented as fever-triggered liver injury, including nine patients that satisfied the criteria of acute liver failure (ALF) in whom c.3596G>A/p.Cys1199Tyr occurred five times. Nine patients had extrahepatic manifestations including short stature, skeletal abnormalities, intellectual disability, ophthalmic abnormalities, low levels of serum immunoglobulins, facial dysmorphism, and cardiac abnormalities. Ten other Chinese patients were collected through a review of published works. Genotype-phenotype analysis in 24 Chinese patients revealed that the percentage of ALF patients with variants in the Sec39 domain was significantly higher than that in the C-terminal (100% vs. 12.5%, P = 0.000), and the percentage of multi-organ/system involvement in patients with variants in the Sec39 domain was significantly lower than that in the C-terminal (40% vs. 100%, P = 0.0128).

CONCLUSIONS

We reported 14 new patients, 10 novel mutations, and a unique recurrent mutation. Correlation analysis indicated that the domain of missense and non-frameshift insertion/deletion mutations in NBAS protein is related to phenotype among Chinese patients.

Severe SOPH syndrome due to a novel NBAS mutation in a 27-year-old woman-Review of this pleiotropic, autosomal recessive disorder: Mystery solved after two decades

Autosomal recessive SOPH syndrome was first described in the Yakuts population of Asia by Maksimova et al. in 2010. It arises from biallelic pathogenic variants in the NBAS gene and is characterized by severe postnatal growth retardation, senile facial appearance, small hands and feet, optic atrophy with loss of visual acuity and color vision, and normal intelligence (OMIM #614800). The presence of Pelger-Hüet anomaly in this disorder led to its name as an acronym for Short stature, Optic nerve atrophy, and Pelger-Hüet anomaly. Recent publications have further contributed to the characterization of this syndrome through additional phenotype-genotype correlations. We review the clinical features described in these publications and report on a 27-year-old woman with dwarfism with osteolysis and multiple skeletal problems, minor anomalies, immunodeficiency, diabetes mellitus, and multiple secondary medical problems. Her condition was considered an unknown autosomal recessive disorder for many years until exome sequencing provided the diagnosis by revealing a founder disease-causing variant that was compound heterozygous with a novel pathogenic variant in NBAS. Based on the major clinical features of this individual and others reported earlier, a revision of the acronym is warranted to facilitate clinical recognition.