Listening to silence and understanding nonsense: exonic mutations that affect splicing

Hereditary coagulation factor VII deficiency caused by novel compound heterozygous mutations c.572-1G>A and c.1037A>C in a Chinese pedigree

Hereditary coagulation factor VII (FVII) deficiency is a rare autosomal recessive bleeding disorder. The aims of this study were to identify and verify the pathogenic mutation sites in a family with hereditary coagulation FVII deficiency, and preliminarily explore the underlying mechanisms. We identified a novel combination of compound heterozygous mutations, c.572-1G>A and c.1037A>C in F7 gene, associated with FVII deficiency. The splice site mutation c.572-1G>A led to a truncation, resulting in the loss of the essential catalytic domain of the FVII protein. The c.1037A>C missense mutation has not been previously reported. Our study revealed that this mutation leads to steric hindrance between residues, causing significant changes in the energy and structure of the FVII protein, ultimately affecting its function. These changes disrupt the normal function of the FVII protein, accelerating the development of inherited FVII deficiency. Moreover, the mRNA expression of the F7 gene and the protein expression of the FVII antigen (FVII: Ag) were significantly lower in the proband, as well as in the proband's parents, compared to the healthy control (P<0.05). Our findings not only elucidate the genetic underpinning of FVII deficiency in the family studied but also contribute a new mutation to the known disease spectrum, potentially assisting in future diagnostic and therapeutic approaches.

ParSE-seq: a calibrated multiplexed assay to facilitate the clinical classification of putative splice-altering variants

Interpreting the clinical significance of putative splice-altering variants outside canonical splice sites remains difficult without time-intensive experimental studies. To address this, we introduce Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed assay to quantify variant effects on RNA splicing. We first apply this technique to study hundreds of variants in the arrhythmia-associated gene SCN5A. Variants are studied in 'minigene' plasmids with molecular barcodes to allow pooled variant effect quantification. We perform experiments in two cell types, including disease-relevant induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). The assay strongly separates known control variants from ClinVar, enabling quantitative calibration of the ParSE-seq assay. Using these evidence strengths and experimental data, we reclassify 29 of 34 variants with conflicting interpretations and 11 of 42 variants of uncertain significance. In addition to intronic variants, we show that many synonymous and missense variants disrupted RNA splicing. Two splice-altering variants in the assay also disrupt splicing and sodium current when introduced into iPSC-CMs by CRISPR-Cas9 editing. ParSE-seq provides high-throughput experimental data for RNA-splicing to support precision medicine efforts and can be readily adopted to study other loss-of-function genotype-phenotype relationships.

The Role of RNA Splicing in Liver Function and Disease: A Focus on Metabolic Dysfunction-Associated Steatotic Liver Disease

RNA splicing is an essential post-transcriptional mechanism that facilitates the excision of introns and the connection of exons to produce mature mRNA, which is essential for gene expression and proteomic diversity. In the liver, precise splicing regulation is critical for maintaining metabolic balance, detoxification, and protein synthesis. This review explores the mechanisms of RNA splicing and the role of splicing factors, particularly in the context of Metabolic Dysfunction-Associated Steatotic Liver Disease (MASLD). This review also highlights how RNA splicing dysregulation can lead to aberrant splicing and impact the progression of liver diseases such as MASLD, with a particular focus on Metabolic Dysfunction-Associated Steatohepatitis (MASH), which represents the advanced stage of MASLD. Recent advances in the clinical application of antisense oligonucleotides (ASOs) to correct splicing errors offer promising therapeutic strategies for restoring normal liver function. Additionally, the dysregulation of splicing observed in liver diseases may serve as a potential diagnostic marker, offering new opportunities for early identification of individuals more susceptible to disease progression. This review provides insights into the molecular mechanisms that govern splicing regulation in the liver, with a particular emphasis on MASLD, and discusses potential therapeutic approaches targeting RNA splicing to treat MASLD and related metabolic disorders.

Splicing-specific transcriptome-wide association uncovers genetic mechanisms for schizophrenia

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

Epigenetic control of skeletal muscle atrophy

Skeletal muscular atrophy is a complex disease involving a large number of gene expression regulatory networks and various biological processes. Despite extensive research on this topic, its underlying mechanisms remain elusive, and effective therapeutic approaches are yet to be established. Recent studies have shown that epigenetics play an important role in regulating skeletal muscle atrophy, influencing the expression of numerous genes associated with this condition through the addition or removal of certain chemical modifications at the molecular level. This review article comprehensively summarizes the different types of modifications to DNA, histones, RNA, and their known regulators. We also discuss how epigenetic modifications change during the process of skeletal muscle atrophy, the molecular mechanisms by which epigenetic regulatory proteins control skeletal muscle atrophy, and assess their translational potential. The role of epigenetics on muscle stem cells is also highlighted. In addition, we propose that alternative splicing interacts with epigenetic mechanisms to regulate skeletal muscle mass, offering a novel perspective that enhances our understanding of epigenetic inheritance's role and the regulatory network governing skeletal muscle atrophy. Collectively, advancements in the understanding of epigenetic mechanisms provide invaluable insights into the study of skeletal muscle atrophy. Moreover, this knowledge paves the way for identifying new avenues for the development of more effective therapeutic strategies and pharmaceutical interventions.

Pathogenic relationship between phenotypes of ARPKD and novel compound heterozygous mutations of PKHD1

Background

To investigate whether the novel mutation of PKHD1 could cause polycystic kidney disease by affecting splicing with a recessive inheritance pattern.

Methods

A nonconsanguineous Chinese couple with two recurrent pregnancies showed fetal enlarged echogenic polycystic kidney and oligoamnios were recruited. Pedigree WES, minigene splicing assay experiment and following bioinformatics analysis were performed to verify the effects, and inheritance pattern of diseasing-causing mutations.

Results

WES revealed that both fetuses were identified as carrying the same novel mutation c.3592_3628 + 45del, p.? and c.11207 T>C, p.(Ile3736Thr) in the PKHD1 gene (NM_138694.4), which inherited from the father and mother respectively. Both bioinformatic method prediction and minigene splicing assay experience results supported the mutation c.3592_3628 + 45del, p.? affects the splicing of the PKHD1 transcript, resulting in exon 31 skipping. Another missense mutation c.11207 T>C, p.(Ile3736Thr) has a low frequency in populations and is predicted to be deleterious by bioinformatic methods.

Conclusion

These findings provide a direct clinical and functional evidence that the truncating mutations of the PKHD1 gene could lead to more severe phenotypes, and cause ARPKD as a homozygous or compound heterozygous pattern. Our study broadens the variant spectrum of the PKHD1 gene and provides a basis for genetic counseling and diagnosis of ARPKD.

A second hotspot for pathogenic exon-skipping variants in CDC45

