GeneSplicer: a new computational method for splice site prediction

MTSplice predicts effects of genetic variants on tissue-specific splicing

We develop the free and open-source model Multi-tissue Splicing (MTSplice) to predict the effects of genetic variants on splicing of cassette exons in 56 human tissues. MTSplice combines MMSplice, which models constitutive regulatory sequences, with a new neural network that models tissue-specific regulatory sequences. MTSplice outperforms MMSplice on predicting tissue-specific variations associated with genetic variants in most tissues of the GTEx dataset, with largest improvements on brain tissues. Furthermore, MTSplice predicts that autism-associated de novo mutations are enriched for variants affecting splicing specifically in the brain. We foresee that MTSplice will aid interpreting variants associated with tissue-specific disorders.

The Impact of Modern Technologies on Molecular Diagnostic Success Rates, with a Focus on Inherited Retinal Dystrophy and Hearing Loss

The identification of pathogenic variants in monogenic diseases has been of interest to researchers and clinicians for several decades. However, for inherited diseases with extremely high genetic heterogeneity, such as hearing loss and retinal dystrophies, establishing a molecular diagnosis requires an enormous effort. In this review, we use these two genetic conditions as examples to describe the initial molecular genetic identification approaches, as performed since the early 90s, and subsequent improvements and refinements introduced over the years. Next, the history of DNA sequencing from conventional Sanger sequencing to high-throughput massive parallel sequencing, a.k.a. next-generation sequencing, is outlined, including their advantages and limitations and their impact on identifying the remaining genetic defects. Moreover, the development of recent technologies, also coined "third-generation" sequencing, is reviewed, which holds the promise to overcome these limitations. Furthermore, we outline the importance and complexity of variant interpretation in clinical diagnostic settings concerning the massive number of different variants identified by these methods. Finally, we briefly mention the development of novel approaches such as optical mapping and multiomics, which can help to further identify genetic defects in the near future.

Transcript analysis for variant classification resolution in a child with primary ciliary dyskinesia

Transcriptional analysis can be utilized to reconcile variants of uncertain significance, particularly those predicted to impact splicing. Laboratory analysis of the predicted mRNA transcript may allow inference of the in vivo impact of the variant and aid prediction of its clinical significance. We present a patient with classical features of primary ciliary dyskinesia (PCD) who was identified to have compound heterozygous variants in the DNAH11 gene (c.10691 + 2T > C, c.13523_13543dup21) via trio whole-exome sequencing in 2013. These variants were originally classified as Mutation and Likely Mutation. However, these variants were downgraded to variants of uncertain significance (VUSs) during reanalysis in 2016 because of uncertainty that they caused a loss of function of the gene. c.10691 + 2T > C is predicted to abrogate the canonical splice site and lead to the skipping of exon 65, but the adjoining of exon 64 and exon 66 in the DNAH11 transcript preserves the reading frame of the resultant protein. c.13523_13543dup21 is located in the last exon of the DNAH11 coding sequence, upstream of the canonical stop codon, which suggests a reduced likelihood to trigger nonsense-mediated decay (NMD). Transcriptional analysis was performed to characterize the impact of the variants, resulting in reclassification of c.10691 + 2T > C to Likely Pathogenic by providing evidence that it results in a deleterious effect and subsequent downstream reclassification of c.13523_13543dup21 to Likely Pathogenic as well. Our case illustrates the potential impact of transcriptional analysis on variant resolution, supporting its usage on variants that exert an unpredictable effect on splicing.

Next-Generation Molecular Investigations in Lysosomal Diseases: Clinical Integration of a Comprehensive Targeted Panel

Diagnosis of lysosomal disorders (LDs) may be hampered by their clinical heterogeneity, phenotypic overlap, and variable age at onset. Conventional biological diagnostic procedures are based on a series of sequential investigations and require multiple sampling. Early diagnosis may allow for timely treatment and prevent clinical complications. In order to improve LDs diagnosis, we developed a capture-based next generation sequencing (NGS) panel allowing the detection of single nucleotide variants (SNVs), small insertions and deletions, and copy number variants (CNVs) in 51 genes related to LDs. The design of the LD panel covered at least coding regions, promoter region, and flanking intronic sequences for 51 genes. The validation of this panel consisted in testing 21 well-characterized samples and evaluating analytical and diagnostic performance metrics. Bioinformatics pipelines have been validated for SNVs, indels and CNVs. The clinical output of this panel was tested in five novel cases. This capture-based NGS panel provides an average coverage depth of 474× which allows the detection of SNVs and CNVs in one comprehensive assay. All the targeted regions were covered above the minimum required depth of 30×. To illustrate the clinical utility, five novel cases have been sequenced using this panel and the identified variants have been confirmed using Sanger sequencing or quantitative multiplex PCR of short fluorescent fragments (QMPSF). The application of NGS as first-line approach to analyze suspected LD cases may speed up the identification of alterations in LD-associated genes. NGS approaches combined with bioinformatics analyses, are a useful and cost-effective tool for identifying the causative variations in LDs.

Intronic variant screening with targeted next-generation sequencing reveals first pseudoexon in LDLR in familial hypercholesterolemia

BACKGROUND AND AIMS

Familial hypercholesterolemia (FH) is caused by pathogenic variants in LDLR, APOB, or PCSK9 genes (designated FH+). However, a significant number of clinical FH patients do not carry these variants (designated FH-). Here, we investigated whether variants in intronic regions of LDLR attribute to FH by affecting pre-mRNA splicing.

METHODS

LDLR introns are partly covered in routine sequencing of clinical FH patients using next-generation sequencing. Deep intronic variants, >20 bp from intron-exon boundary, were considered of interest once (a) present in FH- patients (n = 909) with LDL-C >7 mmol/L (severe FH-) or after in silico analysis in patients with LDL-C >5 mmol/L (moderate FH-) and b) absent in FH + patients (control group). cDNA analysis and co-segregation analysis were performed to assess pathogenicity of the identified variants.

RESULTS

Three unique variants were present in the severe FH- group. One of these was the previously described likely pathogenic variant c.2140+103G>T. Three additional variants were selected based on in silico analyses in the moderate FH- group. One of these variants, c.2141-218G>A, was found to result in a pseudo-exon inclusion, producing a premature stop codon. This variant co-segregated with the hypercholesterolemic phenotype.

CONCLUSIONS

Through a screening approach, we identified a deep intronic variant causal for FH. This finding indicates that filtering intronic variants in FH- patients for the absence in FH + patients might enrich for true FH-causing variants and suggests that intronic regions of LDLR need to be considered for sequencing in FH- patients.

Targeted High-Throughput Sequencing Analysis Results of Osteogenesis Imperfecta Patients from Different Regions of Turkey

Objective: Osteogenesis imperfecta (OI) includes a group of disorders characterized by susceptibility to bone fractures with different severities. The increasing number of genes that may underlie the disorder, along with the broad phenotypic spectrum that overlaps with other skeletal diseases, provided a compelling case for the use of high-throughput sequencing (HTS) technology as an aid to OI diagnoses. The aim of this analysis was to present the data from our 5-year targeted HTS results, that includes the reporting of 9 novel and 24 known mutations, found in OI patients, from 5 different regions of Turkey. Materials and Methods: We performed a retrospective cross-sectional study, reporting the HTS results of 43 patients (23 female and 20 male; mean age: 9.5 years), directed to our center with a suspicion of OI between February 2015 and May 2020. Genetic analyses were also performed for 24 asymptomatic parents to aid the segregation analyses. We utilized an HTS panel targeting the coding regions of 57 genes associated with a reduction, increase, or abnormal development of bone mineralization. In addition, we sequenced the entire coding region of the IFITM5 gene through HTS. Results: Thirty-nine patients had at least one pathogenic/likely pathogenic variation (90.69%) in the COL1A1 (56.41%), COL1A2 (20.51%), FKBP10 (7.7%), P3H1 (5.13%), IFITM5 (5.13%), CTRAP (2.56%), or TMEM38B (2.56%) genes. Nine of the determined pathogenic/likely pathogenic variations were novel. The recurrent pathogenic mutations were c.1081C>T (p.Arg361Ter) (3/43), c.1405C>T (p.Arg469Ter) (2/43), and c.3749del (p.Gly1250AlafsTer81) in COL1A1 gene, along with c.-14C>T variation in the 5'UTR of the IFITM5 gene (2/43) and the c.890_897dup variation in the FKBP10 gene (2/43). Three out of 43 patients were carrying at least one additional variant of unknown significance, highlighting the importance of a multigene panel approach and segregation analyses. Conclusion: We suggest that a targeted HTS panel is a feasible tool for genetic diagnosis of OI in patients.

