An increased specificity score matrix for the prediction of SF2/ASF-specific exonic splicing enhancers

In silico RNA isoform screening to identify potential cancer driver exons with therapeutic applications

Alternative splicing is crucial for cancer progression and can be targeted pharmacologically, yet identifying driver exons genome-wide remains challenging. We propose identifying such exons by associating statistically gene-level cancer dependencies from knockdown viability screens with splicing profiles and gene expression. Our models predict the effects of splicing perturbations on cell proliferation from transcriptomic data, enabling in silico RNA screening and prioritizing targets for splicing-based therapies. We identified 1,073 exons impacting cell proliferation, many from genes not previously linked to cancer. Experimental validation confirms their influence on proliferation, especially in highly proliferative cancer cell lines. Integrating pharmacological screens with splicing dependencies highlights the potential driver exons affecting drug sensitivity. Our models also allow predicting treatment outcomes from tumor transcriptomes, suggesting applications in precision oncology. This study presents an approach to identifying cancer driver exon and their therapeutic potential, emphasizing alternative splicing as a cancer target.

Comprehensive genomic filtering algorithm to expose the cause of skewed X chromosome inactivation. The proof of concept in female haemophilia expression

BACKGROUND

Exploring the expression of X linked disorders like haemophilia A (HA) in females involves understanding the balance achieved through X chromosome inactivation (XCI). Skewed XCI (SXCI) may be involved in symptomatic HA carriers. We aimed to develop an approach for dissecting the specific cause of SXCI and verify its value in HA.

METHODS

A family involving three females (two symptomatic with severe/moderate HA: I.2, the mother, and II.1, the daughter; one asymptomatic: II.2) and two related affected males (I.1, the father and I.3, the maternal uncle) was studied. The genetic analysis included F8 mutational screening, multiplex ligation-dependent probe amplification, SNP microarray, whole exome sequencing (WES) and Sanger sequencing. XCI patterns were assessed in ectoderm/endoderm and mesoderm-derived tissues using AR-based and RP2-based systems.

RESULTS

The comprehensive family analysis identifies I.2 female patient as a heterozygous carrier of F8:p.(Ser1414Ter) excluding copy number variations. A consistent XCI pattern of 99.5% across various tissues was observed. A comprehensive filtering algorithm for WES data was designed, developed and applied to I.2. A Gly58Arg missense variant in VMA21 was revealed as the cause for SXCI.Each step of the variant filtering system takes advantage of publicly available genomic databases, non-SXCI controls and case-specific molecular data, and aligns with established concepts in the theoretical background of SXCI.

CONCLUSION

This study acts as a proof of concept for our genomic filtering algorithm's clinical utility in analysing X linked disorders. Our findings clarify the molecular aspects of SXCI and improve genetic diagnostics and counselling for families with X linked diseases like HA.

Regulating PCCA gene expression by modulation of pseudoexon splicing patterns to rescue enzyme activity in propionic acidemia

Pseudoexons are nonfunctional intronic sequences that can be activated by deep-intronic sequence variation. Activation increases pseudoexon inclusion in mRNA and interferes with normal gene expression. The PCCA c.1285-1416A>G variation activates a pseudoexon and causes the severe metabolic disorder propionic acidemia by deficiency of the propionyl-CoA carboxylase enzyme encoded by PCCA and PCCB. We characterized this pathogenic pseudoexon activation event in detail and identified hnRNP A1 to be important for normal repression. The PCCA c.1285-1416A>G variation disrupts an hnRNP A1-binding splicing silencer and simultaneously creates a splicing enhancer. We demonstrate that blocking this region of regulation with splice-switching antisense oligonucleotides restores normal splicing and rescues enzyme activity in patient fibroblasts and in a cellular model created by CRISPR gene editing. Interestingly, the PCCA pseudoexon offers an unexploited potential to upregulate gene expression because healthy tissues show relatively high inclusion levels. By blocking inclusion of the nonactivated wild-type pseudoexon, we can increase both PCCA and PCCB protein levels, which increases the activity of the heterododecameric enzyme. Surprisingly, we can increase enzyme activity from residual levels in not only patient fibroblasts harboring PCCA missense variants but also those harboring PCCB missense variants. This is a potential treatment strategy for propionic acidemia.

Modulation of Suppressive Activity and Proliferation of Human Regulatory T Cells by Splice-Switching Oligonucleotides Targeting FoxP3 Pre-mRNA

The maturation, development, and function of regulatory T cells (Tregs) are under the control of the crucial transcription factor Forkhead Box Protein 3 (FoxP3). Through alternative splicing, the human FoxP3 gene produces four different splice variants: a full-length variant (FL) and truncated variants with deletions of each of exons 2 (∆2 variant) or 7 (∆7 variant) or a deletion of both exons (∆2∆7 variant). Their involvement in the biology of Tregs as well as their association with autoimmune diseases remains to be clarified. The aim of this work was to induce a single FoxP3 splice variant in human Tregs by splice switching oligonucleotides and to monitor their phenotype and proliferative and suppressive activity. We demonstrated that Tregs from peripheral blood from patients with multiple sclerosis preferentially expressed truncated splice variants, while the FL variant was the major variant in healthy donors. Tregs with induced expression of truncated FoxP3 splice variants demonstrated lower suppressive activity than those expressing FL variants. Reduced suppression was associated with the decreased expression of Treg-associated suppressive surface molecules and the production of cytokines. The deletion of exons 2 and/or 7 also reduced the cell proliferation rate. The results of this study show an association between FoxP3 splice variants and Treg function and proliferation. The modulation of Treg suppressive activity by the induction of the FoxP3 FL variant can become a promising strategy for regenerative immunotherapy.

Multiocular defect in the Old English Sheepdog: A canine form of Stickler syndrome type II associated with a missense variant in the collagen-type gene COL11A1

Multiocular defect has been described in different canine breeds, including the Old English Sheepdog. Affected dogs typically present with multiple and various ocular abnormalities. We carried out whole genome sequencing on an Old English Sheepdog that had been diagnosed with hereditary cataracts at the age of five and then referred to a board-certified veterinary ophthalmologist due to owner-reported visual deterioration. An ophthalmic assessment revealed that there was bilateral vitreal degeneration, macrophthalmos, and spherophakia in addition to cataracts. Follow-up consultations revealed cataract progression, retinal detachment, uveitis and secondary glaucoma. Whole genome sequence filtered variants private to the case, shared with another Old English Sheepdog genome and predicted to be deleterious were genotyped in an initial cohort of six Old English Sheepdogs (three affected by multiocular defect and three control dogs without evidence of inherited eye disease). Only one of the twenty-two variants segregated correctly with multiocular defect. The variant is a single nucleotide substitution, located in the collagen-type gene COL11A1, c.1775T>C, that causes an amino acid change, p.Phe1592Ser. Genotyping of an additional 14 Old English Sheepdogs affected by multiocular defect revealed a dominant mode of inheritance with four cases heterozygous for the variant. Further genotyping of hereditary cataract-affected Old English Sheepdogs revealed segregation of the variant in eight out of nine dogs. In humans, variants in the COL11A1 gene are associated with Stickler syndrome type II, also dominantly inherited.

Insertion/deletions within the bovine FoxO1 gene and their association analysis with growth traits in three Chinese cattle breeds

FOXO1 (FKHR) gene, as a transcription factor, plays a vital role in animal growth and development, participating in many biological processes. The aim of this study was to ascertain Insertion/deletions (Indels) polymorphism within bovine FoxO1 gene in 679 Chinese adult cows and associate them with stature traits. Two Indels (named as Indel-3 and Indel-4, recorded as rs383545622 and rs525318770 in NCBI, respectively) were successfully genotyped by the Once PCR method, which was reliable, rapid and cost effective for simultaneous detection of two or more Indels. Indel-3 and Indel-4 were located at the second intron. All four different haplotypes (H1: D3D4, H2: I3D4, H3: D3I4, H4: I3I4) could be identified, and the D (del-) allele, DD (del-/del-) genotype and D3D4 haplotype retained the highest frequency. However, individuals with DI (D3I3, D4I4 or H1H4/H2H3 genotype) showed significantly better phenotypic traits than those with the other genotypes in Nanyang cattle, showing a hybrid vigor. The results implied that this DI genotype can be applied to early selective breeding to improve the productivity of Nanyang cattle. Our results suggested that these two Indels within the bovine FoxO1 gene might be used as genetic markers for marker-assisted selection (MAS) in cattle breeding and genetics.

A patient-derived organoid-based study identified an ASO targeting SNORD14E for endometrial cancer through reducing aberrant FOXM1 Expression and β-catenin nuclear accumulation

BACKGROUND

Most of the endometrial cancer (EC) patients are diagnosis in early stage with a good prognosis while the patients with locally advanced recurrent or metastatic result in a poor prognosis. Adjuvant therapy could benefit the prognosis of patients with high-risk factors. Unfortunately, the molecular classification of great prognostic value has not yet reached an agreement and need to be further refined. The present study aims to identify new targets that have prognostic value in EC based on the method of EC patient-derived organ-like organs (PDOs), and further investigate their efficacy and mechanism.

