BACKGROUND: Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.
RESULTS: A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5\' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.
CONCLUSIONS: Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

BACKGROUND: Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis.
RESULTS: Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity.
CONCLUSIONS: By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

The human U1 snRNA is encoded by a multigene family consisting of transcribed variants and defective pseudogenes. Many variant U1 (vU1) snRNAs have been demonstrated to not only be transcribed but also processed by the addition of a trimethylated guanosine cap, packaged into snRNPs, and assembled into spliceosomes; however, their capacity to facilitate pre-mRNA splicing has, so far, not been tested. A recent systematic analysis of the human snRNA genes identified 178 U1 snRNA genes that are present in the genome as either tandem arrays or single genes on multiple chromosomes. Of these, 15 were found to be expressed in human tissues and cell lines, although at significantly low levels from their endogenous loci, <0.001% of the canonical U1 snRNA. In this study, we found that placing the variants in the context of the regulatory elements of the RNU1-1 gene improves the expression of many variants to levels comparable to the canonical U1 snRNA. Application of a previously established HeLa cell-based minigene reporter assay to examine the capacity of the vU1 snRNAs to support pre-mRNA splicing revealed that even though the exogenously expressed variant snRNAs were enriched in the nucleus, only a few had a measurable effect on splicing.

Congenital contractural arachnodactyly (CCA) is a rare connective tissue disorder characterized by arachnodactyly, multiple joint contractures, progressive kyphoscoliosis, pectus deformity and abnormal crumpled ears. FBN2 is the only gene currently known to be associated with CCA. In this study, we report on a prenatal case presented with skeletal, cardiac and spinal malformations. And his father had elongated limbs, contractures of the proximal interphalangeal joints, high myopia and scoliosis. We conducted whole exome sequencing (WES) on the fetus-parental trio and a heterozygous variant (hg19 chr5:127,673,685, c.3598 + 4A > G, NM_001999.4) in intron 27 of the FBN2 gene was successfully identified, inherited from the father. Reverse transcriptase-polymerase chain reaction (RT-PCR) was performed to evaluate the potential splicing effect of this variant, which confirmed that the variant caused a deletion of exon 27 (126 bp) by disrupting the splice-donor site and destroyed the 17th calcium-binding epidermal growth factor-like (cbEGF) domain. Our research not only finds the etiology of the disease in affected individuals and expands the mutation spectrum of FBN2 gene, but also provides genetic counseling and fertility guidance for this family.

HLA-DRB1*11:01:01:12N differs from HLA-DRB1*11:01:01:03 by one nucleotide substitution in intron 3 at position c.652+1G>C, hg19.

Duchenne muscular dystrophy is a rare and lethal hereditary disease responsible for progressive muscle wasting due to mutations in the DMD gene. We used the CRISPR-Cas9 Prime editing technology to develop different strategies to correct frameshift mutations in DMD gene carrying the deletion of exon 52 or exons 45 to 52. With optimized epegRNAs, we were able to induce the specific substitution of the GT nucleotides of the splice donor site of exon 53 in up to 32% of HEK293T cells and 28% of patient myoblasts. We also achieved up to 44% and 29% deletion of the G nucleotide of the GT splice site of exon 53, as well as inserted 17% and 5.5% GGG between the GT splice donor site of exon 51 in HEK293T cells and human myoblasts, respectively. The modification of the splice donor site for exon 51 and exon 53 provoke their skipping and allowed exon 50 to connect to exon 53 and allowed exon 44 to connect to exon 54, respectively. These corrections restored the expression of dystrophin as demonstrated by western blot. Thus, Prime editing was used to induce specific substitutions, insertions and deletions in the splice donor sites for exons 51 and 53 to correct the frameshift mutations in DMD gene carrying deletions of exon 52 and exons 45 to 52, respectively.

Predicting the impact of coding and noncoding variants on splicing is challenging, particularly in non-canonical splice sites, leading to missed diagnoses in patients. Existing splice prediction tools are complementary but knowing which to use for each splicing context remains difficult. Here, we describe Introme, which uses machine learning to integrate predictions from several splice detection tools, additional splicing rules, and gene architecture features to comprehensively evaluate the likelihood of a variant impacting splicing. Through extensive benchmarking across 21,000 splice-altering variants, Introme outperformed all tools (auPRC: 0.98) for the detection of clinically significant splice variants. Introme is available at https://github.com/CCICB/introme .

