Functional Splicing Research Articles

Large-scale map of RNA-binding protein interactomes across the mRNA life cycle

mRNAs interact with RNA-binding proteins (RBPs) throughout their processing and maturation. While efforts have assigned RBPs to RNA substrates, less exploration has leveraged protein-protein interactions (PPIs) to study proteins in mRNA life-cycle stages. We generated an RNA-aware, RBP-centric PPI map across the mRNA life cycle in human cells by immunopurification-mass spectrometry (IP-MS) of ∼100 endogenous RBPs with and without RNase, augmented by size exclusion chromatography-mass spectrometry (SEC-MS). We identify 8,742 known and 20,802 unreported interactions between 1,125 proteins and determine that 73% of the IP-MS-identified interactions are RNA regulated. Our interactome links many proteins, some with unknown functions, to specific mRNA life-cycle stages, with nearly half associated with multiple stages. We demonstrate the value of this resource by characterizing the splicing and export functions of enhancer of rudimentary homolog (ERH), and by showing that small nuclear ribonucleoprotein U5 subunit 200 (SNRNP200) interacts with stress granule proteins and binds cytoplasmic RNA differently during stress.

Alternative splicing of Mff regulates AMPK-mediated phosphorylation, mitochondrial fission and antiviral response

Mitochondrial morphology and function change dynamically in response to intracellular signaling and the surrounding environment. The mitochondrial fission factor Mff, which localizes to the outer mitochondrial membrane, mediates not only mitochondrial fission by recruiting the dynamin-related GTPase Drp1 to mitochondrial fission sites but also the double-stranded RNA-induced antiviral response on mitochondria through mitochondrial antiviral signaling (MAVS). Mff is reported to be regulated by AMP-activated protein kinase (AMPK)-mediated protein phosphorylation and alternative pre-mRNA splicing; however, the relationships among RNA splicing, phosphorylation, and multiple functions of Mff have not been fully understood. Here, we showed that mouse Mff has a tissue-specific splicing pattern, and at least eight Mff splice isoforms were expressed in mouse embryonic fibroblasts (MEFs). We introduced single Mff isoforms into Mff knockout MEFs and found that insertion of exon 6 just after the phosphorylation site, by the alternative splicing, reduced its phosphorylation by AMPK and its functions in mitochondrial fission and the antiviral response. In addition, the underlying mechanism repressing these functions was independent of phosphorylation. These results indicate that multiple functions of Mff on mitochondria are regulated by AMPK-mediated phosphorylation and alternative splicing, under the control of energy metabolism and cellular differentiation.

Maize splicing-mediated mRNA surveillance impeded by sugarcane mosaic virus-coded pathogenic protein NIa-Pro.

The eukaryotic mRNA surveillance pathway, a pivotal guardian of mRNA fidelity, stands at the nexus of diverse biological processes, including antiviral immunity. Despite the recognized function of splicing factors on mRNA fate, the intricate interplay shaping the mRNA surveillance pathway remains elusive. We illustrate that the conserved splicing factor U2 snRNP auxiliary factor large subunit B (U2AF65B) modulates splicing of mRNA surveillance complex, contributing to transcriptomic homeostasis in maize. The functionality of the mRNA surveillance pathway requires ZmU2AF65B-mediated normal splicing of upstream frameshift 3 (ZmUPF3) pre-mRNA, encoding a core factor in this pathway. Intriguingly, sugarcane mosaic virus (SCMV)-coded nuclear inclusion protein a protease (NIa-Pro) hinders the splicing function of ZmU2AF65B. Furthermore, NIa-Pro disrupts ZmU2AF65B binding to ZmUPF3 pre-mRNA, leading to dysregulated splicing of ZmUPF3 transcripts and, consequently, impairing mRNA surveillance, thus facilitating viral infection. Together, this study establishes that splicing governs the mRNA surveillance pathway and identifies a pathogenic protein capable of disrupting this regulation to compromise RNA immunity.

Unveiling axolotl transcriptome for tissue regeneration with high-resolution annotation via long-read sequencing

