Panicle architecture is a critical determinant of rice yield and resilience, yet the genetic and environmental factors shaping this trait remain incompletely understood. Here, we applied an integrative genomic approach combining multi-locus association mapping, transcriptome analysis and population genomics to dissect the genetic basis of key panicle traits in rice. We identified robust genetic loci underlying the number of primary branches, panicle length and spikelets per panicle, with many showing sensitivity to temperature, underscoring the importance of gene-environment interactions for yield stability. Notably, we discovered that variation in alternative polyadenylation (APA) of specific transcripts is associated with panicle trait diversity at the population level, suggesting that regulatory mechanisms such as APA are significant contributors to phenotypic plasticity and adaptation. These findings deliver both novel candidate genes in panicle development and mechanistic insights to support the breeding of rice varieties with enhanced productivity and climate resilience.

RNA and proteins can have diverse isoforms due to post-transcriptional and post-translational modifications. A fundamental question is whether these isoforms are mostly beneficial or the result of noisy molecular processes. To assess the plausibility of these explanations, we developed mathematical models depicting different regulatory architectures and investigated isoform evolution under multiple population genetic regimes. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. We found that sub-optimal phenotypes are more likely to evolve when populations are small and/or when the number of cis-loci is large. We also discovered that opposing selection on cis- and trans-acting loci can constrain adaptation, leading to a non-monotonic relationship between effective population size and optimization. More generally, our models provide a quantitative framework for developing statistical tests to analyze empirical data; as a demonstration of this, we analyzed A-to-I RNA editing levels in coleoids and found these to be largely consistent with non-adaptive explanations.

Alternative polyadenylation (APA) plays a crucial role in cancer development and prognosis. However, the molecular characteristics of APA related to non-small cell lung cancer (NSCLC) susceptibility remain understudied, especially in East Asian populations. In this study, we constructed an atlas of APA-regulated 3' untranslated region (3'UTR) and profiled its genetic regulation in 747 lung tissue samples (including tumors and paired normal tissues) from 417 NSCLC Chinese patients. We verified a significant global shortening of 3'UTRs in tumor samples compared to normal samples and underscored the value of APA-regulation as a prognostic marker. The 3'UTR APA quantitative trait loci (3'aQTL) was identified by regressing the percentage of distal poly(A) site usage index (PDUI) value on genetic variants. We found that a significant proportion 3'aQTLs are independent of genetic regulation of expression and are specific in Chinese. We also conducted a 3'UTR APA transcriptome-wide association study (3'aTWAS) by integrating the APA regulation atlas with a genome-wide association study (GWAS) for NSCLC involving 7035 cases and 185,413 cancer-free controls. We identified NSCLC-associated genes, highlighting TUBB, TEAD3, and PPP1R10. Combining the consistent results from colocalization analysis, differential APA analysis, and survival analysis, we provide novel evidence for the role TUBB APA regulation in NSCLC and identified potential upstream regulators. Overall, our study profiled the APA regulation and highlighted the substantial role of APA in NSCLC carcinogenesis and prognosis in East Asian populations.

The endonucleolytic cleavage step of the eukaryotic mRNA 3'-end processing is considered imprecise, which leads to heterogeneity of cleavage site (CS) with hitherto unknown function. Contrary to popular belief, we show that this imprecision in the cleavage is tightly regulated, resulting in the CS heterogeneity (CSH) that controls gene expression in antioxidant response. CSH centres around a primary CS, followed by several subsidiary cleavages determined by CS's positions. Globally and using reporter antioxidant mRNA, we discovered an inverse relationship between the number of CS and the gene expression, with the primary CS exhibiting the highest cleavage efficiency. Strikingly, reducing CSH and increasing primary CS usage induces gene expression. Under oxidative stress (we employ three conditions that induce antioxidant response, tBHQ, H2O2, and NaAsO2) conditions, there is a decrease in the CSH and an increase in the primary CS usage to induce antioxidant gene expression. Key oxidative stress response genes (NQO1, HMOX1, PRDX1, and CAT) also show higher CSH compared to the non-stress response genes and that the number of CSs are reduced to impart cellular response to oxidative stresses. Concomitantly, ectopic expression of one of the key antioxidant response gene (NQO1) driven by the primary CS but not from other subsidiary CSs, or reduction in CSH imparts tolerance to cellular oxidative stresses (H2O2, and NaAsO2). Genome-wide CS analysis of stress response genes also shows a similar result. Compromised CSH or CSH-mediated gene control hampers cellular response to oxidative stress. We establish that oxidative stress induces affinity/strength of cleavage complex assembly, increasing the fidelity of cleavage at the primary CS, thereby reducing CSH inducing antioxidant response. Together, our study reports a novel cleavage imprecision- or CSH-mediated anti-oxidant response mechanism that is distinct and operates downstream but in concert with the transcriptional pathway of oxidative stress induction.

Most eukaryotic genes express more than one mature mRNA, defined as transcript variants. This complex phenomenon arises from various mechanisms, such as using alternative transcription start sites and alternative post-transcriptional processing events. The resulting transcript variants can lead to synthesizing proteins that possess distinct functional domains or may even generate noncoding RNAs, each with unique roles in cellular processes. The generation of these transcript variants is not merely a random occurrence; it is cell-type specific and varies with developmental stages, aging processes, or pathogenesis of diseases. This highlights the biological significance of transcript variants in regulating gene expression and their potential impact on cellular functionality. Despite the biological importance, investigating transcript variants has been hampered by challenges associated with detecting their expression. This review article addresses the advancements in molecular techniques in detecting transcript variants. Traditional methods such as RT-PCR and RT-qPCR can easily detect known transcript variants using primers that target unique exons associated with the variants. Other techniques like RACE-PCR and hybridization-based methods, including Northern blotting, RNase protection assays, and microarrays, have also been utilized to detect transcript variants. Nevertheless, RNA sequencing (RNA-Seq) has emerged as a powerful technique for identifying transcript variants, especially those with previously unknown sequences. The effectiveness of RNA sequencing in transcript variant detection depends on the specific sequencing approach and the precision of data analysis. By understanding the strengths and weaknesses of each laboratory technique, researchers can develop more effective strategies for detecting mRNA transcript variants. This ability will be crucial for our comprehensive understanding of gene regulation and the implications of transcript diversity in various biological contexts.

