Pre-mRNA splicing and human disease

Genome-wide DNA methylation profiles in the raphe nuclei of patients with autism spectrum disorder

AIM

Autism spectrum disorder (ASD) has a strong genetic basis, yet its genetic complexities remain elusive. Current research highlights environmental factors and epigenetic processes, such as DNA methylation, as crucial in ASD development. This exploratory study addresses a gap in understanding epigenetic regulation in the dorsal raphe (DR)-a region regulating multiple neurotransmitters and implicated in ASD-by examining DNA methylation profiles in postmortem ASD and control brains.

METHODS

We comprehensively analyzed genome-wide DNA methylation profiles in the DR brain region (seven controls and five ASD) using the Infinium HumanMethylation450 BeadChip (Illumina). Additionally, quantitative polymerase chain reaction was used to measure messenger RNA levels of differentially methylated genes in ASD (11 controls and six ASD).

RESULTS

We identified differentially methylated regions (DMRs) between ASD and controls. These DMRs were located among various genomic regions, including promoters, gene bodies, and intergenic regions. Notably, we found hypermethylation in genes related to olfaction (e.g. OR2C3), which is regulated by serotonin. Additionally, we observed that the hypomethylation of promoter-associated CpG islands in RABGGTB, a gene related to autophagy and synaptic function, corresponded with its increased expression.

CONCLUSIONS

Our findings reveal extensive DNA methylation changes in critical genomic regions, shedding light on potential mechanisms underlying ASD. The identification of RABGGTB as a novel candidate gene, not listed in the SFARI database, underscores its significance and warrants further research to explore its role in ASD diagnosis. This study enhances our understanding of the epigenetic landscape in ASD, emphasizing the interplay between genetic and environmental factors in its pathophysiology.

Splicing to orchestrate cell fate

Alternative splicing (AS) plays a critical role in gene expression by generating protein diversity from single genes. This review provides an overview of the role of AS in regulating cell fate, focusing on its involvement in processes such as cell proliferation, differentiation, apoptosis, and tumorigenesis. We explore how AS influences the cell cycle, particularly its impact on key stages like G1, S, and G2/M. The review also examines AS in cell differentiation, highlighting its effects on mesenchymal stem cells and neurogenesis, and how it regulates differentiation into adipocytes, osteoblasts, and chondrocytes. Additionally, we discuss the role of AS in programmed cell death, including apoptosis and pyroptosis, and its contribution to cancer progression. Importantly, targeting aberrant splicing mechanisms presents promising therapeutic opportunities for restoring normal cellular function. By synthesizing recent findings, this review provides insights into how AS governs cellular fate and offers directions for future research into splicing regulatory networks.

RNA splicing variants of the novel long non-coding RNA, CyKILR, possess divergent biological functions in non-small cell lung cancer

The CDKN2A gene, responsible for encoding the tumor suppressors p16(INK4A) and p14(ARF), is frequently inactivated in non-small cell lung cancer (NSCLC). Herein, an uncharacterized long non-coding RNA (lncRNA) (ENSG00000267053) on chromosome 19p13.12 was found to be overexpressed in NSCLC cells with an active, wild-type CDKN2A gene. This lncRNA, named cyclin-dependent kinase inhibitor 2A-regulated lncRNA (CyKILR), also correlated with an active WT STK11 gene, which encodes the tumor suppressor, liver kinase B1. CyKILR displayed two splice variants, CyKILRa (exon 3 included) and CyKILRb (exon 3 excluded), which are cooperatively regulated by CDKN2A and STK11 as knockdown of both tumor suppressor genes was required to induce a significant loss of exon 3 inclusion in mature CyKILR RNA. CyKILRa localized to the nucleus, and its downregulation using antisense RNA oligonucleotides enhanced cellular proliferation, migration, clonogenic survival, and tumor incidence. In contrast, CyKILRb localized to the cytoplasm, and its downregulation using small interfering RNA reduced cell proliferation, migration, clonogenic survival, and tumor incidence. Transcriptomics analyses revealed the enhancement of apoptotic pathways with concomitant suppression of key cell-cycle pathways by CyKILRa demonstrating its tumor-suppressive role. CyKILRb inhibited tumor suppressor miRNAs indicating an oncogenic nature. These findings elucidate the intricate roles of lncRNAs in cell signaling and tumorigenesis.

Insights From a Novel Splicing Variant and Recurrent Arginine Variants in the CHD3 Gene Causing Snijders Blok-Campeau Syndrome