Biallelic pathogenic variants in CDC45 are associated with Meier-Gorlin syndrome with craniosynostosis (MGORS type 7), which also includes short stature and absent/hypoplastic patellae. Identified variants act through a hypomorphic loss of function mechanism, to reduce CDC45 activity and impact DNA replication initiation. In addition to missense and premature termination variants, several pathogenic synonymous variants have been identified, most of which cause increased exon skipping of exon 4, which encodes an essential part of the RecJ-orthologue's DHH domain. Here we have identified a second cohort of families segregating CDC45 variants, where patients have craniosynostosis and a reduction in height, alongside common facial dysmorphisms, including thin eyebrows, consistent with MGORS7. Skipping of exon 15 is a consequence of two different variants, including a shared synonymous variant that is enriched in individuals of East Asian ancestry, while other variants in trans are predicted to alter key intramolecular interactions in α/β domain II, or cause retention of an intron within the 3'UTR. Our cohort and functional data confirm exon skipping is a relatively common pathogenic mechanism in CDC45, and highlights the need for alternative splicing events, such as exon skipping, to be especially considered for variants initially predicted to be less likely to cause the phenotype, particularly synonymous variants.

Case report: A novel somatic SDHB variant in a patient with bladder paraganglioma

Background

Paragangliomas (PGL) are rare neuroendocrine tumors derived from the autonomic nervous system paraganglia. Urinary bladder paragangliomas (UBPGL) originate from the sympathetic neurons of the urinary bladder wall and represent 0.7% of all paragangliomas and <0.05% of all bladder tumors. PGL and UBPGL can be associated with SDHB, SDHD, NF1, and VHL gene variants, with the most common germline alterations found in SDHB and VHL.

Case report

We report a case of a 42-year-old woman who presented with menorrhagia/hematuria, uterine leiomyomas, as well as cardiac and bladder masses. The cardiac mass was favored to be a myxoma based on clinical findings, while the bladder mass was diagnosed as UBPGL. A novel SDHB mutation (c.642G>A, p Q214Q), detected in the UBPGL, was proven to be somatic. Although this variant was seemingly synonymous, it was predicted to have a loss of function due to the splice site effect, which was further supported by the immunohistochemical loss of SDHB.

Conclusion

This case highlights the challenges of diagnosing an extremely rare entity, bladder paraganglioma, with an emphasis on the multidisciplinary approach to navigate various clinical and imaging findings that may initially be misleading. In addition, a novel loss of function SDHB variant that could have been overlooked as a synonymous variant is herein reported, while also illustrating the importance of both germline and somatic mutation testing.

Relative impact of THPO mutation causing hereditary thrombocythemia