Axolotls are known for their remarkable regeneration ability. Exploring their transcriptome provides insight into regenerative mechanisms. However, the current annotation of the axolotl transcriptome is limited, leaving the role of unannotated transcripts in regeneration unknown. To discourse this challenge, we exploited long-read sequencing technology, which enables direct observation of full-length RNA transcripts, greatly enhancing the coverage and accuracy of axolotl transcriptome annotation. By utilizing this method, we identified 222 novel gene loci and 4775 novel transcripts, which were quantified using short-read sequencing data. Through the inclusive analysis, we discovered novel homologs, potential functional proteins, noncoding RNAs, and alternative splicing events in key regeneration pathways. In particular, we identified novel transcripts with high protein-coding potential implicated in cell cycle regulation and musculoskeletal development, and regeneration were identified. Interestingly, alternative splice variants were also detected across diverse pathways critical to regeneration. This specifies that these novel transcripts potentially play vital roles underpinning the robust regenerative capacities of axolotls. Single-cell transcriptomic analysis further revealed these isoforms to predominantly exist in axolotl limb chondrocytes and mature tissue cell populations. Overall, the findings significantly advanced consideration of the axolotl transcriptome and provided a new perspective for understanding the mechanisms of regenerative abilities of axolotls.

Identification of Key Prognostic Alternative Splicing Events of Costimulatory Molecule-Related Genes in Colon Cancer.

This study aimed to explore the key alternative splicing events in costimulatory molecule-related genes in colon cancer and to determine their correlation with prognosis. Gene expression RNA-sequencing data, clinical data, and SpliceSeq data of colon cancer were obtained from The Cancer Genome Atlas. Differentially expressed alternative splicing events in genes were identified, Followed by correlation analysis of genes corresponding to differentially expressed alternative splicing events with costimulatory molecule-related genes. Survival analysis was conducted using differentially expressed alternative splicing events in these genes and a prognostic model was constructed. Functional enrichment, proteinprotein interaction network, and splicing factor analyses were performed. In total, 6504 differentially expressed alternative splicing events in 3949 genes were identified between tumor and normal tissues. Correlation analysis revealed 3499 differentially expressed alternative splicing events in 2168 costimulatory molecule-related genes. Moreover, 328 differentially expressed alternative splicing events in 288 costimulatory molecule-related genes were associated with overall survival. The prognostic models constructed using these showed considerable power in predicting survival. The ubiquitin A-52 residue ribosomal protein fusion product 1 and ribosomal protein S9 were the hub nodes in the protein-protein interaction network. Furthermore, one splicing factor, splicing factor proline and glutamine-rich, was significantly associated with patient prognosis. Four splicing factor-alternative splicing pairs were obtained from four alternative splicing events in three genes: TBC1 domain family member 8 B, complement factor H, and mitochondrial fission 1. The identified differentially expressed alternative splicing events of costimulatory molecule-related genes may be used to predict patient prognosis and immunotherapy responses in colon cancer.

U1 snRNA interactions with deep intronic sequences regulate splicing of multiple exons of spinal muscular atrophy genes.