Phenotype Driven Analysis of Whole Genome Sequencing Identifies Deep Intronic Variants that Cause Retinal Dystrophies by Aberrant Exonization

Purpose

To demonstrate the effectiveness of combining retinal phenotyping and focused variant filtering from genome sequencing (GS) in identifying deep intronic disease causing variants in inherited retinal dystrophies.

Methods

Affected members from three pedigrees with classical enhanced S-cone syndrome (ESCS; Pedigree 1), congenital stationary night blindness (CSNB; Pedigree 2), and achromatopsia (ACHM; Pedigree 3), respectively, underwent detailed ophthalmologic evaluation, optical coherence tomography, and electroretinography. The probands underwent panel-based genetic testing followed by GS analysis. Minigene constructs (NR2E3, GPR179 and CNGB3) and patient-derived cDNA experiments (NR2E3 and GPR179) were performed to assess the functional effect of the deep intronic variants.

Results

The electrophysiological findings confirmed the clinical diagnosis of ESCS, CSNB, and ACHM in the respective pedigrees. Panel-based testing revealed heterozygous pathogenic variants in NR2E3 (NM_014249.3; c.119-2A>C; Pedigree 1) and CNGB3 (NM_019098.4; c.1148delC/p.Thr383Ilefs*13; Pedigree 3). The GS revealed heterozygous deep intronic variants in Pedigrees 1 (NR2E3; c.1100+1124G>A) and 3 (CNGB3; c.852+4751A>T), and a homozygous GPR179 variant in Pedigree 2 (NM_001004334.3; c.903+343G>A). The identified variants segregated with the phenotype in all pedigrees. All deep intronic variants were predicted to generate a splice acceptor gain causing aberrant exonization in NR2E3 [89 base pairs (bp)], GPR179 (197 bp), and CNGB3 (73 bp); splicing defects were validated through patient-derived cDNA experiments and/or minigene constructs and rescued by antisense oligonucleotide treatment.

Conclusions

Deep intronic mutations contribute to missing heritability in retinal dystrophies. Combining results from phenotype-directed gene panel testing, GS, and in silico splice prediction tools can help identify these difficult-to-detect pathogenic deep intronic variants.

Machine learning based CRISPR gRNA design for therapeutic exon skipping

Restoring gene function by the induced skipping of deleterious exons has been shown to be effective for treating genetic disorders. However, many of the clinically successful therapies for exon skipping are transient oligonucleotide-based treatments that require frequent dosing. CRISPR-Cas9 based genome editing that causes exon skipping is a promising therapeutic modality that may offer permanent alleviation of genetic disease. We show that machine learning can select Cas9 guide RNAs that disrupt splice acceptors and cause the skipping of targeted exons. We experimentally measured the exon skipping frequencies of a diverse genome-integrated library of 791 splice sequences targeted by 1,063 guide RNAs in mouse embryonic stem cells. We found that our method, SkipGuide, is able to identify effective guide RNAs with a precision of 0.68 (50% threshold predicted exon skipping frequency) and 0.93 (70% threshold predicted exon skipping frequency). We anticipate that SkipGuide will be useful for selecting guide RNA candidates for evaluation of CRISPR-Cas9-mediated exon skipping therapy.

Pharmacogenomics, How to Deal with Different Types of Variants in Next Generation Sequencing Data in the Personalized Medicine Area

Pharmacogenomics (PGx) is the knowledge of diverse drug responses and effects in people, based on their genomic profiles. Such information is considered as one of the main directions to reach personalized medicine in future clinical practices. Since the start of applying next generation sequencing (NGS) methods in drug related clinical investigations, many common medicines found their genetic data for the related metabolizing/shipping proteins in the human body. Yet, the employing of technology is accompanied by big obtained data, which most of them have no clear guidelines for consideration in routine treatment decisions for patients. This review article talks about different types of NGS derived PGx variants in clinical studies and try to display the current and newly developed approaches to deal with pharmacogenetic data with/without clear guidelines for considering in clinical settings.

Generalized Hailey-Hailey disease: Novel splice-site mutations of ATP2C1 gene in Chinese population and a literature review

BACKGROUND

Hailey-Hailey disease (HHD; OMIM: 169600) is an autosomal dominate genodermatosis, characterized by recurrent blisters and erosions clinically and remarkable acantholysis pathologically. The underlying pathogenic factor is the mutation of ATP2C1 gene (OMIM: 604384), which encodes secretory pathway Ca²⁺ /Mn²⁺ -ATPase (SPCA1). Skin folds are the predilection site of HHD. Atypical cases with a generalized pattern have rarely been reported, making it prone to misdiagnosis.

METHODS

In this study, we presented three Chinese pedigrees of Hailey-Hailey disease with generalized skin lesions. ATP2C1 mutations were screened by DNA sequencing and their transcripts were further confirmed by minigene assay. We also performed a literature review of previously published generalized HHD over past two decades together with our cases.

RESULTS

Three splice-site mutations were identified: c.2487+1G>A, c.2126+1G>A, and c.1891-2A>G, which resulted in an exon 25-truncated transcript, two exon 22-truncated transcripts, and two exon 21-truncated transcripts, respectively. The c.2487+1G>A and the c.1891-2A>G mutations are novel mutations which have not been reported before. No clustered mutations of ATP2C1 gene were found in generalized HHD patients in literature along with our novel mutations.

CONCLUSION

We found no hot spot mutations in ATP2C1 correlated with the generalized pattern of HHD. Our study expanded the spectrum of ATP2C1 mutations, which would be useful for disease diagnosis and genetic counseling.

Bi-Allelic Pathogenic Variations in MERTK Including Deletions Are Associated with an Early Onset Progressive Form of Retinitis Pigmentosa

Bi-allelic pathogenic variants in MERTK cause retinitis pigmentosa (RP). Since deletions of more than one exon have been reported repeatedly for MERTK, CNV (copy number variation) analysis of next-generation sequencing (NGS) data has proven important in molecular genetic diagnostics of MERTK. CNV analysis was performed on NGS data of 677 individuals with inherited retinal diseases (IRD) and confirmed by quantitative RT-PCR analysis. Clinical evaluation was based on retrospective records. Clinical re-examination included visual field examination, dark adaption, scotopic and photopic full-field electroretinograms (ffERG), multifocal ERG (mfERG) and optic coherence tomography (OCT). Fourteen variants were detected in MERTK in six individuals, three of which were deletions of more than one exon. Clinical examinations of five out of six individuals revealed a severe phenotype with early-onset generalized retinal dystrophy with night blindness and progressive visual field loss; however, one individual had a milder phenotype. Three individuals had hearing impairments. We show that deletions represent a substantial part of the causative variants in MERTK and emphasize that CNV analysis should be included in the molecular genetic diagnostics of IRDs.

Novel de novo ARCN1 intronic variant causes rhizomelic short stature with microretrognathia and developmental delay

