Performance evaluation of differential splicing analysis methods and splicing analytics platform construction

Identification of deleterious variants in nine polycystic kidney disease affected families

Polycystic kidney disease (PKD) is common genetic renal disorder. In present study, we performed WES to identify pathogenic variant in nine families including 26 patients with PKD and 19 unaffected members. The eight pathogenic variants were identified in known PKD associated genes including PKD1 (n = 6), PKD2 (n = 1), and OFD1 (n = 1) in eight families. There is one missense, one stopgain, two non-frameshifts, two canonical splicing variants, three frameshift variants and one potential non-canonical splicing variant (NCSV) in 8 families. The six variants were novel variants and not reported in ClinVar database. In addition, the compound heterozygous variants in PKHD1 were identified including one frameshift variants (PKHD1: NM_138694.4, c.9841del, p.S3281Lfs*4) and one non-canonical splicing variant (PKHD1: NM_138694.4, c.6332 + 40A > G) which were defined as deleterious variant by four splicing prediction tools (CADD-splice, SpliceAI, Spliceogen, Squirl). We used the minigene method to validate whether the prioritized potential NSCVs disrupt the typical mRNA splicing process and found abnormally larger PCR production of minigene carrying potential NCSV comparing to wild-type minigene. Sanger sequencing confirmed the 39-bp insertion of intron 38 between exon 38 and exon 39, which results in non-frameshift and 13 amino acid insertions. In conclusion, our study expands the variant spectrum and highlight the important role of non-canonical splicing variant in PKD.

Achieving an optimal pregnancy outcome through the combined utilization of micro-TESE and ICSI in cryptorchidism associated with a non-canonical splicing variant in RXFP2

PURPOSE

To identify the genetic cause of a cryptorchidism patient carrying a non-canonical splicing variant highlighted by SPCards platform in RXFP2 and to provide a comprehensive overview of RXFP2 variants with cryptorchidism correlation.

METHODS

We identified a homozygous non-canonical splicing variant by whole-exome sequencing and Sanger sequencing in a case with cryptorchidism and non-obstructive azoospermia (NOA). As the pathogenicity of this non-canonical splicing variant remained unclear, we initially utilized the SPCards platform to predict its pathogenicity. Subsequently, we employed a minigene splicing assay to further evaluate the influence of the identified splicing variant. Microdissection testicular sperm extraction (micro-TESE) combined with intracytoplasmic sperm injection (ICSI) was performed. PubMed and Human Genome Variant Database (HGMD) were queried to search for RXFP2 variants.

RESULTS

We identified a homozygous non-canonical splicing variant (NM_130806: c.1376-12A > G) in RXFP2, and confirmed this variant caused aberrant splicing of exons 15 and 16 of the RXFP2 gene: 11 bases were added in front of exon 16, leading to an abnormal transcript initiation and a frameshift. Fortunately, the patient successfully obtained his biological offspring through micro-TESE combined with ICSI. Four cryptorchidism-associated variants in RXFP2 from 90 patients with cryptorchidism were identified through a literature search in PubMed and HGMD, with different inheritance patterns.

CONCLUSION

This is the first cryptorchidism case carrying a novel causative non-canonical splicing RXFP2 variant. The combined approach of micro-TESE and ICSI contributed to an optimal pregnancy outcome. Our literature review demonstrated that RXFP2 variants caused cryptorchidism in a recessive inheritance pattern, rather than a dominant pattern.

VarCards2: an integrated genetic and clinical database for ACMG-AMP variant-interpretation guidelines in the human whole genome

VarCards, an online database, combines comprehensive variant- and gene-level annotation data to streamline genetic counselling for coding variants. Recognising the increasing clinical relevance of non-coding variations, there has been an accelerated development of bioinformatics tools dedicated to interpreting non-coding variations, including single-nucleotide variants and copy number variations. Regrettably, most tools remain as either locally installed databases or command-line tools dispersed across diverse online platforms. Such a landscape poses inconveniences and challenges for genetic counsellors seeking to utilise these resources without advanced bioinformatics expertise. Consequently, we developed VarCards2, which incorporates nearly nine billion artificially generated single-nucleotide variants (including those from mitochondrial DNA) and compiles vital annotation information for genetic counselling based on ACMG-AMP variant-interpretation guidelines. These annotations include (I) functional effects; (II) minor allele frequencies; (III) comprehensive function and pathogenicity predictions covering all potential variants, such as non-synonymous substitutions, non-canonical splicing variants, and non-coding variations and (IV) gene-level information. Furthermore, VarCards2 incorporates 368 820 266 documented short insertions and deletions and 2 773 555 documented copy number variations, complemented by their corresponding annotation and prediction tools. In conclusion, VarCards2, by integrating over 150 variant- and gene-level annotation sources, significantly enhances the efficiency of genetic counselling and can be freely accessed at http://www.genemed.tech/varcards2/.