Snijders Blok-Campeau syndrome (SNIBCPS, OMIM#618205) is an autosomal dominant neurodevelopmental disorder attributed to pathogenic variants in the chromodomain helicase DNA binding protein 3 (CHD3) gene. To date, more than 100 individuals have been diagnosed with SNIBCPS. The syndrome is characterized by intellectual disability, global developmental delay, speech or language impediments, and dysmorphic features associated with macrocephaly. Additionally, affected individuals may exhibit behavioral issues, hypotonia, and autistic traits. A novel splicing variant (c.5590+1G > T) in the C-terminal 2 region of the CHD3 gene was identified in a patient predominantly exhibiting autistic characteristics. In vitro minigene splicing experiments conducted in HEK293 cells revealed that aberrant splicing resulted in the formation of a cryptic site 46 nucleotides downstream of the 5' splice site. This alteration was predicted to disrupt the reading frame by eliminating the physiological stop codon, consequently causing an extension in protein translation. Furthermore, an additional patient presenting with hypotonia, dysmorphic features, and global developmental delay was documented. This patient harbored a missense variant in the helicase C-terminal domain, c.3505C > T (p. Arg1169Trp). The pathogenic variant was anticipated to impact chromatin remodeling capacity and enzyme activity. Given the high prevalence of arginine residue pathogenic variants in the CHD3 protein and its notable propensity for binding and storing ATP molecules, intriguing insights into the potential effects of arginine residue pathogenic variants on phenotypes are provided. These findings contribute to a more comprehensive understanding of the genetic landscape of SNIBCPS while elucidating potential molecular mechanisms underlying the syndrome.

Alternative Splicing as a Modulator of the Interferon-Gamma Pathway

Interferon-gamma (IFN-γ) is a critical cytokine that plays a pivotal role in immune system regulation. It is a key mediator of both cellular defense mechanisms and antitumor immunity. As the sole member of the type II interferon family, IFN-γ modulates immune responses by activating macrophages, enhancing natural killer cell function, and regulating gene expression across multiple cellular processes. Alternative splicing is a post-transcriptional gene expression regulatory mechanism that generates multiple mature messenger RNAs from a single gene, dramatically increasing proteome diversity without the need of a proportional genome expansion. This process occurs in 90-95% of human genes, with alternative splicing events allowing for the production of diverse protein isoforms that can have distinct-or even opposing-functional properties. Alternative splicing plays a crucial role in cancer immunology, potentially generating tumor neoepitopes and modulating immune responses. However, how alternative splicing affects IFN-γ's activity is still poorly understood. This review explores how alternative splicing regulates the expression and function of both upstream regulators and downstream effectors of IFN-γ, revealing complex mechanisms of gene expression and immune response modulation. Key transcription factors and signaling molecules of the IFN-γ pathway are alternatively spliced, and alternative splicing can dramatically alter IFN-γ signaling, immune cell function, and response to environmental cues. Specific splice variants can enhance or inhibit IFN-γ-mediated immune responses, potentially influencing cancer immunotherapy, autoimmune conditions, and infectious disease outcomes. The emerging understanding of these splicing events offers promising therapeutic strategies for manipulating immune responses through targeted molecular interventions.

Role of silent mutations in KRAS -mutant tumors

ABSTRACT

Silent mutations within the RAS  gene have garnered increasing attention for their potential roles in tumorigenesis and therapeutic strategies. Kirsten-RAS ( KRAS ) mutations, predominantly oncogenic, are pivotal drivers in various cancers. While extensive research has elucidated the molecular mechanisms and biological consequences of active KRAS  mutations, the functional significance of silent mutations remains relatively understudied. This review synthesizes current knowledge on KRAS  silent mutations, highlighting their impact on cancer development. Silent mutations, which do not alter protein sequences but can affect RNA stability and translational efficiency, pose intriguing questions regarding their contribution to tumor biology. Understanding these mutations is crucial for comprehensively unraveling KRAS -driven oncogenesis and exploring novel therapeutic avenues. Moreover, investigations into the clinical implications of silent mutations in KRAS -mutant tumors suggest potential diagnostic and therapeutic strategies. Despite being in early stages, research on KRAS  silent mutations holds promise for uncovering novel insights that could inform personalized cancer treatments. In conclusion, this review underscores the evolving landscape of KRAS  silent mutations, advocating for further exploration to bridge fundamental biology with clinical applications in oncology.

Whole Blood DNA Methylation Analysis Reveals Epigenetic Changes Associated with ARSACS

Autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) is a rare inherited condition described worldwide and characterized by a wide spectrum of heterogeneity in terms of genotype and phenotype. How sacsin loss leads to neurodegeneration is still unclear, and current knowledge indicates that sacsin is involved in multiple functional mechanisms. We hence hypothesized the existence of epigenetic factors, in particular alterations in methylation patterns, that could contribute to ARSACS pathogenesis and explain the pleiotropic effects of SACS further than pathogenic mutations. To investigate this issue, we recruited eight patients affected by ARSACS, four characterized by early onset of the disease and four with late onset. We performed Whole Genome Bisulfite Sequencing using DNA from peripheral blood to define the methylome of patients and compared them with a control group. Our analysis showed that patients with ARSACS exhibit an altered methylation pattern and that the observed differences exist also among affected individuals with different age of onset. Our study provides valuable insights for employing epigenetic biomarkers to assess the severity and progression of this disorder and propels further investigations into the role of epigenetic processes in ARSACS pathogenesis.

UBL5 and Its Role in Viral Infections

Unlike other ubiquitin-like family members, UBL5 is structurally and functionally atypical, and a novel role in various biological processes and diseases has been discovered. UBL5 can stabilize the structure of the spliceosome, can promote post-transcriptional processing, and has been implicated in both DNA damage repair and protein unfolding reactions, as well as cellular mechanisms that are frequently exploited by viruses for their own proliferation during viral infections. In addition, UBL5 can inhibit viral infection by binding to the non-structural protein 3 of rice stripe virus and mediating its degradation. Therefore, UBL5 is an important link between viral infections and immunity, and its study will be beneficial for the prevention and treatment of viral infections in the future. However, a review of the current findings on the role of UBL5 in viral infection has not been undertaken. Therefore, in this review, we summarize the recent progress in understanding the functions of UBL5 and discuss its putative role in viral infections.

One-Size-Fits-Many: Antisense oligonucleotides for rescuing splicing mutations in hotspot exons

Mutations that impact splicing play a significant role in disease etiology but are not fully understood. To characterize the impact of exonic variants on splicing in 71 clinically-actionable disease genes in asymptomatic people, we analyzed 32,112 exonic mutations from ClinVar and Geisinger MyCode using a minigene reporter assay. We identify 1,733 splice-disrupting mutations, of which the most extreme 1-2% of variants are likely to be deleterious. We report that these variants are not distributed evenly across exons but are mostly concentrated in the ∼8% of exons that are most susceptible to splicing mutations (i.e. hotspot exons). We demonstrate that splicing defects in these exons can be reverted by ASOs targeting the splice sites of either their upstream or downstream flanking exons. This finding supports the feasibility of developing single therapeutic ASOs that could revert all splice-altering variants localized to a particular exon.

p12 isoform-2 is a regulatory subunit of human DNA polymerase delta and is dysregulated in various cancers

Dysregulation of human DNA polymerase delta (Polδ) subunits is associated with genome instability and pathological disorders. Genome databases suggest the expression of several spliced variants of subunits which may alter Polδ function. Here, we analyzed the protein-encoding variants of the Polδ subunit p12 and their association with cancer. p12 isoform-2 (p12*) encodes a 79 aa protein with a C-terminal tail distinct from the previously characterized p12. Like p12, p12* dimerizes and interacts with p125 and p50 subunits and is thus an integral component of Polδ. Further, we observed dysregulated p12* expression in low-grade glioma, renal, thyroid, and pancreatic carcinomas. This study identifies a previously unrecognized Polδ complex and highlights a possible regulatory role of p12 variants in cellular phenotypes.

RNA splicing variants of the novel long non-coding RNA, CyKILR, possess divergent biological functions in non-small cell lung cancer

The CDKN2A gene, responsible for encoding the tumor suppressors p16(INK4A) and p14(ARF), is frequently inactivated in non-small cell lung cancer (NSCLC). Herein, an uncharacterized long non-coding RNA (lncRNA) (ENSG00000267053) on chromosome 19p13.12 was found to be overexpressed in NSCLC cells with an active, wild-type CDKN2A gene. This lncRNA, named Cy clin-Dependent K inase I nhibitor 2A-regulated l nc R NA (CyKILR), also correlated with an active WT STK11 gene, which encodes the tumor suppressor, Liver kinase B1. CyKILR displayed two splice variants, CyKILRa (exon 3 included) and CyKILRb (exon 3 excluded), which are cooperatively regulated by CDKN2A and STK11 as knockdown of both tumor suppressor genes was required to induce a significant loss of exon 3 inclusion in mature CyKILR RNA. CyKILRa localized to the nucleus, and its downregulation using antisense RNA oligonucleotides enhanced cellular proliferation, migration, clonogenic survival, and tumor incidence. In contrast, CyKILRb localized to the cytoplasm, and its downregulation using siRNA reduced cell proliferation, migration, clonogenic survival, and tumor incidence. Transcriptomics analyses revealed enhancement of apoptotic pathways with concomitant suppression of key cell cycle pathways by CyKILRa demonstrating its tumor-suppressive role. CyKILRb inhibited tumor suppressor microRNAs indicating an oncogenic nature. These findings elucidate the intricate roles of lncRNAs in cell signaling and tumorigenesis.