The archain 1 (ARCN1) gene encodes the coatomer subunit delta protein and is a component of the COPI coatomer complex, which is involved in retrograde vesical trafficking from the Golgi complex to the endoplasmic reticulum. Variants in ARCN1 have recently been associated with rhizomelic short stature with microcephaly, microretrognathia, and developmental delay. Here we report a 3.5-yr-old boy with microcephaly, global developmental delay, and multiple congenital abnormalities and the ARCN1-related syndrome caused by a novel de novo intronic variant. Whole-exome sequencing of the proband and his parents was utilized to determine the genetic origin of the patient's disorder and identified a de novo variant, NM_001655.5:c.654-15A > G, in the ARCN1 gene. Follow-up functional characterization of mRNA from the patient demonstrated that this variant creates a splicing defect of the ARCN1 mRNA. ARCN1-related syndrome represents an emerging disorder of developmental delay, and this report represents the sixth described patient. Despite the few instances reported in literature, the phenotype is consistent between our patient and previously reported individuals.

Usher syndrome and Nebulin-associated myopathy in a single patient due to variants in MYO7A and NEB

In a patient with Usher syndrome and atypical muscle complaints, we have identified two separate variants in MYO7A andNEB genes by exome sequencing. The homozygous variants in these two recessive genes could explain the full phenotype of our patient.

A CRISPR and high-content imaging assay compliant with ACMG/AMP guidelines for clinical variant interpretation in ciliopathies

Ciliopathies are a broad range of inherited developmental and degenerative diseases associated with structural or functional defects in motile or primary non-motile cilia. There are around 200 known ciliopathy disease genes and whilst genetic testing can provide an accurate diagnosis, 24–60% of ciliopathy patients who undergo genetic testing do not receive a genetic diagnosis. This is partly because following current guidelines from the American College of Medical Genetics and the Association for Molecular Pathology, it is difficult to provide a confident clinical diagnosis of disease caused by missense or non-coding variants, which account for more than one-third of cases of disease. Mutations in PRPF31 are the second most common cause of the degenerative retinal ciliopathy autosomal dominant retinitis pigmentosa. Here, we present a high-throughput high-content imaging assay providing quantitative measure of effect of missense variants in PRPF31 which meets the recently published criteria for a baseline standard in vitro test for clinical variant interpretation. This assay utilizes a new PRPF31^+/– human retinal cell line generated using CRISPR gene editing to provide a stable cell line with significantly fewer cilia in which novel missense variants are expressed and characterised. We show that high-content imaging of cells expressing missense variants in a ciliopathy gene on a null background can allow characterisation of variants according to the cilia phenotype. We hope that this will be a useful tool for clinical characterisation of PRPF31 variants of uncertain significance, and can be extended to variant classification in other ciliopathies.

Improving Genetic Testing in Hereditary Cancer by RNA Analysis: Tools to Prioritize Splicing Studies and Challenges in Applying American College of Medical Genetics and Genomics Guidelines

RNA analyses are a potent tool to identify spliceogenic effects of DNA variants, although they are time-consuming and cannot always be performed. We present splicing assays of 20 variants that represent a variety of mutation types in 10 hereditary cancer genes and attempt to incorporate these results into American College of Medical Genetics and Genomics (ACMG) classification guidelines. Sixteen single-nucleotide variants, 3 exon duplications, and 1 single-exon deletion were selected and prioritized by in silico algorithms. RNA was extracted from short-term lymphocyte cultures to perform RT-PCR and Sanger sequencing, and allele-specific expression was assessed whenever possible. Aberrant transcripts were detected in 14 variants (70%). Variant interpretation was difficult, especially comparing old classification standards to generic ACMG guidelines and a proposal was devised to weigh functional analyses at RNA level. According to the ACMG guidelines, only 12 variants were reclassified as pathogenic/likely pathogenic because the other two variants did not gather enough evidence. This study highlights the importance of RNA studies to improve variant classification. However, it also indicates the challenge of incorporating these results into generic ACMG guidelines and the need to refine these criteria gene specifically. Nevertheless, 60% of variants were reclassified, thus improving genetic counseling and surveillance for carriers of these variants.

Detailed Phenotyping and Therapeutic Strategies for Intronic ABCA4 Variants in Stargardt Disease

Stargardt disease is a progressive retinal disorder caused by bi-allelic mutations in the ABCA4 gene that encodes the ATP-binding cassette, subfamily A, member 4 transporter protein. Over the past few years, we and others have identified several pathogenic variants that reside within the introns of ABCA4, including a recurrent variant in intron 36 (c.5196+1137G>A) of which the pathogenicity so far remained controversial. Detailed clinical characterization of this variant confirmed its pathogenic nature, and classified it as an allele of intermediate severity. Moreover, we discovered several additional ABCA4 variants clustering in intron 36. Several of these variants resulted in aberrant splicing of ABCA4, i.e., the inclusion of pseudoexons, while the splicing defects caused by the recurrent c.5196+1137G>A variant strongly increased upon differentiation of patient-derived induced pluripotent stem cells into retina-like cells. Finally, all splicing defects could be rescued by the administration of antisense oligonucleotides that were designed to specifically block the pseudoexon insertion, including rescue in 3D retinal organoids harboring the c.5196+1137G>A variant. Our data illustrate the importance of intronic variants in ABCA4 and expand the therapeutic possibilities for overcoming splicing defects in Stargardt disease.

EDeepSSP: Explainable deep neural networks for exact splice sites prediction

Splice site prediction is crucial for understanding underlying gene regulation, gene function for better genome annotation. Many computational methods exist for recognizing the splice sites. Although most of the methods achieve a competent performance, their interpretability remains challenging. Moreover, all traditional machine learning methods manually extract features, which is tedious job. To address these challenges, we propose a deep learning-based approach (EDeepSSP) that employs convolutional neural networks (CNNs) architecture for automatic feature extraction and effectively predicts splice sites. Our model, EDeepSSP, divulges the opaque nature of CNN by extracting significant motifs and explains why these motifs are vital for predicting splice sites. In this study, experiments have been conducted on six benchmark acceptors and donor datasets of humans, cress, and fly. The results show that EDeepSSP has outperformed many state-of-the-art approaches. EDeepSSP achieves the highest area under the receiver operating characteristic curve (AUC_ROC) and area under the precision-recall curve (AUC_PR) of 99.32% and 99.26% on human donor datasets, respectively. We also analyze various filter activities, feature activations, and extracted significant motifs responsible for the splice site prediction. Further, we validate the learned motifs of our model against known motifs of JASPAR splice site database.

Next-Generation Sequencing of CYP2C19 in Stent Thrombosis: Implications for Clopidogrel Pharmacogenomics

PURPOSE

Describe CYP2C19 sequencing results in the largest series of clopidogrel-treated cases with stent thrombosis (ST), the closest clinical phenotype to clopidogrel resistance. Evaluate the impact of CYP2C19 genetic variation detected by next-generation sequencing (NGS) with comprehensive annotation and functional studies.

METHODS

Seventy ST cases on clopidogrel identified from the PLATO trial (n = 58) and Mayo Clinic biorepository (n = 12) were matched 1:1 with controls for age, race, sex, diabetes mellitus, presentation, and stent type. NGS was performed to cover the entire CYP2C19 gene. Assessment of exonic variants involved measuring in vitro protein expression levels. Intronic variants were evaluated for potential splicing motif variations.

RESULTS

Poor metabolizers (n = 4) and rare CYP2C19*8, CYP2C19*15, and CYP2C19*11 alleles were identified only in ST cases. CYP2C19*17 heterozygote carriers were observed more frequently in cases (n = 29) than controls (n = 18). Functional studies of CYP2C19 exonic variants (n = 11) revealed 3 cases and only 1 control carrying a deleterious variant as determined by in vitro protein expression studies. Greater intronic variation unique to ST cases (n = 169) compared with controls (n = 84) was observed with predictions revealing 13 allele candidates that may lead to a potential disruption of splicing and a loss-of-function effect of CYP2C19 in ST cases.

CONCLUSION

NGS detected CYP2C19 poor metabolizers and paradoxically greater number of so-called rapid metabolizers in ST cases. Rare deleterious exonic variation occurs in 4%, and potentially disruptive intronic alleles occur in 16% of ST cases. Additional studies are required to evaluate the role of these variants in platelet aggregation and clopidogrel metabolism.