Novel MEIOB pathogenic variants including a homozygous non-canonical splicing variant, cause meiotic arrest and human non-obstructive azoospermia

Non-obstructive azoospermia (NOA) is the most severe form of human male infertility, and the genetic causes of NOA with meiotic arrest remain largely unclear. In this study, we identified novel compound heterozygous MEIOB variants (c.814C > T: p.R272X and c.976G > A: p.A326T) and a previously undescribed homozygous non-canonical splicing variant of MEIOB (c.528 + 3A > C) in two NOA-affected individuals from two irrelevant Chinese families. MEIOB missense variant (p.A326T) significantly reduced protein abundance and nonsense variant (p.R272X) produced a truncated protein. Both of two variants impaired the MEIOB-SPATA22 interaction. The MEIOB non-canonical splicing variant resulted in whole Exon 6 skipping by minigene assay, which was predicted to produce a frameshift truncated protein (p.S111Rfs*32). Histological and immunostaining analysis indicated that both patients exhibited a similar phenotype as we previously reported in Meiob mutant mice, that is, absence of spermatids in seminiferous tubules and meiotic arrest. Our study identified three novel pathogenic variants of MEIOB in NOA patients, extending the mutation spectrum of the MEIOB and highlighting the contribution of meiotic recombination related genes in human fertility.

Prioritizing de novo potential non-canonical splicing variants in neurodevelopmental disorders

BACKGROUND

Genomic variants outside of the canonical splicing site (±2) may generate abnormal mRNA splicing, which are defined as non-canonical splicing variants (NCSVs). However, the clinical interpretation of NCSVs in neurodevelopmental disorders (NDDs) is largely unknown.

METHODS

We investigated the contribution of NCSVs to NDDs from 345,787 de novo variants (DNVs) in 47,574 patients with NDDs. We performed functional enrichment and protein-protein interaction analysis to assess the association between genes carrying prioritised NCSVs and NDDs. Minigene was used to validate the impact of NCSVs on mRNA splicing.

FINDINGS

We observed significantly more NCSVs (p = 0.02, odds ratio [OR] = 2.05) among patients with NDD than in controls. Both canonical splicing variants (CSVs) and NCSVs contributed to an equal proportion of patients with NDD (0.76% vs. 0.82%). The candidate genes carrying NCSVs were associated with glutamatergic synapse and chromatin remodelling. Minigene successfully validated 59 of 79 (74.68%) NCSVs that led to abnormal splicing in 40 candidate genes, and 9 of the genes (ARID1B, KAT6B, TCF4, SMARCA2, SHANK3, PDHA1, WDR45, SCN2A, SYNGAP1) harboured recurrent NCSVs with the same variant present in more than two unrelated patients with NDD. Moreover, 36 of 59 (61.02%) NCSVs are novel clinically relevant variants, including 34 unreported and 2 clinically conflicting interpretations or of uncertain significance NCSVs in the ClinVar database.

INTERPRETATION

This study highlights the common pathology and clinical importance of NCSVs in unsolved patients with NDD.

FUNDING

The present study was funded by grants from the National Natural Science Foundation of China, China Postdoctoral Science Foundation, the Hunan Youth Science and Technology Innovation Talent Project, the Provincial Natural Science Foundation of Hunan, The Scientific Research Program of FuRong laboratory, and the Natural Science Project of the University of Anhui Province.

FexSplice: A LightGBM-Based Model for Predicting the Splicing Effect of a Single Nucleotide Variant Affecting the First Nucleotide G of an Exon

Single nucleotide variants (SNVs) affecting the first nucleotide G of an exon (Fex-SNVs) identified in various diseases are mostly recognized as missense or nonsense variants. Their effect on pre-mRNA splicing has been seldom analyzed, and no curated database is available. We previously reported that Fex-SNVs affect splicing when the length of the polypyrimidine tract is short or degenerate. However, we cannot readily predict the splicing effects of Fex-SNVs. We here scrutinized the available literature and identified 106 splicing-affecting Fex-SNVs based on experimental evidence. We similarly identified 106 neutral Fex-SNVs in the dbSNP database with a global minor allele frequency (MAF) of more than 0.01 and less than 0.50. We extracted 115 features representing the strength of splicing cis-elements and developed machine-learning models with support vector machine, random forest, and gradient boosting to discriminate splicing-affecting and neutral Fex-SNVs. Gradient boosting-based LightGBM outperformed the other two models, and the length and nucleotide compositions of the polypyrimidine tract played critical roles in the discrimination. Recursive feature elimination showed that the LightGBM model using 15 features achieved the best performance with an accuracy of 0.80 ± 0.12 (mean and SD), a Matthews Correlation Coefficient (MCC) of 0.57 ± 0.15, an area under the curve of the receiver operating characteristics curve (AUROC) of 0.86 ± 0.08, and an area under the curve of the precision-recall curve (AUPRC) of 0.87 ± 0.09 using a 10-fold cross-validation. We developed a web service program, named FexSplice that accepts a genomic coordinate either on GRCh37/hg19 or GRCh38/hg38 and returns a predicted probability of aberrant splicing of A, C, and T variants.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