Reduced transcriptome analysis in zebrafish uncovers disruptors of spliceosome and ribosome biosynthesis

Abnormal biosynthesis of spliceosomes and ribosomes can lead to their dysfunction, which in turn disrupts protein synthesis and results in various diseases. While genetic factors have been extensively studied, our understanding of how environmental compounds interfere with spliceosome and ribosome biosynthesis remains limited. In the present study, we employed a Reduced Transcriptome Analysis (RTA) approach, integrating large-scale transcriptome data sets of zebrafish and compiling a specific zebrafish gene panel focusing on the spliceosome and ribosome, to elucidate the potential disruptors targeting their biosynthesis. Transcriptomic data sets for 118 environmental substances and 1400 related gene expression profiles were integrated resulting in 513 exposure signatures. Among these substances, several categories including PCB126, transition metals Lanthanum (La) and praseodymium (Pr), heavy metals Cd2+ and AgNO3 and atrazine were highlighted for inducing the significant transcriptional alterations. Furthermore, we found that the transcriptional patterns were distinct between categories, yet overlapping patterns were generally observed within each group. For instance, over 82 % differentially expressed ribosomal genes were shared between La and Pr within the equivalent concentration range. Additionally, transcriptional complexities were also evident across various organs and developmental stages of zebrafish, with notable differences in the inhibition of the transcription of various spliceosome subunits. Overall, our results provide novel insights into the understanding of the adverse effects of environmental compounds, thereby contributing to their environmental risk assessments.

The interplay between epitranscriptomic RNA modifications and neurodegenerative disorders: Mechanistic insights and potential therapeutic strategies

Neurodegenerative disorders encompass a group of age-related conditions characterized by the gradual decline in both the structure and functionality of the central nervous system (CNS). RNA modifications, arising from the epitranscriptome or RNA-modifying protein mutations, have recently been observed to contribute significantly to neurodegenerative disorders. Specific modifications like N6-methyladenine (m6A), N1-methyladenine (m1A), 5-methylcytosine (m5C), pseudouridine and adenosine-to-inosine (A-to-I) play key roles, with their regulators serving as crucial therapeutic targets. These epitranscriptomic changes intricately control gene expression, influencing cellular functions and contributing to disease pathology. Dysregulation of RNA metabolism, affecting mRNA processing and noncoding RNA biogenesis, is a central factor in these diseases. This review underscores the complex relationship between RNA modifications and neurodegenerative disorders, emphasizing the influence of RNA modification and the epitranscriptome, exploring the function of RNA modification enzymes in neurodegenerative processes, investigating the functional consequences of RNA modifications within neurodegenerative pathways, and evaluating the potential therapeutic advancements derived from assessing the epitranscriptome.

Profiling genetically driven alternative splicing across the Indonesian archipelago

One of the regulatory mechanisms influencing the functional capacity of genes is alternative splicing (AS). Previous studies exploring the splicing landscape of human tissues have shown that AS has contributed to human biology, especially in disease progression and the immune response. Nonetheless, this phenomenon remains poorly characterized across human populations, and it is unclear how genetic and environmental variation contribute to AS. Here, we examine a set of 115 Indonesian samples from three traditional island populations spanning the genetic ancestry cline that characterizes Island Southeast Asia. We conduct a global AS analysis between islands to ascertain the degree of functionally significant AS events and their consequences. Using an event-based statistical model, we detected over 1,500 significant differential AS events across all comparisons. Additionally, we identify over 6,000 genetic variants associated with changes in splicing (splicing quantitative trait loci [sQTLs]), some of which are driven by Papuan-like genetic ancestry, and only show partial overlap with other publicly available sQTL datasets derived from other populations. Computational predictions of RNA binding activity reveal that a fraction of these sQTLs directly modulate the binding propensity of proteins involved in the splicing regulation of immune genes. Overall, these results contribute toward elucidating the role of genetic variation in shaping gene regulation in one of the most diverse regions in the world.

Alternative splicing in the genome of HPV and its regulation

Persistent infection with high-risk human papillomavirus (HR-HPV) is the main cause of cervical cancer. These chronic infections are characterized by high expression of the HPV E6 and E7 oncogenes and the absence of the L1 and L2 capsid proteins. The regulation of HPV gene expression plays a crucial role in both the viral life cycle and rare oncogenic events. Alternative splicing of HPV mRNA is a key mechanism in post-transcriptional regulation. Through alternative splicing, HPV mRNA is diversified into various splice isoforms with distinct coding potentials, encoding multiple proteins and influencing the expression of HPV genes. The spliced mRNAs derived from a donor splicing site within the E6 ORF and one of the different acceptor sites located in the early mRNA contain E6 truncated mRNAs, named E6*. E6* is one of the extensively studied splicing isoforms. However, the role of E6* proteins in cancer progression remains controversial. Here, we reviewed and compared the alternative splicing events occurring in the genomes of HR-HPV and LR-HPV. Recently, new HPV alternative splicing regulatory proteins have been continuously discovered, and we have updated the regulation of HPV alternative splicing. In addition, we summarized the functions of known splice isoforms from three aspects: anti-tumorigenic, tumorigenic, and other cancer-related functions, including not only E6*, but also E6^E7, E8^E2, and so on. Comprehending their contributions to cancer development enhances insights into the carcinogenic mechanisms of HPV and explores the potential utility of alternative splicing in the diagnosis and treatment of cervical cancer.

OpenASO: RNA Rescue - designing splice-modulating antisense oligonucleotides through community science

Splice-modulating antisense oligonucleotides (ASOs) are precision RNA-based drugs that are becoming an established modality to treat human disease. Previously, we reported the discovery of ASOs that target a novel, putative intronic RNA structure to rescue splicing of multiple pathogenic variants of F8 exon 16 that cause hemophilia A. However, the conventional approach to discovering splice-modulating ASOs is both laborious and expensive. Here, we describe an alternative paradigm that integrates data-driven RNA structure prediction and community science to discover splice-modulating ASOs. Using a splicing-deficient pathogenic variant of F8 exon 16 as a model, we show that 25% of the top-scoring molecules designed in the Eterna OpenASO challenge have a statistically significant impact on enhancing exon 16 splicing. Additionally, we show that a distinct combination of ASOs designed by Eterna players can additively enhance the inclusion of the splicing-deficient exon 16 variant. Together, our data suggests that crowdsourcing designs from a community of citizen scientists may accelerate the discovery of splice-modulating ASOs with potential to treat human disease.