Phenotype-genotype correlations in a pseudodominant Stargardt disease pedigree due to a novel ABCA4 deletion-insertion variant causing a splicing defect

BACKGROUND

Deletion-insertion (delins) variants in the retina-specific ATP-binding cassette transporter gene, subfamily A, member 4 (ABCA4) accounts for <1% in Stargardt disease. The consequences of these delins variants on splicing cannot be predicted with certainty without supporting in vitro data.

METHODS

Candidate ABCA4 variants were revealed by genetic and segregation analysis of a family with pseudodominant Stargardt disease using a commercial panel and Sanger sequencing. RNA extracted from patient-derived fibroblasts was analyzed by RT-PCR to evaluate splicing behavior of the ABCA4 variants.

RESULTS

Affected members carrying the novel c.6031_6044delinsAGTATTTAACCAATATTT variant in exon 44 presented with contrasting phenotypes; from early-onset cone-rod dystrophy to late-onset macular dystrophy. This variant resulted in a 56-nucleotide deletion in the mutant allele by activation of a cryptic splice acceptor site which disrupts the reading frame and results in a premature termination codon (p.Ile2003LeufsTer41). If translated, the crucial functional domains near the C-terminus would be truncated from the ABCA4 protein.

CONCLUSION

This work demonstrates the intrafamilial phenotypic variability in a pseudodominant Stargardt disease pedigree and the use of patient-derived fibroblasts to evaluate the effect of a novel ABCA4 delins variant on splicing to complement in silico pathogenicity assessment.

Spliceogen: an integrative, scalable tool for the discovery of splice-altering variants

MOTIVATION

In silico prediction tools are essential for identifying variants which create or disrupt cis-splicing motifs. However, there are limited options for genome-scale discovery of splice-altering variants.

RESULTS

We have developed Spliceogen, a highly scalable pipeline integrating predictions from some of the individually best performing models for splice motif prediction: MaxEntScan, GeneSplicer, ESRseq and Branchpointer.

AVAILABILITY AND IMPLEMENTATION

