Splicing-accessible coding 3'UTRs control protein stability and interaction networks

Alternative Splicing of BnABF4L Mediates Response to Abiotic Stresses in Rapeseed (Brassica napus L.)

ABRE BINDING FACTOR 4 (ABF4) is a pivotal regulatory gene in the abscisic acid (ABA) signaling pathway, and changes in its expression levels can modulate the plant's stress resistance. To further explore the specific regulatory mechanisms of alternative splicing (AS) in the ABA signaling pathway and to identify new breakthroughs for breeding high stress-resistant varieties of Brassica napus, we identified 17 homologous genes of ABF4 in the genome. Utilizing bioinformatics techniques, we analyzed their motifs, conserved domains, and cis-acting elements of their promoters. Through transcriptome data from the stress-tolerant dwarf strain ndf2 and its parental line 3529, we uncovered a significantly differentially expressed ABF4 gene, which we named BnABF4L. Subsequently, we analyzed the AS events of BnABF4L under normal growth conditions and different abiotic stresses, as well as the impact of different transcript variants' 5' untranslated region (5'UTR) on gene translation. BnABF4L undergoes alternative 3' splice site (A3SS) selection to produce three transcripts (V1-V3) with divergent 5'UTRs. While V1 translation is suppressed by upstream ORFs (uORFs), V2/V3 exhibit enhanced translational efficiency. Under stress, ndf2 shifts splicing toward V3, circumventing uORF-mediated repression to upregulate stress-adapted isoforms. We validated the inhibitory effect of upstream open reading frames (uORFs) on protein-coding open reading frame (pORFs) and, based on the collective experimental results, proposed the flexible regulatory mechanism of AS events of BnABF4L in response to stress. Our findings provide new insights for future studies on stress resistance in rapeseed as well as for research on the regulation of alternative splicing mechanisms in the ABA signaling pathway.

Immediate early splicing controls translation in activated T-cells and is mediated by hnRNPC2 phosphorylation

The fast and transient induction of immediate early genes orchestrates the cellular response to various stimuli. These stimuli trigger phosphorylation cascades that promote immediate early gene transcription independent of de novo protein synthesis. Here we show that the same phosphorylation cascades also target the splicing machinery, inducing an analogous splicing switch that we call immediate early splicing (IES). We characterize hnRNPC2-controlled IES, which depends on the MEK-ERK pathway and the T cell-specific kinase PKCθ. This splicing switch mainly targets components of the translation machinery, such as mRNAs encoding ribosomal proteins and eIF5A. Inducing the eIF5A IES protein variant is by itself sufficient to reduce global translation, and consistently, we observe reduced de novo protein synthesis early after T cell activation. We suggest that immediate early splicing and the ensuing transient decrease in translation efficiency help to coordinate the extensive changes in gene expression during T cell activation. Together, these findings set a paradigm for fast and transient alternative splicing in the immediate cellular response to activation, and provide evidence for its functional relevance during T-cell stimulation.

Clusters of mammalian conserved RNA structures in UTRs associate with RBP binding sites

RNA secondary structures play essential roles in the formation of the tertiary structure and function of a transcript. Recent genome-wide studies highlight significant potential for RNA structures in the mammalian genome. However, a major challenge is assigning functional roles to these structured RNAs. In this study, we conduct a guilt-by-association analysis of clusters of computationally predicted conserved RNA structure (CRSs) in human untranslated regions (UTRs) to associate them with gene functions. We filtered a broad pool of ∼500 000 human CRSs for UTR overlap, resulting in 4734 and 24 754 CRSs from the 5' and 3' UTR of protein-coding genes, respectively. We separately clustered these CRSs for both sets using RNAscClust, obtaining 793 and 2403 clusters, each containing an average of five CRSs per cluster. We identified overrepresented binding sites for 60 and 43 RNA-binding proteins co-localizing with the clustered CRSs. Furthermore, 104 and 441 clusters from the 5' and 3' UTRs, respectively, showed enrichment for various Gene Ontologies, including biological processes such as 'signal transduction', 'nervous system development', molecular functions like 'transferase activity' and the cellular components such as 'synapse' among others. Our study shows that significant functional insights can be gained by clustering RNA structures based on their structural characteristics.

Identification of Noncoding Functional Regulatory Variants of STIM1 Gene in Idiopathic Pulmonary Arterial Hypertension