The U1 small nuclear RNA (snRNA) forms ribonucleoprotein particles (RNPs) such as U1 snRNP and U1-TAF15 snRNP. U1 snRNP is one of the most studied RNPs due to its critical role in pre-mRNA splicing in defining the 5' splice site (5'ss) of every exon through direct interactions with sequences at exon/intron junctions. Recent reports support the role of U1 snRNP in all steps of transcription, namely initiation, elongation, and termination. Functions of U1-TAF15 snRNP are less understood, though it associates with the transcription machinery and may modulate pre-mRNA splicing by interacting with the 5'ss and/or 5'ss-like sequences within the pre-mRNA. An anti-U1 antisense oligonucleotide (ASO) that sequesters the 5' end of U1 snRNA inhibits the functions of U1 snRNP, including transcription and splicing. However, it is not known if the inhibition of U1 snRNP influences post-transcriptional regulation of pre-mRNA splicing through deep intronic sequences. We examined the effect of an anti-U1 ASO that sequesters the 5' end of U1 snRNA on transcription and splicing of all internal exons of the spinal muscular atrophy (SMA) genes, SMN1 and SMN2. Our study was enabled by the employment of a multi-exon-skipping detection assay (MESDA) that discriminates against prematurely terminated transcripts. We employed an SMN2 super minigene to determine if anti-U1 ASO differently affects splicing in the context of truncated introns. We observed substantial skipping of multiple internal exons of SMN1 and SMN2 triggered by anti-U1 treatment. Suggesting a role for U1 snRNP in interacting with deep intronic sequences, early exons of the SMN2 super minigene with truncated introns were resistant to anti-U1 induced skipping. Consistently, overexpression of engineered U1 snRNAs targeting the 5'ss of early SMN1 and SMN2 exons did not prevent exon skipping caused by anti-U1 treatment. Our results uncover a unique role of the U1 snRNA-associated RNPs in splicing regulation executed through deep intronic sequences. Findings are significant for developing novel therapies for SMA based on deep intronic targets.

CUTS RNA Biosensor for the Real-Time Detection of TDP-43 Loss-of-Function.

Given the mounting evidence implicating TDP-43 dysfunction in several neurodegenerative diseases, there is a pressing need to establish accessible tools to sense and quantify TDP-43 loss-of-function (LOF). These tools are crucial for assessing potential disease contributors and exploring therapeutic candidates in TDP-43 proteinopathies. Here, we develop a sensitive and accurate real-time sensor for TDP-43 LOF: the CUTS (CFTR UNC13A TDP-43 Loss-of-Function) system. This system combines previously reported cryptic exons regulated by TDP-43 with a reporter, enabling the tracking of TDP-43 LOF through live microscopy and RNA/protein-based assays. We demonstrate CUTS' effectiveness in detecting LOF caused by TDP-43 mislocalization and RNA binding dysfunction, and pathological aggregation. Our results highlight the sensitivity and accuracy of the CUTS system in detecting and quantifying TDP-43 LOF, opening avenues to explore unknown TDP-43 interactions that regulate its function. In addition, by replacing the fluorescent tag in the CUTS system with the coding sequence for TDP-43, we show significant recovery of its function under TDP-43 LOF conditions, highlighting CUTS' potential for self-regulating gene therapy applications. In summary, CUTS represents a versatile platform for evaluating TDP-43 LOF in real-time and advancing gene-replacement therapies in neurodegenerative diseases associated with TDP-43 dysfunction.

Two homozygous adjacent novel missense mutations in DYSF gene caused dysferlinopathy due to splicing abnormalities.

Background: Dysferlinopathy is an autosomal recessive disorder caused by mutations in the DYSF gene. This study reported two homozygous adjacent missense mutations in the DYSF gene, presenting clinically with bilateral lower limb weakness and calf swelling. Two homozygous adjacent missense mutations in the DYSF gene may be associated with the development of dysferlinopathy, but the exact mechanism needs further investigation. Methods: A retrospective analysis of clinical data from a dysferlinopathy-affected family was conducted. Peripheral blood samples were collected from members of this family for whole-exome sequencing (WES) and copy number variation analysis. Sanger sequencing was employed to confirm potential pathogenic variants. The Human Splicing Finder, SpliceAI, and varSEAK database were used to predict the effect of mutations on splicing function. The pathogenic mechanism of aberrant splicing in dysferlinopathy due to two homozygous adjacent missense mutations in the DYSF gene was determined by an in vivo splicing assay and an in vitro minigene assay. Results: The proband was a 42-year-old woman who presented with weakness of the lower limbs for 2 years and edema of the lower leg. Two homozygous DYSF variants, c.5628C>A p. D1876E and c.5633A>T p. Y1878F, were identified in the proband. Bioinformatics databases suggested that the mutation c.5628C>A of DYSF had no significant impact on splicing signals. Human Splicing Finder Version 2.4.1 suggested that the c.5633A>T of DYSF mutation caused alteration of auxiliary sequences and significant alteration of the ESE/ESS motif ratio. VarSEAK and SpliceAI suggested that the c.5633A>T of DYSF mutation had no splicing effect. Both an in vivo splicing assay and an in vitro minigene assay showed two adjacent mutations: c.5628C>A p. D1876E and c.5633A>T p. Y1878F in the DYSF gene leading to an Exon50 jump that resulted in a 32-aa amino acid deletion within the protein. Point mutation c.5628C>A p. D1876E in the DYSF gene affected splicing in vitro, while point mutation c.5633A>T p. Y1878F in the DYSF gene did not affect splicing function. Conclusion: This study confirmed for the first time that two homozygous mutations of DYSF were associated with the occurrence of dysferlinopathy. The c.5628C>A p. D1876E mutation in DYSF affected the splicing function and may be one of the contributing factors to the pathogenicity.