Spliceogen is available as a command line tool which accepts VCF/BED inputs and handles both single nucleotide variants (SNVs) and indels (https://github.com/VCCRI/Spliceogen). SNV databases with prediction scores are also available, covering all possible SNVs at all genomic positions within all Gencode-annotated multi-exon transcripts.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Epigenome-based splicing prediction using a recurrent neural network

Alternative RNA splicing provides an important means to expand metazoan transcriptome diversity. Contrary to what was accepted previously, splicing is now thought to predominantly take place during transcription. Motivated by emerging data showing the physical proximity of the spliceosome to Pol II, we surveyed the effect of epigenetic context on co-transcriptional splicing. In particular, we observed that splicing factors were not necessarily enriched at exon junctions and that most epigenetic signatures had a distinctly asymmetric profile around known splice sites. Given this, we tried to build an interpretable model that mimics the physical layout of splicing regulation where the chromatin context progressively changes as the Pol II moves along the guide DNA. We used a recurrent-neural-network architecture to predict the inclusion of a spliced exon based on adjacent epigenetic signals, and we showed that distinct spatio-temporal features of these signals were key determinants of model outcome, in addition to the actual nucleotide sequence of the guide DNA strand. After the model had been trained and tested (with >80% precision-recall curve metric), we explored the derived weights of the latent factors, finding they highlight the importance of the asymmetric time-direction of chromatin context during transcription.

Deleterious mis-splicing of STK11 caused by a novel single-nucleotide substitution in the 3' polypyrimidine tract of intron five

Background

Pathogenic variants in STK11, also designated as LKB1, cause Peutz–Jeghers syndrome, which is a rare autosomal dominant disorder characterized by mucocutaneous pigmentation changes, polyposis, and a high risk of cancer.

Methods

A male meeting the clinical diagnostic criteria for Peutz–Jeghers syndrome underwent next‐generation sequencing. To validate the predicted splicing impact of a detected STK11 variant, we performed RNA‐Seq on mRNA extracted from patient‐derived Epstein‐Barr virus‐transformed lymphocytes treated with cycloheximide to inhibit nonsense‐mediated decay ex vivo.

Results

Blood testing identified a novel single‐nucleotide substitution, NM_000455.4:c.735‐10C>A, at the end of the 3′ polypyrimidine tract of intron five in STK11. RNA‐Seq confirmed a predicted eight base pair insertion in the mRNA transcript. Following inhibition of nonsense‐mediated decay, the out‐of‐frame insertion was detected in 50% of all RNA‐Seq reads. This confirmed a strong, deleterious splicing impact of the variant.

Conclusion

We characterized a novel likely pathogenic germline variant in intron five of STK11 associated with Peutz–Jeghers syndrome. The study highlights RNA‐Seq as a useful supplement in hereditary cancer predisposition testing.

A novel missense variant in MYO3A is associated with autosomal dominant high-frequency hearing loss in a German family

Background

MYO3A, encoding the myosin IIIA protein, is associated with autosomal recessive and autosomal dominant nonsyndromic hearing loss. To date, only two missense variants located in the motor‐head domain of MYO3A have been described in autosomal dominant families with progressive, mild‐to‐profound sensorineural hearing loss. These variants alter the ATPase activity of myosin IIIA.

Methods

Exome sequencing of a proband from a three‐generation German family with prelingual, moderate‐to‐profound, high‐frequency hearing loss was performed. Segregation analysis confirmed a dominant inheritance pattern. Regression analysis of mean hearing level thresholds per individual and ear was performed at high‐, mid‐, and low‐frequencies.

Results

A novel heterozygous missense variant c.716T>C, p.(Leu239Pro) in the kinase domain of MYO3A was identified that is predicted in silico as disease causing. High‐frequency, progressive hearing loss was identified.

Conclusion

Correlation analysis of pure‐tone hearing thresholds revealed progressive hearing loss, especially in the high‐frequencies. In the present study, we report the first dominant likely pathogenic variant in MYO3A in a European family and further support MYO3A as an autosomal dominant hearing loss gene.

HRness in Breast and Ovarian Cancers

Ovarian and breast cancers are currently defined by the main pathways involved in the tumorigenesis. The majority are carcinomas, originating from epithelial cells that are in constant division and subjected to cyclical variations of the estrogen stimulus during the female hormonal cycle, therefore being vulnerable to DNA damage. A portion of breast and ovarian carcinomas arises in the context of DNA repair defects, in which genetic instability is the backdrop for cancer initiation and progression. For these tumors, DNA repair deficiency is now increasingly recognized as a target for therapeutics. In hereditary breast/ovarian cancers (HBOC), tumors with BRCA1/2 mutations present an impairment of DNA repair by homologous recombination (HR). For many years, BRCA1/2 mutations were only screened on germline DNA, but now they are also searched at the tumor level to personalize treatment. The reason of the inactivation of this pathway remains uncertain for most cases, even in the presence of a HR-deficient signature. Evidence indicates that identifying the mechanism of HR inactivation should improve both genetic counseling and therapeutic response, since they can be useful as new biomarkers of response.

Splice2Deep: An ensemble of deep convolutional neural networks for improved splice site prediction in genomic DNA

Background

The accurate identification of the exon/intron boundaries is critical for the correct annotation of genes with multiple exons. Donor and acceptor splice sites (SS) demarcate these boundaries. Therefore, deriving accurate computational models to predict the SS are useful for functional annotation of genes and genomes, and for finding alternative SS associated with different diseases. Although various models have been proposed for the in silico prediction of SS, improving their accuracy is required for reliable annotation. Moreover, models are often derived and tested using the same genome, providing no evidence of broad application, i.e. to other poorly studied genomes.

Results

With this in mind, we developed the Splice2Deep models for SS detection. Each model is an ensemble of deep convolutional neural networks. We evaluated the performance of the models based on the ability to detect SS in Homo sapiens, Oryza sativa japonica, Arabidopsis thaliana, Drosophila melanogaster, and Caenorhabditis elegans. Results demonstrate that the models efficiently detect SS in other organisms not considered during the training of the models. Compared to the state-of-the-art tools, Splice2Deep models achieved significantly reduced average error rates of 41.97% and 28.51% for acceptor and donor SS, respectively. Moreover, the Splice2Deep cross-organism validation demonstrates that models correctly identify conserved genomic elements enabling annotation of SS in new genomes by choosing the taxonomically closest model.

Conclusions

The results of our study demonstrated that Splice2Deep both achieved a considerably reduced error rate compared to other state-of-the-art models and the ability to accurately recognize SS in other organisms for which the model was not trained, enabling annotation of poorly studied or newly sequenced genomes. Splice2Deep models are implemented in Python using Keras API; the models and the data are available at https://github.com/SomayahAlbaradei/Splice_Deep.git.

The difficult management of persistent, non-focal congenital hyperinsulinism: A retrospective review from a single, tertiary center

BACKGROUND/OBJECTIVE

Congenital hyperinsulinism (CHI) is a rare, heterogeneous disease with transient or persistent hypoglycemia. Histologically, focal, diffuse, and atypical forms of CHI exist, and at least 11 disease-causing genes have been identified.

METHODS

We retrospectively evaluated the treatment and outcome of a cohort of 40 patients with non-focal, persistent CHI admitted to the International Hyperinsulinism Center, Denmark, from January 2000 to May 2017.

RESULTS

Twenty-two patients (55%) could not be managed with medical monotherapy (diazoxide or octreotide) and six (15%) patients developed severe potential side effects to medication. Surgery was performed in 17 (43%) patients with resection of 66% to 98% of the pancreas. Surgically treated patients had more frequently K_ATP -channel gene mutations (surgical treatment 12/17 vs conservative treatment 6/23, P = .013), highly severe disease (15/17 vs 13/23, P = .025) and clinical onset <30 days of age (15/17 vs 10/23, P = .004). At last follow-up at median 5.3 (range: 0.3-31.3) years of age, 31/40 (78%) patients still received medical treatment, including 12/17 (71%) after surgery. One patient developed diabetes after a 98% pancreatic resection. Problematic treatment status was seen in 7/40 (18%). Only 8 (20%) had clinical remission (three spontaneous, five after pancreatic surgery). Neurodevelopmental impairment (n = 12, 30%) was marginally associated with disease severity (P = .059).

CONCLUSIONS

Persistent, non-focal CHI remains difficult to manage. Neurological impairment in 30% suggests a frequent failure of prompt and adequate treatment. A high rate of problematic treatment status at follow-up demonstrates an urgent need for new medical treatment modalities.

Genetic variants of uncertain significance: How to match scientific rigour and standard of proof in sudden cardiac death?

In many SCD cases, in particular in pediatric age, autopsy can be completely negative and then a post-mortem genetic testing (molecular autopsy) is indicated. In NGS era finding new/rare variants is extremely frequent and, when only variants of unknown significance are found, molecular autopsy fails to find a cause of death. We describe the emblematic case of the sudden death of a 7-year-old girl. We performed a full-body micro-CT analysis, an accurate autopsy, a serum tryptase test and toxicological tests. Since the only macroscopic abnormality we found was a myocardial bridging (length: 1,1 cm, thickness: 0,5 cm) of the left anterior descending coronary artery, a molecular autopsy has been performed. NGS analysis on victim DNA detected rare variants in DPP6, MYH7, SCN2B and NOTCH1 and segregation analysis was then achieved. On the basis of ACMG/AMP (clinical) guidelines, all the found variants were classified as of unknown significance. In other words, both the macroscopic and genetic anomalies we found were of uncertain significance and then the autopsy failed to find the cause of the death. Our case raises three main discussion points: (a) economical, ethical and legal limitations of genetic investigation; (b) risk that genetic testing does not succeed in finding a certain cause of the death; (c) absence of specific guidelines to face the problem of VUS in forensic cases.

In or Out? New Insights on Exon Recognition through Splice-Site Interdependency

Noncanonical splice-site mutations are an important cause of inherited diseases. Based on in vitro and stem-cell-based studies, some splice-site variants show a stronger splice defect than expected based on their predicted effects, suggesting that other sequence motifs influence the outcome. We investigated whether splice defects due to human-inherited-disease-associated variants in noncanonical splice-site sequences in ABCA4, DMD, and TMC1 could be rescued by strengthening the splice site on the other side of the exon. Noncanonical 5′- and 3′-splice-site variants were selected. Rescue variants were introduced based on an increase in predicted splice-site strength, and the effects of these variants were analyzed using in vitro splice assays in HEK293T cells. Exon skipping due to five variants in noncanonical splice sites of exons in ABCA4, DMD, and TMC1 could be partially or completely rescued by increasing the predicted strengths of the other splice site of the same exon. We named this mechanism “splicing interdependency”, and it is likely based on exon recognition by splicing machinery. Awareness of this interdependency is of importance in the classification of noncanonical splice-site variants associated with disease and may open new opportunities for treatments.

Genetic landscape of preterm birth due to cervical insufficiency: Comprehensive gene analysis and patient next-generation sequencing data interpretation

Preterm delivery is both a traumatizing experience for the patient and a burden on the healthcare system. A condition distinguishable by its phenotype in prematurity is cervical insufficiency, where certain cases exhibit a strong genetic component. Despite genomic advancements, little is known about the genetics of human cervix remodeling during pregnancy. Using selected gene approaches, a few studies have demonstrated an association of common gene variants with cervical insufficiency. However, until now, no study has employed comprehensive methods to investigate this important subject matter. In this study, we asked: i) are there genes reliably linked to cervical insufficiency and, if so, what are their roles? and ii) what is the proportion of cases of non-syndromic cervical insufficiency attributable to these genetic variations? We performed next-generation sequencing on 21 patients with a clinical presentation of cervical insufficiency. To assist the sequencing data interpretation, we retrieved all known genes implicated in cervical functioning through a systematic literature analysis and additional gene searches. These genes were then classified according to their relation to the questions being posed by the study. Patients' sequence variants were filtered for pathogenicity and assigned a likelihood of being contributive to phenotype development. Gene extraction and analysis revealed 12 genes primarily linked to cervical insufficiency, the majority of which are known to cause collagenopathies. Ten patients carried disruptive variants potentially contributive to the development of non-syndromic cervical insufficiency. Pathway enrichment analysis of variant genes from our cohort revealed an increased variation burden in genes playing roles in tissue mechanical and biomechanical properties, i.e. collagen biosynthesis and cell-extracellular matrix communications. Consequently, the proposed idea of cervical insufficiency being a subtle form of collagenopathy, now strengthened by our genetic findings, might open up new opportunities for improved patient evaluation and management.

AG-exclusion zone revisited: Lessons to learn from 91 intronic NF1 3' splice site mutations outside the canonical AG-dinucleotides

Uncovering frequent motives of action by which variants impair 3′ splice site (3′ss) recognition and selection is essential to improve our understanding of this complex process. Through several mini‐gene experiments, we demonstrate that the pyrimidine (Y) to purine (R) transversion NM_000267.3(NF1):c.1722‐11T>G, although expected to weaken the polypyrimidine tract, causes exon skipping primarily by introducing a novel AG in the AG‐exclusion zone (AGEZ) between the authentic 3′ss AG and the branch point. Evaluation of 90 additional noncanonical intronic NF1 3′ss mutations confirmed that 63% of all mutations and 89% (49/55) of the single‐nucleotide variants upstream of positions ‐3 interrupt the AGEZ. Of these AGEZ‐interrupting mutations, 24/49 lead to exon skipping suggesting that absence of AG in this region is necessary for accurate 3′ss selection already in the initial steps of splicing. The analysis of 91 noncanonical NF1 3′ss mutations also shows that 90% either introduce a novel AG in the AGEZ, cause a Y>R transversion at position ‐3 or remove ≥2 Ys in the AGEZ. We confirm in a validation cohort that these three motives distinguish spliceogenic from splice‐neutral variants with 85% accuracy and, therefore, are generally applicable to select among variants of unknown significance those likely to affect splicing.