Germline mutations of THPO were reported as causes of hereditary thrombocythemia. Six previously reported distinct sites of the mutation were clustered at the 5`-untranslated region or the exon 3 splicing donor site of the THPO gene. Each mutation was identified in an independent pedigree, and the differences between the mutations were not compared. We cloned six distinct THPO mutations (THPO c.-47delG, THPO c.-31G>T, THPO c.13G>A, THPO c.13+1G>A, THPO c.13+2T>C, and THPO c.13+5G>A) and compared the molecular mechanisms that underlie the increased production of THPO protein. At the transcript level, all of the mutations except THPO c.-47delG showed an exon 3 skipping transcript, including two mutations (THPO c.-31G>T and THPO c.13+5G>A) that were distant from the splicing donor site. THPO c.-47delG showed the same full-length transcript as that of the wild-type transcript. At the protein level, all mutations resulted in a higher level of production of thrombopoietin (THPO) protein compared with wild-type THPO. There are only two distinct patterns of mechanisms for increased production of THPO: exon 3 skipping that deleted upstream suppressive open reading frame (ORF)7 and one base deletion that shifted ORF7 to connect to the initial codon of THPO in-frame. The common mechanisms of hereditary thrombocytosis due to THPO mutations are unleashed THPO translations, which are usually suppressed by upstream out-of-frame ORF7.

Three exonic variants in the PHEX gene cause aberrant splicing in a minigene assay

Background: X-linked hypophosphatemia (XLH, OMIM 307800) is a rare phosphorus metabolism disorder caused by PHEX gene variants. Many variants simply classified as missense or nonsense variants were only analyzed at the DNA level. However, growing evidence indicates that some of these variants may alter pre-mRNA splicing, causing diseases. Therefore, this study aimed to use bioinformatics tools and a minigene assay to ascertain the effects of PHEX variations on pre-mRNA splicing. Methods: We analyzed 174 variants in the PHEX gene described as missense or nonsense variants. Finally, we selected eight candidate variants using bioinformatics tools to evaluate their effects on pre-mRNA splicing using a minigene assay system. The complementary DNA (cDNA) sequence for the PHEX gene (RefSeq NM_000444.6) serves as the basis for DNA variant numbering. Results: Of the eight candidate variants, three were found to cause abnormal splicing. Variants c.617T>G p.(Leu206Trp) and c.621T>A p.(Tyr207*) in exon 5 altered the splicing of pre-mRNA, owing to the activation of a cryptic splice site in exon 5, which produced an aberrant transcript lacking a part of exon 5, whereas variant c.1700G>C p.(Arg567Pro) in exon 16 led to the activation of a cryptic splice site in intron 16, resulting in a partial inclusion of intron 16. Conclusion: Our study employed a minigene system, which has a great degree of flexibility to assess abnormal splicing patterns under the circumstances of patient mRNA samples that are not available, to explore the impact of the exonic variants on pre-mRNA splicing. Based on the aforementioned experimental findings, we demonstrated the importance of analyzing exonic variants at the mRNA level.

Pathway-based, reaction-specific annotation of disease variants for elucidation of molecular phenotypes

Germline and somatic mutations can give rise to proteins with altered activity, including both gain and loss-of-function. The effects of these variants can be captured in disease-specific reactions and pathways that highlight the resulting changes to normal biology. A disease reaction is defined as an aberrant reaction in which a variant protein participates. A disease pathway is defined as a pathway that contains a disease reaction. Annotation of disease variants as participants of disease reactions and disease pathways can provide a standardized overview of molecular phenotypes of pathogenic variants that is amenable to computational mining and mathematical modeling. Reactome (https://reactome.org/), an open source, manually curated, peer-reviewed database of human biological pathways, in addition to providing annotations for >11 000 unique human proteins in the context of ∼15 000 wild-type reactions within more than 2000 wild-type pathways, also provides annotations for >4000 disease variants of close to 400 genes as participants of ∼800 disease reactions in the context of ∼400 disease pathways. Functional annotation of disease variants proceeds from normal gene functions, described in wild-type reactions and pathways, through disease variants whose divergence from normal molecular behaviors has been experimentally verified, to extrapolation from molecular phenotypes of characterized variants to variants of unknown significance using criteria of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Reactome's data model enables mapping of disease variant datasets to specific disease reactions within disease pathways, providing a platform to infer pathway output impacts of numerous human disease variants and model organism orthologs, complementing computational predictions of variant pathogenicity. Database URL: https://reactome.org/.

Genetic analysis of PALB2 gene WD40 domain in canine mammary tumour patients

BACKGROUND

DNA repair mechanisms are essential for tumorigenesis and disruption of HR mechanism is an important predisposing factor of human breast cancers (BC). PALB2 is an important part of the HR. There are similarities between canine mammary tumours (CMT) and BCs. As its human counterpart, PALB2 mutations could be a predisposing factor of CMT.

OBJECTIVES

In this study, we aimed to investigate the impacts of PALB2 variants on tumorigenesis and canine mammary tumor (CMT) malignancy.

METHODS

We performed Sanger sequencing to detect germline mutations in the WD40 domain of the canine PALB2 gene in CMT patients. We conducted in silico analysis to investigate the variants, and compared the germline PALB2 mutations in humans that cause breast cancer (BC) with the variants detected in dogs with CMT.

RESULTS

We identified an intronic (c.3096+8C>G) variant, two exonic (p.A1050V and p.R1354R) variants, and a 3' UTR variant (c.4071T>C). Of these, p.R1354R and c.4071T>C novel variants were identified for the first time in this study. We found that the p.A1050V mutation had a significant effect. However, we could not determine sufficient similarity due to the differences in nucleotide/amino acid sequences between two species. Nonetheless, possible variants of human sequences in the exact location as their dog counterparts are associated with several cancer types, implying that the variants could be crucial for tumorigenesis in dogs. Our results did not show any effect of the variants on tumor malignancy.

CONCLUSIONS

The current project is the first study investigating the relationship between the PALB2 gene WD40 domain and CMTs. Our findings will contribute to a better understanding of the pathogenic mechanism of the PALB2 gene in CMTs. In humans, variant positions in canines have been linked to cancer-related phenotypes such as familial BC, endometrial tumor, and hereditary cancer predisposition syndrome. The results of bioinformatics analyses should be investigated through functional tests or case-control studies.

Association between Genetic Polymorphism of SCN1A, GABRA1 and ABCB1 and Drug Responsiveness in Vietnamese Epileptic Children

Background and Objectives: Drug resistant epilepsy (DRE) is a major hurdle in epilepsy, which hinders clinical care, patients' management and treatment outcomes. DRE may partially result from genetic variants that alter proteins responsible for drug targets and drug transporters in the brain. We aimed to examine the relationship between SCN1A, GABRA1 and ABCB1 polymorphism and drug response in epilepsy children in Vietnam. Materials and Methods: In total, 213 children diagnosed with epilepsy were recruited in this study (101 were drug responsive and 112 were drug resistant). Sanger sequencing had been performed in order to detect six single nucleotide polymorphisms (SNPs) belonging to SCN1A (rs2298771, rs3812718, rs10188577), GABRA1 (rs2279020) and ABCB1 (rs1128503, rs1045642) in study group. The link between SNPs and drug response status was examined by the Chi-squared test or the Fisher's exact test. Results: Among six investigated SNPs, two SNPs showed significant difference between the responsive and the resistant group. Among those, heterozygous genotype of SCN1A rs2298771 (AG) were at higher frequency in the resistant patients compared with responsive patients, playing as risk factor of refractory epilepsy. Conversely, the heterozygous genotype of SCN1A rs3812718 (CT) was significantly lower in the resistant compared with the responsive group. No significant association was found between the remaining four SNPs and drug response. Conclusions: Our study demonstrated a significant association between the SCN1A genetic polymorphism which increased risk of drug-resistant epilepsy in Vietnamese epileptic children. This important finding further supports the underlying molecular mechanisms of SCN1A genetic variants in the pathogenesis of drug-resistant epilepsy in children.

Quantifying negative selection on synonymous variants

Widespread adoption of DNA sequencing has resulted in large numbers of genetic variants, whose contribution to disease is not easily determined. Although many types of variation are known to disrupt cellular processes in predictable ways, for some categories of variants, the effects may not be directly detectable. A particular example is synonymous variants, that is, those single-nucleotide variants that create a codon substitution, such that the produced amino acid sequence is unaffected. Contrary to the original theory suggesting that synonymous variants are benign, there is a growing volume of research showing that, despite their "silent" mechanism of action, some synonymous variation may be deleterious. Here, we studied the extent of the negative selective pressure acting on different classes of synonymous variants by analyzing the relative enrichment of synonymous singleton variants in the human exomes provided by gnomAD. Using a modification of the mutability-adjusted proportion of singletons (MAPS) metric as a measure of purifying selection, we found that some classes of synonymous variants are subject to stronger negative selection than others. For instance, variants that reduce codon optimality undergo stronger selection than optimality-increasing variants. Besides, selection affects synonymous variants implicated in splice-site-loss or splice-site-gain events. To understand what drives this negative selection, we tested a number of predictors in the aim to explain the variability in the selection scores. Our findings provide insights into the effects of synonymous variants at the population level, highlighting the specifics of the role that these variants play in health and disease.

Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Pleiotropic effects of different exonic nucleotide changes at the same position contribute to hemophilia B phenotypic variation

BACKGROUND

The disease-causing effects of genetic variations often depend on their location within a gene. Exonic changes generally lead to alterations in protein production, secretion, activity, or clearance. However, owing to the overlap between proteins and splicing codes, missense variants can also affect messenger RNA splicing, thus adding a layer of complexity and influencing disease phenotypes.

OBJECTIVES

To extensively characterize a panel of 13 exonic variants in the F9 gene occurring at 6 different factor IX positions and associated with varying severities of hemophilia B (HB).

METHODS

Computational predictions, splicing analysis, and recombinant factor IX assays were exploited to characterize F9 variants.

RESULTS

We demonstrated that 5 (38%) of 13 selected F9 exonic variants have pleiotropic effects. Although bioinformatic approaches accurately classified effects, extensive experimental assays were required to elucidate and deepen the molecular mechanisms underlying the pleiotropic effects. Importantly, their characterization was instrumental in developing tailored RNA therapeutics based on engineered U7 small nuclear RNA to mask cryptic splice sites and compensatory U1 small nuclear RNA to enhance exon definition.

CONCLUSION

Overall, albeit a multitool bioinformatic approach suggested the molecular effects of multiple HB variants, the deep investigation of molecular mechanisms revealed insights into the HB phenotype-genotype relationship, enabling accurate classification of HB variants. Importantly, knowledge of molecular mechanisms allowed the development of tailored RNA therapeutics, which can also be translated to other genetic diseases.

RNA-binding proteins potentially regulate the alternative splicing of apoptotic genes during knee osteoarthritis progression

BACKGROUND

Alternative splicing (AS) is a principal mode of genetic regulation and one of the most widely used mechanisms to generate structurally and functionally distinct mRNA and protein variants. Dysregulation of AS may result in aberrant transcription and protein products, leading to the emergence of human diseases. Although considered important for regulating gene expression, genome-wide AS dysregulation, underlying mechanisms, and clinical relevance in knee osteoarthritis (OA) remain unelucidated. Therefore, in this study, we elucidated and validated AS events and their regulatory mechanisms during OA progression.

RESULTS

In this study, we identified differentially expressed genes between human OA and healthy meniscus samples. Among them, the OA-associated genes were primarily enriched in biological pathways such as extracellular matrix organization and ossification. The predominant OA-associated regulated AS (RAS) events were found to be involved in apoptosis during OA development. The expression of the apoptosis-related gene BCL2L13, XAF1, and NF2 were significantly different between OA and healthy meniscus samples. The construction of a covariation network of RNA-binding proteins (RBPs) and RAS genes revealed that differentially expressed RBP genes LAMA2 and CUL4B may regulate the apoptotic genes XAF1 and BCL2L13 to undergo AS events during OA progression. Finally, RT-qPCR revealed that CUL4B expression was significantly higher in OA meniscus samples than in normal controls and that the AS ratio of XAF1 was significantly different between control and OA samples; these findings were consistent with their expected expression and regulatory relationships.

CONCLUSIONS

Differentially expressed RBPs may regulate the AS of apoptotic genes during knee OA progression. XAF1 and its regulator, CUL4B, may serve as novel biomarkers and potential therapeutic targets for this disease.

Detecting and understanding meaningful cancerous mutations based on computational models of mRNA splicing

Cancer research has long relied on non-silent mutations. Yet, it has become overwhelmingly clear that silent mutations can affect gene expression and cancer cell fitness. One fundamental mechanism that apparently silent mutations can severely disrupt is alternative splicing. Here we introduce Oncosplice, a tool that scores mutations based on models of proteomes generated using aberrant splicing predictions. Oncosplice leverages a highly accurate neural network that predicts splice sites within arbitrary mRNA sequences, a greedy transcript constructor that considers alternate arrangements of splicing blueprints, and an algorithm that grades the functional divergence between proteins based on evolutionary conservation. By applying this tool to 12M somatic mutations we identify 8K deleterious variants that are significantly depleted within the healthy population; we demonstrate the tool's ability to identify clinically validated pathogenic variants with a positive predictive value of 94%; we show strong enrichment of predicted deleterious mutations across pan-cancer drivers. We also achieve improved patient survival estimation using a proposed set of novel cancer-involved genes. Ultimately, this pipeline enables accelerated insight-gathering of sequence-specific consequences for a class of understudied mutations and provides an efficient way of filtering through massive variant datasets - functionalities with immediate experimental and clinical applications.

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants

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.

Three exonic variants in the COL4A5 gene alter RNA splicing in a minigene assay

BACKGROUND

X-linked Alport syndrome (XLAS) is an inherited renal disease caused by rare variants of COL4A5 on chromosome Xq22. Many studies have indicated that single nucleotide variants (SNVs) in exons can disrupt normal splicing process of the pre-mRNA by altering various splicing regulatory signals. The male patients with XLAS have a strong genotype-phenotype correlation. Confirming the effect of variants on splicing can help to predict kidney prognosis. This study aimed to investigate whether single nucleotide substitutions, located within three bases at the 5' end of the exons or internal position of the exons in COL4A5 gene, cause aberrant splicing process.

METHODS

We analyzed 401 SNVs previously presumed missense and nonsense variants in COL4A5 gene by bioinformatics programs and identified candidate variants that may affect the splicing of pre-mRNA via minigene assays.

RESULTS

Our study indicated three of eight candidate variants induced complete or partial exon skipping. Variants c.2678G>C and c.2918G>A probably disturb classic splice sites leading to corresponding exon skipping. Variant c.3700C>T may disrupt splicing enhancer motifs accompanying with generation of splicing silencer sequences resulting in the skipping of exon 41.

CONCLUSION

Our study revealed that two missense variants positioned the first nucleotides of the 5' end of COL4A5 exons and one internal exonic nonsense variant caused aberrant splicing. Importantly, this study emphasized the necessity of assessing the effects of SNVs at the mRNA level.

RNA elements that control human papillomavirus mRNA splicing-targets for therapy?

Human papillomaviruses (HPVs) cause more than 4.5% of all cancer in the world and more than half of these cases are attributed to human papillomavirus type 16 (HPV16). Prophylactic vaccines are available but antiviral drugs are not. Novel targets for therapy are urgently needed. Alternative RNA splicing is extensively used by HPVs to express all their genes and HPV16 is no exception. This process must function to perfection since mis-splicing could perturb the HPV gene expression program by altering mRNA levels or by generating dysfunctional mRNAs. Cis-acting RNA elements on the viral mRNAs and their cognate cellular trans-acting factors control papillomavirus RNA splicing. The precise but delicate nature of the splicing process renders splicing sensitive to interference. As such, papillomavirus RNA splicing is a potential target for therapy. Here we summarize our current understanding of cis-acting HPV16 RNA elements that control HPV16 mRNA splicing via cellular proteins and discuss how they may be exploited as targets for therapy to papillomavirus infections and cancer.

Synonymous edits in the Escherichia coli genome have substantial and condition-dependent effects on fitness

CRISPR-Cas-based genome editing is widely used in bacteria at scales ranging from construction of individual mutants to massively parallel libraries. This procedure relies on guide RNA-directed cleavage of the genome followed by repair with a template that introduces a desired mutation along with synonymous "immunizing" mutations to prevent re-cleavage of the genome after editing. Because the immunizing mutations do not change the protein sequence, they are often assumed to be neutral. However, synonymous mutations can change mRNA structures in ways that alter levels of the encoded proteins. We have tested the assumption that immunizing mutations are neutral by constructing a library of over 50,000 edits that consist of only synonymous mutations in Escherichia coli. Thousands of edits had substantial effects on fitness during growth of E. coli on acetate, a poor carbon source that is toxic at high concentrations. The percentage of high-impact edits varied considerably between genes and at different positions within genes. We reconstructed clones with high-impact edits and found that 69% indeed had significant effects on growth in acetate. Interestingly, fewer edits affected fitness during growth in glucose, a preferred carbon source, suggesting that changes in protein expression caused by synonymous mutations may be most important when an organism encounters challenging conditions. Finally, we showed that synonymous edits can have widespread effects; a synonymous edit at the 5' end of ptsI altered expression of hundreds of genes. Our results suggest that the synonymous immunizing edits introduced during CRISPR-Cas-based genome editing should not be assumed to be innocuous.

Broken silence: 22,841 predicted deleterious synonymous variants identified in the human exome through computational analysis

Synonymous single nucleotide variants (sSNVs) do not alter the primary structure of a protein, thus it was previously accepted that they were neutral. Recently, several studies demonstrated their significance to a range of diseases. Still, variant prioritization strategies lack focus on sSNVs. Here, we identified 22,841 deleterious synonymous variants in 125,748 human exomes using two in silico predictors (SilVA and CADD). While 98.2% of synonymous variants are classified as neutral, 1.8% are predicted to be deleterious, yielding an average of 9.82 neutral and 0.18 deleterious sSNVs per exome. Further investigation of prediction features via Heterogeneous Ensemble Feature Selection revealed that impact on amino acid sequence and conservation carry the most weight for a deleterious prediction. Thirty nine detrimental sSNVs are not rare and are located on disease associated genes. Ten distinct putatively non-deleterious sSNVs are likely to be under positive selection in the North-Western European and East Asian populations. Taken together our analysis gives voice to the so-called silent mutations as we propose a robust framework for evaluating the deleteriousness of sSNVs in variant prioritization studies.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

A founder COL4A3 pathogenic variant resulting in Alport syndrome and thin basement membrane disease: a case report series

Alport syndrome is a rare genetic condition characterized by kidney disease, hearing impairment, and ocular abnormalities. It exhibits various inheritance patterns involving pathogenic variants in COL4A3, COL4A4, and COL4A5 genes. The phenotypes can range from isolated hematuria with a non-progressive or very slowly progressive course to progressive kidney disease with extrarenal abnormalities. Timely diagnosis of Alport syndrome facilitates the early and effective implementation of treatment, as well as genetic counseling. Here, we report the COL4A3 c.765G > A, p.((=)) mutation in three ethnically Azerbaijani, apparently unrelated, consanguineous families from the village of Algeti in the Marneuli region of Georgia. We speculate that this variant could represent a founder mutation within this population and recommend offering genetic testing to Algeti village residents with persistent hematuria.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays

BACKGROUND

Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes.

RESULTS

We benchmark eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compare experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, is lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieve the best overall performance at distinguishing disruptive and neutral variants, and controlling for overall call rate genome-wide, SpliceAI and Pangolin have superior sensitivity. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues.

CONCLUSION

SpliceAI and Pangolin show the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

Adenine base editor-mediated splicing remodeling activates noncanonical splice sites

Adenine base editors (ABEs) are genome-editing tools that have been harnessed to introduce precise A•T to G•C conversion. The discovery of split genes revealed that all introns contain two highly conserved dinucleotides, canonical "AG" (acceptor) and "GT" (donor) splice sites. ABE can directly edit splice acceptor sites of the adenine (A) base, leading to aberrant gene splicing, which may be further adopted to remodel splicing. However, spliced isoforms triggered with ABE have not been well explored. To address it, we initially generated a cell line harboring C-terminal enhanced GFP (eGFP)-tagged β-actin (ACTB), in which the eGFP signal can track endogenous β-actin expression. Expectedly, after the editing of splice acceptor sites, we observed a dramatical decrease in the percentage of eGFP-positive cells and generation of splicing products with the noncanonical splice site. Furthermore, we manipulated Peroxidasin in mouse embryos with ABE, in which a noncanonical acceptor was activated to remodel splicing, successfully generating a mouse disease model of anophthalmia and severely malformed microphthalmia. Collectively, we demonstrate that ABE-mediated splicing remodeling can activate a noncanonical acceptor to manipulate human and mouse genomes, which will facilitate the investigation of basic and translational medicine studies.

Polymorphism and expression level of the FADS3 gene and associated with the growth traits in Hu sheep

Growth traits are the economically important traits of sheep, and screening for genes related to growth and development is helpful for the genetic improvement of ovine growth traits. The fatty acid desaturase 3 (FADS3) is one of the important genes affecting the synthesis and accumulation of polyunsaturated fatty acids in animals. In this study, the expression levels of the FADS3 gene and polymorphism of the FADS3 gene associated with growth traits in Hu sheep were detected using quantitative real-time PCR (qRT-PCR), Sanger sequencing, and KAspar assay. The result showed that the expression levels of the FADS3 gene were widely expressed in all tissues, and the expression level of FADS3 in the lung was significantly higher than in other tissues (p < .05). Then, the polymorphism locus g. 2918 A > C was detected in intron 2 of the FADS3 gene, and associated analysis showed that the mutation in the FADS3 gene was associated significantly with growth traits (including body weight, body height, body length, and chest circumference, p < .05). Therefore, individuals with AA genotype showed significantly better growth traits than those with CC genotype, and FADS3 gene could be a candidate gene for improving growth traits in Hu sheep.

Efficient exon skipping by base-editor-mediated abrogation of exonic splicing enhancers

Duchenne muscular dystrophy (DMD) is a severe genetic disease caused by the loss of the dystrophin protein. Exon skipping is a promising strategy to treat DMD by restoring truncated dystrophin. Here, we demonstrate that base editors (e.g., targeted AID-mediated mutagenesis [TAM]) are able to efficiently induce exon skipping by disrupting functional redundant exonic splicing enhancers (ESEs). By developing an unbiased and high-throughput screening to interrogate exonic sequences, we successfully identify novel ESEs in DMD exons 51 and 53. TAM-CBE (cytidine base editor) induces near-complete skipping of the respective exons by targeting these ESEs in patients' induced pluripotent stem cell (iPSC)-derived cardiomyocytes. Combined with strategies to disrupt splice sites, we identify suitable single guide RNAs (sgRNAs) with TAM-CBE to efficiently skip most DMD hotspot exons without substantial double-stranded breaks. Our study thus expands the repertoire of potential targets for CBE-mediated exon skipping in treating DMD and other RNA mis-splicing diseases.

Characterization of pre-mRNA Splicing Defects Caused by CLCN5 and OCRL Mutations and Identification of Novel Variants Associated with Dent Disease

Dent disease (DD) is an X-linked renal tubulopathy characterized by low-molecular-weight proteinuria, hypercalciuria, nephrocalcinosis, nephrolithiasis and progressive renal failure. Two-thirds of cases are associated with inactivating variants in the CLCN5 gene (Dent disease 1, DD1) and a few present variants in the OCRL gene (Dent disease 2, DD2). The aim of the present study was to test the effect on the pre-mRNA splicing process of DD variants, described here or in the literature, and describe the clinical and genotypic features of thirteen unrelated patients with suspected DD. All patients presented tubular proteinuria, ten presented hypercalciuria and five had nephrolithiasis or nephrocalcinosis. CLCN5 and OCRL genes were analyzed by Sanger sequencing. Nine patients showed variants in CLCN5 and four in OCRL; eight of these were new. Bioinformatics tools were used to select fifteen variants with a potential effect on pre-mRNA splicing from our patients' group and from the literature, and were experimentally tested using minigene assays. Results showed that three exonic missense mutations and two intronic variants affect the mRNA splicing process. Our findings widen the genotypic spectrum of DD and provide insight into the impact of variants causing DD.

Pathway-based, reaction-specific annotation of disease variants for elucidation of molecular phenotypes

Disease variant annotation in the context of biological reactions and pathways can provide a standardized overview of molecular phenotypes of pathogenic mutations that is amenable to computational mining and mathematical modeling. Reactome, an open source, manually curated, peer-reviewed database of human biological pathways, provides annotations for over 4000 disease variants of close to 400 genes in the context of ∼800 disease reactions constituting ∼400 disease pathways. Functional annotation of disease variants proceeds from normal gene functions, through disease variants whose divergence from normal molecular behaviors has been experimentally verified, to extrapolation from molecular phenotypes of characterized variants to variants of unknown significance using criteria of the American College of Medical Genetics and Genomics (ACMG). Reactome's pathway-based, reaction-specific disease variant dataset and data model provide a platform to infer pathway output impacts of numerous human disease variants and model organism orthologs, complementing computational predictions of variant pathogenicity.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

Identification of two short peptide motifs from serine/arginine-rich protein ribonucleic acid recognition motif-1 domain acting as splicing regulators

Background

Serine/arginine-rich (SR) proteins regulate pre-mRNA splicing. However, structurally similar proteins often behave differently in splicing regulation and the underlying mechanisms are largely unknown. Here, using SMN1/2 minigenes we extensively analyzed four SR proteins, SRSF1/5/6/9.

Methods

In this study, the effects of these proteins on SMN1/2 exon 7 splicing when tethered at either intron 6 or 7 were evaluated using an MS2-tethering assay. Deletion analysis in four SR proteins and co-overexpression analysis were performed.

Results

Splicing outcomes varied among all four SR proteins, SRSF1 and SRSF5 function the same at the two sites, acting as repressor and stimulator, respectively; while SRSF6 and SRSF9 promote exon 7 inclusion at only one site. Further, the key domains of each SR proteins were investigated, which identified a potent inhibitory nonapeptide in the C-terminus of SRSF1/9 ribonucleic acid recognition motif-1 (RRM1) and a potent stimulatory heptapeptide at the N-terminus of SRSF5/6 RRM1.

Conclusion

The insight of the four SR proteins and their domains in affecting SMN gene splicing brings a new perspective on the modes of action of SR proteins; and the functional peptides obtained here offers new ideas for developing splice switching-related therapies.

Genetic Factors of Teeth Impaction: Polymorphic and Haplotype Variants of PAX9, MSX1, AXIN2, and IRF6 Genes

In recent research, there has been a growing awareness of the role of genetic factors in the positioning and eruption of teeth in the maxilla and mandible. This study aimed to evaluate the potential of specific polymorphic markers of single nucleotide polymorphisms (SNPs) located within the PAX9, MSX1, AXIN2, and IRF6 genes to determine the predisposition to tooth impaction. The study participants were divided into two groups: the first group consisted of individuals with at least one impacted secondary tooth. In contrast, the second group (control group) had no impacted teeth in their jaws. To analyze the genes, real-time PCR (polymerase chain reaction) and TaqMan probes were utilized to detect the selected polymorphisms. The findings suggest that disruptions in the structure and function of the mentioned genetic factors such as polymorphic and haplotype variants of PAX9, MSX1, AXIN2, and IRF6 genes, which play a direct role in tooth and periodontal tissue development, might be significant factors in tooth impaction in individuals with genetic variations. Therefore, it is reasonable to hypothesize that tooth impaction may be influenced, at least in part, by the presence of specific genetic markers, including different allelic variants of the PAX9, AXIN2, and IRF6 genes, and especially MSX1.

ParSE-seq: A Calibrated Multiplexed Assay to Facilitate the Clinical Classification of Putative Splice-altering Variants

Background

Interpreting the clinical significance of putative splice-altering variants outside 2-base pair canonical splice sites remains difficult without functional studies.

Methods

We developed Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed minigene-based assay, to test variant effects on RNA splicing quantified by high-throughput sequencing. We studied variants in SCN5A, an arrhythmia-associated gene which encodes the major cardiac voltage-gated sodium channel. We used the computational tool SpliceAI to prioritize exonic and intronic candidate splice variants, and ClinVar to select benign and pathogenic control variants. We generated a pool of 284 barcoded minigene plasmids, transfected them into Human Embryonic Kidney (HEK293) cells and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), sequenced the resulting pools of splicing products, and calibrated the assay to the American College of Medical Genetics and Genomics scheme. Variants were interpreted using the calibrated functional data, and experimental data were compared to SpliceAI predictions. We further studied some splice-altering missense variants by cDNA-based automated patch clamping (APC) in HEK cells and assessed splicing and sodium channel function in CRISPR-edited iPSC-CMs.

Results

ParSE-seq revealed the splicing effect of 224 SCN5A variants in iPSC-CMs and 244 variants in HEK293 cells. The scores between the cell types were highly correlated (R2=0.84). In iPSCs, the assay had concordant scores for 21/22 benign/likely benign and 24/25 pathogenic/likely pathogenic control variants from ClinVar. 43/112 exonic variants and 35/70 intronic variants with determinate scores disrupted splicing. 11 of 42 variants of uncertain significance were reclassified, and 29 of 34 variants with conflicting interpretations were reclassified using the functional data. SpliceAI computational predictions correlated well with experimental data (AUC = 0.96). We identified 20 unique SCN5A missense variants that disrupted splicing, and 2 clinically observed splice-altering missense variants of uncertain significance had normal function when tested with the cDNA-based APC assay. A splice-altering intronic variant detected by ParSE-seq, c.1891-5C>G, also disrupted splicing and sodium current when introduced into iPSC-CMs at the endogenous locus by CRISPR editing.

Conclusions

ParSE-seq is a calibrated, multiplexed, high-throughput assay to facilitate the classification of candidate splice-altering variants.

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.

Genome editing in the treatment of ocular diseases

Genome-editing technologies have ushered in a new era in gene therapy, providing novel therapeutic strategies for a wide range of diseases, including both genetic and nongenetic ocular diseases. These technologies offer new hope for patients suffering from previously untreatable conditions. The unique anatomical and physiological features of the eye, including its immune-privileged status, size, and compartmentalized structure, provide an optimal environment for the application of these cutting-edge technologies. Moreover, the development of various delivery methods has facilitated the efficient and targeted administration of genome engineering tools designed to correct specific ocular tissues. Additionally, advancements in noninvasive ocular imaging techniques and electroretinography have enabled real-time monitoring of therapeutic efficacy and safety. Herein, we discuss the discovery and development of genome-editing technologies, their application to ocular diseases from the anterior segment to the posterior segment, current limitations encountered in translating these technologies into clinical practice, and ongoing research endeavors aimed at overcoming these challenges.

The novel compound heterozygous variants identified in a Chinese family with glucose phosphate isomerase deficiency and pathogenicity analysis

BACKGROUND AND AIMS

Glucose phosphate isomerase (GPI) deficiency is an extremely rare autosomal recessive disorder caused by mutations in the GPI gene. In this research, the proband displaying typical manifestations of haemolytic anaemia and his family members were recruited to analyse the pathogenicity of the detected variants.

METHODS

Peripheral blood samples were collected from the family members and genomic DNA was extracted and targeted for capture and sequencing. The effect of the candidate pathogenic variants on splicing was further investigated using the minigene splicing system. The computer simulation was also used for further analysis of the detected data.

RESULTS

The proband carried the compound heterozygous variants c.633 + 3 A > G and c.295G > T in the GPI gene, which have never been reported before. In the genealogy, co-segregation of the mutant genotype with the phenotype was established. The minigene study showed that intronic mutations resulted in abnormal pre-mRNA splicing. Specifically, the two aberrant transcripts: r.546_633del and r.633 + 1_633 + 2insGT were transcribed by the minigene plasmid expressing the c.633 + 3 A > G variant. The missense mutation c.295G > T in exon 3 resulted in altering glycine at codon 87 to cysteine which was predicted to be pathogenic in an in silico analysis. Deeper analyses revealed that the Gly87Cys missense mutation led to steric hindrance. Compared to the wild-type, the mutation G87C led to a significant increase in intermolecular forces.

CONCLUSION

Overall, the novel compound heterozygous variants in the GPI gene contributed to the etiology of the disease. Genetic testing can assist in the diagnosis. The novel gene variants identified in the present study has further expanded the mutational spectrum of GPI deficiency, which can better guide family counselling.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

Stripes and loss of color in ball pythons (Python regius) are associated with variants affecting endothelin signaling

Color patterns in nonavian reptiles are beautifully diverse, but little is known about the genetics and development of these patterns. Here, we investigated color patterning in pet ball pythons (Python regius), which have been bred to show color phenotypes that differ dramatically from the wildtype form. We report that several color phenotypes in pet animals are associated with putative loss-of-function variants in the gene encoding endothelin receptor EDNRB1: (1) frameshift variants in EDNRB1 are associated with conversion of the normal mottled color pattern to skin that is almost fully white, (2) missense variants affecting conserved sites of the EDNRB1 protein are associated with dorsal, longitudinal stripes, and (3) substitutions at EDNRB1 splice donors are associated with subtle changes in patterning compared to wildtype. We propose that these phenotypes are caused by loss of specialized color cells (chromatophores), with loss ranging from severe (fully white) to moderate (dorsal striping) to mild (subtle changes in patterning). Our study is the first to describe variants affecting endothelin signaling in a nonavian reptile and suggests that reductions in endothelin signaling in ball pythons can produce a variety of color phenotypes, depending on the degree of color cell loss.

Analysis of 200 unrelated individuals with a constitutional NF1 deep intronic pathogenic variant reveals that variants flanking the alternatively spliced NF1 exon 31 [23a] cause a classical neurofibromatosis type 1 phenotype while altering predominantly NF1 isoform type II

Neurofibromatosis type 1 results from loss-of-function NF1 pathogenic variants (PVs). Up to 30% of all NF1 PVs disrupt mRNA splicing, including deep intronic variants. Here, we retrospectively investigated the spectrum of NF1 deep intronic PVs in a cohort of 8,090 unrelated individuals from the University of Alabama at Birmingham (UAB) dataset with a molecularly confirmed neurofibromatosis type 1. All variants were identified through their effect on the NF1 transcript, followed by variant characterization at the DNA-level. A total of 68 distinct variants, which were ≥ 20 nucleotides away from the closest exon-intron junction, were identified in 2.5% unrelated individuals with NF1 (200/8,090). Nine different pathogenic splice variants, identified in 20 probands, led to exonization of different parts of intron 30 [23.2] or 31 [23a]. The two major NF1 transcript isoforms, distinguished by the absence (type I) or presence (type II) of the alternatively spliced cassette exon 31 [23a], are equally expressed in blood in control individuals without NF1 or NF1-affected individuals carrying their PV not in the introns flanking exon 31 [23a]. By fragment and cloning analysis we demonstrated that the exonization of intron 31 [23a] sequences due to deep intronic PV predominantly affects the NF1 isoform II. Seven additional (likely) pathogenic NF1 deep intronic variants not observed in the UAB dataset were found by classification of 36 variants identified by a literature search. Hence, the unique list of these 75 deep intronic (likely) PVs should be included in any comprehensive NF1 testing strategy.

A novel PKHD1 splicing variant identified in a fetus with autosomal recessive polycystic kidney disease

Objective: Variants of the polycystic kidney and hepatic disease 1 (PKHD1) gene are associated with autosomal recessive polycystic kidney disease (ARPKD). This study aimed to identify the genetic causes in a Chinese pedigree with ARPKD and design a minigene construct of the PKHD1 gene to investigate the impact of its variants on splicing. Methods: Umbilical cord samples from the proband and peripheral blood samples from his parents were collected, and genomic DNA was extracted for whole-exome sequencing (WES). Bioinformatic analysis was used to identify potential genetic causes, and Sanger sequencing confirmed the existence of variants within the pedigree. A minigene assay was performed to validate the effects of an intronic variant on mRNA splicing. Results: Two variants, c.9455del (p.N3152Tfs*10) and c.2408-13C>G, were identified in the PKHD1 gene (NM_138694.4) by WES; the latter has not been previously reported. In silico analysis predicted that this intronic variant is potentially pathogenic. Bioinformatic splice prediction tools revealed that the variant is likely to strongly impact splice site function. An in vitro minigene assay revealed that c.2408-13C>G can cause aberrant splicing, resulting in the retention of 12 bp of intron 23. Conclusion: A novel pathogenic variant of PKHD1, c.2408-13C>G, was found in a fetus with ARPKD, which enriches the variant spectrum of the PKHD1 gene and provides a basis for genetic counseling and the diagnosis of ARPKD. Minigenes are optimal to determine whether intron variants can cause aberrant splicing.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Base editors: development and applications in biomedicine

Base editor (BE) is a gene-editing tool developed by combining the CRISPR/Cas system with an individual deaminase, enabling precise single-base substitution in DNA or RNA without generating a DNA double-strand break (DSB) or requiring donor DNA templates in living cells. Base editors offer more precise and secure genome-editing effects than other conventional artificial nuclease systems, such as CRISPR/Cas9, as the DSB induced by Cas9 will cause severe damage to the genome. Thus, base editors have important applications in the field of biomedicine, including gene function investigation, directed protein evolution, genetic lineage tracing, disease modeling, and gene therapy. Since the development of the two main base editors, cytosine base editors (CBEs) and adenine base editors (ABEs), scientists have developed more than 100 optimized base editors with improved editing efficiency, precision, specificity, targeting scope, and capacity to be delivered in vivo, greatly enhancing their application potential in biomedicine. Here, we review the recent development of base editors, summarize their applications in the biomedical field, and discuss future perspectives and challenges for therapeutic applications.

A wave of deep intronic mutations in X-linked Alport Syndrome

X-linked Alport syndrome (XLAS) is an inherited kidney disease caused exclusively by pathogenic variants in the COL4A5 gene. In 10-20% of cases, DNA sequencing of COL4A5 exons or flanking regions cannot identify molecular causes. Here, our objective was to use a transcriptomic approach to identify causative events in a group of 19 patients with XLAS without identified mutation by Alport gene panel sequencing. Bulk RNAseq and/or targeted RNAseq using a capture panel of kidney genes was performed. Alternative splicing events were compared to those of 15 controls by a developed bioinformatic score. When using targeted RNAseq, COL4A5 coverage was found to be 23-fold higher than with bulk RNASeq and revealed 30 significant alternative splicing events in 17 of the 19 patients. After computational scoring, a pathogenic transcript was found in all patients. A causative variant affecting COL4A5 splicing and absent in the general population was identified in all cases. Altogether, we developed a simple and robust method for identification of aberrant transcripts due to pathogenic deep-intronic COL4A5 variants. Thus, these variants, potentially targetable by specific antisense oligonucleotide therapies, were found in a high percentage of patients with XLAS in whom pathogenic variants were missed by conventional DNA sequencing.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays

Background

Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes.

Results

We benchmarked eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compared experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, was lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieved the best overall performance at distinguishing disruptive and neutral variants. Controlling for overall call rate genome-wide, SpliceAI and Pangolin also showed superior overall sensitivity for identifying SDVs. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues.

Conclusion

SpliceAI and Pangolin showed the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

Targeting splicing factors for cancer therapy

Alternative splicing (AS) of mRNAs is an essential regulatory mechanism in eukaryotic gene expression. AS misregulation, caused by either dysregulation or mutation of splicing factors, has been shown to be involved in cancer development and progression, making splicing factors suitable targets for cancer therapy. In recent years, various types of pharmacological modulators, such as small molecules and oligonucleotides, targeting distinct components of the splicing machinery, have been under development to treat multiple disorders. Although these approaches have promise, targeting the core spliceosome components disrupts the early stages of spliceosome assembly and can lead to nonspecific and toxic effects. New research directions have been focused on targeting specific splicing factors for a more precise effect. In this Perspective, we will highlight several approaches for targeting splicing factors and their functions and suggest ways to improve their specificity.

Detection of unintended on-target effects in CRISPR genome editing by DNA donors carrying diagnostic substitutions

CRISPR nucleases can introduce double-stranded DNA breaks in genomes at positions specified by guide RNAs. When repaired by the cell, this may result in the introduction of insertions and deletions or nucleotide substitutions provided by exogenous DNA donors. However, cellular repair can also result in unintended on-target effects, primarily larger deletions and loss of heterozygosity due to gene conversion. Here we present a strategy that allows easy and reliable detection of unintended on-target effects as well as the generation of control cells that carry wild-type alleles but have demonstratively undergone genome editing at the target site. Our 'sequence-ascertained favorable editing' (SAFE) donor approach relies on the use of DNA donor mixtures containing the desired nucleotide substitutions or the wild-type alleles together with combinations of additional 'diagnostic' substitutions unlikely to have any effects. Sequencing of the target sites then results in that two different sequences are seen when both chromosomes are edited with 'SAFE' donors containing different sets of substitutions, while a single sequence indicates unintended effects such as deletions or gene conversion. We analyzed more than 850 human embryonic stem cell clones edited with 'SAFE' donors and detect all copy number changes and almost all clones with gene conversion.

Identification and in silico characterization of CSRP3 synonymous variants in dilated cardiomyopathy

BACKGROUND

Synonymous variations have always been ignored while studying the underlying genetic mechanisms for most of the human diseases. However, recent studies have suggested that these silent changes in the genome can alter the protein expression and folding.

METHODS AND RESULTS

CSRP3, which is a well-known candidate gene associated with dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM), was screened for 100 idiopathic DCM cases and 100 controls. Three synonymous variations were identified viz., c.96G > A, p.K32=; c.336G > A, p.A112=; c.354G > A, p.E118=. A comprehensive in silico analysis was performed using various web based widely accepted tools, Mfold, Codon Usage, HSF3.1 and RNA22. Mfold predicted structural changes in all the variants except c.96 G > A (p.K32=), however it predicted changes in the stability of mRNA due to all the synonymous variants. Codon bias was observed as evident by the Relative Synonymous Codon Usage and Log Ratio of Codon Usage Frequencies. The Human Splicing Finder also predicted remarkable changes in the regulatory elements in the variants c.336G > A and c.354 G > A. The miRNA target prediction using varied modes available in RNA22 revealed that 70.6% of the target sites of miRNAs in CSRP3 were altered due to variant c.336G > A while 29.41% sites were completely lost.

CONCLUSION

Findings of the present study suggest that synonymous variants revealed striking deviations in the structural conformation of mRNA, stability of mRNA, relative synonymous codon usage, splicing and miRNA binding sites from the wild type suggesting their possible role in the pathogenesis of DCM, either by destabilizing the mRNA structure, or codon usage bias or else altering the cis-acting regulatory elements during splicing.

Individual Genomic Loci and mRNA Levels of Immune Biomarkers Associated with Pneumonia Susceptibility in Baladi Goats

The effectiveness of breeding for inherent disease resistance in animals could be considerably increased by identifying the genes and mutations that cause diversity in disease resistance. One hundred and twenty adult female Baladi goats (sixty pneumonic and sixty apparently healthy) were used in this study. DNA and RNA were extracted from blood samples collected from the jugular vein of each goat. SLC11A1, CD-14, CCL2, TLR1, TLR7, TLR8, TLR9, β defensin, SP110, SPP1, BP1, A2M, ADORA3, CARD15, IRF3, and SCART1 SNPs that have been previously found to be associated with pneumonia resistance/susceptibility were identified via PCR-DNA sequencing. The pneumonic and healthy goats differed significantly, according to a Chi-square analysis of the discovered SNPs. The mRNA levels of the studied immune markers were noticeably greater in the pneumonic goats than in the healthy ones. The findings could support the significance of the use of immune gene expression profiles and nucleotide variations as biomarkers for the susceptibility/resistance to pneumonia and provide a practical management technique for Baladi goats. These results also suggest a potential strategy for lowering pneumonia in goats by employing genetic markers linked to an animal’s ability to fend off infection in selective breeding.