BACKGROUND

STIM1 (stromal interaction molecule 1) regulates store-operated calcium entry and is involved in pulmonary artery vasoconstriction and pulmonary artery smooth muscle cell proliferation, leading to pulmonary arterial hypertension (PAH).

METHODS

Bioinformatics analysis and a 2-stage matched case-control study were conducted to screen for noncoding variants that may potentially affect STIM1 transcriptional regulation in 242 patients with idiopathic PAH and 414 healthy controls. Luciferase reporter assay, real-time quantitative polymerase chain reaction, western blot, 5-ethynyl-2'-deoxyuridine (EdU) assay, and intracellular Ca2+ measurement were performed to study the mechanistic roles of those STIM1 noncoding variants in PAH.

RESULTS

Five noncoding variants (rs3794050, rs7934581, rs3750996, rs1561876, and rs3750994) were identified and genotyped using Sanger sequencing. Rs3794050, rs7934581, and rs1561876 were associated with idiopathic PAH (recessive model, all P<0.05). Bioinformatics analysis showed that these 3 noncoding variants possibly affect the enhancer function of STIM1 or the microRNA (miRNA) binding to STIM1. Functional validation performed in HEK293 and pulmonary artery smooth muscle cells demonstrated that the noncoding variant rs1561876-G (STIM1 mutant) had significantly stronger transcriptional activity than the wild-type counterpart, rs1561876-A, by affecting the transcriptional regulatory function of both hsa-miRNA-3140-5p and hsa-miRNA-4766-5p. rs1561876-G enhanced intracellular Ca2+ signaling in human pulmonary artery smooth muscle cells secondary to calcium-sensing receptor activation and promoted proliferation of pulmonary artery smooth muscle cells under both normoxia and hypoxia conditions, suggesting a possible contribution to PAH development.

CONCLUSIONS

The potential clinical implications of the 3 noncoding variants of STIM1, rs3794050, rs7934581, and rs1561876, are 2-fold, as they may help predict the risk and prognosis of idiopathic PAH and guide investigations on novel therapeutic pathway(s).

Thermoregulated transcriptomics: the molecular basis and biological significance of temperature-dependent alternative splicing

Temperature-dependent alternative splicing (AS) is a crucial mechanism for organisms to adapt to varying environmental temperatures. In mammals, even slight fluctuations in body temperature are sufficient to drive significant AS changes in a concerted manner. This dynamic regulation allows organisms to finely tune gene expression and protein isoform diversity in response to temperature cues, ensuring proper cellular function and physiological adaptation. Understanding the molecular mechanisms underlying temperature-dependent AS thus provides valuable insights into the intricate interplay between environmental stimuli and gene expression regulation. In this review, we provide an overview of recent advances in understanding temperature-regulated AS across various biological processes and systems. We will discuss the machinery sensing and translating temperature cues into changed AS patterns, the adaptation of the splicing regulatory machinery to extreme temperatures, the role of temperature-dependent AS in shaping the transcriptome, functional implications and the development of potential therapeutics targeting temperature-sensitive AS pathways.

RNA structure profiling at single-cell resolution reveals new determinants of cell identity

RNA structure is critical for multiple steps in gene regulation. However, how the structures of transcripts differ both within and between individual cells is unknown. Here we develop a SHAPE-inspired method called single-cell structure probing of RNA transcripts that enables simultaneous determination of transcript secondary structure and abundance at single-cell resolution. We apply single-cell structure probing of RNA transcripts to human embryonic stem cells and differentiating neurons. Remarkably, RNA structure is more homogeneous in human embryonic stem cells compared with neurons, with the greatest homogeneity found in coding regions. More extensive heterogeneity is found within 3' untranslated regions and is determined by specific RNA-binding proteins. Overall RNA structure profiles better discriminate cell type identity and differentiation stage than gene expression profiles alone. We further discover a cell-type variable region of 18S ribosomal RNA that is associated with cell cycle and translation control. Our method opens the door to the systematic characterization of RNA structure-function relationships at single-cell resolution.

A biofertilizing fungal endophyte of cranberry plants suppresses the plant pathogen Diaporthe