Mortality in a neonate with molybdenum cofactor deficiency illustrates the need for a comprehensive rapid precision medicine system

Neonatal encephalopathy with seizures is a presentation in which rapid whole-genome sequencing (rWGS) has shown clinical utility and improved outcomes. We report a neonate who presented on the third day of life with seizures refractory to antiepileptic medications and neurologic and computerized tomographic findings consistent with severe generalized brain swelling. rWGS revealed compound heterozygous variants in the molybdenum cofactor synthesis gene, type 1A (MOCS1 c.*7 + 5G > A and c.377G > A); a provisional diagnosis of molybdenum cofactor deficiency on day of life 4. An emergency investigational new drug application for intravenous replacement of the MOCS1 product, cyclic pyranopterin monophosphate, was considered, but felt unsuitable in light of the severity of disease and delay in the start of treatment. The patient died on day of life 9 despite having a precise molecular diagnosis within the first week of life. This case illustrates that an rWGS-based molecular diagnosis within the first week of life may be insufficient to improve outcomes. However, it did inform clinical decision-making with regard to resuscitation and predicted long-term outcome. We suggest that to achieve optimal reductions in morbidity and mortality, rWGS must be implemented within a comprehensive rapid precision medicine system (CRPM). Akin to newborn screening (NBS), CRPM will have onboarding, diagnosis, and precision medicine implementation components developed in response to patient and parental needs. Education of health-care providers in a learning model in which ongoing data analyses informs system improvement will be essential for optimal effectiveness of CRPM.

Novel Loss-of-Function Variants in CDC14A are Associated with Recessive Sensorineural Hearing Loss in Iranian and Pakistani Patients

CDC14A encodes the Cell Division Cycle 14A protein and has been associated with autosomal recessive non-syndromic hearing loss (DFNB32), as well as hearing impairment and infertile male syndrome (HIIMS) since 2016. To date, only nine variants have been associated in patients whose initial symptoms included moderate-to-profound hearing impairment. Exome analysis of Iranian and Pakistani probands who both showed bilateral, sensorineural hearing loss revealed a novel splice site variant (c.1421+2T>C, p.?) that disrupts the splice donor site and a novel frameshift variant (c.1041dup, p.Ser348Glnfs*2) in the gene CDC14A, respectively. To evaluate the pathogenicity of both loss-of-function variants, we analyzed the effects of both variants on the RNA-level. The splice variant was characterized using a minigene assay. Altered expression levels due to the c.1041dup variant were assessed using RT-qPCR. In summary, cDNA analysis confirmed that the c.1421+2T>C variant activates a cryptic splice site, resulting in a truncated transcript (c.1414_1421del, p.Val472Leufs*20) and the c.1041dup variant results in a defective transcript that is likely degraded by nonsense-mediated mRNA decay. The present study functionally characterizes two variants and provides further confirmatory evidence that CDC14A is associated with a rare form of hereditary hearing loss.

SpliceFinder: ab initio prediction of splice sites using convolutional neural network

Background

Identifying splice sites is a necessary step to analyze the location and structure of genes. Two dinucleotides, GT and AG, are highly frequent on splice sites, and many other patterns are also on splice sites with important biological functions. Meanwhile, the dinucleotides occur frequently at the sequences without splice sites, which makes the prediction prone to generate false positives. Most existing tools select all the sequences with the two dimers and then focus on distinguishing the true splice sites from those pseudo ones. Such an approach will lead to a decrease in false positives; however, it will result in non-canonical splice sites missing.

Result

We have designed SpliceFinder based on convolutional neural network (CNN) to predict splice sites. To achieve the ab initio prediction, we used human genomic data to train our neural network. An iterative approach is adopted to reconstruct the dataset, which tackles the data unbalance problem and forces the model to learn more features of splice sites. The proposed CNN obtains the classification accuracy of 90.25%, which is 10% higher than the existing algorithms. The method outperforms other existing methods in terms of area under receiver operating characteristics (AUC), recall, precision, and F1 score. Furthermore, SpliceFinder can find the exact position of splice sites on long genomic sequences with a sliding window. Compared with other state-of-the-art splice site prediction tools, SpliceFinder generates results in about half lower false positive while keeping recall higher than 0.8. Also, SpliceFinder captures the non-canonical splice sites. In addition, SpliceFinder performs well on the genomic sequences of Drosophila melanogaster, Mus musculus, Rattus, and Danio rerio without retraining.

Conclusion

Based on CNN, we have proposed a new ab initio splice site prediction tool, SpliceFinder, which generates less false positives and can detect non-canonical splice sites. Additionally, SpliceFinder is transferable to other species without retraining. The source code and additional materials are available at https://gitlab.deepomics.org/wangruohan/SpliceFinder.

Cylindrical spirals in two families: Clinical and genetic investigations

Cylindrical spirals are a rare ultrastructural finding on muscle biopsy, with fewer than 20 reported cases since its first description in 1979. These structures are sometimes observed with tubular aggregates and are thought to comprise longitudinal sarcoplasmic reticulum. While mutations in genes encoding key components of Ca²⁺ handling (ORAI1 and STIM1) underlie tubular aggregate myopathy, no causative genes have been associated with cylindrical spirals. Here we describe two families with cylindrical spirals on muscle biopsy with a suspected genetic cause. In one family we identified a known truncating variant in EBF3, previously associated with a neurodevelopmental disorder. The affected individuals in this family present with clinical features overlapping with those described for EBF3 disease. An isolated proband in the second family harbours bi-allelic truncating variants in TTN and her clinical course and other features on biopsy are highly concordant for titinopathy. From experimental studies, EBF3 is known to be involved in Ca²⁺ regulation in muscle, thus EBF3 dysregulation may represent a novel mechanism of impaired Ca²⁺ handling leading to cylindrical spirals. Additional cases of EBF3 disease or titinopathy with cylindrical spirals need to be identified to support the involvement of these genes in the pathogenesis of cylindrical spirals.

Identification of splice defects due to noncanonical splice site or deep-intronic variants in ABCA4

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

Machine Learning Approaches for the Prioritization of Genomic Variants Impacting Pre-mRNA Splicing

Defects in pre-mRNA splicing are frequently a cause of Mendelian disease. Despite the advent of next-generation sequencing, allowing a deeper insight into a patient's variant landscape, the ability to characterize variants causing splicing defects has not progressed with the same speed. To address this, recent years have seen a sharp spike in the number of splice prediction tools leveraging machine learning approaches, leaving clinical geneticists with a plethora of choices for in silico analysis. In this review, some basic principles of machine learning are introduced in the context of genomics and splicing analysis. A critical comparative approach is then used to describe seven recent machine learning-based splice prediction tools, revealing highly diverse approaches and common caveats. We find that, although great progress has been made in producing specific and sensitive tools, there is still much scope for personalized approaches to prediction of variant impact on splicing. Such approaches may increase diagnostic yields and underpin improvements to patient care.

Using the Chou's 5-steps rule to predict splice junctions with interpretable bidirectional long short-term memory networks

Neural models have been able to obtain state-of-the-art performances on several genome sequence-based prediction tasks. Such models take only nucleotide sequences as input and learn relevant features on their own. However, extracting the interpretable motifs from the model remains a challenge. This work explores various existing visualization techniques in their ability to infer relevant sequence information learnt by a recurrent neural network (RNN) on the task of splice junction identification. The visualization techniques have been modulated to suit the genome sequences as input. The visualizations inspect genomic regions at the level of a single nucleotide as well as a span of consecutive nucleotides. This inspection is performed based on the modification of input sequences (perturbation based) or the embedding space (back-propagation based). We infer features pertaining to both canonical and non-canonical splicing from a single neural model. Results indicate that the visualization techniques produce comparable performances for branchpoint detection. However, in the case of canonical donor and acceptor junction motifs, perturbation based visualizations perform better than back-propagation based visualizations, and vice-versa for non-canonical motifs. The source code of our stand-alone SpliceVisuL tool is available at https://github.com/aaiitggrp/SpliceVisuL.