Investigation of cryptic JAG1 splice variants as a cause of Alagille syndrome and performance evaluation of splice predictor tools

Haploinsufficiency of JAG1 is the primary cause of Alagille syndrome (ALGS), a rare, multisystem disorder. The identification of JAG1 intronic variants outside of the canonical splice region as well as missense variants, both of which lead to uncertain associations with disease, confuses diagnostics. Strategies to determine whether these variants affect splicing include the study of patient RNA or minigene constructs, which are not always available or can be laborious to design, as well as the utilization of computational splice prediction tools. These tools, including SpliceAI and Pangolin, use algorithms to calculate the probability that a variant results in a splice alteration, expressed as a Δ score, with higher Δ scores (>0.2 on a 0-1 scale) positively correlated with aberrant splicing. We studied the consequence of 10 putative splice variants in ALGS patient samples through RNA analysis and compared this to SpliceAI and Pangolin predictions. We identified eight variants with aberrant splicing, seven of which had not been previously validated. Combining these data with non-canonical and missense splice variants reported in the literature, we identified a predictive threshold for SpliceAI and Pangolin with high sensitivity (Δ score >0.6). Moreover, we showed reduced specificity for variants with low Δ scores (<0.2), highlighting a limitation of these tools that results in the misidentification of true splice variants. These results improve genomic diagnostics for ALGS by confirming splice effects for seven variants and suggest that the integration of splice prediction tools with RNA analysis is important to ensure accurate clinical variant classifications.

Network Modeling and Control of Dynamic Disease Pathways, Review and Perspectives

Dynamic disease pathways are a combination of complex dynamical processes among bio-molecules in a cell that leads to diseases. Network modeling of disease pathways considers disease-related bio-molecules (e.g. DNA, RNA, transcription factors, enzymes, proteins, and metabolites) and their interaction (e.g. DNA methylation, histone modification, alternative splicing, and protein modification) to study disease progression and predict therapeutic responses. These bio-molecules and their interactions are the basic elements in the study of the misregulation in the disease-related gene expression that lead to abnormal cellular responses. Gene regulatory networks, cell signaling networks, and metabolic networks are the three major types of intracellular networks for the study of the cellular responses elicited from extracellular signals. The disease-related cellular responses can be prevented or regulated by designing control strategies to manipulate these extracellular or other intracellular signals. The paper reviews the regulatory mechanisms, the dynamic models, and the control strategies for each intracellular network. The applications, limitations and the prospective for modeling and control are also discussed.

A reliable and quick method for screening alternative splicing variants for low-abundance genes

Alternative splicing (AS) is a universal phenomenon in eukaryotes, and it is still challenging to identify AS events. Several methods have been developed to identify AS events, such as expressed sequence tags (EST), microarrays and RNA-seq. However, EST has limitations in identifying low-abundance genes, while microarray and RNA-seq are high-throughput technologies, and PCR-based technology is needed for validation. To overcome the limitations of EST and shortcomings of high-throughput technologies, we established a method to identify AS events, especially for low-abundance genes, by reverse transcription (RT) PCR with gene-specific primers (GSPs) followed by nested PCR. This process includes two major steps: 1) the use of GSPs to amplify as long as the specific gene segment and 2) multiple rounds of nested PCR to screen the AS and confirm the unknown splicing variants. With this method, we successfully identified three new splicing variants, namely, GenBank Accession No. HM623886 for the bdnf gene (GenBank GeneID: 12064), GenBank Accession No. JF417977 for the trkc gene (GenBank GeneID: 18213) and GenBank Accession No. HM623888 for the glb-18 gene (GenBank GeneID: 172485). In addition to its reliability and simplicity, the method is also cost-effective and labor-intensive. In conclusion, we developed an RT-nested PCR method using gene-specific primers to efficiently identify known and novel AS variants. This approach overcomes the limitations of existing methods for detecting rare transcripts. By enabling the discovery of new isoforms, especially for low-abundance genes, this technique can aid research into aberrant splicing in disease. Future studies can apply this method to uncover AS variants involved in cancer, neurodegeneration, and other splicing-related disorders.

A systematic screen identifies Saf5 as a link between splicing and transcription in fission yeast

Splicing is an important step of gene expression regulation in eukaryotes, as there are many mRNA precursors that can be alternatively spliced in different tissues, at different cell cycle phases or under different external stimuli. We have developed several integrated fluorescence-based in vivo splicing reporter constructs that allow the quantification of fission yeast splicing in vivo on intact cells, and we have compared their splicing efficiency in a wild type strain and in a prp2-1 (U2AF65) genetic background, showing a clear dependency between Prp2 and a consensus signal at 5' splicing site (5'SS). To isolate novel genes involved in regulated splicing, we have crossed the reporter showing more intron retention with the Schizosaccharomyces pombe knock out collection. Among the candidate genes involved in the regulation of splicing, we have detected strong splicing defects in two of the mutants -Δcwf12, a member of the NineTeen Complex (NTC) and Δsaf5, a methylosome subunit that acts together with the survival motor neuron (SMN) complex in small nuclear ribonucleoproteins (snRNP) biogenesis. We have identified that strains with mutations in cwf12 have inefficient splicing, mainly when the 5'SS differs from the consensus. However, although Δsaf5 cells also have some dependency on 5'SS sequence, we noticed that when one intron of a given pre-mRNA was affected, the rest of the introns of the same pre-mRNA had high probabilities of being also affected. This observation points Saf5 as a link between transcription rate and splicing.

Nematode histone H2A variant evolution reveals diverse histories of retention and loss and evidence for conserved core-like variant histone genes

Histone variants are paralogs that replace canonical histones in nucleosomes, often imparting novel functions. However, how histone variants arise and evolve is poorly understood. Reconstruction of histone protein evolution is challenging due to large differences in evolutionary rates across gene lineages and sites. Here we used intron position data from 108 nematode genomes in combination with amino acid sequence data to find disparate evolutionary histories of the three H2A variants found in Caenorhabditis elegans: the ancient H2A.ZHTZ-1, the sperm-specific HTAS-1, and HIS-35, which differs from the canonical S-phase H2A by a single glycine-to-alanine C-terminal change. Although the H2A.ZHTZ-1 protein sequence is highly conserved, its gene exhibits recurrent intron gain and loss. This pattern suggests that specific intron sequences or positions may not be important to H2A.Z functionality. For HTAS-1 and HIS-35, we find variant-specific intron positions that are conserved across species. Patterns of intron position conservation indicate that the sperm-specific variant HTAS-1 arose more recently in the ancestor of a subset of Caenorhabditis species, while HIS-35 arose in the ancestor of Caenorhabditis and its sister group, including the genus Diploscapter. HIS-35 exhibits gene retention in some descendent lineages but gene loss in others, suggesting that histone variant use or functionality can be highly flexible. Surprisingly, we find the single amino acid differentiating HIS-35 from core H2A is ancestral and common across canonical Caenorhabditis H2A sequences. Thus, we speculate that the role of HIS-35 lies not in encoding a functionally distinct protein, but instead in enabling H2A expression across the cell cycle or in distinct tissues. This work illustrates how genes encoding such partially-redundant functions may be advantageous yet relatively replaceable over evolutionary timescales, consistent with the patchwork pattern of retention and loss of both genes. Our study shows the utility of intron positions for reconstructing evolutionary histories of gene families, particularly those undergoing idiosyncratic sequence evolution.