Transcriptome-wide alternative polyadenylation (APA) is involved in both innate and adaptive immune responses of immune cells. Downregulation of the CPSF6 protein, one of the 3' end-processing factors, mediates APA in macrophages with responses to virus infection and plays an important role in its anti-virus effect. However, the signaling pathway and molecular mechanism underlying the downregulation of the CPSF6 protein remain elusive. Here, we found that MAVS triggers the nuclear import of the E3 ligase SYVN1 mediated by NUP153 in response to vesicular stomatitis virus infection. Then, SYVN1 catalyzes K48-linked polyubiquitination of CPSF6, resulting in degradation of CPSF6 via the proteasome and then transcriptome-wide APA and anti-virus effects. Our results identify an antiviral mechanism via APA regulation based on ubiquitination modification of the CPSF6 protein, which may serve as a target for developing immune interventions.

Conserved proteins are characterized by their functions remaining nearly constant throughout evolutionary history, both vertically through time and horizontally across species. In this review, we focus on a class of conserved proteins known as the Muscleblind-like (MBNL) family. As RNA-binding proteins, MBNL family members interact with pre-mRNAs through evolutionarily conserved tandem zinc finger domains and play critical roles in various RNA metabolic processes, including alternative splicing, mRNA stability, trafficking, regulation of subcellular localization, and alternative polyadenylation. Dysregulation of MBNL proteins can lead to severe consequences. Initially, research primarily associated MBNL proteins with myotonic dystrophy. However, recent studies have revealed their involvement in a broad spectrum of physiological and pathological processes, such as embryonic tissue differentiation and circulatory disorders. Furthermore, the emerging role of MBNL proteins in cancer sheds light on a novel aspect of these evolutionarily ancient proteins. This review provides a comprehensive overview of the MBNL family, emphasizing its structure, the mechanisms underlying its biological functions, and its roles in various diseases.Subject terms: Muscleblind-like-like protein, RNA-binding proteins, Alternative splicing, Tumor, Myotonic dystrophy.

Epigenes are defined as proteins that perform post-translational modification of histones or DNA, reading of post-translational modifications, form complexes with epigenetic factors or changing the general structure of chromatin. This specialized group of proteins is responsible for controlling the organization of genomic DNA in a cell-type specific fashion, controlling normal development in a spatial and temporal fashion. Moreover, mutations in epigenes have been implicated as causal in germline pediatric disorders and as driver mutations in cancer. Despite their importance to human disease, to date, there has not been a systematic analysis of the sources of functional diversity for epigenes at large. Epigenes' unique functions that require the assembly of pools within the nucleus suggest that their structure and amino acid composition would have been enriched for features that enable efficient assembly of chromatin and DNA for transcription, splicing, and post-translational modifications.

In this study, we assess the functional diversity stemming from gene structure, isoforms, protein domains, and multiprotein complex formation that drive the functions of established epigenes. We found that there are specific structural features that enable epigenes to perform their variable roles depending on the cellular and environmental context. First, epigenes are significantly larger and have more exons compared with non-epigenes which contributes to increased isoform diversity. Second epigenes participate in more multimeric complexes than non-epigenes. Thirdly, given their proposed importance in membraneless organelles, we show epigenes are enriched for substantially larger intrinsically disordered regions (IDRs). Additionally, we assessed the specificity of their expression profiles and showed epigenes are more ubiquitously expressed consistent with their enrichment in pediatric syndromes with intellectual disability, multiorgan dysfunction, and developmental delay. Finally, in the L1000 dataset, we identify drugs that can potentially be used to modulate expression of these genes.

Here we identify significant differences in isoform usage, disordered domain content, and variable binding partners between human epigenes and non-epigenes using various functional genomics datasets from Ensembl, ENCODE, GTEx, HPO, LINCS L1000, and BrainSpan. Our results contribute new knowledge to the growing field focused on developing targeted therapies for diseases caused by epigene mutations, such as chromatinopathies and cancers.

Neurodegenerative diseases (NDDs) are one of the prevailing conditions characterized by progressive neuronal loss. Polyadenylation (PA) and alternative polyadenylation (APA) are the two main post-transcriptional events that regulate neuronal gene expression and protein production. This systematic review analyzed the available literature on the role of PA and APA in NDDs, with an emphasis on their contributions to disease development. A comprehensive literature search was performed using the PubMed, Scopus, Cochrane, Google Scholar, Embase, Web of Science, and ProQuest databases. The search strategy was developed based on the framework introduced by Arksey and O'Malley and supplemented by the inclusion and exclusion criteria. The study selection was performed by two independent reviewers. Extraction and data organization were performed in accordance with the predefined variables. Subsequently, quantitative and qualitative analyses were performed. Forty-seven studies were included, related to a variety of NDDs, namely Alzheimer's disease, Parkinson's disease, Huntington's disease, and amyotrophic lateral sclerosis. Disease induction was performed using different models, including human tissues, animal models, and cultured cells. Most investigations were related to PA, although some were related to APA or both. Amyloid precursor protein (APP), Tau, SNCA, and STMN2 were the major genes identified; most of the altered PA patterns were related to mRNA stability and translation efficiency. This review particularly underscores the key roles of PA and APA in the pathogenesis of NDDs through their mechanisms that contribute to gene expression dysregulation, protein aggregation, and neuronal dysfunction. Insights into these mechanisms may lead to new therapeutic strategies focused on the modulation of PA and APA activities. Further research is required to investigate the translational potential of targeting these pathways for NDD treatment.