Parsing the roles of DExD-box proteins DDX39A and DDX39B in alternative RNA splicing.

DExD-box RNA proteins DDX39A and DDX39B are highly homologous paralogs that are conserved in vertebrates. They are required for energy-driven reactions involved in RNA processing. Although we have some understanding of how their functions overlap in RNA nuclear export, our knowledge of whether or not these proteins have specific or redundant functions in RNA splicing is limited. Our previous work has shown that DDX39B is responsible for regulating the splicing of important immune transcripts IL7R and FOXP3. In this study, we aimed to investigate whether DDX39A, a highly homologous paralog of DDX39B, plays a similar role in regulating alternative RNA splicing. We find that DDX39A and DDX39B have significant redundancy in their gene targets, but there are targets that uniquely require one or the other paralog. For instance, DDX39A is incapable of complementing defective splicing of IL7R exon 6 when DDX39B is depleted. This exon and other cassette exons that specifically depend on DDX39B have U-poor/C-rich polypyrimidine tracts in the upstream intron and this variant polypyrimidine tract is required for DDX39B dependency. This study provides evidence that despite a high degree of functional redundancy, DDX39A and DDX39B are selectively required for the splicing ofspecific pre-mRNAs.

Assessing the Impact of Novel BRCA1 Exon 11 Variants on Pre-mRNA Splicing.

Our study focused on assessing the effects of three newly identified BRCA1 exon 11 variants (c.1019T>C, c.2363T>G, and c.3192T>C) on breast cancer susceptibility. Using computational predictions and experimental splicing assays, we evaluated their potential as pathogenic mutations. Our in silico analyses suggested that the c.2363T>G and c.3192T>C variants could impact both splicing and protein function, resulting in the V340A and V788G mutations, respectively. We further examined their splicing effects using minigene assays in MCF7 and SKBR3 breast cancer cell lines. Interestingly, we found that the c.2363T>G variant significantly altered splicing patterns in MCF7 cells but not in SKBR3 cells. This finding suggests a potential influence of cellular context on the variant's effects. While attempts to correlate in silico predictions with RNA binding factors were inconclusive, this observation underscores the complexity of splicing regulation. Splicing is governed by various factors, including cellular contexts and protein interactions, making it challenging to predict outcomes accurately. Further research is needed to fully understand the functional consequences of the c.2363T>G variant in breast cancer pathogenesis. Integrating computational predictions with experimental data will provide valuable insights into the role of alternative splicing regulation in different breast cancer types and stages.

GEN1 as a risk factor for human congenital anomalies of the kidney and urinary tract