RNA-Seq Perspectives to Improve Clinical Diagnosis

In recent years, high-throughput next-generation sequencing technology has allowed a rapid increase in diagnostic capacity and precision through different bioinformatics processing algorithms, tools, and pipelines. The identification, annotation, and classification of sequence variants within different target regions are now considered a gold standard in clinical genetic diagnosis. However, this procedure lacks the ability to link regulatory events such as differential splicing to diseases. RNA-seq is necessary in clinical routine in order to interpret and detect among others splicing events and splicing variants, as it would increase the diagnostic rate by up to 10-35%. The transcriptome has a very dynamic nature, varying according to tissue type, cellular conditions, and environmental factors that may affect regulatory events such as splicing and the expression of genes or their isoforms. RNA-seq offers a robust technical analysis of this complexity, but it requires a profound knowledge of computational/statistical tools that may need to be adjusted depending on the disease under study. In this article we will cover RNA-seq analyses best practices applied to clinical routine, bioinformatics procedures, and present challenges of this approach.

High diagnostic yield of direct Sanger sequencing in the diagnosis of neuronal ceroid lipofuscinoses

BACKGROUND

Neuronal ceroid lipofuscinoses are neurodegenerative disorders. To investigate the diagnostic yield of direct Sanger sequencing of the CLN genes, we reviewed Molecular Genetics Laboratory Database for molecular genetic test results of the CLN genes from a single clinical molecular diagnostic laboratory.

METHODS

We reviewed electronic patient charts. We used consent forms and Research Electronic Data Capture questionnaires for the patients from outside of our Institution. We reclassified all variants in the CLN genes.

RESULTS

Six hundred and ninety three individuals underwent the direct Sanger sequencing of the CLN genes for the diagnosis of neuronal ceroid lipofuscinoses. There were 343 symptomatic patients and 350 family members. Ninety-one symptomatic patients had molecular genetic diagnosis of neuronal ceroid lipofuscinoses including CLN1 (PPT1) (n = 10), CLN2 (TPP1) (n = 33), CLN3 (n = 17), CLN5 (n = 7), CLN6 (n = 10), CLN7 (MFSD8) (n = 10), and CLN8 (n = 4) diseases. The diagnostic yield of direct Sanger sequencing of CLN genes was 27% in symptomatic patients. We report detailed clinical and investigation results of 33 NCL patients. Juvenile onset CLN1 (PPT1) and adult onset CLN6 diseases were nonclassical phenotypes.

CONCLUSION

In our study, the diagnostic yield of direct Sanger sequencing was close to diagnostic yield of whole exome sequencing. Developmental regression, cognitive decline, visual impairment and cerebral and/or cerebellar atrophy in brain MRI are significant clinical and neuroimaging denominators to include NCL in the differential diagnosis.

Prevalence of ABCA4 Deep-Intronic Variants and Related Phenotype in An Unsolved "One-Hit" Cohort with Stargardt Disease

We investigated the prevalence of reported deep-intronic variants in a French cohort of 70 patients with Stargardt disease harboring a monoallelic pathogenic variant on the exonic regions of ABCA4. Direct Sanger sequencing of selected intronic regions of ABCA4 was conducted. Complete phenotypic analysis and correlation with the genotype was performed in case a known intronic pathogenic variant was identified. All other variants found on the analyzed sequences were queried for minor allele frequency and possible pathogenicity by in silico predictions. The second mutated allele was found in 14 (20%) subjects. The three known deep-intronic variants found were c.5196+1137G>A in intron 36 (6 subjects), c.4539+2064C>T in intron 30 (4 subjects) and c.4253+43G>A in intron 28 (4 subjects). Even though the phenotype depends on the compound effect of the biallelic variants, a genotype-phenotype correlation suggests that the c.5196+1137G>A was mostly associated with a mild phenotype and the c.4539+2064C>T with a more severe one. A variable effect was instead associated with the variant c.4253+43G>A. In addition, two novel variants, c.768+508A>G and c.859-245_859-243delinsTGA never associated with Stargardt disease before, were identified and a possible splice defect was predicted in silico. Our study calls for a larger cohort analysis including targeted locus sequencing and 3D protein modeling to better understand phenotype-genotype correlations associated with deep-intronic changes and patients’ selection for clinical trials.

Identification and characterization of an IgG sequence variant with an 11 kDa heavy chain C-terminal extension using a combination of mass spectrometry and high-throughput sequencing analysis

Protein primary structure is a potential critical quality attribute for biotherapeutics. Identifying and characterizing any sequence variants present is essential for product development. A sequence variant ~11 kDa larger than the expected IgG mass was observed by size-exclusion chromatography and two-dimensional liquid chromatography coupled with online mass spectrometry. Further characterization indicated that the 11 kDa was added to the heavy chain (HC) Fc domain. Despite the relatively large mass addition, only one unknown peptide was detected by peptide mapping. To decipher the sequence, the transcriptome of the manufacturing cell line was characterized by Illumina RNA-seq. Transcriptome reconstruction detected an aberrant fusion transcript, where the light chain (LC) constant domain sequence was fused to the 3ʹ end of the HC transcript. Translation of this fusion transcript generated an extended peptide sequence at the HC C-terminus corresponding to the observed 11 kDa mass addition. Nanopore-based genome sequencing showed multiple copies of the plasmid had integrated in tandem with one copy missing the 5ʹ end of the plasmid, deleting the LC variable domain. The fusion transcript was due to read-through of the HC terminator sequence into the adjacent partial LC gene and an unexpected splicing event between a cryptic splice-donor site at the 3ʹ end of the HC and the splice acceptor site at the 5ʹ end of the LC constant domain. Our study demonstrates that combining protein physicochemical characterization with genomic and transcriptomic analysis of the manufacturing cell line greatly improves the identification of sequence variants and understanding of the underlying molecular mechanisms.

Pathogenic Abnormal Splicing Due to Intronic Deletions that Induce Biophysical Space Constraint for Spliceosome Assembly