Gene expression or gene regulation studies often assume one gene expresses one mRNA. However, contrary to the conventional idea, a single gene in mammalian cells can express multiple transcript variants translated into several different proteins. The transcript variants are generated through transcription from alternative start sites and alternative post-transcriptional processing of the precursor mRNA (pre-mRNA). In addition, gene mutations and RNA editing further enhance the diversity of the transcript variants. The transcript variants can encode proteins with various domains, expanding the functional repertoire of a single gene. Some transcript variants may not encode proteins but function as non-coding RNAs and regulate gene expression. The expression level of the transcript variants may vary between cell types or within the same cells under different biological conditions. Transcript variants are characteristic of cell differentiation in a particular tissue, and the variants may play a key role in normal development and aging. Studies also reported that some transcript variants may have roles in disease pathogenesis. The biological significances urge studying the complexity of gene expression at the transcript level. This article updates the molecular basis of transcript variants in mammalian cells, including the formation mechanisms and potential roles in host biology. Gaining insight into the transcript variants will not only identify novel mechanisms of gene regulation but also unravel the role of the variants in health and disease.

Alternative polyadenylation (APA) is a major mechanism of post-transcriptional regulation that affects mRNA stability, localization and translation efficiency. Previous pan-cancer studies have revealed that APA is frequently disrupted in cancer and is associated with patient outcomes. Yet, little is known about cancer type-specific APA alterations. Here, we integrated RNA-sequencing data from a Korean cohort (GEO: GSE40419) and The Cancer Genome Atlas (TCGA) to comprehensively analyze APA alterations in lung adenocarcinomas (LUADs). Comparing expression levels of core genes involved in polyadenylation, we find that overall, the set of 28 of 31 genes are upregulated, with CSTF2 particularly upregulated. We observed broad and recurrent APA changes in LUAD growth-promoting genes. In addition, we find enrichment of APA events in genes associated with known LUAD pathways and an increased heterogeneity in polyadenylation (polyA) site usage of proliferation-associated genes. Upon further investigation, we report smoking-specific APA changes are also highly relevant to LUAD development. Overall, our in-depth analysis reveals APA as an important driver for the molecular and clinical features of lung adenocarcinoma.

RNA-binding proteins (RBPs) are essential modulators in the regulation of mRNA processing. The binding patterns, interactions, and functions of most RBPs are not well-characterized. Previous studies have shown that motif context is an important contributor to RBP binding specificity, but its precise role remains unclear. Despite recent computational advances to predict RBP binding, existing methods are challenging to interpret and largely lack a categorical focus on RBP motif contexts and RBP-RBP interactions. There remains a need for interpretable predictive models to disambiguate the contextual determinants of RBP binding specificity in vivo . Here, we present a novel and comprehensive pipeline to address these knowledge gaps. We devise a Natural Language Processing-based decomposition method to deconstruct sequences into entities consisting of a central target k -mer and its flanking regions, then use this representation to formulate the RBP binding prediction task as a weakly supervised Multiple Instance Learning problem. To interpret our predictions, we introduce a deterministic motif discovery algorithm designed to handle our data structure, recapitulating the established motifs of numerous RBPs as validation. Importantly, we characterize the binding motifs and binding contexts for 71 RBPs, with many of them being novel. Finally, through feature integration, transitive inference, and a new cross-prediction approach, we propose novel cooperative and competitive RBP-RBP interaction partners and hypothesize their potential regulatory functions. In summary, we present a complete computational strategy for investigating the contextual determinants of specific RBP binding, and we demonstrate the significance of our findings in delineating RBP binding patterns, interactions, and functions.

Non-syndromic orofacial clefts (NSOFCs) are complex congenital abnormalities involving both environmental and genetic factors involved in orofacial development. This study aimed to investigate the genetic association of specific genetic variants at different CYRIA gene loci with the development of NSOFCs in Kuwait.

Four genetic variants (rs7552, rs3758249, rs3821949, and rs3917201) at four selected gene loci (CYRIA, FOXE1, MSX1, and TGFB3) were genotyped in a total of 240 DNA samples (patients (n = 114) and random controls (n = 126)) employing TaqMan® allele discrimination assay. For each variant and its genotype, the frequencies were determined and tested for Hardy-Weinberg Equilibrium. Genotype frequencies was compared between patients and controls using Pearson's test. Logistic regression analyses were employed to test for the associations of the four selected variants with the occurrence of NSOFCSs.

Significant differences in the distribution of genotypes between cases and controls, rs7552, rs3821949, and rs3917201 were found to have a positive association with NSOFCs. After adjusting for gender, the GG genotype of the rs7552 variant, the AG genotype of the rs3821949 variant, and the CC genotype of the rs3917201 variant showed nearly a two-fold increased risk of NSOFC (p < 0.05).

This study reports significant findings on the contribution and modest effect of CYRIA rs7552, MSX1 rs3821949, and TGFB3 rs3917201 in the development of NSOFCs. Our findings provide further evidence on the molecular mechanism and the role of the selected genes in NSOFCs.

Although rare non-coding variants (RVs) play crucial roles in complex traits and diseases, understanding their mechanisms and identifying disease-associated RVs continue to be major challenges. Here we constructed a comprehensive atlas of alternative polyadenylation (APA) outliers (aOutliers), including 1334 3' UTR and 200 intronic aOutliers, from 15,201 samples across 49 human tissues. These aOutliers exhibit unique characteristics from transcription or splicing outliers, with a pronounced RV enrichment. Mechanistically, aOutlier-RVs alter poly(A) signals and splicing sites, and perturbation indeed triggers APA events. Furthermore, we developed a Bayesian-based APA RV prediction model, which successfully pinpointed a specific set of 1799 RVs impacting 278 genes with significantly large disease effect sizes. Notably, we observed a convergence effect between rare and common cancer variants, exemplified by regulation in the DDX18 gene. Together, this study introduced an APA-enhanced framework for genome annotation, underscoring APA's role in uncovering functional RVs linked to complex traits and diseases.