BackgroundCongenital anomalies of the kidney and urinary tract (CAKUT) are prevalent birth defects. Although pathogenic CAKUT genes are known, they are insufficient to reveal the causes for all patients. Our previous studies indicated GEN1 as a pathogenic gene of CAKUT in mice, and this study further investigated the correlation between GEN1 and human CAKUT.MethodsIn this study, DNA from 910 individuals with CAKUT was collected; 26 GEN1 rare variants were identified, and two GEN1 (missense) variants in a non-CAKUT group were found. Mainly due to the stability results of the predicted mutant on the website, in vitro, 10 variants (eight CAKUT, two non-CAKUT) were selected to verify mutant protein stability. In addition, mainly based on the division of the mutation site located in the functional region of the GEN1 protein, 8 variants (six CAKUT, two non-CAKUT) were selected to verify enzymatic hydrolysis, and the splice variant GEN1 (c.1071 + 3(IVS10) A > G) was selected to verify shear ability. Based on the results of in vitro experiments and higher frequency, three sites with the most significant functional change were selected to build mouse models.ResultsProtein stability changed in six variants in the CAKUT group. Based on electrophoretic mobility shift assay of eight variants (six CAKUT, two non-CAKUT), the enzymatic hydrolysis and DNA-binding abilities of mutant proteins were impaired in the CAKUT group. The most serious functional damage was observed in the Gen1 variant that produced a truncated protein. A mini-gene splicing assay showed that the variant GEN1 (c.1071 + 3(IVS10) A > G) in the CAKUT group significantly affected splicing function. An abnormal exon10 was detected in the mini-gene splicing assay. Point-mutant mouse strains were constructed (Gen1: c.1068 + 3 A > G, p.R400X, and p.T105R) based on the variant frequency in the CAKUT group and functional impairment in vitro study and CAKUT phenotypes were replicated in each.ConclusionOverall, our findings indicated GEN1 as a risk factor for human CAKUT.

Identification of a novel intronic variant of ATP6V0A2 in a Han-Chinese family with cutis laxa.

Cutis laxa is a connective tissue disease caused by abnormal synthesis or secretion of skin elastic fibers, leading to skin flabby and saggy in various body parts. It can be divided into congenital cutis laxa and acquired cutis laxa, and inherited cutis laxa syndromes is more common in clinic. In this study, we reported a case of a Han-Chinese male newborn with ATP6V0A2 gene variant leading to cutis laxa. The proband was identified by whole-exome sequencing to determine the novel variant, and their parents were verified by Sanger sequencing. Bioinformatics analysis and minigene assay were used to verify the effect of this variant on splicing function. The main manifestations of the proband are skin laxity, abnormal facial features, and enlargement of the anterior fontanelle. Whole-exome sequencing showed that the newborn carried a non-canonical splicing-site variant c.117 + 5G > T, p. (?) in ATP6V0A2 gene. Sanger sequencing showed that both parents of the proband carried the heterozygous variant. The results of bioinformatics analysis and minigene assay displayed that the variant site affected the splicing function of pre-mRNA of the ATP6V0A2 gene. In this study, it was identified that ATP6V0A2 gene c. 117 + 5G > T may be the cause of the disease. The non-canonical splicing variants of ATP6V0A2 gene were rarely reported in the past, and this variant expanded the variants spectrum of the gene. The functional study of minigene assay plays a certain role in improving the level of evidence for the pathogenicity of splicing variants, which lays a foundation for prenatal counseling and follow-up gene therapy.

MALAT1: A Long Non-Coding RNA with Multiple Functions and Its Role in Processes Associated with Fat Deposition.

Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) belongs to the lncRNA molecules, which are involved in transcriptional and epigenetic regulation and the control of gene expression, including the mechanism of chromatin remodeling. MALAT1 was first discovered during carcinogenesis in lung adenocarcinoma, hence its name. In humans, 66 of its isoforms have been identified, and in pigs, only 2 are predicted, for which information is available in Ensembl databases (Ensembl Release 111). MALAT1 is expressed in numerous tissues, including adipose, adrenal gland, heart, kidney, liver, ovary, pancreas, sigmoid colon, small intestine, spleen, and testis. MALAT1, as an lncRNA, shows a wide range of functions. It is involved in the regulation of the cell cycle, where it has pro-proliferative effects and high cellular levels during the G1/S and mitotic (M) phases. Moreover, it is involved in invasion, metastasis, and angiogenesis, and it has a crucial function in alternative splicing during carcinogenesis. In addition, MALAT1 plays a significant role in the processes of fat deposition and adipogenesis. The human adipose tissue stem cells, during differentiation into adipocytes, secrete MALAT1 as one the most abundant lncRNAs in the exosomes. MALAT1 expression in fat tissue is positively correlated with adipogenic FABP4 and LPL. This lncRNA is involved in the regulation of PPARγ at the transcription stage, fatty acid metabolism, and insulin signaling. The wide range of MALAT1 functions makes it an interesting target in studies searching for drugs to prevent obesity development in humans. In turn, in farm animals, it can be a source of selection markers to control the fat tissue content.