A precise genetic diagnosis is the single most important step for families with genetic disorders to enable personalized and preventative medicine. In addition to genetic variants in coding regions (exons) that can change a protein sequence, abnormal pre-mRNA splicing can be devastating for the encoded protein, inducing a frameshift or in-frame deletion/insertion of multiple residues. Non-coding variants that disrupt splicing are extremely challenging to identify. Stemming from an initial clinical discovery in two index Australian families, we define 25 families with genetic disorders caused by a class of pathogenic non-coding splice variant due to intronic deletions. These pathogenic intronic deletions spare all consensus splice motifs, though they critically shorten the minimal distance between the 5' splice-site (5'SS) and branchpoint. The mechanistic basis for abnormal splicing is due to biophysical constraint precluding U1/U2 spliceosome assembly, which stalls in A-complexes (that bridge the 5'SS and branchpoint). Substitution of deleted nucleotides with non-specific sequences restores spliceosome assembly and normal splicing, arguing against loss of an intronic element as the primary causal basis. Incremental lengthening of 5'SS-branchpoint length in our index EMD case subject defines 45-47 nt as the critical elongation enabling (inefficient) spliceosome assembly for EMD intron 5. The 5'SS-branchpoint space constraint mechanism, not currently factored by genomic informatics pipelines, is relevant to diagnosis and precision medicine across the breadth of Mendelian disorders and cancer genomics.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5′ and 3′ consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Tailoring the American College of Medical Genetics and Genomics and the Association for Molecular Pathology Guidelines for the Interpretation of Sequenced Variants in the FBN1 Gene for Marfan Syndrome: Proposal for a Disease- and Gene-Specific Guideline

BACKGROUND

The introduction of next-generation sequencing techniques has substantially increased the identification of new genetic variants and hence the necessity of accurate variant interpretation. In 2015, the American College of Medical Genetics and Genomics and the Association for Molecular Pathology proposed new variant interpretation guidelines. Gene-specific characteristics were, however, not considered, sometimes leading to inconsistent variant interpretation.

METHODS

To allow a more uniform interpretation of variants in the FBN1 (fibrillin-1) gene, causing Marfan syndrome, we tailored these guidelines to this gene and disease. We adapted 15 of the 28 general criteria and classified 713 FBN1 variants previously identified in our laboratory as causal mutation or variant of uncertain significance according to these adapted guidelines. We then compared the agreement between previous methods and the adapted American College of Medical Genetics and Genomics and the Association for Molecular Pathology criteria.

RESULTS

Agreement between the methods was 86.4% (K-alpha, 0.6). Application of the tailored guidelines resulted in an increased number of variants of uncertain significance (14.5% to 24.2%). Of the 85 variants that were downscaled to likely benign or variant of uncertain significance, 59.7% were missense variants outside a well-established functional site. Available clinical- or segregation data, necessary to further classify these types of variants, were in many cases insufficient to aid the classification.

CONCLUSIONS

Our study shows that classification of variants remains challenging and may change over time. Currently, a higher level of evidence is necessary to classify a variant as pathogenic. Gene-specific guidelines may be useful to allow a more precise and uniform interpretation of the variants to accurately support clinical decision-making.

A deep intronic splice mutation of STAT3 underlies hyper IgE syndrome by negative dominance

Heterozygous in-frame mutations in coding regions of human STAT3 underlie the only known autosomal dominant form of hyper IgE syndrome (AD HIES). About 5% of familial cases remain unexplained. The mutant proteins are loss-of-function and dominant-negative when tested following overproduction in recipient cells. However, the production of mutant proteins has not been detected and quantified in the cells of heterozygous patients. We report a deep intronic heterozygous STAT3 mutation, c.1282-89C>T, in 7 relatives with AD HIES. This mutation creates a new exon in the STAT3 complementary DNA, which, when overexpressed, generates a mutant STAT3 protein (D427ins17) that is loss-of-function and dominant-negative in terms of tyrosine phosphorylation, DNA binding, and transcriptional activity. In immortalized B cells from these patients, the D427ins17 protein was 2 kDa larger and 4-fold less abundant than wild-type STAT3, on mass spectrometry. The patients' primary B and T lymphocytes responded poorly to STAT3-dependent cytokines. These findings are reminiscent of the impaired responses of leukocytes from other patients with AD HIES due to typical STAT3 coding mutations, providing further evidence for the dominance of the mutant intronic allele. These findings highlight the importance of sequencing STAT3 introns in patients with HIES without candidate variants in coding regions and essential splice sites. They also show that AD HIES-causing STAT3 mutant alleles can be dominant-negative even if the encoded protein is produced in significantly smaller amounts than wild-type STAT3.

AGT haplotype in ITGA4 gene is related to antibody-mediated rejection in heart transplant patients

Introduction

One of the main problems involved in heart transplantation (HT) is antibody-mediated rejection (AMR). Many aspects of AMR are still unresolved, including its etiology, diagnosis and treatment. In this project, we hypothesize that variants in genes involved in B-cell biology in HT patients can yield diagnostic and prognostic information about AMR.

Methods

Genetic variants in 61 genes related to B-cell biology were analyzed by next generation sequencing in 46 HT patients, 23 with and 23 without AMR.

Results

We identified 3 single nucleotide polymorphisms in ITGA4 gene (c.1845G>A, c.2633A>G, and c.2883C>T) that conformed the haplotype AGT-ITGA4. This haplotype is associated with the development of AMR. Moreover, AMR patients with the haplotype AGT-ITGA4 present lower levels of integrin α-4 in serum samples compared to the reference GAC haplotype in control patients.

Conclusion

We can conclude that polymorphisms in genes related to the biology of B-cells could have an important role in the development of AMR. In fact, the AGT haplotype in ITGA4 gene could potentially increase the risk of AMR.

Comprehensive elaboration of database resources utilized in next-generation sequencing-based tumor somatic mutation detection

The rapid evolution of next-generation sequencing (NGS)-based tumor genomic profile detection and the emergence of molecularly targeted therapies have enabled precision oncology. In NGS-based analysis, various types of databases have been developed to perform different functions. However, many problems still exist when using these public databases. Therefore, it is important to better understand the characteristics and limitations of each database and have them complement each other to provide useful clinical evidence for NGS testing. In this review, we elaborate on the important role of databases and their concrete applications in NGS-based somatic mutation detection. We introduce the typically used databases for sequence alignment, variant filtration, and variant interpretation, and compare the differences between the databases with similar functions. Subsequently, we determine the limitations of each database and provide the corresponding solutions. Furthermore, we present an overview diagram to clearly illustrate the database used in the entire NGS-based somatic mutation detection pipeline.

Cost-effective molecular inversion probe-based ABCA4 sequencing reveals deep-intronic variants in Stargardt disease

PURPOSE

Stargardt disease (STGD1) is caused by biallelic mutations in ABCA4, but many patients are genetically unsolved due to insensitive mutation-scanning methods. We aimed to develop a cost-effective sequencing method for ABCA4 exons and regions carrying known causal deep-intronic variants.

METHODS

Fifty exons and 12 regions containing 14 deep-intronic variants of ABCA4 were sequenced using double-tiled single molecule Molecular Inversion Probe (smMIP)-based next-generation sequencing. DNAs of 16 STGD1 cases carrying 29 ABCA4 alleles and of four healthy persons were sequenced using 483 smMIPs. Thereafter, DNAs of 411 STGD1 cases with one or no ABCA4 variant were sequenced. The effect of novel noncoding variants on splicing was analyzed using in vitro splice assays.

RESULTS

Thirty-four ABCA4 variants previously identified in 16 STGD1 cases were reliably identified. In 155/411 probands (38%), two causal variants were identified. We identified 11 deep-intronic variants present in 62 alleles. Two known and two new noncanonical splice site variants showed splice defects, and one novel deep-intronic variant (c.4539+2065C>G) resulted in a 170-nt mRNA pseudoexon insertion (p.[Arg1514Lysfs*35,=]).

CONCLUSIONS

smMIPs-based sequence analysis of coding and selected noncoding regions of ABCA4 enabled cost-effective mutation detection in STGD1 cases in previously unsolved cases.

A novel phenotype with splicing mutation identified in a Chinese family with desminopathy

BACKGROUND

Desminopathy, a hereditary myofibrillar myopathy, mainly results from the desmin gene (DES) mutations. Desminopathy involves various phenotypes, mainly including different cardiomyopathies, skeletal myopathy, and arrhythmia. Combined with genotype, it helps us precisely diagnose and treat for desminopathy.

METHODS

Sanger sequencing was used to characterize DES variation, and then a minigene assay was used to verify the effect of splice-site mutation on pre-mRNA splicing. Phenotypes were analyzed based on clinical characteristics associated with desminopathy.

RESULTS

A splicing mutation (c.735+1G>T) in DES was detected in the proband. A minigene assay revealed skipping of the whole exon 3 and transcription of abnormal pre-mRNA lacking 32 codons. Another affected family member who carried the identical mutation, was identified with a novel phenotype of desminopathy, non-compaction of ventricular myocardium. There were 2 different phenotypes varied in cardiomyopathy and skeletal myopathy among the 2 patients, but no significant correlation between genotype and phenotype was identified.

CONCLUSIONS

We reported a novel phenotype with a splicing mutation in DES, enlarging the spectrum of phenotype in desminopathy. Molecular studies of desminopathy should promote our understanding of its pathogenesis and provide a precise molecular diagnosis of this disorder, facilitating clinical prevention and treatment at an early stage.