METHODS

The Cancer Genome Atlas (TCGA) database was used to determine SNORD14E expression. The effects of SNORD14E were investigated using CCK8, Transwell, wound-healing assays, and a xenograft model experiment; apoptosis was measured by flow cytometry. Antisense oligonucleotide (ASO) targeting SNORD14E was designed and patient-derived organoids (PDO) models in EC patients was established. A xenograft mouse and PDO model were employed to evaluate the effects of ASO targeting SNORD14E. RNA-seq, Nm-seq, and RNA immunoprecipitation (RIP) experiments were employed to confirm the alternative splicing (AS) and modification induced by SNORD14E. A minigene reporter gene assay was conducted to confirm AS and splicing factors on a variable exon. Actinomycin-d (Act-D) and Reverse Transcription at Low deoxy-ribonucleoside triphosphate concentrations followed by PCR (RTL-P) were utilized to confirm the effects of 2'-O methylation modification on FOXM1.

RESULTS

We found that SNORD14E was overexpressed in EC tissues and patients with high expressed SNORD14E were distributed in the TCGA biomolecular classification subgroups without difference. Further, SNORD14E could reduce disease-free survival (DFS) and recurrence free survival (RFS) of EC patients. SNORD14E promoted proliferation, migration, and invasion and inhibited the apoptosis of EC cells in vitro. ASOs targeting SNORD14E inhibited cell proliferation, migration, invasion while promoted cell apoptosis. ASOs targeting SNORD14E inhibited tumor growth in the xenograft mouse model. TCGA-UCEC database showed that the proportion of patients with high expression of SNORD14E in middle-high risk and high-risk patients recommended by EMSO-ESGO-ESTRO guidelines for adjuvant therapy is more than 50%. Next, we enrolled 8 cases of high-risk and high-risk EC patients according to EMSO-ESGO-ESTRO guidelines and successfully constructed EC-PDOs. ASOs targeting SNORD14E inhibited the EC-PDO growth. Mechanistically, SNORD14E could recognize the mRNA of FOXM1 and recruit SRSF1 to promote the shearing of the variable exon VIIa of FOXM1, resulting in the overexpression of the FOXM1 malignant subtypes FOXM1b and FOXM1c. In addition, SNORD14E modified FOXM1 mRNA with 2`-O-methylation, which prolonged the half-life of FOXM1 mRNA. The nucleus accumulation of β-catenin caused by aberrant FOXM1 expression led to EC progression.

CONCLUSIONS

ASO targeting SNORD14E can be an effective treatment for EC.

Precision medicine using whole genome sequencing in a cat identifies a novel COL5A1 variant for classical Ehlers-Danlos syndrome

BACKGROUND

Ehlers-Danlos syndromes (EDS) are a heterogeneous group of heritable connective tissue disorders occurring in both human and veterinary patients. The genetics of these disorders are poorly described in small animal patients.

HYPOTHESIS/OBJECTIVES

Define the clinical manifestations and genetic cause of a suspected form of EDS in a cat.

ANIMALS

A 14-week-old male domestic medium hair cat was presented with skin hyperextensibility and fragility. The classic tragic facial expression was observed as well as chronic pruritus and mild hyperesthesia.

METHODS

Blood samples and a skin biopsy sample were collected from the affected cat. Clinical examinations, histology, electron microscopy and whole genome sequencing were conducted to characterize the clinical presentation and identify possible pathogenic DNA variants to support a diagnosis. Criteria defining variant pathogenicity were examined including human disease variant databases.

RESULTS

Histology showed sparse, disorganized collagen and an increase in cutaneous mast cells. Electron microscopy identified ultrastructural defects commonly seen in collagen type V alpha 1 chain (COL5A1) variants including flower-like collagen fibrils in cross-section. Whole genome sequencing and comparison with 413 cats in the 99 Lives Cat Genome Sequencing Consortium database identified a novel splice acceptor site variant at exon 4 in COL5A1 (c.501-2A>C).

CONCLUSIONS AND CLINICAL IMPORTANCE

Our report broadens the current understanding of EDS in veterinary patients and supports the use of precision medicine techniques in clinical veterinary practice. The classification of variants for pathogenicity should be considered in companion animals.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Deletion induced splicing in RIC3 drives nicotinic acetylcholine receptor regulation with implications for endoplasmic reticulum stress in human astrocytes

Nicotinic acetylcholine receptor (nAChR) dysregulation in astrocytes is reported in neurodegenerative disorders. Modulation of nAChRs through agonists confers protection to astrocytes from stress but regulation of chaperones involved in proteostasis with pathological implications is unclear. Resistance to inhibitors of cholinesterase 3 (RIC3), a potential chaperone of nAChRs is poorly studied in humans. We characterized RIC3 in astrocytes derived from an isogenic wild-type and Cas9 edited "del" human iPSC line harboring a 25 bp homozygous deletion in exon2. Altered RIC3 transcript ratio due to deletion induced splicing and an unexpected gain of α7nAChR expression were observed in "del" astrocytes. Transcriptome analysis showed higher expression of neurotransmitter/G-protein coupled receptors mediated by cAMP and calcium/calmodulin-dependent kinase signaling with increased cytokines/glutamate secretion. Functional implications examined using tunicamycin induced ER stress in wild-type astrocyte stress model showed cell cycle arrest, RIC3 upregulation, reduction in α7nAChR membrane levels but increased α4nAChR membrane expression. Conversely, tunicamycin-treated "del" astrocytes showed a comparatively higher α4nAChR membrane expression and upsurged cAMP signaling. Furthermore, reduced expression of stress markers CHOP, phospho-PERK and lowered XBP1 splicing in western blot and qPCR, validated by proteome-based pathway analysis indicated lowered disease severity. Findings indicate (i) a complex RNA regulatory mechanism via exonic deletion induced splicing; (ii) RIC-3 as a disordered protein having contrasting effects on co-expressed nAChR subtypes under basal/stress conditions; and (iii) RIC3 as a potential drug target against ER stress in astrocytes for neurodegenerative/nicotine-related brain disorders. Cellular rescue mechanism through deletion induced exon skipping may encourage ASO-based therapies for tauopathies.

Molecular mechanisms underlying the defects of two novel mutations in the HSD17B3 gene found in the Tunisian population

17β-hydroxysteroid dehydrogenase type 3 (17β-HSD3) converts Δ4-androstene-3,17-dione (androstenedione) to testosterone. It is expressed almost exclusively in the testes and is essential for appropriate male sexual development. More than 70 mutations in the HSD17B3 gene that cause 17β-HSD3 deficiency and result in 46,XY Disorders of Sex Development (46,XY DSD) have been reported. This study describes three novel Tunisian cases with mutations in HSD17B3. The first patient is homozygous for the previously reported mutation p.C206X. The inheritance of this mutation seemed to be independent of consanguineous marriage, which can be explained by its high frequency in the Tunisian population. The second patient has a novel splice site mutation in intron 6 at position c.490 -6 T > C. A splicing assay revealed a complete omission of exon 7 in the resulting HSD17B3 mRNA transcript. Skipping of exon 7 in HSD17B3 is predicted to cause a frame shift in exon 8 that affects the catalytic site and results in a truncation in exon 9, leading to an inactive enzyme. The third patient is homozygous for the novel missense mutation p.K202M, representing the first mutation identified in the catalytic tetrad of 17β-HSD3. Site-directed mutagenesis and enzyme activity measurements revealed a completely abolished 17β-HSD3 activity of the p.K202M mutant, despite unaffected protein expression, compared to the wild-type enzyme. Furthermore, the present study emphasizes the importance of genetic counselling, detabooization of 46,XY DSD, and a sensitization of the Tunisian population for the risks of consanguineous marriage.

Activated PI3K delta syndrome 1 mutations cause neutrophilia in zebrafish larvae

People with activated PI3 kinase delta syndrome 1 (APDS1) suffer from immune deficiency and severe bronchiectasis. APDS1 is caused by dominant activating mutations of the PIK3CD gene that encodes the PI3 kinase delta (PI3Kδ) catalytic subunit. Despite the importance of innate immunity defects in bronchiectasis, there has been limited investigation of neutrophils or macrophages in APDS1 patients or mouse models. Zebrafish embryos provide an ideal system to study neutrophils and macrophages. We used CRISPR-Cas9 and CRISPR-Cpf1, with oligonucleotide-directed homologous repair, to engineer zebrafish equivalents of the two most prevalent human APDS1 disease mutations. These zebrafish pik3cd alleles dominantly caused excessive neutrophilic inflammation in a tail-fin injury model. They also resulted in total body neutrophilia in the absence of any inflammatory stimulus but normal numbers of macrophages. Exposure of zebrafish to the PI3Kδ inhibitor CAL-101 reversed the total body neutrophilia. There was no apparent defect in neutrophil maturation or migration, and tail-fin regeneration was unimpaired. Overall, the finding is of enhanced granulopoeisis, in the absence of notable phenotypic change in neutrophils and macrophages.

Polymorphism of glutathione S-transferase P1 of dogs with mammary tumours

Mammary tumours constitute more than half of neoplasms in female dogs from different countries. Genome sequences are associated with cancer susceptibility but there is little information available about genetic polymorphisms of glutathione S-transferase P1 (GSTP1) in canine cancers. The aim of this study was to find single nucleotide polymorphisms (SNPs) in GSTP1 of dogs (Canis lupus familiaris) with mammary tumours compared to healthy dogs and to determine the association between GSTP1 polymorphisms and the occurrence of these tumours. The study population included 36 client-owned female dogs with mammary tumours and 12 healthy female dogs, with no previous diagnosis of cancer. DNA was extracted from blood and amplified by PCR assay. PCR-products were sequenced by Sanger method and analysed manually. The 33 polymorphisms were found in GSTP1: 1 coding SNP (exon 4), 24 non-coding SNPs (9 in exon 1), 7 deletions and 1 insertion. The 17 polymorphisms have been found in introns 1, 4, 5 and 6. The dogs with mammary tumours have significant difference from healthy in SNPs I4 c.1018 + 123 T > C (OR 13.412, 95%CI 1.574-114.267, P = .001), I5 c.1487 + 27 T > C (OR 10.737, 95%CI 1.260-91.477, P = .004), I5 c.1487 + 842 G > C (OR 4.714, 95% CI 1.086-20.472, P = .046) and I6 c.2481 + 50 A > G (OR 12.000, 95% CI 1.409-102.207, P = .002). SNP E5 c.1487 T > C and I5 c.1487 + 829 delG also differed significantly (P = .03) but not to the confidence interval. The study, for the first time, showed a positive association of SNPs in GSTP1 with mammary tumours of dogs, that can possibly be used to predict the occurrence of this pathology.

Effects of splicing-regulatory polymorphisms in ABCC2, ABCG2, and ABCB1 on methotrexate exposure in Chinese children with acute lymphoblastic leukemia

PURPOSE

Adenosine triphosphate (ATP)-binding cassette (ABC) transporters play an important role in the response to methotrexate (MTX). In this study, we investigated the frequency distribution of three splicing-regulatory polymorphisms in ABC transporters (ABCC2 rs2273697 G>A, ABCG2 rs2231142 G>T, and ABCB1 rs1128503 A>G) and their effects on MTX concentrations and the clinical outcome in a Chinese pediatric population with acute lymphoblastic leukemia (ALL).

METHODS

A fluorescence polarization immunoassay was used to measure the serum MTX concentrations in 24 h (C_24h) and 42 h (C_42h). The Sequenom MassARRAY system was used for single-nucleotide polymorphism (SNP) genotyping.

RESULTS

The study population had significantly lower frequencies of ABCC2 rs2273697 A, ABCG2 rs2231142 G, and ABCB1 rs1128503 G than African and European samples (P < 0.05). The dose-normalized MTX concentrations after 24 h and the proportion of C_42h > 0.5 µmol/L were significantly lower in patients with the ABCG2 rs2231142 GG genotype than in patients with the GT or TT genotype (P = 0.01 and 0.006, respectively). No significant effects on MTX pharmacokinetics were observed for ABCC2 rs2273697 and ABCB1 rs1128503 polymorphisms. Bioinformatics analysis suggested that the three SNPs overlapped with the putative binding sites of several splicing factors.

CONCLUSION

In conclusion, our study confirmed the ethnicity-based differences in the distribution of the three investigated SNPs. The ABCG2 rs2231142 polymorphism exerted a significant effect on the level of MTX exposure. These findings may help explain the variability in MTX responses and optimize MTX treatment in pediatric patients with ALL.

Correction of the Splicing Defect Caused by a Recurrent Variant in ABCA4 (c.769-784C>T) That Underlies Stargardt Disease

Stargardt disease is an inherited retinal disease caused by biallelic mutations in the ABCA4 gene, many of which affect ABCA4 splicing. In this study, nine antisense oligonucleotides (AONs) were designed to correct pseudoexon (PE) inclusion caused by a recurrent deep-intronic variant in ABCA4 (c.769-784C>T). First, the ability of AONs to skip the PE from the final ABCA4 mRNA transcript was assessed in two cellular models carrying the c.769-784C>T variant: a midigene assay using HEK293T cells and patient-derived fibroblasts. Based on the splicing-correcting ability of each individual AON, the three most efficacious AONs targeting independent regions of the PE were selected for a final assessment in photoreceptor precursor cells (PPCs). The final analysis in the PPC model confirmed high efficacy of AON2, -5, and -7 in promoting PE exclusion. Among the three AONs, AON2 is chosen as the lead candidate for further optimization, hereby showcasing the high potential of AONs to correct aberrant splicing events driven by deep-intronic variants.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

Background

Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.

Results

Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in “cryptic” splicing regions beyond canonical donor or acceptor splice sites.

Conclusion

Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

Diagnostic reasoning in neurogenetics: a general approach

Establishing the definitive diagnosis of a neurogenetic disease is usually a complex task. However, like any type of clinical diagnostic reasoning, an organized process of development and consideration of diagnostic hypotheses may guide neurologists and medical geneticists to solve this difficult task. The aim of the present review is to propose a general method for diagnostic reasoning in neurogenetics, with the definition of the main neurological syndrome and its associated topographical diagnosis, followed by the identification of major and secondary neurological syndromes, extraneurological findings, and inheritance pattern. We also discuss general rules and knowledge requirements of the ordering physician to request genetic testing and information on how to interpret genetic variants in a genetic report. By guiding the requests for genetic testing according to an organized model of diagnostic reasoning and with the availability of specific treatments, clinicians may find greater resoluteness and efficacy in the diagnostic investigation, shortening the struggle of patients for a definitive diagnosis.

Alternative splicing of HSPA12A pre-RNA by SRSF11 contributes to metastasis potential of colorectal cancer

BACKGROUND

Dysregulation of alternative splicing (AS) induced by serine/arginine-rich proteins has recently been linked to cancer metastasis. Nonetheless, as a member of the serine/arginine-rich protein family, the involvement of SRSF11 in colorectal cancer (CRC) is unknown.

METHODS

The TCGA dataset and clinical samples were used to assess SRSF11 expression levels in CRC. For SRSF11, functional experiments were conducted both in vitro and in vivo. RNA-seq technology was used to analyze and screen SRSF11-triggered AS events, which were then confirmed by in vivo UV crosslinking and immunoprecipitation (CLIP) and mini-gene reporter assays. Jalview software was used to determine the preferential binding motif with relation to exon skipping (ES) events. Furthermore, coimmunoprecipitation (Co-IP) and Phospho-tag SDS-PAGE experiments were used to investigate PAK5-mediated phosphorylation regulation on SRSF11, and in vitro kinase experiments validated the interaction.

RESULTS

In CRC, SRSF11 was discovered to be overexpressed and associated with a poor prognosis. And SRSF11 played a pro-metastatic role in vitro and in vivo. By screening SRSF11-regulated AS events, we identified the binding motif of SRSF11-triggered splicing-switching of HSPA12A AS, which specifically regulated HSPA12A AS by directly binding to a motif in exon 2. Mechanistically, the HSPA12A transcript with exon 2 retention increased N-cadherin expression by promoting RNA stability. Furthermore, the oncogenic kinase PAK5 phosphorylated SRSF11 at serine 287, protecting it from ubiquitination degradation.

CONCLUSIONS

SRSF11 exerts pro-metastatic effects in CRC by inhibiting the AS of HSPA12A pre-RNA. Our findings point to SRSF11-regulated HSPA12A splicing as a novel relationship between SRSF11-regulated splicing and CRC metastasis and suggest a PAK5/SRSF11/HSPA12A axis as a potential therapeutic target and prognostic biomarker in CRC.

HvGBSSI mutation at the splicing receptor site affected RNA splicing and decreased amylose content in barley

Granule-bound starch synthase I (HvGBSSI) is encoded by the barley waxy (Wx-1) gene and is the sole enzyme in the synthesis of amylose. Here, a Wx-1 mutant was identified from an ethyl methane sulfonate (EMS)-mutagenized barley population. There were two single-base mutations G1086A and A2424G in Wx-1 in the mutant (M2-1105). The G1086A mutation is located at the 3' splicing receptor (AG) site of the fourth intron, resulting in an abnormal RNA splicing. The A2424G mutation was a synonymous mutation in the ninth intron. The pre-mRNA of Wx-1 was incorrectly spliced and transcribed into two abnormal transcripts. The type I transcript had a 6 bp deletion in the 5' of fifth exon, leading to a translated HvGBSSI protein lacking two amino acids with a decreased starch-binding capacity. In the type II transcript, the fourth intron was incorrectly cleaved and retained, resulting in the premature termination of the barley Wx-1 gene. The mutations in the Wx-1 decreased the enzymatic activity of the HvGBSSI enzyme and resulted in a decreased level in amylose content. This work sheds light on a new Wx-1 gene inaction mechanism.

lncRNA260 siRNA Accelerates M2 Macrophage Polarization and Alleviates Oxidative Stress via Inhibiting IL28RA Gene Alternative Splicing

The macrophage transformation of inflammatory M1 to anti-inflammatory M2 could be promoted by activating PI3K/AKT signaling pathway. In our previous study, it was found that downregulation of lncRNA260 could ameliorate hypoxic cardiomyocyte injury by regulating IL28RA through the activation of PI3K/AKT signaling pathways. It was suggested that lncRNA260 siRNA could promote the macrophages toward M2 polarization by regulating IL28RA. In this study, lncRNA260 siRNA was used to observe its effect on the polarization of murine bone marrow-derived macrophages (BMDM) and investigate its related mechanisms. lncRNA 260 specific siRNA were designed and synthesized which were transfected into murine BMDM with liposomes. The experiment was divided into three groups: Hypoxia group, Hypoxia+lncRNA 260-specific siRNA transfection group, and Normoxia group. The CD206-APC/CD11b-FITC or CD206-FITC/CD107b (Mac-3) double positive proportions were used to compare the M2 polarization proportions in the hypoxia process by using the immunofluorescence staining method. The p-AKT, Arg 1, PI3KCG, IL28RAV1, and IL28RAV2 protein expression changes were observed by using the western blot method. Compared with the Normoxia group, the M2 proportions were significantly decreased in the Hypoxia group (P < 0.05). Compared with the hypoxia group, the M2 proportions were significantly increased in the Hypoxia+lncRNA260 siRNA transfection group (P < 0.05). In the Hypoxia group, the ratios of Arg 1/β-Actin, p-AKT/β-Actin, PI3KCG/β-Actin, and IL28RAV1/β-Actin were significantly lower than those in the Normoxia group (P < 0.05). After transfection with lncRNA260 siRNA, the ratios of Arg1/β-Actin, p-AKT/β-Actin, PI3KCG/β-Actin, and IL28RAV1/β-Actin were significantly higher than those in the Hypoxia group (P < 0.05). Compared with the Normoxia group, the IL28RAV2/β-Actin in the Hypoxia group was significantly increased (P < 0.05). After transfection with lncRNA260 siRNA, the ratio of IL28RAV2/β-Actin was significantly decreased than that in the Hypoxia group (P < 0.05). lncRNA260 siRNA could promote the M2 polarization of the hypoxia macrophages by reducing the IL28RAV2 alternative splicing variant, which might be related to the activation of the JAK-STAT and PI3K/AKT signaling pathways. It will provide a new strategy for the anti-inflammation, antioxidative stress therapy, and cardiac remodeling after AMI.

Paediatric BCOR ‐associated sarcomas with a novel long spliced internal tandem duplication of BCOR exon 15

Clear cell sarcoma of the kidney (CCSK) and primitive myxoid mesenchymal tumour of infancy (PMMTI) are paediatric sarcomas that most commonly harbour internal tandem duplications (ITDs) of exon 15 of the BCOR gene, in the range of 87–114 base pairs (bp). Some cases, instead, have BCOR‐CCNB3 or YWHAE‐NUTM2 gene fusions. About 10% of cases lack any of these genetic alterations when tested by standard methods. Two cases of CCSK and one PMMTI lacking the aforementioned mutations were analysed using Archer FusionPlex technology. Two related BCOR exon 15 RNA transcripts with ITDs of lengths 388 and 96 bp were detected in each case; only the 388 bp transcript was identified when genomic DNA was sequenced. In silico analysis of this transcript revealed acceptor and donor splice sites indicating that, at the RNA level, the 388‐bp transcript was likely spliced to form the 96‐bp transcript. The results were confirmed by Sanger sequencing using primers targeting the ITD breakpoint. This novel and unusually long ITD segment is difficult to identify by DNA sequencing using typical primer design strategies flanking entire duplicated segments because it exceeds the typical read lengths of most sequencing platforms as well as the usual fragment lengths obtained from formalin‐fixed paraffin‐embedded material. As diagnosis of CCSK and PMMTI may be challenging by morphology and immunohistochemistry alone, it is important to identify mutations in these cases. Knowledge of this novel BCOR ITD is important in relation to primer design for detection by sequencing, and using RNA versus DNA for sequencing.

A Study of Combined Genotype Effects of SHCBP1 on Wool Quality Traits in Chinese Merino

SHCBP1 (Shc SH2-domain binding protein 1) is a member of the Src and collagen homolog (Shc) protein family and is closely associated with multiple signaling pathways that play important roles during hair follicle induction, morphogenesis, and cycling. The purpose of this study was to investigate SHCBP1 gene expression, polymorphisms, and the association between SHCBP1 and wool quality traits in Chinese Merino sheep. The SHCBP1 gene was shown, by qPCR, to be ubiquitously expressed in sheep tissues and differentially expressed in the adult skin of Chinese Merino and Suffolk sheep. Four SNPs (termed SHCBP1SNPs 1-4) were identified by Sanger sequencing and were located in exon 2, intron 9, intron 12, and exon 13 of the sheep SHCBP1 gene, respectively. SHCBP1SNPs 3 and 4 (rs411176240 and rs160910635) were significantly associated with wool crimp (P < 0.05). The combined polymorphism (SHCBP1SNP3-SHCBP1SNP4) was significantly associated with wool crimp (P < 0.05). Bioinformatics analysis showed that the SNPs associated with wool crimp (SHCBP1SNPs 3 and 4) might affect the pre-mRNA splicing by creating binding sites for serine-arginine-rich proteins and that SHCBP1SNP4 might alter the SHCBP1 mRNA and protein secondary structure. Our results suggest that SHCBP1 influences wool crimp and SHCBP1SNPs 3 and 4 might be useful markers for marker-assisted selection and sheep breeding.

Alternatively spliced CSF3R isoforms in SRSF2 P95H mutated myeloid neoplasms

Alternatively spliced colony stimulating factor 3 receptor (CSF3R) isoforms Class III and Class IV are observed in myelodysplastic syndromes (MDS), but their roles in disease remain unclear. We report that the MDS-associated splicing factor SRSF2 affects the expression of Class III and Class IV isoforms and perturbs granulopoiesis. Add-back of the Class IV isoform in Csf3r-null mouse progenitor cells increased granulocyte progenitors with impaired neutrophil differentiation, while add-back of the Class III produced dysmorphic neutrophils in fewer numbers. These CSF3R isoforms were elevated in patients with myeloid neoplasms harboring SRSF2 mutations. Using in vitro splicing assays, we confirmed increased Class III and Class IV transcripts when SRSF2 P95 mutations were co-expressed with the CSF3R minigene in K562 cells. Since SRSF2 regulates splicing partly by recognizing exonic splicing enhancer (ESE) sequences on pre-mRNA, deletion of either ESE motifs within CSF3R exon 17 decreased Class IV transcript levels without affecting Class III. CD34+ cells expressing SRSF2 P95H showed impaired neutrophil differentiation in response to G-CSF and was accompanied by increased levels of Class IV. Our findings suggest that SRSF2 P95H promotes Class IV splicing by binding to key ESE sequences in CSF3R exon 17, and that SRSF2, when mutated, contributes to dysgranulopoiesis.

A novel missense variant in endothelin-2 (EDN2) causes a growth and respiratory lethal syndrome in bovine

The high level of fragmentation of the Spanish Lidia cattle breed, divided into lineages called ‘castas’ and into herds within lineages based on reproductive isolation, increases the risk of homozygosity and the outbreak of recessive genetic defects. Since 2004, an increasing number of calves have been identified in a Lidia herd with signs of severe growth retardation, respiratory alterations and juvenile lethality, which constitutes a novel inherited syndrome in cattle and was subsequently termed growth and respiratory lethal syndrome. We performed a genome‐wide association study on a cohort of 13 affected calves and 24 putative non‐carrier parents, mapping the disease to a wide 6 cM region on bovine chromosome 3 (p < 10⁻⁷). Whole genome re‐sequencing of three affected calves and three putative non‐carrier parents identified a novel missense variant (c.149G>A|p.Cys50Tyr) in exon 2 of the endothelin 2 (EDN2) gene. Bioinformatic analyses of p.Cys50Tyr effects predicted them to be damaging for both the structure and the function of the edn2 protein, and to create a new site of splicing that may also affect the pattern of pre‐mRNA splicing and exon definition. Sanger sequencing of this variant on the rest of the sample set confirmed the segregation pattern obtained with whole genome re‐sequencing. The identification of the causative variant and the development of a diagnostic genetic test enable the efficient design of matings to keep the effective population size as high as possible, as well as providing insights into the first EDN2‐associated hereditary disease in cattle or other species.

A novel leaky splice variant in centromere protein J (CENPJ)-associated Seckel syndrome

Primary microcephaly and Seckel syndrome are rare genetically and clinically heterogenous brain development disorders. Several exonic/splicing mutations are reported for these disorders to date, but ∼40% of all cases remain unexplained. We aimed to uncover the genetic correlate(s) in a family of multiple siblings with microcephaly. A novel homozygous intronic variant (NC_000013.10:g.25459823T>C) in CENPJ (13q12) segregating with all four affected male siblings was identified by exome sequencing and validated by targeted linkage approach (logarithm of the odds score 1.8 at θ 0.0). RT-PCR of CENPJ in affected siblings using their EBV derived cell lines showed aberrant transcripts suggestive of exon skipping confirmed by Sanger sequencing. Significantly reduced wild type transcript/protein in the affected siblings having the splice variant indicates a leaky gene expression of pathological relevance. Based on known CENPJ function, assessing for mitotic alterations revealed defect in centrosome duplication causing mono/multicentrosome(s) at prophase, delayed metaphase, and unequal chromosomal segregation in patient cells. Clinical features witnessed in this study expand the spectrum of CENPJ-associated primary microcephaly and Seckel syndrome. Furthermore, besides the importance of regulatory variants in classical monogenic disorders these findings provide new insights into splice site biology with possible implications for ASO-based therapies.

Antisense Oligonucleotide Induction of the hnRNPA1b Isoform Affects Pre-mRNA Splicing of SMN2 in SMA Type I Fibroblasts

Spinal muscular atrophy (SMA) is a severe, debilitating neuromuscular condition characterised by loss of motor neurons and progressive muscle wasting. SMA is caused by a loss of expression of SMN1 that encodes the survival motor neuron (SMN) protein necessary for the survival of motor neurons. Restoration of SMN expression through increased inclusion of SMN2 exon 7 is known to ameliorate symptoms in SMA patients. As a consequence, regulation of pre-mRNA splicing of SMN2 could provide a potential molecular therapy for SMA. In this study, we explored if splice switching antisense oligonucleotides could redirect the splicing repressor hnRNPA1 to the hnRNPA1b isoform and restore SMN expression in fibroblasts from a type I SMA patient. Antisense oligonucleotides (AOs) were designed to promote exon 7b retention in the mature mRNA and induce the hnRNPA1b isoform. RT-PCR and western blot analysis were used to assess and monitor the efficiency of different AO combinations. A combination of AOs targeting multiple silencing motifs in hnRNPA1 pre-mRNA led to robust hnRNPA1b induction, which, in turn, significantly increased expression of full-length SMN (FL-SMN) protein. A combination of PMOs targeting the same motifs also strongly induced hnRNPA1b isoform, but surprisingly SMN2 exon 5 skipping was detected, and the PMO cocktail did not lead to a significant increase in expression of FL-SMN protein. We further performed RNA sequencing to assess the genome-wide effects of hnRNPA1b induction. Some 3244 genes were differentially expressed between the hnRNPA1b-induced and untreated SMA fibroblasts, which are functionally enriched in cell cycle and chromosome segregation processes. RT-PCR analysis demonstrated that expression of the master regulator of these enrichment pathways, MYBL2 and FOXM1B, were reduced in response to PMO treatment. These findings suggested that induction of hnRNPA1b can promote SMN protein expression, but not at sufficient levels to be clinically relevant.

Silent mutations reveal therapeutic vulnerability in RAS Q61 cancers

RAS family members are the most frequently mutated oncogenes in human cancers. Although KRAS(G12C)-specific inhibitors show clinical activity in patients with cancer^1-3, there are no direct inhibitors of NRAS, HRAS or non-G12C KRAS variants. Here we uncover the requirement of the silent KRAS^G60G mutation for cells to produce a functional KRAS(Q61K). In the absence of this G60G mutation in KRAS^Q61K, a cryptic splice donor site is formed, promoting alternative splicing and premature protein termination. A G60G silent mutation eliminates the splice donor site, yielding a functional KRAS(Q61K) variant. We detected a concordance of KRAS^Q61K and a G60G/A59A silent mutation in three independent pan-cancer cohorts. The region around RAS Q61 is enriched in exonic splicing enhancer (ESE) motifs and we designed mutant-specific oligonucleotides to interfere with ESE-mediated splicing, rendering the RAS(Q61) protein non-functional in a mutant-selective manner. The induction of aberrant splicing by antisense oligonucleotides demonstrated therapeutic effects in vitro and in vivo. By studying the splicing necessary for a functional KRAS(Q61K), we uncover a mutant-selective treatment strategy for RAS^Q61 cancer and expose a mutant-specific vulnerability, which could potentially be exploited for therapy in other genetic contexts.

Cytoprotective Activity of Polyamines Is Associated with the Alternative Splicing of RAD51A Pre-mRNA in Normal Human CD4+ T Lymphocytes

Physiological polyamines are ubiquitous polycations with pleiotropic biochemical activities, including regulation of gene expression and cell proliferation as well as modulation of cell signaling. They can also decrease DNA damage and promote cell survival. In the present study, we demonstrated that polyamines have cytoprotective effects on normal human CD4+ T lymphocytes but not on cancer Jurkat or K562 cells. Pretreatment of lymphocytes with polyamines resulted in a significant reduction in cells with DNA damage induced by doxorubicin, cisplatin, or irinotecan, leading to an increase in cell survival and viability. The induction of RAD51A expression was in response to DNA damage in both cancer and normal cells. However, in normal cells, putrescin pretreatment resulted in alternative splicing of RAD51A and the switch of the predominant expression from the splice variant with the deletion of exon 4 to the full-length variant. Induction of RAD51A alternative splicing by splice-switching oligonucleotides resulted in a decrease in DNA damage and cell protection against cisplatin-induced apoptosis. The results of this study suggest that the cytoprotective activity of polyamines is associated with the alternative splicing of RAD51A pre-mRNA in normal human CD4+ T lymphocytes. The difference in the sensitivity of normal and cancer cells to polyamines may become the basis for the use of these compounds to protect normal lymphocytes during lymphoblastic chemotherapy.

PDE4B Proposed as a High Myopia Susceptibility Gene in Chinese Population

Myopia is the most common cause of refractive error worldwide. High myopia is a severe type of myopia, which usually accompanies pathological changes in the fundus. To identify high myopia susceptibility genes, DNA-pooling based genome-wide association analysis was used to search for a correlation between single nucleotide polymorphisms and high myopia in a Han Chinese cohort (cases vs. controls in discovery stage: 507 vs. 294; replication stage 1: 991 vs. 1,025; replication stage 2: 1,021 vs. 52,708). Three variants (rs10889602T/G, rs2193015T/C, rs9676191A/C) were identified as being significantly associated with high myopia in the discovery, and replication stage. rs10889602T/G is located at the third intron of phosphodiesterase 4B (PDE4B), whose functional assays were performed by comparing the effects of rs10889602T/T deletion of this risk allele on PDE4B and COL1A1 gene and protein expression levels in the rs10889602T/T^del/del, rs10889602T/T^del/wt, and normal control A549 cell lines. The declines in the PDE4B and COL1A1 gene expression levels were larger in the rs10889602T/T deleted A549 cells than in the normal control A549 cells (one-way ANOVA, p < 0.001). The knockdown of PDE4B by siRNA in human scleral fibroblasts led to downregulation of COL1A1. This correspondence between the declines in rs10889602 of the PDE4B gene, PDE4B knockdown, and COL1A1 protein expression levels suggest that PDE4B may be a novel high myopia susceptibility gene, which regulates myopia progression through controlling scleral collagen I expression levels. More studies are needed to determine if there is a correlation between PDE4B and high myopia in other larger sample sized cohorts.

Associations of MCM8 rs3761873 and rs16991617 variants with abnormal uterine bleeding induced by copper intrauterine device

AIM

Intrauterine device (IUD) is a commonly used contraceptive method worldwide. Abnormal uterine bleeding (AUB) is one of the most common side effects of Cu-IUDs. Since AUB varies among Cu-IUD users, changes in the bleeding-related genetic factors may contribute to AUB. This study aimed to determine the genetic risk factors of AUB after Cu-IUD insertion.

METHODS

We conducted a case-control study on women who experienced AUB after Cu-IUD insertion (case:control = 62:59). Six candidate variants were genotyped using the Sequenom MassARRAY. Genotype and allele frequencies were analyzed using SHEsisPlus. We performed Pearson's Chi-squared test to analyze categorical data, and ESEfinder to predict the impact on splicing regulation.

RESULTS

MCM8 coding sequence variants: rs3761873-A>C was in Exon 7 and rs16991617 A>G was in Exon 12 of all 19 exons, both of which were significantly different between cases and controls (p_allele  = 0.039 and p_genotype  = 0.092). rs6022 and rs6029 in F5 gene and rs3761873 and rs16991617 in the MCM8 gene showed strong linkage disequilibrium (R²  > 0.8). ESEfinder indicated that the variants of MCM8 may affect the splicing regulation.

CONCLUSIONS

MCM8 rs376187 and rs16991617 were associated with AUB in Cu-IUDs users. MCM8 may play a role in AUB by regulating functions of reproductive organs and primary ovarian insufficiency. Our findings may improve the understanding of the genetic basis of AUB caused by Cu-IUDs.

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

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

Cas9-specific immune responses compromise local and systemic AAV CRISPR therapy in multiple dystrophic canine models

Adeno-associated virus (AAV)-mediated CRISPR-Cas9 editing holds promise to treat many diseases. The immune response to bacterial-derived Cas9 has been speculated as a hurdle for AAV-CRISPR therapy. However, immunological consequences of AAV-mediated Cas9 expression have thus far not been thoroughly investigated in large mammals. We evaluate Cas9-specific immune responses in canine models of Duchenne muscular dystrophy (DMD) following intramuscular and intravenous AAV-CRISPR therapy. Treatment results initially in robust dystrophin restoration in affected dogs but also induces muscle inflammation, and Cas9-specific humoral and cytotoxic T-lymphocyte (CTL) responses that are not prevented by the muscle-specific promoter and transient prednisolone immune suppression. In normal dogs, AAV-mediated Cas9 expression induces similar, though milder, immune responses. In contrast, other therapeutic (micro-dystrophin and SERCA2a) and reporter (alkaline phosphatase, AP) vectors result in persistent expression without inducing muscle inflammation. Our results suggest Cas9 immunity may represent a critical barrier for AAV-CRISPR therapy in large mammals.

Functional pre-therapeutic evaluation by genome editing of variants of uncertain significance of essential tumor suppressor genes

BACKGROUND

Targeted therapies in oncology are promising but variants of uncertain significance (VUS) limit their use for clinical management and necessitate functional testing in vitro. Using BRCA1 and BRCA2 variants, which have consequences on PARP inhibitor sensitivity, and POLE variants, potential biomarkers of immunotherapy response, we developed a rapid functional assay based on CRISPR-Cas9 genome editing to determine the functional consequences of these variants having potentially direct implications on patients' access to targeted therapies.

METHODS

We first evaluated the functional impact of 26 BRCA1 and 7 BRCA2 variants by editing and comparing NGS results between the variant of interest and a silent control variant. Ten of these variants had already been classified as benign or pathogenic and were used as controls. Finally, we extended this method to the characterization of POLE VUS.

RESULTS

For the 23 variants that were unclassified or for which conflicting interpretations had been reported, 15 were classified as functionally normal and 6 as functionally abnormal. Another two variants were found to have intermediate consequences, both with potential impacts on splicing. We then compared these scores to the patients' responses to PARP inhibitors when possible. Finally, to prove the application of our method to the classification of variants from other tumor suppressor genes, we exemplified with three POLE VUS. Among them, two were classified with an intermediate functional impact and one was functionally abnormal. Eventually, four POLE variants previously classified in databases were also evaluated. However, we found evidence of a discordance with the classification, variant p.Leu424Val being found here functionally normal.

CONCLUSIONS

Our new rapid functional assay can be used to characterize the functional implication of BRCA1 and BRCA2 variants, giving patients whose variants were evaluated as functionally abnormal access to PARP inhibitor treatment. Retrospective analysis of patients' responses to PARP inhibitors, where accessible, was consistent with our functional score evaluation and confirmed the accuracy of our protocol. This method could potentially be extended to the classification of VUS from all essential tumor suppressor genes and can be performed within a timeframe compatible with clinical applications, thereby having a direct theranostic impact.

Nucleotides in both donor and acceptor splice sites are responsible for choice in NAGNAG tandem splice sites

Among alternative splicing events in the human transcriptome, tandem NAGNAG acceptor splice sites represent an appreciable proportion. Both proximal and distal NAG can be used to produce two splicing isoforms differing by three nucleotides. In some cases, the upstream exon can be alternatively spliced as well, which further increases the number of possible transcripts. In this study, we showed that NAG choice in tandem splice site depends considerably not only on the concerned acceptor, but also on the upstream donor splice site sequence. Using an extensive set of experiments with systematically modified two-exonic minigene systems of AFAP1L2 or CSTD gene, we recognized the third and fifth intronic upstream donor splice site position and the tandem acceptor splice site region spanning from -10 to +2, including NAGNAG itself, as the main drivers. In addition, competition between different branch points and their composition were also shown to play a significant role in NAG choice. All these nucleotide effects appeared almost additive, which explained the high variability in proximal versus distal NAG usage.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

A new non-aggregative splicing isoform of human Tau is decreased in Alzheimer's disease

Tauopathies, including Alzheimer's disease (AD) and frontotemporal lobar degeneration with Tau pathology (FTLD-tau), are a group of neurodegenerative disorders characterized by Tau hyperphosphorylation. Post-translational modifications of Tau such as phosphorylation and truncation have been demonstrated to be an essential step in the molecular pathogenesis of these tauopathies. In this work, we demonstrate the existence of a new, human-specific truncated form of Tau generated by intron 12 retention in human neuroblastoma cells and, to a higher extent, in human RNA brain samples, using qPCR and further confirming the results on a larger database of human RNA-seq samples. Diminished protein levels of this new Tau isoform are found by Westernblotting in Alzheimer's patients' brains (Braak I n = 3; Braak II n = 6, Braak III n = 3, Braak IV n = 1, and Braak V n = 10, Braak VI n = 8) with respect to non-demented control subjects (n = 9), suggesting that the lack of this truncated isoform may play an important role in the pathology. This new Tau isoform exhibits similar post-transcriptional modifications by phosphorylation and affinity for microtubule binding, but more interestingly, is less prone to aggregate than other Tau isoforms. Finally, we present evidence suggesting this new Tau isoform could be linked to the inhibition of GSK3β, which would mediate intron 12 retention by modulating the serine/arginine rich splicing factor 2 (SRSF2). Our results show the existence of an important new isoform of Tau and suggest that further research on this less aggregation-prone Tau may help to develop future therapies for Alzheimer's disease and other tauopathies.

Diversity in STK4 Deficiency and Review of the Literature

BACKGROUND

Serine-threonine kinase-4 (STK4) deficiency is an autosomal recessive (AR) combined immunodeficiency (CID).

OBJECTIVE

We aimed to define characteristic clinical and laboratory features to aid the differential diagnosis and determine the most suitable therapy.

METHODS

In addition to nine patients diagnosed, we reviewed 15 patients from medical literature. We compared B lymphocyte subgroups of our cohort with age-matched healthy controls.

RESULTS

In our cohort, the median age at symptom onset and age of diagnosis are 6years-8months (mo)(6-248mo) and 7years-5mo (6-260mo), respectively. The main clinical findings were infections (9/9), autoimmune/inflammatory diseases (7/9), and atopy (4/9). CD4 lymphopenia (9/9), lymphopenia (7/9), intermittent eosinophilia (4/9), transient neutropenia (3/9), low immunoglobulin (Ig) M (4/9), and high IgE (4/9) were common. Decreased recent thymic emigrants, naive and central memory T cells, albeit increased effector memory T cells were present. The increase in plasmablasts (p=0.003) and the decrease in switched memory B cells (p=0.022) were significant. Out of a total of 24 patients, cutaneous viral infections (n=20), recurrent pneumonia (n=18), Epstein Barr Virus (EBV)-associated lymphoproliferation (n=11), atopic dermatitis (n=10), autoimmune cytopenia (n=7), and lymphoma (n=6) were frequently seen. Lymphopenia, CD4 lymphopenia, high Ig G, A, and E were the most common laboratory characteristics.

CONCLUSION

The differential diagnosis with AR-hyperimmunoglobulin E syndrome is crucial as atopy and CD4 lymphopenia are prominent in both diseases. Immunoglobulins and antibacterial/antiviral prophylaxis are the mainstays of treatment. Clinicians may use immunomodulatory therapies during inflammatory/autoimmune complications. However, more data is needed to recommend hematopoietic stem cell transplantation (HSCT) as a safe therapy.

Characterization of Rheumatoid Arthritis Risk-Associated SNPs and Identification of Novel Therapeutic Sites Using an In-Silico Approach

Simple Summary

Rheumatoid arthritis (RA) is a complex disease resulting from multiple genetic and environmental pathogenic factors. The genetic factors include single-nucleotide polymorphisms (SNPs), which have been reported to be associated with RA, but their specific role in the pathogenesis of RA remains unexplained. This study explains the potential role of RA risk-associated SNPs in its pathogenesis in order to provide a basis for understanding the genetic complexity of RA. Several roles of these SNPs are described in this study, and may also aid in the design of a therapeutic strategy for RA. Furthermore, novel potential therapeutic sites have also been researched, resulting in the identification of three novel therapeutic targets. The therapeutic strategies for the treatment of RA include inflammatory pathway-targeting drugs, which alleviate inflammation in joints. There is always a need for novel therapeutic targets that can play a role in alleviating inflammation in autoimmune diseases including RA. Therefore, these novel therapeutic sites are very important, and further experimental studies are required.

Abstract

Single-nucleotide polymorphisms (SNPs) are reported to be associated with many diseases, including autoimmune diseases. In rheumatoid arthritis (RA), about 152 SNPs are reported to account for ~15% of its heritability. These SNPs may result in the alteration of gene expression and may also affect the stability of mRNA, resulting in diseased protein. Therefore, in order to predict the underlying mechanism of these SNPs and identify novel therapeutic sites for the treatment of RA, several bioinformatics tools were used. The damaging effect of 23 non-synonymous SNPs on proteins using different tools suggested four SNPs, including rs2476601 in PTPN22, rs5029941 and rs2230926 in TNFAIP3, and rs34536443 in TYK2, to be the most damaging. In total, 42 of 76 RA-associated intronic SNPs were predicted to create or abolish potential splice sites. Moreover, the analysis of 11 RA-associated UTR SNPs indicated that only one SNP, rs1128334, located in 3′UTR of ETS1, caused functional pattern changes in BRD-BOX. For the identification of novel therapeutics sites to treat RA, extensive gene–gene interaction network interactive pathways were established, with the identification of 13 potential target sites for the development of RA drugs, including three novel target genes. The anticipated effect of these findings on RA pathogenesis may be further validated in both in vivo and in vitro studies.

SRSF1 promotes the inclusion of exon 3 of SRA1 and the invasion of hepatocellular carcinoma cells by interacting with exon 3 of SRA1pre-mRNA

Steroid receptor RNA activator 1 (SRA1) has been described as a novel transcriptional co-activator that affects the migration of cancer cells. Through RT-PCR, we identified that skipping exon 3 of SRA1 produces two isoforms, including the truncated short isoform, SRA1-S, and the long isoform, SRA1-L. However, the effect of these two isomers on the migration of HCC cells, as well as the specific mechanism of exon 3 skipping remain unclear. In this study, we found up regulated expression of SRSF1 and SRA1-L in highly metastatic HCCLM3, as well as in HCCs with SRSF1 demonstrating the strongest correlation with SRA1-L. In contrast, we observed a constitutively low expression of SRA1-S and SRSF1 in lowly metastatic HepG2 cells. Overexpression of SRSF1 or SRA1-L promoted migration and invasion by increasing the expression of CD44, while SRA1-S reversed the effect of SRSF1 and SRA1-L in vitro. In addition, lung metastasis in mice revealed that, knockdown of SRSF1 or SRA1-L inhibited the migration of HCC cells, while SRA1-L overexpression abolished the effect of SRSF1 knockout and instead promoted HCC cells migration in vivo. More importantly, RNA immunoprecipitation and Cross-link immunoprecipitation analyses showed that SRSF1 interacts with exon 3 of SRA1 to up regulate the expression of SRA1-L in HCC cells. RNA pull-down results indicated that SRSF1 could also bind to exon 3 of SRA1 in vitro. Finally, minigene -MS2 mutation experiments showed that mutation of the SRA1 exon 3 binding site for SRSF1 prevented the binding of SRA1 pre-mRNA. In summary, our results provide experimental evidence that SRA1 exon 3 inclusion is up regulated by SRSF1 to promote tumor invasion and metastasis in hepatocellular carcinoma.

Antisense oligonucleotide-based treatment of retinitis pigmentosa caused by USH2A exon 13 mutations

Mutations in USH2A are among the most common causes of syndromic and non-syndromic retinitis pigmentosa (RP). The two most recurrent mutations in USH2A, c.2299delG and c.2276G > T, both reside in exon 13. Skipping exon 13 from the USH2A transcript presents a potential treatment modality in which the resulting transcript is predicted to encode a slightly shortened usherin protein. Morpholino-induced skipping of ush2a exon 13 in zebrafish ush2a^rmc1 mutants resulted in the production of usherinΔexon 13 protein and a completely restored retinal function. Antisense oligonucleotides were investigated for their potential to selectively induce human USH2A exon 13 skipping. Lead candidate QR-421a induced a concentration-dependent exon 13 skipping in induced pluripotent stem cell (iPSC)-derived photoreceptor precursors from an Usher syndrome patient homozygous for the c.2299delG mutation. Mouse surrogate mQR-421a reached the retinal outer nuclear layer after a single intravitreal injection and induced a detectable level of exon skipping until at least 6 months post-injection. In conclusion, QR-421a-induced exon skipping proves to be a highly promising treatment option for RP caused by mutations in USH2A exon 13.

Widespread Alternative Splicing Changes in Metastatic Breast Cancer Cells

Aberrant alternative splicing (AS) is a hallmark of cancer and a potential target for novel anti-cancer therapeutics. Breast cancer-associated AS events are known to be linked to disease progression, metastasis, and survival of breast cancer patients. To identify altered AS programs occurring in metastatic breast cancer, we perform a global analysis of AS events by using RNA-mediated oligonucleotide annealing, selection, and ligation coupled with next-generation sequencing (RASL-seq). We demonstrate that, relative to low-metastatic, high-metastatic breast cancer cells show different AS choices in genes related to cancer progression. Supporting a global reshape of cancer-related splicing profiles in metastatic breast cancer we found an enrichment of RNA-binding motifs recognized by several splicing regulators, which have aberrant expression levels or activity during breast cancer progression, including SRSF1. Among SRSF1-regulated targets we found DCUN1D5, a gene for which skipping of exon 4 in its pre-mRNA introduces a premature termination codon (PTC), thus generating an unstable transcript degraded by nonsense-mediated mRNA decay (NMD). Significantly, distinct breast cancer subtypes show different DCUN1D5 isoform ratios with metastatic breast cancer expressing the highest level of the NMD-insensitive DCUN1D5 mRNA, thus showing high DCUN1D5 expression levels, which are ultimately associated with poor overall and relapse-free survival in breast cancer patients. Collectively, our results reveal global AS features of metastatic breast tumors, which open new possibilities for the treatment of these aggressive tumor types.

Investigating the ACE2 polymorphisms in COVID-19 susceptibility: An in silico analysis

Background

Novel coronavirus (SARS‐CoV‐2) became an epidemic disease and lead to a pneumonia outbreak first in December 2019 in Wuhan, China. The symptoms related to coronavirus disease‐19 (COVID‐19) were different ranging from mild to severe lung injury and multi‐organ failure symptoms, eventually leading to death, especially in older patients with other co‐morbidities. The receptor of this virus in the human cell is angiotensin‐converting enzyme 2 (ACE2).

Methods

In this paper, we aimed to perform an in silico analysis of the frequently studied variants of the ACE2 gene and determine the effects of the variants in mRNA secondary structure and binding affinity of cellular factors. Fourteen single‐nucleotide polymorphisms were selected based on previous studies and investigated.

Results

All of the variants were analyzed in the RNAsnp database and three revealed a significant p‐value. The spliceAid2 database prediction showed that 7 out of 14 SNPs caused an alteration in a way that only the wild or mutated form was able to bind to proteins. The latter database also reported that three SNPs produces a dual form in which different specific proteins can bind to the sequence in a specific form (either wild or mutated form).

Conclusion

Altogether, these estimations revealed the potential of variants in manipulation of the final stable form of ACE2 that can lead to different COVID‐19 susceptibility.

A Common 3'UTR Variant of the PHOX2B Gene Is Associated With Infant Life-Threatening and Sudden Death Events in the Italian Population

Heterozygous mutations in the Paired like homeobox 2b (PHOX2B) gene are causative of congenital central hypoventilation syndrome (CCHS), a rare monogenic disorder belonging to the family of neurocristopathies and due to a defective development of the autonomic nervous system. Most patients manifest sudden symptoms within 1 year of birth, mainly represented by central apnea and cyanosis episodes. The sudden appearance of hypoxic manifestations in CCHS and their occurrence during sleep resemble two other unexplained perinatal disorders, apparent life-threatening event (ALTE) and sudden and unexpected infant death (SUID), among which the vast majority is represented by sudden infant death syndrome (SIDS). Differently from CCHS, characterized by Mendelian autosomal dominant inheritance, ALTE and SIDS are complex traits, where common genetic variants, together with external factors, may exert an additive effect with symptoms likely manifesting only over a “threshold.” Given the similarities observed among the three abovementioned perinatal disorders, in this work, we have analyzed the frequency of PHOX2B common variants in two groups of Italian idiopathic ALTE (IALTE) and SUIDs/SIDS patients. Here, we report that the c^*161G>A (rs114290493) SNP of the 3′UTR PHOX2B (i) became overrepresented in the two sets of patients compared to population matched healthy controls, and (ii) associated with decreased PHOX2B gene expression, likely mediated by miR-204, a microRNA already known to bind the 3′UTR of the PHOX2B gene. Overall, these results suggest that, at least in the Italian population, the SNP c^*161G>A (rs114290493) does contribute, presumably in association with others mutations or polymorphisms, to confer susceptibility to sudden unexplained perinatal life-threatening or fatal disorders by increasing the effect of miR-204 in inducing PHOX2B expression down-regulation. However, these are preliminary observations that need to be confirmed on larger cohorts to achieve a clinical relevance.

Gene variants of coagulation related proteins that interact with SARS-CoV-2

Thrombosis is a recognized complication of Coronavirus disease of 2019 (COVID-19) and is often associated with poor prognosis. There is a well-recognized link between coagulation and inflammation, however, the extent of thrombotic events associated with COVID-19 warrants further investigation. Poly(A) Binding Protein Cytoplasmic 4 (PABPC4), Serine/Cysteine Proteinase Inhibitor Clade G Member 1 (SERPING1) and Vitamin K epOxide Reductase Complex subunit 1 (VKORC1), which are all proteins linked to coagulation, have been shown to interact with SARS proteins. We computationally examined the interaction of these with SARS-CoV-2 proteins and, in the case of VKORC1, we describe its binding to ORF7a in detail. We examined the occurrence of variants of each of these proteins across populations and interrogated their potential contribution to COVID-19 severity. Potential mechanisms, by which some of these variants may contribute to disease, are proposed. Some of these variants are prevalent in minority groups that are disproportionally affected by severe COVID-19. Therefore, we are proposing that further investigation around these variants may lead to better understanding of disease pathogenesis in minority groups and more informed therapeutic approaches.

Evidence for the loss and recovery of SLAMF9 during human evolution: implications on Dollo's law

Signaling lymphocyte activation molecule family member 9 (SLAMF9) is a cell surface protein of the CD2/SLAM family of leukocyte surface receptors. It is conserved throughout mammals and has roles in the initiation of inflammatory responses and regulation of plasmacytoid dendritic cell function. Through comparison of reference sequences encoding SLAMF9 in human, mouse, and primate sequences, we have determined that the SLAMF9 gene underwent successive mutation events, resulting in the loss of the protein and subsequent recovery of a less stable version. The mutations included a single base pair deletion in the second exon and a change in the splice acceptor site of that same exon. These changes would have had the effect of creating and later repairing a frameshift in the coding sequence. These events took place since the divergence of the human lineage from the chimpanzee-human last common ancestor and represent the first known case of the functional loss and recovery of a gene within the human lineage.

Structure of SRSF1 RRM1 bound to RNA reveals an unexpected bimodal mode of interaction and explains its involvement in SMN1 exon7 splicing

The human prototypical SR protein SRSF1 is an oncoprotein that contains two RRMs and plays a pivotal role in RNA metabolism. We determined the structure of the RRM1 bound to RNA and found that the domain binds preferentially to a CN motif (N is for any nucleotide). Based on this solution structure, we engineered a protein containing a single glutamate to asparagine mutation (E87N), which gains the ability to bind to uridines and thereby activates SMN exon7 inclusion, a strategy that is used to cure spinal muscular atrophy. Finally, we revealed that the flexible inter-RRM linker of SRSF1 allows RRM1 to bind RNA on both sides of RRM2 binding site. Besides revealing an unexpected bimodal mode of interaction of SRSF1 with RNA, which will be of interest to design new therapeutic strategies, this study brings a new perspective on the mode of action of SRSF1 in cells.

SRSF1 is an oncoprotein that plays important roles in RNA metabolism. We reveal the structure of the human SRSF1 RRM1 bound to RNA, and propose a bimodal mode of interaction of the protein with RNA. A single mutation in RRM1 changed SRSF1 specificity for RNA and made it active on SMN2 exon7 splicing.

Development of antisense-mediated myostatin knockdown for the treatment of insulin resistance

Myostatin is a negative regulator of muscle mass and its inhibition represents a promising strategy for the treatment of muscle disorders and type 2 diabetes. However, there is currently no clinically effective myostatin inhibitor, and therefore novel methods are required. We evaluated the use of antisense phosphorodiamidate morpholino oligomers (PMO) to reduce myostatin expression in skeletal muscle and measured their effects on muscle mass and glucose uptake. C57/Bl6 mice received intramuscular or intravenous injections of anti-myostatin PMOs. Repeated intramuscular administration lead to a reduction in myostatin transcript levels (~ 20–40%), and an increase in muscle mass in chow and high-fat diet (HFD)-fed mice, but insulin-stimulated glucose uptake was reduced in PMO-treated muscles of HFD-fed mice. Five weekly intravenous administrations of 100 nmol PMO did not reduce myostatin expression, and therefore had no significant physiological effects. Unexpectedly, exon skipping levels were higher after intramuscular administration of PMO in HFD- than chow-fed mice. These results suggest that a modest PMO-induced reduction in myostatin transcript levels is sufficient to induce an increase in muscle mass, but that a greater degree of inhibition may be required to improve muscle glucose uptake.

Novel mutations in Uridyl-diphosphate-glucuronosyl-transferase 1A1 (UGT1A1) gene in Tunisian patients with unconjugated hyperbilirubinemia

INTRODUCTION

Unconjugated hyperbilirubinemia (UCB) is a feature of Gilbert's syndrome (GS) and Crigler-Najjar's syndrome (CNS), which are two hereditary defects in bilirubin metabolism. Both syndromes are linked to mutations in the UGT1A1 gene, which cause either the decrease or the absence of the UGT1A1 enzymatic activity. Here, we investigated the molecular basis of the UGT1A1 gene in Tunisian patients presenting with unconjugated hyperbilirubinemia.

METHODS

Twenty-four patients with UCB were investigated. The screening protocol for hemoglobinopathies, enzymopathies, and membrane defects was executed in all patients. Afterward, the molecular analysis of the entire UGT1A1 gene was performed by DNA Sanger sequencing. Several bioinformatic tools were used to explore the effects of novel mutations.

RESULTS

Fifteen different UGT1A1 variations were identified, among which four are described here for the first time. In exon 5, the c.1412C > G; p.(Ala471Gly) and c.1589C > T; p.(Ser530Phe) mutations were detected in patients presenting with CNS type I and GS, respectively. In the 3'UTR region of UGT1A1, the c.*90C > T mutation was detected in 3 patients with CNS type I. In the same region, the c.*388C > T defect was found in a GS patient. A deleterious and damaging effect on the UGT1A1 protein were predicted for both exonic mutations. Furthermore, novel microRNAs were identified as targetting the mutated sequences for the 3'UTR mutations.

CONCLUSION

Our study provides novel data on UCB among Tunisians. Furthermore, we report four novel mutations associated with both GS and CNS. The identification of these mutations increases the spectrum of the UGT1A1 mutations and contributes to an understanding of the molecular abnormalities associated with unconjugated hyperbilirubinemia.

Implications of Splicing Alterations in the Onset and Phenotypic Variability of a Family with Subclinical Manifestation of Peutz-Jeghers Syndrome: Bioinformatic and Molecular Evidence

Peutz–Jeghers Syndrome (PJS) is an autosomal dominant pre-cancerous disorder caused in 80–90% of cases by germline mutations in the tumor suppressor gene STK11. We performed a genetic test of the STK11 gene in two Italian young sisters suspected of PJS, since they showed pathognomonic café au lait spots in absence of other symptoms and familiarity. Sequencing of all exons of STK11 gene and other 8 genes, suggested to be involved in hamartomatous syndromes, (PTEN, BMPR1A, SDHB, SDHD, SMAD4, AKT1, ENG, PIK3CA) led to the identification in both the probands of a novel germline silent mutation named c.597 G>A, hitting the last nucleotide of exon 4. Interestingly, genetic testing of the two probands’ parents showed that their unaffected father was carrier of this mutation. Moreover, he carried a second intronic substitution named c.465-51 T>C (rs2075606) which was not inherited by his daughters. We also observed that all the family members carrying the c.597 G>A mutation presented an aberrant splice variant of STK11 mRNA lacking exon 4. Furthermore, in silico analysis of c.465-51 T>C substitution showed that it may activate an Enhancer Splicing Element. Finally, qRT-PCR analysis of STK11 expression levels showed a slight downregulation of the wild type allele in the father and a 2-fold downregulation in the probands compared to the unaffected mother. Our results have led the hypothesis that the c.465-51 T>C intronic variant, which segregates with the wild type allele, could increase the splicing effectiveness of STK11 wild-type allele and compensate the side effect of the c.597 G>A splicing mutation, being responsible for the phenotypic variability observed within this family. This finding highlight the importance of RNA analysis in genetic testing, remarking that silent DNA variant can often be splicing variant involved in disease onset and progression. The identification of these variants has a crucial role to ensure an appropriate follow-up and cancer prevention in at-risk individuals.

Splicing Characteristics of Dystrophin Pseudoexons and Identification of a Novel Pathogenic Intronic Variant in the DMD Gene

Pseudoexon (PE) inclusion has been implicated in various dystrophinopathies; however, its splicing characteristics have not been fully investigated. This study aims to analyze the splicing characteristics of dystrophin PEs and compare them with those of dystrophin canonical exons (CEs). Forty-two reported dystrophin PEs were divided into a splice site (ss) group and a splicing regulatory element (SRE) group. Five dystrophin PEs with characteristics of poison exons were identified and categorized as the possible poison exon group. The comparative analysis of each essential splicing signal among different groups of dystrophin PEs and dystrophin CEs revealed that the possible poison exon group had a stronger 3′ ss compared to any other group. As for auxiliary SREs, different groups of dystrophin PEs were found to have a smaller density of diverse types of exonic splicing enhancers and a higher density of several types of exonic splicing silencers compared to dystrophin CEs. In addition, the possible poison exon group had a smaller density of 3′ ss intronic splicing silencers compared to dystrophin CEs. To our knowledge, our findings indicate for the first time that poison exons might exist in DMD (the dystrophin gene) and present with different splicing characteristics than other dystrophin PEs and CEs.