Phosphorylation of the compartmentalized PKA substrate TAF15 regulates RNA-protein interactions.

Spatiotemporal-controlled second messengers alter molecular interactions of central signaling nodes for ensuring physiological signal transmission. One prototypical second messenger molecule which modulates kinase signal transmission is the cyclic-adenosine monophosphate (cAMP). The main proteinogenic cellular effectors of cAMP are compartmentalized protein kinase A (PKA) complexes. Their cell-type specific compositions precisely coordinate substrate phosphorylation and proper signal propagation which is indispensable for numerous cell-type specific functions. Here we present evidence that TAF15, which is implicated in the etiology of amyotrophic lateral sclerosis, represents a novel nuclear PKA substrate. In cross-linking and immunoprecipitation experiments (iCLIP) we showed that TAF15 phosphorylation alters the binding to target transcripts related to mRNA maturation, splicing and protein-binding related functions. TAF15 appears to be one of multiple PKA substrates that undergo RNA-binding dynamics upon phosphorylation. We observed that the activation of the cAMP-PKA signalingaxis caused a change in the composition of a collection of RNA species that interact with TAF15. This observation appears to be a broader principle in the regulation of molecular interactions, as we identified a significant enrichment of RNA-binding proteins within endogenous PKA complexes. We assume that phosphorylation of RNA-binding domains adds another layer of regulation to binary protein-RNAs interactions with consequences to RNA features including binding specificities, localization, abundance and composition.

Exploring salt tolerance mechanisms using machine learning for transcriptomic insights: case study in Spartina alterniflora.

Salt stress poses a significant threat to global cereal crop production, emphasizing the need for a comprehensive understanding of salt tolerance mechanisms. Accurate functional annotations of differentially expressed genes are crucial for gaining insights into the salt tolerance mechanism. The challenge of predicting gene functions in under-studied species, especially when excluding infrequent GO terms, persists. Therefore, we proposed the use of NetGO 3.0, a machine learning-based annotation method that does not rely on homology information between species, to predict the functions of differentially expressed genes under salt stress. Spartina alterniflora, a halophyte with salt glands, exhibits remarkable salt tolerance, making it an excellent candidate for in-depth transcriptomic analysis. However, current research on the S. alterniflora transcriptome under salt stress is limited. In this study we used S. alterniflora as an example to investigate its transcriptional responses to various salt concentrations, with a focus on understanding its salt tolerance mechanisms. Transcriptomic analysis revealed substantial changes impacting key pathways, such as gene transcription, ion transport, and ROS metabolism. Notably, we identified a member of the SWEET gene family in S. alterniflora, SA_12G129900.m1, showing convergent selection with the rice ortholog SWEET15. Additionally, our genome-wide analyses explored alternative splicing responses to salt stress, providing insights into the parallel functions of alternative splicing and transcriptional regulation in enhancing salt tolerance in S. alterniflora. Surprisingly, there was minimal overlap between differentially expressed and differentially spliced genes following salt exposure. This innovative approach, combining transcriptomic analysis with machine learning-based annotation, avoids the reliance on homology information and facilitates the discovery of unknown gene functions, and is applicable across all sequenced species.

ARGLU1 enhances promoter-proximal pausing of RNA polymerase II and stimulates DNA damage repair.