Time-resolved profiling of RNA binding proteins throughout the mRNA life cycle

mRNAs continually change their protein partners throughout their lifetimes, yet our understanding of mRNA-protein complex (mRNP) remodeling is limited by a lack of temporal data. Here, we present time-resolved mRNA interactome data by performing pulse metabolic labeling with photoactivatable ribonucleoside in human cells, UVA crosslinking, poly(A)+ RNA isolation, and mass spectrometry. This longitudinal approach allowed the quantification of over 700 RNA binding proteins (RBPs) across ten time points. Overall, the sequential order of mRNA binding aligns well with known functions, subcellular locations, and molecular interactions. However, we also observed RBPs with unexpected dynamics: the transcription-export (TREX) complex recruited posttranscriptionally after nuclear export factor 1 (NXF1) binding, challenging the current view of transcription-coupled mRNA export, and stress granule proteins prevalent in aged mRNPs, indicating roles in late stages of the mRNA life cycle. To systematically identify mRBPs with unknown functions, we employed machine learning to compare mRNA binding dynamics with Gene Ontology (GO) annotations. Our data can be explored at chronology.rna.snu.ac.kr.

Housekeeping U1 snRNA facilitates antiviral innate immunity by promoting TRIM25-mediated RIG-I activation

U1 small nuclear RNA (snRNA) is an abundant and evolutionarily conserved 164-nucleotide RNA species that functions in pre-mRNA splicing, and it is considered to be a housekeeping non-coding RNA. However, the role of U1 snRNA in regulating host antiviral immunity remains largely unexplored. Here, we find that RNVU1-18, a U1 pseudogene, is significantly upregulated in the host infected with RNA viruses, including influenza and respiratory syncytial virus. Overexpression of U1 snRNA protects cells against RNA viruses, while knockdown of U1 snRNA leads to more viral burden in vitro and in vivo. Knockout of RNVU1-18 is sufficient to impair the type I interferon-dependent antiviral innate immunity. U1 snRNA is required to fully activate the retinoic acid-inducible gene I (RIG-I)-dependent antiviral signaling, since it interacts with tripartite motif 25 (TRIM25) and enhances the RIG-I-TRIM25 interaction to trigger K63-linked ubiquitination of RIG-I. Our study reveals the important role of housekeeping U1 snRNA in regulating host antiviral innate immunity and restricting RNA virus infection.

HOXA1 3'UTR Methylation Is a Potential Prognostic Biomarker in Oral Squamous cell Carcinoma

BACKGROUND

HOXA1 is a prognostic marker and a potential predictive biomarker for radioresistance in head and neck tumors. Its overexpression has been associated with promoter methylation and a worse prognosis in oral squamous cell carcinoma (OSCC) patients. However, opposite outcomes are also described. The effect of the methylation of this gene on different gene regions, other than the promoter, remains uncertain. We investigated the methylation profile at different genomic regions of HOXA1 in OSCC and correlated differentially methylated CpG sites with clinicopathological data.

METHODS

The HOXA1 DNA methylation status was evaluated by analyzing data from The Cancer Genome Atlas and three Gene Expression Omnibus datasets. Significant differentially methylated CpG sites were considered with a |∆β| ≥ 0.10 and a Bonferroni-corrected p-value < 0.01. Differentially methylated CpGs were validated by pyrosequencing using two independent cohorts of 15 and 47 OSCC patients, respectively.

RESULTS

Compared to normal tissues, we found significantly higher DNA methylation levels in the 3'UTR region of HOXA1 in OSCC. Higher methylation levels in tumor samples were positively correlated with smoking habits and patients' overall survival.

CONCLUSIONS

Our findings suggest that HOXA1 gene body methylation is a promising prognostic biomarker for OSCC with potential clinical applications in patient monitoring.

CRISPR activation to characterize splice-altering variants in easily accessible cells

Genetic variants that affect mRNA splicing are a major cause of hereditary disorders, but the spliceogenicity of variants is challenging to predict. RNA diagnostics of clinically accessible tissues enable rapid functional characterization of splice-altering variants within their natural genetic context. However, this analysis cannot be offered to all individuals as one in five human disease genes are not expressed in easily accessible cell types. To overcome this problem, we have used CRISPR activation (CRISPRa) based on a dCas9-VPR mRNA-based delivery platform to induce expression of the gene of interest in skin fibroblasts from individuals with suspected monogenic disorders. Using this ex vivo splicing assay, we characterized the splicing patterns associated with germline variants in the myelin protein zero gene (MPZ), which is exclusively expressed in Schwann cells of the peripheral nerves, and the spastin gene (SPAST), which is predominantly expressed in the central nervous system. After overnight incubation, CRISPRa strongly upregulated MPZ and SPAST transcription in skin fibroblasts, which enabled splice variant profiling using reverse transcription polymerase chain reaction, next-generation sequencing, and long-read sequencing. Our investigations show proof of principle of a promising genetic diagnostic tool that involves CRISPRa to activate gene expression in easily accessible cells to study the functional impact of genetic variants. The procedure is easy to perform in a diagnostic laboratory with equipment and reagents all readily available.

Sox9 regulates alternative splicing and pancreatic beta cell function

Despite significant research, mechanisms underlying the failure of islet beta cells that result in type 2 diabetes (T2D) are still under investigation. Here, we report that Sox9, a transcriptional regulator of pancreas development, also functions in mature beta cells. Our results show that Sox9-depleted rodent beta cells have defective insulin secretion, and aging animals develop glucose intolerance, mimicking the progressive degeneration observed in T2D. Using genome editing in human stem cells, we show that beta cells lacking SOX9 have stunted first-phase insulin secretion. In human and rodent cells, loss of Sox9 disrupts alternative splicing and triggers accumulation of non-functional isoforms of genes with key roles in beta cell function. Sox9 depletion reduces expression of protein-coding splice variants of the serine-rich splicing factor arginine SRSF5, a major splicing enhancer that regulates alternative splicing. Our data highlight the role of SOX9 as a regulator of alternative splicing in mature beta cell function.