Plasmodium falciparum is an obligate human parasite of the phylum Apicomplexa and is the causative agent of the most lethal form of human malaria. Although N6-methyladenosine modification is thought to be one of the major post-transcriptional regulatory mechanisms for stage-specific gene expression in apicomplexan parasites, the precise base position of m6A in mRNAs or noncoding RNAs in these parasites remains unknown. Here, we report global nucleotide-resolution mapping of m6A residues in P. falciparum using DART-seq technology, which quantitatively displayed a stage-specific, dynamic distribution pattern with enrichment near mRNA 3' ends. In this process we identified 894, 788, and 1,762 m6A-modified genes in Ring, Trophozoite and Schizont stages respectively, with an average of 5-7 m6A sites per-transcript at the individual gene level. Notably, several genes involved in malaria pathophysiology, such as KAHRP, ETRAMPs, SERA and stress response genes, such as members of Heat Shock Protein (HSP) family are highly enriched in m6A and therefore could be regulated by this RNA modification. Since we observed preferential methylation at the 3' ends of P. falciparum transcripts and because malaria polyadenylation specificity factor PfCPSF30 harbors an m6A reader 'YTH' domain, we reasoned that m6A might play an important role in 3'-end processing of malaria mRNAs. To investigate this, we used two complementary high-throughput RNA 3'-end mapping approaches, which provided an initial framework to explore potential roles of m6A in the regulation of alternative polyadenylation (APA) during malaria development in human hosts.

Autophagy is an essential recycling and quality control pathway which preserves cellular and organismal homeostasis. As a catabolic process, autophagy degrades damaged and aged intracellular components in response to conditions of stress, including nutrient deprivation, oxidative and genotoxic stress. Autophagy is a highly adaptive and dynamic process which requires an intricately coordinated molecular control. Here we provide an overview of how autophagy is regulated post-transcriptionally, through RNA processing events, epitranscriptomic modifications and non-coding RNAs. We further discuss newly revealed RNA-binding properties of core autophagy machinery proteins and review recent indications of autophagy's ability to impact cellular RNA homeostasis. From a physiological perspective, we examine the biological implications of these emerging regulatory layers of autophagy, particularly in the context of nutrient deprivation and tumorigenesis.

Aberrant gene expression lies at the heart of many pathologies. This review will point out how 3' end processing, the final mRNA-maturation step in the transcription cycle, is surprisingly prone to regulated as well as stochastic variations with a wide range of consequences. Whereas smaller variations contribute to the plasticity of gene expression, larger alternations to 3' end processing and coupled transcription termination can lead to pathological consequences. These can be caused by the local mutation of one gene or affect larger numbers of genes systematically, if aspects of the mechanisms of 3' end processing and transcription termination are altered.

Single-cell analysis of human peripheral blood cells provides insights into innate and adaptive immune systems. However, robust protocols are essential to ensuring single-cell sequencing data quality and cell viability. Here, we present a protocol for acquiring high-quality single-cell multi-omics data from human peripheral blood mononuclear cells (PBMCs). We describe steps for collecting human blood followed by single-cell sequencing, whole-genome sequencing, and metabolome and proteome analysis of PBMCs using modified multi-omics sample processing.

The African sharptooth catfish (Clarias gariepinus) assumes significance in aquaculture, given its role as a farmed freshwater species with modified gill structures functioning as an air-breathing organ (ABO). To provide a scientific basis for further elucidating the air-breathing formation mechanism and deeply utilizing the genetic resources of Clarias gariepinus, we utilized the PacBio sequencing platform to acquire a comprehensive full-length transcriptome from five juvenile developmental stages and various adult tissues, including the ABO, gills, liver, skin, and muscle. We generated 25,766,688 high-quality reads, with an average length of 2,006 bp and an N50 of 2,241 bp. Following rigorous quality control, 34,890 (97.7 %) of the high-quality isoforms were mapped to the reference genome for gene and transcript annotation, yielding 387 novel isoforms and 14,614 new isoforms. Additionally, we identified 28,582 open reading frames, 48 SNPs, 5,464 variable splices, and 6,141 variable polyadenylation sites, along with 475 long non-coding RNAs. Many DEGs were involved with low oxygen GO terms and KEGG pathways, such as response to stimulus, biological regulation and catalytic activities. Furthermore, it was found that transcription factors such as zf-C2H2, Homeobox, bHLH, and MYB could underpin the African sharptooth catfish's developmental plasticity and its capacity to adapt its morphology and function to its environment. Through the comprehensive analysis of its genomic characteristics, it was found that the African sharptooth catfish has developed a series of unique respiratory adaptive mechanisms during the evolutionary process, These results not only advances the understanding of genetic adaptations to hypoxia in Clarias fish but also provides a valuable framework for future studies aimed at improving aquaculture practices,besides provide important references and inspirations for the evolution of aquatic organisms.

Temperature is a pivotal ecological factor in the regulation of insect survival and reproduction [...].

Aging is a natural physiological process involving biological and genetic pathways. Growing evidence suggests that alterations in the epigenome during aging result in transcriptional changes, which play a significant role in the onset of age-related diseases, including cancer, cardiovascular disease, diabetes, and neurodegenerative disorders. For this reason, the epigenetic alterations in aging and age-related diseases have been reviewed, and the major extrinsic factors influencing these epigenetic alterations have been identified. In addition, the role of the gut microbiome and its metabolites as epigenetic modifiers has been addressed.