Arginine and glutamate rich 1 (ARGLU1) is a poorly understood cellular protein with functions in RNA splicing and transcription. Computational prediction suggests that ARGLU1 contains intrinsically disordered regions and lacks any known structural or functional domains. We used adenovirus Early protein 1A (E1A) to probe for critical regulators of important cellular pathways and identified ARGLU1 as a significant player in transcription and the DNA damage response pathway. Transcriptional effects induced by ARGLU1 occur via enhancement of promoter-proximal RNA polymerase II pausing, likely by inhibiting the interaction between JMJD6 and BRD4. When overexpressed, ARGLU1 increases the growth rate of cancer cells, while its knockdown leads to growth arrest. Significantly, overexpression of ARGLU1 increased cancer cell resistance to genotoxic drugs and promoted DNA damage repair. These results identify new roles for ARGLU1 in cancer cell survival and the DNA damage repair pathway, with potential clinical implications for chemotherapy resistance.

Abstract 1693: Mitochondrial-associated DNAJA3 variant and its role in NASH-driven hepatocellular carcinoma

Abstract Nonalcoholic fatty liver disease (NAFLD) is increasingly recognized as a major risk factor for hepatocellular carcinoma (HCC). Despite the known progression from NAFLD to nonalcoholic steatohepatitis (NASH) and eventually HCC, the genetic underpinnings driving this transition remain underexplored. A comprehensive genomics assessment was conducted, encompassing a genome-wide study of variants linked to body fat distribution in 344,369 individuals. Correlating these, significant variants associated with NASH and HCC were discerned in a cohort of 1,009 participants from the NCI-UMD study. Candidate variants underwent further eQTL analysis. Our investigation unveiled the rs3747579-TT variant emerged as significantly associated with NASH-related HCC, serving as an eQTL for mitochondrial DNAJA3. HCC patients harboring this variant presented with diminished DNAJA3 expression, correlating with an unfavorable prognosis. Exploration into the 3D genome architecture has identified potential chromatin topology sites that may influence the regulation of DNAJA3 by rs3747579. Specifically, enhancer loops exhibit allele-specific interactions with rs3747579 in the regulation of DNAJA3. Notably, the rs3747579-CC variant was found to interact with RBFOX2. Beyond its known RNA splicing functions, RBFOX2 might have implications in the chromatin landscape, influencing chromatin dynamics and gene transcription. RBFOX2 knockdown experiments showcased significantly reduced DNAJA3 expression, and luciferase assays highlighted increased activity for the rs3747579-CC allele versus its TT counterpart. In summary, our findings spotlight the rs3747579-TT variant as a potential oncogenic factor in DNAJA3, elucidating its role in the genesis of NASH and progression to HCC. This understanding could pave the way for novel therapeutic interventions in HCC stemming from NAFLD/NASH. Citation Format: Yuto Shiode, Ching-Wen Chang, Xin Wei Wang. Mitochondrial-associated DNAJA3 variant and its role in NASH-driven hepatocellular carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1693.

Identification of splicing factors signature predicting prognosis risk and the mechanistic roles of novel oncogenes in HNSCC

Head and neck squamous cell carcinoma (HNSCC) is the most frequent subtype of head and neck cancer, generally with a poor prognosis and limited therapeutic options due to its highly heterogeneous malignancy. In this study, we screened functional splicing regulatory RNA binding proteins (RBPs) that were closely related with the prognosis of HNSCC patients and showed significant expression differences between HNSCC tumors and normal tissues. Based on this finding, we chose six candidate genes (HNRNPC, ZCRB1, RBM12B, SF3A2, SF3B3, and SRSF11) to generate a prognostic prediction model and validated the accuracy of the prognostic model for predicting patient survival outcomes. We found that the risk score predicted by our model can serve as an independent prognostic predictor. Notably, HNSCC tumors showing higher expression of SF3B3, HNRNPC, or ZCRB1 possessed higher risk scores in the discovered prediction model. The investigation of the underlying mechanism validated that knockdown of SF3B3, HNRNPC, and ZCRB1 separately induced a substantial impairment of HNSCC cell survival. Conversely, overexpression of each of the three genes promoted tumor cellular proliferation. High throughput RNA sequencing analysis revealed that changes in the expression of SF3B3 and HNRNPC remarkably affected alternative splicing of genes related to cell cycle regulation, whereas the depletion of ZCRB1 contributed to aberrant splicing events involving in DNA damage response. In addition, the prognostic prediction model's risk score was demonstrated to be related with the immune infiltration score. Particularly, SF3B3 has a negative correlation with CD8A expression. Therefore, our findings provide promising prognosis predictors and potential therapeutic targets for better treatment efficacy of HNSCC.