Fungi colonizing plants are gaining attention because of their ability to promote plant growth and suppress pathogens. While most studies focus on endosymbionts from grasses and legumes, the large and diverse group of ericaceous plants has been much neglected. We recently described one of the very few fungal endophytes promoting the growth of the Ericaceae Vaccinium macrocarpon (American cranberry), notably the Codinaeella isolate EC4. Here, we show that EC4 also suppresses fungal pathogens, which makes it a promising endophyte for sustainable cranberry cultivation. By dual-culture assays on agar plates, we tested the potential growth suppression (or biocontrol) of EC4 on other microbes, notably 12 pathogenic fungi and one oomycete reported to infect not only cranberry but also blueberry, strawberry, tomato plants, rose bushes and olive trees. Under greenhouse conditions, EC4 protects cranberry plantlets infected with one of the most notorious cranberry-plant pathogens, Diaporthe vaccinii, known to cause upright dieback and berry rot. The nuclear genome sequence of EC4 revealed a large arsenal of genes potentially involved in biocontrol. About ∼60 distinct clusters of genes are homologs of secondary metabolite gene clusters, some of which were shown in other fungi to synthesize nonribosomal peptides and polyketides, but in most cases, the exact compounds these clusters may produce are unknown. The EC4 genome also encodes numerous homologs of hydrolytic enzymes known to degrade fungal cell walls. About half of the nearly 250 distinct glucanases and chitinases are likely involved in biocontrol because they are predicted to be secreted outside the cell. Transcriptome analysis shows that the expression of about a quarter of the predicted secondary-metabolite gene clusters and glucan and chitin-degrading genes of EC4 is stimulated when it is co-cultured with D. vaccinii. Some of the differentially expressed EC4 genes are alternatively spliced exclusively in the presence of the pathogen, altering the proteins' domain content and subcellular localization signal, thus adding a second level of proteome adaptation in response to habitat competition. To our knowledge, this is the first report of Diaporthe-induced alternative splicing of biocontrol genes.

Srsf1 and Elavl1 act antagonistically on neuronal fate choice in the developing neocortex by controlling TrkC receptor isoform expression

The seat of higher-order cognitive abilities in mammals, the neocortex, is a complex structure, organized in several layers. The different subtypes of principal neurons are distributed in precise ratios and at specific positions in these layers and are generated by the same neural progenitor cells (NPCs), steered by a spatially and temporally specified combination of molecular cues that are incompletely understood. Recently, we discovered that an alternatively spliced isoform of the TrkC receptor lacking the kinase domain, TrkC-T1, is a determinant of the corticofugal projection neuron (CFuPN) fate. Here, we show that the finely tuned balance between TrkC-T1 and the better known, kinase domain-containing isoform, TrkC-TK+, is cell type-specific in the developing cortex and established through the antagonistic actions of two RNA-binding proteins, Srsf1 and Elavl1. Moreover, our data show that Srsf1 promotes the CFuPN fate and Elavl1 promotes the callosal projection neuron (CPN) fate in vivo via regulating the distinct ratios of TrkC-T1 to TrkC-TK+. Taken together, we connect spatio-temporal expression of Srsf1 and Elavl1 in the developing neocortex with the regulation of TrkC alternative splicing and transcript stability and neuronal fate choice, thus adding to the mechanistic and functional understanding of alternative splicing in vivo.

A KLHL40 3' UTR splice-altering variant causes milder NEM8, an under-appreciated disease mechanism

Nemaline myopathy 8 (NEM8) is typically a severe autosomal recessive disorder associated with variants in the kelch-like family member 40 gene (KLHL40). Common features include fetal akinesia, fractures, contractures, dysphagia, respiratory failure and neonatal death. Here, we describe a 26-year-old man with relatively mild NEM8. He presented with hypotonia and bilateral femur fractures at birth, later developing bilateral Achilles' contractures, scoliosis, and elbow and knee contractures. He had walking difficulties throughout childhood and became wheelchair bound from age 13 after prolonged immobilization. Muscle magnetic resonance imaging at age 13 indicated prominent fat replacement in his pelvic girdle, posterior compartments of thighs and vastus intermedius. Muscle biopsy revealed nemaline bodies and intranuclear rods. RNA sequencing and western blotting of patient skeletal muscle indicated significant reduction in KLHL40 mRNA and protein, respectively. Using gene panel screening, exome sequencing and RNA sequencing, we identified compound heterozygous variants in KLHL40; a truncating 10.9 kb deletion in trans with a likely pathogenic variant (c.*152G > T) in the 3' untranslated region (UTR). Computational tools SpliceAI and Introme predicted the c.*152G > T variant created a cryptic donor splice site. RNA-seq and in vitro analyses indicated that the c.*152G > T variant induces multiple de novo splicing events that likely provoke nonsense mediated decay of KLHL40 mRNA explaining the loss of mRNA expression and protein abundance in the patient. Analysis of 3' UTR variants in ClinVar suggests variants that introduce aberrant 3' UTR splicing may be underrecognized in Mendelian disease. We encourage consideration of this mechanism during variant curation.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