SpliceProt 2.0: A Sequence Repository of Human, Mouse, and Rat Proteoforms

SpliceProt 2.0 is a public proteogenomics database that aims to list the sequence of known proteins and potential new proteoforms in human, mouse, and rat proteomes. This updated repository provides an even broader range of computationally translated proteins and serves, for example, to aid with proteomic validation of splice variants absent from the reference UniProtKB/SwissProt database. We demonstrate the value of SpliceProt 2.0 to predict orthologous proteins between humans and murines based on transcript reconstruction, sequence annotation and detection at the transcriptome and proteome levels. In this release, the annotation data used in the reconstruction of transcripts based on the methodology of ternary matrices were acquired from new databases such as Ensembl, UniProt, and APPRIS. Another innovation implemented in the pipeline is the exclusion of transcripts predicted to be susceptible to degradation through the NMD pathway. Taken together, our repository and its applications represent a valuable resource for the proteogenomics community.

ZCCHC17 Modulates Neuronal RNA Splicing and Supports Cognitive Resilience in Alzheimer's Disease

ZCCHC17 is a putative master regulator of synaptic gene dysfunction in Alzheimer's disease (AD), and ZCCHC17 protein declines early in AD brain tissue, before significant gliosis or neuronal loss. Here, we investigate the function of ZCCHC17 and its role in AD pathogenesis using data from human autopsy tissue (consisting of males and females) and female human cell lines. Co-immunoprecipitation (co-IP) of ZCCHC17 followed by mass spectrometry analysis in human iPSC-derived neurons reveals that ZCCHC17's binding partners are enriched for RNA-splicing proteins. ZCCHC17 knockdown results in widespread RNA-splicing changes that significantly overlap with splicing changes found in AD brain tissue, with synaptic genes commonly affected. ZCCHC17 expression correlates with cognitive resilience in AD patients, and we uncover an APOE4-dependent negative correlation of ZCCHC17 expression with tangle burden. Furthermore, a majority of ZCCHC17 interactors also co-IP with known tau interactors, and we find a significant overlap between alternatively spliced genes in ZCCHC17 knockdown and tau overexpression neurons. These results demonstrate ZCCHC17's role in neuronal RNA processing and its interaction with pathology and cognitive resilience in AD, and suggest that the maintenance of ZCCHC17 function may be a therapeutic strategy for preserving cognitive function in the setting of AD pathology.

Long-read single-cell sequencing reveals expressions of hypermutation clusters of isoforms in human liver cancer cells

The protein diversity of mammalian cells is determined by arrays of isoforms from genes. Genetic mutation is essential in species evolution and cancer development. Accurate long-read transcriptome sequencing at single-cell level is required to decipher the spectrum of protein expressions in mammalian organisms. In this report, we developed a synthetic long-read single-cell sequencing technology based on LOOPSeq technique. We applied this technology to analyze 447 transcriptomes of hepatocellular carcinoma (HCC) and benign liver from an individual. Through Uniform Manifold Approximation and Projection analysis, we identified a panel of mutation mRNA isoforms highly specific to HCC cells. The evolution pathways that led to the hyper-mutation clusters in single human leukocyte antigen molecules were identified. Novel fusion transcripts were detected. The combination of gene expressions, fusion gene transcripts, and mutation gene expressions significantly improved the classification of liver cancer cells versus benign hepatocytes. In conclusion, LOOPSeq single-cell technology may hold promise to provide a new level of precision analysis on the mammalian transcriptome.

The splicing factor DHX38 enables retinal development through safeguarding genome integrity

DEAH-Box Helicase 38 (DHX38) is a pre-mRNA splicing factor and also a disease-causing gene of autosomal recessive retinitis pigmentosa (arRP). The role of DHX38 in the development and maintenance of the retina remains largely unknown. In this study, by using the dhx38 knockout zebrafish model, we demonstrated that Dhx38 deficiency causes severe differentiation defects and apoptosis of retinal progenitor cells (RPCs) through disrupted mitosis and increased DNA damage. Furthermore, we found a significant accumulation of R-loops in the dhx38-deficient RPCs and human cell lines. Finally, we found that DNA replication stress is the prerequisite for R-loop-induced DNA damage in the DHX38 knockdown cells. Taken together, our study demonstrates a necessary role of DHX38 in the development of retina and reveals a DHX38/R-loop/replication stress/DNA damage regulatory axis that is relatively independent of the known functions of DHX38 in mitosis control.

Heterogeneous splicing patterns resulting from KIF5A variants associated with amyotrophic lateral sclerosis

Single-nucleotide variants (SNVs) in the gene encoding Kinesin Family Member 5A (KIF5A), a neuronal motor protein involved in anterograde transport along microtubules, have been associated with amyotrophic lateral sclerosis (ALS). ALS is a rapidly progressive and fatal neurodegenerative disease that primarily affects the motor neurons. Numerous ALS-associated KIF5A SNVs are clustered near the splice-site junctions of the penultimate exon 27 and are predicted to alter the carboxy-terminal (C-term) cargo-binding domain of KIF5A. Mis-splicing of exon 27, resulting in exon exclusion, is proposed to be the mechanism by which these SNVs cause ALS. Whether all SNVs proximal to exon 27 result in exon exclusion is unclear. To address this question, we designed an in vitro minigene splicing assay in human embryonic kidney 293 cells, which revealed heterogeneous site-specific effects on splicing: only 5' splice-site (5'ss) SNVs resulted in exon skipping. We also quantified splicing in select clustered, regularly interspaced, short palindromic repeats-edited human stem cells, differentiated to motor neurons, and in neuronal tissues from a 5'ss SNV knock-in mouse, which showed the same result. Moreover, the survival of representative 3' splice site, 5'ss, and truncated C-term variant KIF5A (v-KIF5A) motor neurons was severely reduced compared with wild-type motor neurons, and overt morphological changes were apparent. While the total KIF5A mRNA levels were comparable across the cell lines, the total KIF5A protein levels were decreased for v-KIF5A lines, suggesting an impairment of protein synthesis or stability. Thus, despite the heterogeneous effect on ribonucleic acid splicing, KIF5A SNVs similarly reduce the availability of the KIF5A protein, leading to axonal transport defects and motor neuron pathology.

A partial form of AIRE deficiency underlies a mild form of autoimmune polyendocrine syndrome type 1