Characterization of the chloroplast genome of Gleditsia species and comparative analysis

The genus Gleditsia has significant medicinal and economic value, but information about the chloroplast genomic characteristics of Gleditsia species has been limited. Using the Illumina sequencing, we assembled and annotated the whole chloroplast genomes of seven Gleditsia species (Gleditsia sinensis, Gleditsia japonica var. delavayi (G. delavayi), G. fera, G. japonica, G. microphylla, Fructus Gleditsiae Abnormalis (Zhū Yá Zào), G. microphylla mutant). The assembled genomes revealed that Gleditsia species have a typical circular tetrad structure, with genome sizes ranging from 162,746 to 170,907 bp. Comparative genomic analysis showed that most (65.8–75.8%) of the abundant simple sequence repeats in Gleditsia and Gymnocladus species were located in the large single copy region. The Gleditsia chloroplast genome prefer T/A-ending codons and avoid C/G-ending codons, positive selection was acting on the rpoA, rpl20, atpB, ndhA and ycf4 genes, most of the chloroplast genes of Gleditsia species underwent purifying selection. Expansion and contraction of the inverted repeat (IR)/single copy (SC) region showed similar patterns within the Gleditsia genus. Polymorphism analysis revealed that coding regions were more conserved than non-coding regions, and the IR region was more conserved than the SC region. Mutational hotspots were mostly found in intergenic regions such as “rps16-trnQ”, “trnT-trnL”, “ndhG-ndhI”, and "rpl32-trnL” in Gleditsia. Phylogenetic analysis showed that G. fera is most closely related to G. sinensis,G. japonica and G. delavayi are relatively closely related. Zhū Yá Zào can be considered a bud mutation of the G. sinensis. The albino phenotype of G. microphylla mutant is not caused by variations in the chloroplast genome, and that the occurrence of the albino phenotype may be due to mutations in chloroplast-related genes involved in splicing or localization functions. This study will help us enhance our exploration of the genetic evolution and geographical origins of the Gleditsia genus.

Temperature-sensitive splicing defects in Arabidopsis mitochondria caused by mutations in the ROOT PRIMORDIUM DEFECTIVE 1 gene.

Group II introns in plant organelles have lost splicing autonomy and require the assistance of nuclear-encoded trans-factors whose roles remain to be elucidated. These factors can be mono- or poly-specific with respect to the number of introns whose splicing they facilitate. Poly-acting splicing factors are often essential and their genetic identification may benefit from the use of conditional mutations. Temperature-sensitive (TS) mutations in the ROOT PRIMORDIUM DEFECTIVE 1 (RPD1) gene were initially selected for their inhibitory effect on root formation in Arabidopsis. Further analysis revealed that RPD1 encodes a mitochondria-targeted RNA-binding protein family member, suggesting a role in mitochondrial gene expression and making its role in root formation enigmatic. We analysed the function of RPD1 and found that it is required for the removal of 9 mitochondrial group II introns and that the identified TS mutations affect the splicing function of RPD1. These results support that the inhibition of adventitious root formation at non-permissive temperature results from a reduction in RPD1 activity and thus mitochondrial activity. We further show that RPD1 physically associates in vivo with the introns whose splicing it facilitates. Preliminary mapping indicates that RPD1 may not bind to the same regions within all of its intron targets, suggesting potential variability in its influence on splicing activation.