SRC homology 3 domains: multifaceted binding modules

The assembly of complexes following the detection of extracellular signals is often controlled by signaling proteins comprising multiple peptide binding modules. The SRC homology (SH)3 family represents the archetypical modular protein interaction module, with ~300 annotated SH3 domains in humans that regulate an impressive array of signaling processes. We review recent findings regarding the allosteric contributions of SH3 domains host protein context, their phosphoregulation, and their roles in phase separation that challenge the simple model in which SH3s are considered to be portable domains binding to specific proline-rich peptide motifs.

YBX1 mediates alternative splicing and maternal mRNA decay during pre-implantation development

Background

In mammals, maternal gene products decay and zygotic genome activation (ZGA) during maternal to zygotic transition (MZT) is critical for the early embryogenesis. Y-box binding protein YBX1 plays vital roles in RNA stabilization and transcriptional regulation, but its roles remain to be elucidated during pre-implantation development.

Methods

In the present study, we re-analyzed transcriptional level of YBX1 in mice, human, bovine, and goat embryos using public RNA-seq datasets. We further performed siRNA microinjection to knock down the expression of YBX1, and RNA sequencing of the 8-cell stage embryos in the control and YBX1 knockdown group. To reveal the regulation mechanisms of YBX1, we conducted differentially expression analysis, alternative splicing (AS) analysis, enrichment analysis, and 5-EU staining using DESeq2, rMATs, clusterProfiler, and immunofluorescence technique, respectively.

Results

The expression of YBX1 was increased during MZT in goat, bovine, human, and mice, but significantly decreased in YBX1 knockdown embryos compared with the controls, suggesting successfully knockdown of YBX1. The percentage of blastocyst was decreased, while embryos blocked at the 2- and 4-cell stage were increased in YBX1 knockdown embryos compared to the controls. Using RNA-seq, we identified 1623 up-regulated and 3531 down-regulated genes in the 8-cell stage YBX1 knockdown embryos. Of note, the down-regulated genes were enriched in regulation of RNA/mRNA stability and spliceosome, suggesting that YBX1 might medicate RNA stability and AS. To this end, we identified 3284 differential AS events and 1322 differentially expressed maternal mRNAs at the 8-cell stage YBX1 knockdown embryos. Meanwhile, the splicing factors and mRNA decay-related genes showed aberrant expression, and the transcriptional activity during ZGA in goat and mice was compromised when YBX1 was knocked down.

Conclusion

YBX1 serves an important role in maternal mRNA decay, alternative splicing, and the transcriptional activity required for early embryogenesis, which will broaden the current understanding of YBX1 functions during the stochastic reprogramming events.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00743-4.

Application of deep learning algorithm on whole genome sequencing data uncovers structural variants associated with multiple mental disorders in African American patients

Mental disorders present a global health concern, while the diagnosis of mental disorders can be challenging. The diagnosis is even harder for patients who have more than one type of mental disorder, especially for young toddlers who are not able to complete questionnaires or standardized rating scales for diagnosis. In the past decade, multiple genomic association signals have been reported for mental disorders, some of which present attractive drug targets. Concurrently, machine learning algorithms, especially deep learning algorithms, have been successful in the diagnosis and/or labeling of complex diseases, such as attention deficit hyperactivity disorder (ADHD) or cancer. In this study, we focused on eight common mental disorders, including ADHD, depression, anxiety, autism, intellectual disabilities, speech/language disorder, delays in developments, and oppositional defiant disorder in the ethnic minority of African Americans. Blood-derived whole genome sequencing data from 4179 individuals were generated, including 1384 patients with the diagnosis of at least one mental disorder. The burden of genomic variants in coding/non-coding regions was applied as feature vectors in the deep learning algorithm. Our model showed ~65% accuracy in differentiating patients from controls. Ability to label patients with multiple disorders was similarly successful, with a hamming loss score less than 0.3, while exact diagnostic matches are around 10%. Genes in genomic regions with the highest weights showed enrichment of biological pathways involved in immune responses, antigen/nucleic acid binding, chemokine signaling pathway, and G-protein receptor activities. A noticeable fact is that variants in non-coding regions (e.g., ncRNA, intronic, and intergenic) performed equally well as variants in coding regions; however, unlike coding region variants, variants in non-coding regions do not express genomic hotspots whereas they carry much more narrow standard deviations, indicating they probably serve as alternative markers.