Primary ciliary dyskinesia (PCD) is a rare hereditary orphan condition that results in variable phenotypes, including infertility. About 50 gene variants are reported in the scientific literature to cause PCD, and among them, dynein axonemal assembly factor 4 ( DNAAF4 ) has been recently reported. DNAAF4 has been implicated in the preassembly of a multiunit dynein protein essential for the normal function of locomotory cilia as well as flagella. In the current study, a single patient belonging to a Chinese family was recruited, having been diagnosed with PCD and asthenoteratozoospermia. The affected individual was a 32-year-old male from a nonconsanguineous family. He also had abnormal spine structure and spinal cord bends at angles diagnosed with scoliosis. Medical reports, laboratory results, and imaging data were investigated. Whole-exome sequencing, Sanger sequencing, immunofluorescence analysis, hematoxylin-eosin staining, and in silico functional analysis, including protein modeling and docking studies, were used. The results identified DNAAF4 disease-related variants and confirmed their pathogenicity. Genetic analysis through whole-exome sequencing identified two pathogenic biallelic variants in the affected individual. The identified variants were a hemizygous splice site c.784-1G>A and heterozygous 20.1 Kb deletion at the DNAAF4 locus, resulting in a truncated and functionless DNAAF4 protein. Immunofluorescence analysis indicated that the inner dynein arm was not present in the sperm flagellum, and sperm morphological analysis revealed small sperm with twisted and curved flagella or lacking flagella. The current study found novel biallelic variants causing PCD and asthenoteratozoospermia, extending the range of DNAAF4 pathogenic variants in PCD and associated with the etiology of asthenoteratozoospermia. These findings will improve our understanding of the etiology of PCD.

X-linked dominant hypophosphatemia (XLH), the most common form of hereditary hypophosphatemic rickets/osteomalacia, is caused by loss-of-function phosphate-regulating endopeptidase homolog X-linked gene (PHEX) variants. However, synonymous PHEX variants are rare in XLH. We report a 7-year-old boy with hypophosphatemia, short stature, and lower limb deformity. Whole-exome sequencing, reverse transcription-polymerase chain reaction, and Sanger sequencing were performed to identify the pathogenicity of the variant. A novel synonymous PHEX variant (NM_000444.4:c.1530 C>T, p.Arg510Arg) was detected in the proband. Further analysis revealed a 58-bp deletion at the 5\' site of exon 14 during splicing. This study extends the genetic spectrum of XLH and confirms the rarity and significance of synonymous PHEX variants.

Alternative splicing accounts for a considerable portion of transcriptomic diversity, as most protein-coding genes are spliced into multiple mRNA isoforms. However, errors in splicing patterns can give rise to mis-splicing with pathological consequences, such as the congenital diseases familial dysautonomia, Duchenne muscular dystrophy, and spinal muscular atrophy. Small nuclear RNA (snRNA) components of the U snRNP family have been proposed as a therapeutic modality for the treatment of mis-splicing. U1 snRNAs offer great promise, with prior studies demonstrating in vivo efficacy, suggesting additional preclinical development is merited. Improvements in enabling technologies, including screening methodologies, gene delivery vectors, and relevant considerations from gene editing approaches justify further advancement of U1 snRNA as a therapeutic and research tool. With the goal of providing a user-friendly protocol, we compile and demonstrate a methodological toolkit for sequence-specific targeted perturbation of alternatively spliced pre-mRNA with engineered U1 snRNAs. We observe robust modulation of endogenous pre-mRNA transcripts targeted in two contrasting splicing contexts, SMN2 exon 7 and FAS exon 6, exhibiting the utility and applicability of engineered U1 snRNA to both inclusion and exclusion of targeted exons. We anticipate that these demonstrations will contribute to the usability of U1 snRNA in investigating splicing modulation in eukaryotic cells, increasing accessibility to the broader research community.