Long-term exposure to extrinsic factors such as air pollution, diet, drug use, environmental chemicals, microbial infections, physical activity, radiation, and stress provoke epigenetic changes in the host through several endocrine and immune pathways, potentially accelerating the aging process. Diverse studies have reported that the gut microbiome plays a critical role in regulating brain cell functions through DNA methylation and histone modifications. The interaction between genes and the gut microbiome serves as a source of adaptive variation, contributing to phenotypic plasticity. However, the molecular mechanisms and signaling pathways driving this process are still not fully understood.

Extrinsic factors are potential inducers of epigenetic alterations, which may have important implications for longevity. The gut microbiome serves as an epigenetic effector influencing host gene expression through histone and DNA modifications, while bidirectional interactions with the host and the underexplored roles of microbial metabolites and non-bacterial microorganisms such as fungi and viruses highlight the need for further research.

RNA splicing and 3'-cleavage and polyadenylation (CPA) are essential processes for the maturation of RNA. There have been extensive independent studies of these regulated processing events, including alternative splicing (AS) and alternative polyadenylation (APA). However, growing evidence suggests potential crosstalk between splicing and 3'-end processing in regulating AS or APA. Here, we first provide a brief overview of the molecular machines involved in splicing and 3'-end processing events, and then review recent studies on the functions and mechanisms of the crosstalk between the two processes. On the one hand, 3'-end processing can affect splicing, as 3'-end processing factors and CPA-generated polyA tail promote the splicing of the last intron. Beyond that, 3'-end processing factors can also influence the splicing of internal and terminal exons. Those 3'-end processing factors can also interact with different RNA-binding proteins (RBPs) to exert their effects on AS. The length of 3' untranslated region (3' UTR) can affect the splicing of upstream exons. On the other hand, splicing and CPA may compete within introns in generating different products. Furthermore, splicing within the 3' UTR is a significant factor contributing to 3' UTR diversity. Splicing also influences 3'-end processing through the actions of certain splicing factors. Interestingly, some classical RBPs play dual roles in both splicing and 3'-end processing. Finally, we discuss how long-read sequencing technologies aid in understanding the coordination of AS-APA events and envision that these findings may potentially promote the development of new strategies for disease diagnosis and treatment.

A-to-I RNA editing is an RNA modification that alters the RNA sequence relative to the its genomic blueprint. It is catalyzed by double-stranded RNA-specific adenosine deaminase (ADAR) enzymes, and contributes to the complexity and diversification of the proteome. Advancement in the study of A-to-I RNA editing has been facilitated by computational approaches for accurate mapping and quantification of A-to-I RNA editing based on sequencing data. In this chapter we review some of the main computational approaches currently used, describe potential hurdles, challenges and pitfalls, and discuss possible ways to mitigate them.

Alternative polyadenylation (APA) is an important driver of transcriptome diversity that generates messenger RNA isoforms with distinct 3' ends. The rapid development of single-cell and spatial transcriptomic technologies opened up new opportunities for exploring APA data to discover hidden cell subpopulations invisible in conventional gene expression analysis. However, conventional gene-level analysis tools are not fully applicable to APA data, and commonly used unsupervised dimensionality reduction methods often disregard experimentally derived annotations such as cell type identities. Here, we proposed a supervised analytical framework termed spvAPA, specifically used for APA analysis from both single-cell and spatial transcriptomics data. First, an iterative imputation method based on weighted nearest neighbor was designed to recover missing APA signatures, by integrating both gene expression and APA modalities. Second, a supervised feature selection method based on sparse partial least squares discriminant analysis was devised to identify APA features distinguishing cell types or spatial morphologies. Additionally, spvAPA improves the visualization of high-dimensional data for discovering novel cell subtypes, which considers APA features and dual modalities of gene expression and APA. Evaluations across nine single-cell and spatial transcriptomics datasets demonstrate the effectiveness and applicability of spvAPA. spvAPA is available at https://github.com/BMILAB/spvAPA.

Long-read sequencing (LRS) technologies have the potential to revolutionize scientific discoveries in RNA biology through the comprehensive identification and quantification of full-length mRNA isoforms. Despite great promise, challenges remain in the widespread implementation of LRS technologies for RNA-based applications, including concerns about low coverage, high sequencing error, and robust computational pipelines. Although much focus has been placed on defining mRNA exon composition and structure with LRS data, less careful characterization has been done of the ability to assess the terminal ends of isoforms, specifically, transcription start and end sites. Such characterization is crucial for completely delineating full mRNA molecules and regulatory consequences. However, there are substantial inconsistencies in both start and end coordinates of LRS reads spanning a gene, such that LRS reads often fail to accurately recapitulate annotated or empirically derived terminal ends of mRNA molecules. Here, we describe the specific challenges of identifying and quantifying mRNA terminal ends with LRS technologies and how these issues influence biological interpretations of LRS data. We then review recent experimental and computational advances designed to alleviate these problems, with ideal use cases for each approach. Finally, we outline anticipated developments and necessary improvements for the characterization of terminal ends from LRS data.

The transcriptome serves as a bridge that links genomic variation to phenotypic diversity. A vast number of studies using next-generation RNA sequencing (RNA-seq) over the last 2 decades have emphasized the essential roles of the plant transcriptome in response to developmental and environmental conditions, providing numerous insights into the dynamic changes, evolutionary traces, and elaborate regulation of the plant transcriptome. With substantial improvement in accuracy and throughput, direct RNA sequencing (DRS) has emerged as a new and powerful sequencing platform for precise detection of native and full-length transcripts, overcoming many limitations such as read length and PCR bias that are inherent to short-read RNA-seq. Here, we review recent advances in dissecting the complexity and diversity of plant transcriptomes using DRS as the main technological approach, covering many aspects of RNA metabolism, including novel isoforms, poly(A) tails, and RNA modification, and we propose a comprehensive workflow for processing of plant DRS data. Many challenges to the application of DRS in plants, such as the need for machine learning tools tailored to plant transcriptomes, remain to be overcome, and together we outline future biological questions that can be addressed by DRS, such as allele-specific RNA modification. This technology provides convenient support on which the connection of distinct RNA features is tightly built, sustainably refining our understanding of the biological functions of the plant transcriptome.

Batten disease is a group of rare inherited neurodegenerative diseases. Juvenile CLN3 disease is the most prevalent type, and the most common pathogenic variant shared by most patients is the "1-kb" deletion which removes two internal coding exons (7 and 8) in CLN3. Previously, we identified two transcripts in patient fibroblasts homozygous for the 1-kb deletion: the 'major' and 'minor' transcripts. To understand the full variety of disease transcripts and their role in disease pathogenesis, it is necessary to first investigate CLN3 transcription in "healthy" samples without juvenile CLN3 disease.

We leveraged PacBio long-read RNA sequencing datasets from ENCODE to investigate the full range of CLN3 transcripts across various tissues and cell types in human control samples. Then we sought to validate their existence using data from different sources.

We found that a readthrough gene affects the quantification and annotation of CLN3. After taking this into account, we detected over 100 novel CLN3 transcripts, with no dominantly expressed CLN3 transcript. The most abundant transcript has median usage of 42.9%. Surprisingly, the known disease-associated 'major' transcripts are detected. Together, they have median usage of 1.5% across 22 samples. Furthermore, we identified 48 CLN3 ORFs, of which 26 are novel. The predominant ORF that encodes the canonical CLN3 protein isoform has median usage of 66.7%, meaning around one-third of CLN3 transcripts encode protein isoforms with different stretches of amino acids. The same ORFs could be found with alternative UTRs. Moreover, we were able to validate the translational potential of certain transcripts using public mass spectrometry data.

Overall, these findings provide valuable insights into the complexity of CLN3 transcription, highlighting the importance of studying both canonical and non-canonical CLN3 protein isoforms as well as the regulatory role of UTRs to fully comprehend the regulation and function(s) of CLN3. This knowledge is essential for investigating the impact of the 1-kb deletion and rare pathogenic variants on CLN3 transcription and disease pathogenesis.

Glioblastoma multiforme (GBM) is a highly aggressive brain tumor with a poor prognosis. The cleavage factor Im 25 (CFIm25), a crucial component of the CFIm complex, plays a key role in regulating the length of the mRNA 3'-UTR and has been implicated in various cancers, including GBM. This study sought to investigate the regulatory influence of specific microRNAs (miRNAs) on CFIm25 expression in GBM, a highly aggressive brain tumor. Bioinformatics analysis identified miRNA candidates targeting CFIm25 mRNA, and gene expression profiles from the NCBI database (GSE90603) were used for further analysis. Expression levels of CFIm25 and selected miRNAs were assessed using qRT-PCR in GBM clinical samples (n = 20) and non-malignant brain tissues (n = 5). Additionally, the MTT assay was performed to examine the effect of miRNA overexpression on U251 cell viability. Lentivectors expressing the identified miRNAs were employed to experimentally validate their regulatory role on CFIm25 in U251 cell lines, and Western blot analysis was conducted to determine CFIm25 protein levels. We observed significantly increased levels of miR-23, miR-24, and miR-27 expression, associated with a marked reduction in CFIm25 expression in GBM samples compared to non-malignant brain tissues. In particular, overexpression of miR-23, miR-24, and miR-27 in U251 cells resulted in CFIm25 downregulation at both the mRNA and protein levels, while their inhibition increased CFIm25 and reduced cell proliferation. These observations strongly implicate miR-23, miR-24, and miR-27 in regulating CFIm25 expression in GBM, emphasizing their potential as promising therapeutic targets for enhancing treatment responses in glioblastoma.

Alternative cleavage and polyadenylation within introns (intronic APA) generate shorter mRNA isoforms; however, their physiological significance remains elusive. In this study, we developed a comprehensive workflow to analyze intronic APA profiles using the mammalian target of rapamycin (mTOR)-regulated transcriptome as a model system. Our investigation revealed two contrasting effects within the transcriptome in response to fluctuations in cellular mTOR activity: an increase in intronic APA for a subset of genes and a decrease for another subset of genes. The application of this workflow to RNA-seq data from The Cancer Genome Atlas demonstrated that this dichotomous intronic APA pattern is a consistent feature in transcriptomes across both normal tissues and various cancer types. Notably, our analyses of protein length changes resulting from intronic APA events revealed two distinct phenomena in proteome programming: a loss of functional domains due to significant changes in protein length or minimal alterations in C-terminal protein sequences within unstructured regions. Focusing on conserved intronic APA events across 10 different cancer types highlighted the prevalence of the latter cases in cancer transcriptomes, whereas the former cases were relatively enriched in normal tissue transcriptomes. These observations suggest potential, yet distinct, roles for intronic APA events during pathogenic processes and emphasize the abundance of protein isoforms with similar lengths in the cancer proteome. Furthermore, our investigation into the isoform-specific functions of JMJD6 intronic APA events supported the hypothesis that alterations in unstructured C-terminal protein regions lead to functional differences. Collectively, our findings underscore intronic APA events as a discrete molecular signature present in both normal tissues and cancer transcriptomes, highlighting the contribution of APA to the multifaceted functionality of the cancer proteome.

Functional variations in coding and noncoding RNAs are crucial in tumorigenesis, with cancer-specific alterations often resulting from chemical modifications and posttranscriptional processes mediated by enzymes. These RNA variations have been linked to tumor cell proliferation, growth, metastasis, and drug resistance and are valuable for identifying diagnostic or prognostic cancer biomarkers. The diversity of posttranscriptional RNA modifications, such as splicing, polyadenylation, methylation, and editing, is particularly significant due to their prevalence and impact on cancer progression. Additionally, other modifications, including RNA acetylation, circularization, miRNA isomerization, and pseudouridination, are recognized as key contributors to cancer development. Understanding the mechanisms underlying these RNA modifications in cancer can enhance our knowledge of cancer biology and facilitate the development of innovative therapeutic strategies. Targeting these RNA modifications and their regulatory enzymes may pave the way for novel RNA-based therapies, enabling tailored interventions for specific cancer subtypes. This review provides a comprehensive overview of the roles and mechanisms of various coding and noncoding RNA modifications in cancer progression and highlights recent advancements in RNA-based therapeutic applications.

Systemic sclerosis (SSc) is a life-threatening autoimmune disease characterized by widespread fibrosis in the skin and several internal organs. Nudix Hydrolase 21 (NUDT2 or CFIm25) downregulation in fibroblasts is known to play detrimental roles in both skin and lung fibrosis. This study aims to investigate the upstream mechanisms that lead to NUDT21 repression in skin fibrosis. We identified transforming growth factor β (TGFβ1) as the primary cytokine that downregulated NUDT21 in normal skin fibroblasts. In the bleomycin-induced dermal fibrosis model, consistent with the peak activation of TGFβ1 at the late fibrotic stage, NUDT21 was downregulated at this stage, and delayed NUDT21 knockdown during this fibrotic phase led to enhanced fibrotic response to bleomycin. Further investigation suggested TGFβ downregulated NUDT21 through microRNA (miRNA) 181a and 181b induction. Both miR-181a and miR-181b were elevated in bleomycin-induced skin fibrosis in mice and primary fibroblasts isolated from SSc patients, and they directly targeted NUDT21 and led to its downregulation in skin fibroblasts. Functional studies demonstrated that miR-181a and miR-181b inhibitors attenuated bleomycin-induced skin fibrosis in mice in association with decreased NUDT21 expression, while miR-181a and miR-181b mimics promoted bleomycin-induced fibrosis. Overall, these findings suggest a novel role for miR-181a/b in SSc pathogenesis by repressing NUDT21 expression.

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

Polyglutamine spinocerebellar ataxias (PolyQ SCAs) represent a group of monogenetic diseases in which the expanded polyglutamine repeats give rise to a mutated protein. The abnormally expanded polyglutamine protein produces aggregates and toxic species, causing neuronal dysfunction and neuronal death. The main symptoms of these disorders include progressive ataxia, motor dysfunction, oculomotor impairment, and swallowing problems. Nowadays, the current treatments are restricted to symptomatic alleviation, and no existing therapeutic strategies can reduce or stop the disease progression. Even though the origin of these disorders has been associated with polyglutamine-induced toxicity, RNA toxicity has recently gained relevance in polyQ SCAs molecular pathogenesis. Therefore, the research's focus on RNA metabolism has been increasing, especially on RNA-binding proteins (RBPs). The present review summarizes RNA metabolism, exposing the different processes and the main RBPs involved. We also explore the mechanisms by which RBPs are dysregulated in PolyQ SCAs. Finally, possible therapies targeting the RNA metabolism are presented as strategies to reverse neuropathological anomalies and mitigate physical symptoms.

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis-acting loci, trans-acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis-acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis-loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis- and trans-acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

Recent development of RNA velocity uses master equations to establish the kinetics of the life cycle of RNAs from unspliced RNA to spliced RNA (i.e., mature RNA) to degradation. To feed this kinetic analysis, simultaneous measurement of unspliced RNA and spliced RNA in single cells is greatly desired. However, the majority of single-cell RNA-seq chemistry primarily captures mature RNA species to measure gene expressions. Here, we develop a one-step total-RNA chemistry-based single-cell RNA-seq method: snapTotal-seq. We benchmark this method with multiple single-cell RNA-seq assays in their performance in kinetic analysis of cell cycle by RNA velocity. Next, with LASSO regression between transcription factors, we identify the critical regulatory hubs mediating the cell cycle dynamics. We also apply snapTotal-seq to profile the oncogene-induced senescence and identify the key regulatory hubs governing the entry of senescence. Furthermore, from the comparative analysis of unspliced RNA and spliced RNA, we identify a significant portion of genes whose expression changes occur in spliced RNA but not to the same degree in unspliced RNA, indicating these gene expression changes are mainly controlled by post-transcriptional regulation. Overall, we demonstrate that snapTotal-seq can provide enriched information about gene regulation, especially during the transition between cell states.

Microscopic understanding of protein-RNA interactions is important for different biological activities, such as RNA transport, translation, splicing, silencing, etc. Polyadenine (Poly(A)) binding proteins (PABPs) make up a class of regulatory proteins that play critical roles in protecting the poly(A) tails of cellular mRNAs from nuclease degradation. In this work, we performed molecular dynamics simulations to investigate the conformational modifications of human PABP protein and poly(A) RNA that occur during complexation. It is demonstrated that the intermediate linker domain of the protein transforms from a disordered coil-like structure to a helical form during the recognition process, leading to the formation of the complex. On the other hand, disordered collapsed coil-like RNA on complexation has been found to transform into a rigid extended conformation. Importantly, the binding free energy calculation showed that the thermodynamic stability of the complex is primarily guided by favorable hydrophobic interactions between the protein and the RNA.

To explore the association between apaQTL/eQTL-SNPs and the susceptibility to silicosis. A silicosis-related GWAS was initially conducted to screen for single nucleotide polymorphisms (SNPs) associated with the risk of silicosis. Candidate SNPs with apaQTL and eQTL functions were then obtained from the 3'aQTL-atlas and GTEx databases. Subsequently, additional case-control studies were performed to validate the relationship between the candidate apaQTL/eQTL-SNPs and the risk of silicosis. Finally, experiments were conducted to illustrate APA events occurring at different alleles of the identified apaQTL/eQTL-SNPs. The combined results of the GWAS and iMLDR validations indicate that the variant T allele of the rs2974341 located on SMIM19 (additive model: OR = 0.66, the 95% CI = 0.53-0.84, P = 0.001) and the variant T allele of the rs2390488 located on TMTC4 (additive model: OR = 0.72, 95% CI = 0.57-0.90, P = 0.005) were significantly associated with decreased risk of developing silicosis susceptibility. Furthermore, 3'RACE experiments verified the presence of two poly (A) sites (proximal and distal) in SMIM19, rs2974341 may remotely regulate the binding between miRNA-3646 and SMIM19 with its high LD locus rs2974353 to affect the expression level of SMIM19. The rs2974341 variant T allele may contribute to the generation of the shorter 3'UTR transcript of SMIM19 and affect the binding of miRNA-3646 to the target gene SMIM19. The apaQTL/eQTL-SNPs may provide new perspectives for evaluating the regulatory function of SNPs in the development of silicosis.

As an important mechanism in the post-transcriptional regulation of eukaryotic gene expression, alternative polyadenylation (APA) plays a key role in biological processes such as cell proliferation and differentiation. However, the role and dynamic pattern of APA during Litopenaeus vannamei metamorphosis are poorly understood. Here, RNA-seq data covering from the embryo to the maturation (16 time points) of L. vannamei were utilized. We identified 247 differentially expressed APA events between early and adult stages, and through fuzzy mean clustering analysis, we discovered five dynamic APA patterns. Among them, the gradual elongation of the 3'UTR is the major APA pattern that changes over time, and its genes are enriched in the pathways of protein and energy metabolism. Finally, we constructed mRNA-miRNA and PPI networks and detected several central miRNAs that may regulate L. vannamei development. Our results revealed the complex APA mechanisms in L. vannamei metamorphosis, shedding new light on post-transcriptional regulation of crustacean metamorphosis.

Our ability to edit genomes lags behind our capacity to sequence them, but the growing understanding of CRISPR biology and its application to genome, epigenome and transcriptome engineering is narrowing this gap. In this Review, we discuss recent developments of various CRISPR-based systems that can transiently or permanently modify the genome and the transcriptome. The discovery of further CRISPR enzymes and systems through functional metagenomics has meaningfully broadened the applicability of CRISPR-based editing. Engineered Cas variants offer diverse capabilities such as base editing, prime editing, gene insertion and gene regulation, thereby providing a panoply of tools for the scientific community. We highlight the strengths and weaknesses of current CRISPR tools, considering their efficiency, precision, specificity, reliance on cellular DNA repair mechanisms and their applications in both fundamental biology and therapeutics. Finally, we discuss ongoing clinical trials that illustrate the potential impact of CRISPR systems on human health.

RNA binding proteins (RBPs) play a central in the post-transcriptional regulation of gene expression, which can account for up to 50% of all variations in protein expression within a cell. Following their binding to target RNAs, RBPs most typically confer changes in gene expression through modulation of alternative spicing, RNA stabilization/degradation, or ribosome loading/translation rate. All of these post-transcriptional regulatory processes have been shown to play a functional role in pathological cardiac remodeling, and a growing body of evidence is beginning to identify the mechanistic contribution of individual RBPs and their cardiac RNA targets. This review highlights the mechanisms of RBP-dependent post-transcriptional gene regulation in cardiomyocytes and fibroblasts and our current understanding of how RNA binding proteins functionally contribute to pathological cardiac remodeling.

The length of 3' untranslated regions (3'UTR) is highly regulated during many transitions in cell state, including T cell activation, through the process of alternative polyadenylation (APA). However, the regulatory mechanisms and functional consequences of APA remain largely unexplored. Here we present a detailed analysis of the temporal and condition-specific regulation of APA following activation of primary human CD4+ T cells. We find that global APA changes are regulated temporally and CD28 costimulatory signals enhance a subset of these changes. Most APA changes upon T cell activation involve 3'UTR shortening, although a set of genes enriched for function in the mTOR pathway exhibit 3'UTR lengthening. While upregulation of the core polyadenylation machinery likely induces 3'UTR shortening following prolonged T cell stimulation; a significant program of APA changes occur prior to cellular proliferation or upregulation of the APA machinery. Motif analysis suggests that at least a subset of these early changes in APA are driven by upregulation of RBM3, an RNA-binding protein which competes with the APA machinery for binding. Together this work expands our understanding of the impact and mechanisms of APA in response to T cell activation and suggests new mechanisms by which APA may be regulated.

Messenger RNA precursors (pre-mRNA) generally undergo 3' end processing by cleavage and polyadenylation (CPA), which is specified by a polyadenylation site (PAS) and adjacent RNA sequences and regulated by a large variety of core and auxiliary CPA factors. To date, most of the human CPA factors have been discovered through biochemical and proteomic studies. However, genetic identification of the human CPA factors has been hampered by the lack of a reliable genome-wide screening method. We describe here a dual fluorescence readthrough reporter system with a PAS inserted between two fluorescent reporters. This system enables measurement of the efficiency of 3' end processing in living cells. Using this system in combination with a human genome-wide CRISPR/Cas9 library, we conducted a screen for CPA factors. The screens identified most components of the known core CPA complexes and other known CPA factors. The screens also identified CCNK/CDK12 as a potential core CPA factor, and RPRD1B as a CPA factor that binds RNA and regulates the release of RNA polymerase II at the 3' ends of genes. Thus, this dual fluorescence reporter coupled with CRISPR/Cas9 screens reliably identifies bona fide CPA factors and provides a platform for investigating the requirements for CPA in various contexts.