Autoimmune polyendocrine syndrome type 1 (APS-1) is caused by mutations in the autoimmune regulator (AIRE) gene. Most patients present with severe chronic mucocutaneous candidiasis and organ-specific autoimmunity from early childhood, but the clinical picture is highly variable. AIRE is crucial for negative selection of T cells, and scrutiny of different patient mutations has previously highlighted many of its molecular mechanisms. In patients with a milder adult-onset phenotype sharing a mutation in the canonical donor splice site of intron 7 (c.879+1G>A), both the predicted altered splicing pattern with loss of exon 7 (AireEx7-/-) and normal full-length AIRE mRNA were found, indicating leaky rather than abolished mRNA splicing. Analysis of a corresponding mouse model demonstrated that the AireEx7-/- mutant had dramatically impaired transcriptional capacity of tissue-specific antigens in medullary thymic epithelial cells but still retained some ability to induce gene expression compared with the complete loss-of-function AireC313X-/- mutant. Our data illustrate an association between AIRE activity and the severity of autoimmune disease, with implications for more common autoimmune diseases associated with AIRE variants, such as primary adrenal insufficiency, pernicious anemia, type 1 diabetes, and rheumatoid arthritis.

Pathogenicity analysis and splicing rescue of a classical splice site variant (c.1343+1G>T) of CNOT1 gene associated with neurodevelopmental disorders

Mutations in the CNOT1 gene lead to an incurable rare neurological disorder mainly manifested as a clinical spectrum of intellectual disability, developmental delay, seizures, and behavioral problems. In this study, we investigated a classical splice site variant of CNOT1 (c.1343+1G>T) associated with neurodevelopmental disorders, which was a master regulator, orchestrating gene expression, RNA deadenylation, and protein ubiquitination. To link CNOT1 dysfunction with the neurodevelopmental phenotype observed in a patient, in vitro minigene assay was used to verify the effect of CNOT1 gene splice site variant c.1343+1G>T on mRNA splicing. We also explored the impact of transient transfection introducing modified U1 snRNA on correcting the splicing variant. Through minigene expression in mammalian cells, we demonstrated that the variant induced complete exon 12 skipping, which explained the patient's clinical condition and provided additional genetic diagnosis evidence for the clinical significance of the variant. Moreover, we confirmed that the aberrant splice pattern could be partially corrected by the modified U1 snRNA at the mRNA level, which provided strong evidence for the therapeutic potential of modified U1 snRNA in neutralizing the hazardous effect of incorrect splicing patterns.

Deletions of singular U1 snRNA gene significantly interfere with transcription and 3'-end mRNA formation

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1-2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

Birth and death view of DNA, RNA, and proteins

The potential of cells to guide their genome and configure genes to express at a given time and in response to specific stimuli is pivotal to regulate cellular processes such as tissue differentiation, organogenesis, organismal development, homeostasis, and disease. In this review, we focus on the diverse mechanisms involved in DNA replication and its degradation, mRNA synthesis, and associated regulation such as RNA capping, splicing, tailing, and export. mRNA turnover including Decapping, deadenylation, RNA interference, and Nonsense mediated mRNA decay followed by protein translation, post-translational modification, and protein turnover. We highlight recent advances in understanding the complex series of molecular mechanisms responsible for the remarkable cellular regulatory mechanisms.

Alternative splicing events and function in the tumor microenvironment: New opportunities and challenges

Alternative splicing controls gene expression at the transcriptional level, producing structurally and functionally distinct protein heterodimers. Aberrant alternative splicing greatly affects cell development and plays an important role in the invasion and metastasis of many types of cancer. Recently, it has been shown that alternative splicing can alter the tumor microenvironment and regulate processes such as remodeling, immunity, and inflammation in the tumor microenvironment. However, there is no comprehensive literature review of the complex relationship between alternative splicing and the tumor microenvironment. Therefore, this review aims to collect all the latest data on this topic and provide a new perspective on the therapeutic and potential prognostic markers of cancer.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

Long-read single-cell sequencing reveals expressions of hypermutation clusters of isoforms in human liver cancer cells

The protein diversity of mammalian cells is determined by arrays of isoforms from genes. Genetic mutation is essential in species evolution and cancer development. Accurate Long-read transcriptome sequencing at single-cell level is required to decipher the spectrum of protein expressions in mammalian organisms. In this report, we developed a synthetic long-read single-cell sequencing technology based on LOOPseq technique. We applied this technology to analyze 447 transcriptomes of hepatocellular carcinoma (HCC) and benign liver from an individual. Through Uniform Manifold Approximation and Projection (UMAP) analysis, we identified a panel of mutation mRNA isoforms highly specific to HCC cells. The evolution pathways that led to the hyper-mutation clusters in single human leukocyte antigen (HLA) molecules were identified. Novel fusion transcripts were detected. The combination of gene expressions, fusion gene transcripts, and mutation gene expressions significantly improved the classification of liver cancer cells versus benign hepatocytes. In conclusion, LOOPseq single-cell technology may hold promise to provide a new level of precision analysis on the mammalian transcriptome.

Sex-specific splicing occurs genome-wide during early Drosophila embryogenesis

Sex-specific splicing is an essential process that regulates sex determination and drives sexual dimorphism. Yet, how early in development widespread sex-specific transcript diversity occurs was unknown because it had yet to be studied at the genome-wide level. We use the powerful Drosophila model to show that widespread sex-specific transcript diversity occurs early in development, concurrent with zygotic genome activation. We also present a new pipeline called time2Splice to quantify changes in alternative splicing over time. Furthermore, we determine that one of the consequences of losing an essential maternally deposited pioneer factor called CLAMP (chromatin-linked adapter for MSL proteins) is altered sex-specific splicing of genes involved in diverse biological processes that drive development. Overall, we show that sex-specific differences in transcript diversity exist even at the earliest stages of development..

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans

Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.

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.

The implications of alternative pre-mRNA splicing in cell signal transduction

Cells produce multiple mRNAs through alternative splicing, which ensures proteome diversity. Because most human genes undergo alternative splicing, key components of signal transduction pathways are no exception. Cells regulate various signal transduction pathways, including those associated with cell proliferation, development, differentiation, migration, and apoptosis. Since proteins produced through alternative splicing can exhibit diverse biological functions, splicing regulatory mechanisms affect all signal transduction pathways. Studies have demonstrated that proteins generated by the selective combination of exons encoding important domains can enhance or attenuate signal transduction and can stably and precisely regulate various signal transduction pathways. However, aberrant splicing regulation via genetic mutation or abnormal expression of splicing factors negatively affects signal transduction pathways and is associated with the onset and progression of various diseases, including cancer. In this review, we describe the effects of alternative splicing regulation on major signal transduction pathways and highlight the significance of alternative splicing.

ZCCHC17 modulates neuronal RNA splicing and supports cognitive resilience in Alzheimer's disease