Alternative splicing: a new breakthrough for understanding tumorigenesis and potential clinical applications

BACKGROUND

Alternative splicing (AS) is a post-transcriptional process that produces transcript variants, thus leading to transcriptome complexity. Recently, the scope of AS studies has been greatly expanded toward clinical applications owing to the abundance of RNA sequencing data.

OBJECTIVE

This review consists of two parts. We first summarize bioinformatic resources that are useful for large-scale cancer-related AS studies. We then highlight the research efforts to utilize AS events for predicting clinical outcomes and planning therapeutic strategies.

RESULTS

Computational approaches to interrogate AS events have been reviewed under three categories: (1) databases to provide functional and clinical annotation of AS events, (2) analytical tools to identify cancer-associated AS event, and (3) methods to identify splicing-related DNA variants and splicing-derived neoantigens. We also present the recent progress in exploring the clinical utility of AS under four categories: (1) identification of AS events for cancer prognosis, (2) utilization of AS events in molecular classification of various cancers, (3) regulatory mechanisms of AS underlying drug resistance, and (4) potential use of AS in cancer therapy.

CONCLUSION

This review will be helpful for understanding the biological implications of AS in cancer and facilitate the development of AS markers for cancer prognosis and treatment. We anticipate that future studies will lead to the application of genome-wide AS profiles in cancer precision medicine.

SpliceAI-visual: a free online tool to improve SpliceAI splicing variant interpretation

SpliceAI is an open-source deep learning splicing prediction algorithm that has demonstrated in the past few years its high ability to predict splicing defects caused by DNA variations. However, its outputs present several drawbacks: (1) although the numerical values are very convenient for batch filtering, their precise interpretation can be difficult, (2) the outputs are delta scores which can sometimes mask a severe consequence, and (3) complex delins are most often not handled. We present here SpliceAI-visual, a free online tool based on the SpliceAI algorithm, and show how it complements the traditional SpliceAI analysis. First, SpliceAI-visual manipulates raw scores and not delta scores, as the latter can be misleading in certain circumstances. Second, the outcome of SpliceAI-visual is user-friendly thanks to the graphical presentation. Third, SpliceAI-visual is currently one of the only SpliceAI-derived implementations able to annotate complex variants (e.g., complex delins). We report here the benefits of using SpliceAI-visual and demonstrate its relevance in the assessment/modulation of the PVS1 classification criteria. We also show how SpliceAI-visual can elucidate several complex splicing defects taken from the literature but also from unpublished cases. SpliceAI-visual is available as a Google Colab notebook and has also been fully integrated in a free online variant interpretation tool, MobiDetails ( https://mobidetails.iurc.montp.inserm.fr/MD ).

Computational prediction of human deep intronic variation

BACKGROUND

The adoption of whole-genome sequencing in genetic screens has facilitated the detection of genetic variation in the intronic regions of genes, far from annotated splice sites. However, selecting an appropriate computational tool to discriminate functionally relevant genetic variants from those with no effect is challenging, particularly for deep intronic regions where independent benchmarks are scarce.

RESULTS

In this study, we have provided an overview of the computational methods available and the extent to which they can be used to analyze deep intronic variation. We leveraged diverse datasets to extensively evaluate tool performance across different intronic regions, distinguishing between variants that are expected to disrupt splicing through different molecular mechanisms. Notably, we compared the performance of SpliceAI, a widely used sequence-based deep learning model, with that of more recent methods that extend its original implementation. We observed considerable differences in tool performance depending on the region considered, with variants generating cryptic splice sites being better predicted than those that potentially affect splicing regulatory elements. Finally, we devised a novel quantitative assessment of tool interpretability and found that tools providing mechanistic explanations of their predictions are often correct with respect to the ground - information, but the use of these tools results in decreased predictive power when compared to black box methods.

CONCLUSIONS

Our findings translate into practical recommendations for tool usage and provide a reference framework for applying prediction tools in deep intronic regions, enabling more informed decision-making by practitioners.