Alternative splicing and liver disease

Alternative splicing produces complex and dynamic changes in the protein isoforms that are necessary for the proper biological functioning of the metabolic pathways involved in liver development and hepatocyte homeostasis. Changes in the physiological state of alternatively spliced forms are increasingly linked to liver pathologies. This may occur when the expression or function of the set of proteins controlling the alternative splicing processes are altered by external effectors such as oxidative stress and other environmental variations. Studies addressing these modifications reveal a complex interplay between the expression levels of different proteins that regulate the alternative splicing process as well as the changes in alternative splicing. This interplay results in a cascade of different protein isoforms that correlate with the progression of non-alcoholic fatty liver disease, hepatocellular carcinoma, and alcoholic liver disease. However, research on the detailed molecular mechanism underlying the production of these isoforms is needed. It is imperative to identify the physiological processes affected by the differentially spliced isoforms and confirm their role on the onset and maintenance of the pathology. This is required to design potential therapeutic approaches targeting the key splicing changes to revert the pathological condition as well as identify prognostic markers. In this review, we describe the complexity of the splicing process through an example to encourage researchers to go down this path. Subsequently, rather than a catalog of splicing events we have hand-picked and discuss a few selected studies of specific liver pathologies and suggested ways to focus research on these areas.

Nanopore sequencing reveals endogenous NMD-targeted isoforms in human cells

BACKGROUND

Nonsense-mediated mRNA decay (NMD) is a eukaryotic, translation-dependent degradation pathway that targets mRNAs with premature termination codons and also regulates the expression of some mRNAs that encode full-length proteins. Although many genes express NMD-sensitive transcripts, identifying them based on short-read sequencing data remains a challenge.

RESULTS

To identify and analyze endogenous targets of NMD, we apply cDNA Nanopore sequencing and short-read sequencing to human cells with varying expression levels of NMD factors. Our approach detects full-length NMD substrates that are highly unstable and increase in levels or even only appear when NMD is inhibited. Among the many new NMD-targeted isoforms that our analysis identifies, most derive from alternative exon usage. The isoform-aware analysis reveals many genes with significant changes in splicing but no significant changes in overall expression levels upon NMD knockdown. NMD-sensitive mRNAs have more exons in the 3΄UTR and, for those mRNAs with a termination codon in the last exon, the length of the 3΄UTR per se does not correlate with NMD sensitivity. Analysis of splicing signals reveals isoforms where NMD has been co-opted in the regulation of gene expression, though the main function of NMD seems to be ridding the transcriptome of isoforms resulting from spurious splicing events.

CONCLUSIONS

Long-read sequencing enables the identification of many novel NMD-sensitive mRNAs and reveals both known and unexpected features concerning their biogenesis and their biological role. Our data provide a highly valuable resource of human NMD transcript targets for future genomic and transcriptomic applications.

A Novel, Universally Active C-terminal Protein Degradation Signal Generated by Alternative Splicing

Proteome integrity is crucial for cellular homeostasis and adaptation to stress conditions such as hypoxia. One mechanism for rapid adaptation of the proteome in response to changing environmental signals is alternative splicing. In addition to generating different protein isoforms, alternative splicing is also capable of controlling total protein levels by the regulated synthesis of non-productive mRNA isoforms. The hypoxia-induced isoform E of the tumor suppressor MAX is produced by retention and translation of the last intron. This leads to an alternative C-terminus that harbors a potent C-degron, the isoE degron. Strikingly, the isoE degron represents a universal protein degradation signal that is not only functional in mammalian cells, but also in yeast and even in bacteria. Essential for efficient protein decay is a conserved (F/W)xxW motif. Degradation of isoE tagged proteins is mediated by the proteasome in eukaryotes and Lon protease in bacteria. Thus, the isoE degron is a broadly applicable and highly efficient tool in protein analyses.