ZCCHC17 is a putative master regulator of synaptic gene dysfunction in Alzheimer's Disease (AD), and ZCCHC17 protein declines early in AD brain tissue, before significant gliosis or neuronal loss. Here, we investigate the function of ZCCHC17 and its role in AD pathogenesis. Co-immunoprecipitation of ZCCHC17 followed by mass spectrometry analysis in human iPSC-derived neurons reveals that ZCCHC17's binding partners are enriched for RNA splicing proteins. ZCCHC17 knockdown results in widespread RNA splicing changes that significantly overlap with splicing changes found in AD brain tissue, with synaptic genes commonly affected. ZCCHC17 expression correlates with cognitive resilience in AD patients, and we uncover an APOE4 dependent negative correlation of ZCCHC17 expression with tangle burden. Furthermore, a majority of ZCCHC17 interactors also co-IP with known tau interactors, and we find significant overlap between alternatively spliced genes in ZCCHC17 knockdown and tau overexpression neurons. These results demonstrate ZCCHC17's role in neuronal RNA processing and its interaction with pathology and cognitive resilience in AD, and suggest that maintenance of ZCCHC17 function may be a therapeutic strategy for preserving cognitive function in the setting of AD pathology.

In Vivo Efficacy and Safety Evaluations of Therapeutic Splicing Correction Using U1 snRNA in the Mouse Retina

Efficacy and safety considerations constitute essential steps during development of in vivo gene therapies. Herein, we evaluated efficacy and safety of splice factor-based treatments to correct mutation-induced splice defects in an Opa1 mutant mouse line. We applied adeno-associated viruses to the retina. The viruses transduced retinal cells with an engineered U1 snRNA splice factor designed to correct the Opa1 splice defect. We found the treatment to be efficient in increasing wild-type Opa1 transcripts. Correspondingly, Opa1 protein levels increased significantly in treated eyes. Measurements of retinal morphology and function did not reveal therapy-related side-effects supporting the short-term safety of the treatment. Alterations of potential off-target genes were not detected. Our data suggest that treatments of splice defects applying engineered U1 snRNAs represent a promising in vivo therapeutic approach. The therapy increased wild-type Opa1 transcripts and protein levels without detectable morphological, functional or genetic side-effects in the mouse eye. The U1 snRNA-based therapy can be tailored to specific disease gene mutations, hence, raising the possibility of a wider applicability of this promising technology towards treatment of different inherited retinal diseases.

MBNL-dependent impaired development within the neuromuscular system in myotonic dystrophy type 1

AIMS

Myotonic dystrophy type I (DM1) is one of the most frequent muscular dystrophies in adults. Although DM1 has long been considered mainly a muscle disorder, growing evidence suggests the involvement of peripheral nerves in the pathogenicity of DM1 raising the question of whether motoneurons (MNs) actively contribute to neuromuscular defects in DM1.

METHODS

By using micropatterned 96-well plates as a coculture platform, we generated a functional neuromuscular model combining DM1 and muscleblind protein (MBNL) knock-out human-induced pluripotent stem cells-derived MNs and human healthy skeletal muscle cells.

RESULTS

This approach led to the identification of presynaptic defects which affect the formation or stability of the neuromuscular junction at an early developmental stage. These neuropathological defects could be reproduced by the loss of RNA-binding MBNL proteins, whose loss of function in vivo is associated with muscular defects associated with DM1. These experiments indicate that the functional defects associated with MNs can be directly attributed to MBNL family proteins. Comparative transcriptomic analyses also revealed specific neuronal-related processes regulated by these proteins that are commonly misregulated in DM1.

CONCLUSIONS

Beyond the application to DM1, our approach to generating a robust and reliable human neuromuscular system should facilitate disease modelling studies and drug screening assays.

Nuclear envelope transmembrane proteins involved in genome organization are misregulated in myotonic dystrophy type 1 muscle

Myotonic dystrophy type 1 is a multisystemic disorder with predominant muscle and neurological involvement. Despite a well described pathomechanism, which is primarily a global missplicing due to sequestration of RNA-binding proteins, there are still many unsolved questions. One such question is the disease etiology in the different affected tissues. We observed alterations at the nuclear envelope in primary muscle cell cultures before. This led us to reanalyze a published RNA-sequencing dataset of DM1 and control muscle biopsies regarding the misregulation of NE proteins. We could identify several muscle NE protein encoding genes to be misregulated depending on the severity of the muscle phenotype. Among these misregulated genes were NE transmembrane proteins (NETs) involved in nuclear-cytoskeletal coupling as well as genome organization. For selected genes, we could confirm that observed gene-misregulation led to protein expression changes. Furthermore, we investigated if genes known to be under expression-regulation by genome organization NETs were also misregulated in DM1 biopsies, which revealed that misregulation of two NETs alone is likely responsible for differential expression of about 10% of all genes being differentially expressed in DM1. Notably, the majority of NETs identified here to be misregulated in DM1 muscle are mutated in Emery-Dreifuss muscular dystrophy or clinical similar muscular dystrophies, suggesting a broader similarity on the molecular level for muscular dystrophies than anticipated. This shows not only the importance of muscle NETs in muscle health and disease, but also highlights the importance of the NE in DM1 disease progression.

Gene body methylation in cancer: molecular mechanisms and clinical applications

DNA methylation is an important epigenetic mechanism that regulates gene expression. To date, most DNA methylation studies have focussed on CpG islands in the gene promoter region, and the mechanism of methylation and the regulation of gene expression after methylation have been clearly elucidated. However, genome-wide methylation studies have shown that DNA methylation is widespread not only in promoters but also in gene bodies. Gene body methylation is widely involved in the expression regulation of many genes and is closely related to the occurrence and progression of malignant tumours. This review focusses on the formation of gene body methylation patterns, its regulation of transcription, and its relationship with tumours, providing clues to explore the mechanism of gene body methylation in regulating gene transcription and its significance and application in the field of oncology.

mTOR Contributes to the Proteome Diversity through Transcriptome-Wide Alternative Splicing

The mammalian target of rapamycin (mTOR) pathway is crucial in energy metabolism and cell proliferation. Previously, we reported transcriptome-wide 3′-untranslated region (UTR) shortening by alternative polyadenylation upon mTOR activation and its impact on the proteome. Here, we further interrogated the mTOR-activated transcriptome and found that hyperactivation of mTOR promotes transcriptome-wide exon skipping/exclusion, producing short isoform transcripts from genes. This widespread exon skipping confers multifarious regulations in the mTOR-controlled functional proteomics: AS in coding regions widely affects the protein length and functional domains. They also alter the half-life of proteins and affect the regulatory post-translational modifications. Among the RNA processing factors differentially regulated by mTOR signaling, we found that SRSF3 mechanistically facilitates exon skipping in the mTOR-activated transcriptome. This study reveals a role of mTOR in AS regulation and demonstrates that widespread AS is a multifaceted modulator of the mTOR-regulated functional proteome.