An extensive network of coupling among gene expression machines

eIF4E orchestrates mRNA processing, RNA export and translation to modify specific protein production

The eukaryotic translation initiation factor eIF4E acts as a multifunctional factor that simultaneously influences mRNA processing, export, and translation in many organisms. Its multifactorial effects are derived from its capacity to bind to the methyl-7-guanosine cap on the 5'end of mRNAs and thus can act as a cap chaperone for transcripts in the nucleus and cytoplasm. In this review, we describe the multifactorial roles of eIF4E in major mRNA-processing events including capping, splicing, cleavage and polyadenylation, nuclear export and translation. We discuss the evidence that eIF4E acts at two levels to generate widescale changes to processing, export and ultimately the protein produced. First, eIF4E alters the production of components of the mRNA processing machinery, supporting a widescale reprogramming of multiple mRNA processing events. In this way, eIF4E can modulate mRNA processing without physically interacting with target transcripts. Second, eIF4E also physically interacts with both capped mRNAs and components of the RNA processing or translation machineries. Further, specific mRNAs are sensitive to eIF4E only in particular mRNA processing events. This selectivity is governed by the presence of cis-acting elements within mRNAs known as USER codes that recruit relevant co-factors engaging the appropriate machinery. In all, we describe the molecular bases for eIF4E's multifactorial function and relevant regulatory pathways, discuss the basis for selectivity, present a compendium of ~80 eIF4E-interacting factors which play roles in these activities and provide an overview of the relevance of its functions to its oncogenic potential. Finally, we summarize early-stage clinical studies targeting eIF4E in cancer.

Nanodynamo quantifies subcellular RNA dynamics revealing extensive coupling between steps of the RNA life cycle

The coordinated action of transcriptional and post-transcriptional machineries shapes gene expression programs at steady state and determines their concerted response to perturbations. We have developed Nanodynamo, an experimental and computational workflow for quantifying the kinetic rates of nuclear and cytoplasmic steps of the RNA life cycle. Nanodynamo is based on mathematical modelling following sequencing of native RNA from cellular fractions and polysomes. We have applied this workflow to triple-negative breast cancer cells, revealing widespread post-transcriptional RNA processing that is mutually exclusive with its co-transcriptional counterpart. We used Nanodynamo to unravel the coupling between transcription, processing, export, decay and translation machineries. We have identified a number of coupling interactions within and between the nucleus and cytoplasm that largely contribute to coordinating how cells respond to perturbations that affect gene expression programs. Nanodynamo will be instrumental in unravelling the determinants and regulatory processes involved in the coordination of gene expression responses.

Predicting splicing patterns from the transcription factor binding sites in the promoter with deep learning

BACKGROUND

Alternative splicing is a pivotal mechanism of post-transcriptional modification that contributes to the transcriptome plasticity and proteome diversity in metazoan cells. Although many splicing regulations around the exon/intron regions are known, the relationship between promoter-bound transcription factors and the downstream alternative splicing largely remains unexplored.

RESULTS

In this study, we present computational approaches to unravel the regulatory relationship between promoter-bound transcription factor binding sites (TFBSs) and the splicing patterns. We curated a fine dataset that includes DNase I hypersensitive site sequencing and transcriptomes across fifteen human tissues from ENCODE. Specifically, we proposed different representations of TF binding context and splicing patterns to examine the associations between the promoter and downstream splicing events. While machine learning models demonstrated potential in predicting splicing patterns based on TFBS occupancies, the limitations in the generalization of predicting the splicing forms of singleton genes across diverse tissues was observed with carefully examination using different cross-validation methods. We further investigated the association between alterations in individual TFBS at promoters and shifts in exon splicing efficiency. Our results demonstrate that the convolutional neural network (CNN) models, trained on TF binding changes in the promoters, can predict the changes in splicing patterns. Furthermore, a systemic in silico substitutions analysis on the CNN models highlighted several potential splicing regulators. Notably, using empirical validation using K562 CTCFL shRNA knock-down data, we showed the significant role of CTCFL in splicing regulation.

CONCLUSION

In conclusion, our finding highlights the potential role of promoter-bound TFBSs in influencing the regulation of downstream splicing patterns and provides insights for discovering alternative splicing regulations.

Minor Spliceosomal 65K/RNPC3 Interacts with ANKRD11 and Mediates HDAC3-Regulated Histone Deacetylation and Transcription

RNA splicing is crucial in the multilayer regulatory networks for gene expression, making functional interactions with DNA- and other RNA-processing machineries in the nucleus. However, these established couplings are all major spliceosome-related; whether the minor spliceosome is involved remains unclear. Here, through affinity purification using Drosophila lysates, an interaction is identified between the minor spliceosomal 65K/RNPC3 and ANKRD11, a cofactor of histone deacetylase 3 (HDAC3). Using a CRISPR/Cas9 system, Deletion strains are constructed and found that both Dm65KΔ/Δ and Dmankrd11Δ/Δ mutants have reduced histone deacetylation at Lys9 of histone H3 (H3K9) and Lys5 of histone H4 (H4K5) in their heads, exhibiting various neural-related defects. The 65K-ANKRD11 interaction is also conserved in human cells, and the HsANKRD11 middle-uncharacterized domain mediates Hs65K association with HDAC3. Cleavage under targets and tagmentation (CUT&Tag) assays revealed that HsANKRD11 is a bridging factor, which facilitates the synergistic common chromatin-binding of HDAC3 and Hs65K. Knockdown (KD) of HsANKRD11 simultaneously decreased their common binding, resulting in reduced deacetylation of nearby H3K9. Ultimately, this study demonstrates that expression changes of many genes caused by HsANKRD11-KD are due to the decreased common chromatin-binding of HDAC3 and Hs65K and subsequently reduced deacetylation of H3K9, illustrating a novel and conserved coupling mechanism that links the histone deacetylation with minor spliceosome for the regulation of gene expression.

The RNA-binding protein Nab2 regulates levels of the RhoGEF Trio to govern axon and dendrite morphology

The Drosophila RNA-binding protein (RBP) Nab2 acts in neurons to regulate neurodevelopment and is orthologous to the human intellectual disability-linked RBP, ZC3H14. Nab2 governs axon projection in mushroom body neurons and limits dendritic arborization of class IV sensory neurons in part by regulating splicing events in ∼150 mRNAs. Analysis of the Sex-lethal (Sxl) mRNA revealed that Nab2 promotes an exon-skipping event and regulates m6A methylation on Sxl pre-mRNA by the Mettl3 methyltransferase. Mettl3 heterozygosity broadly rescues Nab2null phenotypes implying that Nab2 acts through similar mechanisms on other RNAs, including unidentified targets involved in neurodevelopment. Here, we show that Nab2 and Mettl3 regulate the removal of a 5'UTR (untranslated region) intron in the trio pre-mRNA. Trio utilizes two GEF domains to balance Rac and RhoGTPase activity. Intriguingly, an isoform of Trio containing only the RhoGEF domain, GEF2, is depleted in Nab2null nervous tissue. Expression of Trio-GEF2 rescues projection defects in Nab2null axons and dendrites, while the GEF1 Rac1-regulatory domain exacerbates these defects, suggesting Nab2-mediated regulation Trio-GEF activities. Collectively, these data indicate that Nab2-regulated processing of trio is critical for balancing Trio-GEF1 and -GEF2 activity and show that Nab2, Mettl3, and Trio function in a common pathway that shapes axon and dendrite morphology.

Arabidopsis mRNA export factor MOS11: molecular interactions and role in abiotic stress responses

Transcription and export (TREX) is a multi-subunit complex that links synthesis, processing and export of mRNAs. It interacts with the RNA helicase UAP56 and export factors such as MOS11 and ALYs to facilitate nucleocytosolic transport of mRNAs. Plant MOS11 is a conserved, but sparsely researched RNA-binding export factor, related to yeast Tho1 and mammalian CIP29/SARNP. Using biochemical approaches, the domains of Arabidopsis thaliana MOS11 required for interaction with UAP56 and RNA-binding were identified. Further analyses revealed marked genetic interactions between MOS11 and ALY genes. Cell fractionation in combination with transcript profiling demonstrated that MOS11 is required for export of a subset of mRNAs that are shorter and more GC-rich than MOS11-independent transcripts. The central α-helical domain of MOS11 proved essential for physical interaction with UAP56 and for RNA-binding. MOS11 is involved in the nucleocytosolic transport of mRNAs that are upregulated under stress conditions and accordingly mos11 mutant plants turned out to be sensitive to elevated NaCl concentrations and heat stress. Collectively, our analyses identify functional interaction domains of MOS11. In addition, the results establish that mRNA export is critically involved in the plant response to stress conditions and that MOS11 plays a prominent role at this.

HPV and RNA Binding Proteins: What We Know and What Remains to Be Discovered

Papillomavirus gene regulation is largely post-transcriptional due to overlapping open reading frames and the use of alternative polyadenylation and alternative splicing to produce the full suite of viral mRNAs. These processes are controlled by a wide range of cellular RNA binding proteins (RPBs), including constitutive splicing factors and cleavage and polyadenylation machinery, but also factors that regulate these processes, for example, SR and hnRNP proteins. Like cellular RNAs, papillomavirus RNAs have been shown to bind many such proteins. The life cycle of papillomaviruses is intimately linked to differentiation of the epithelial tissues the virus infects. For example, viral late mRNAs and proteins are expressed only in the most differentiated epithelial layers to avoid recognition by the host immune response. Papillomavirus genome replication is linked to the DNA damage response and viral chromatin conformation, processes which also link to RNA processing. Challenges with respect to elucidating how RBPs regulate the viral life cycle include consideration of the orchestrated spatial aspect of viral gene expression in an infected epithelium and the epigenetic nature of the viral episomal genome. This review discusses RBPs that control viral gene expression, and how the connectivity of various nuclear processes might contribute to viral mRNA production.

Structured and disordered regions of Ataxin-2 contribute differently to the specificity and efficiency of mRNP granule formation

Ataxin-2 (ATXN2) is a gene implicated in spinocerebellar ataxia type II (SCA2), amyotrophic lateral sclerosis (ALS) and Parkinsonism. The encoded protein is a therapeutic target for ALS and related conditions. ATXN2 (or Atx2 in insects) can function in translational activation, translational repression, mRNA stability and in the assembly of mRNP-granules, a process mediated by intrinsically disordered regions (IDRs). Previous work has shown that the LSm (Like-Sm) domain of Atx2, which can help stimulate mRNA translation, antagonizes mRNP-granule assembly. Here we advance these findings through a series of experiments on Drosophila and human Ataxin-2 proteins. Results of Targets of RNA Binding Proteins Identified by Editing (TRIBE), co-localization and immunoprecipitation experiments indicate that a polyA-binding protein (PABP) interacting, PAM2 motif of Ataxin-2 may be a major determinant of the mRNA and protein content of Ataxin-2 mRNP granules. Experiments with transgenic Drosophila indicate that while the Atx2-LSm domain may protect against neurodegeneration, structured PAM2- and unstructured IDR- interactions both support Atx2-induced cytotoxicity. Taken together, the data lead to a proposal for how Ataxin-2 interactions are remodelled during translational control and how structured and non-structured interactions contribute differently to the specificity and efficiency of RNP granule condensation as well as to neurodegeneration.

Nuclear mRNA decay: regulatory networks that control gene expression

Proper regulation of mRNA production in the nucleus is critical for the maintenance of cellular homoeostasis during adaptation to internal and environmental cues. Over the past 25 years, it has become clear that the nuclear machineries governing gene transcription, pre-mRNA processing, pre-mRNA and mRNA decay, and mRNA export to the cytoplasm are inextricably linked to control the quality and quantity of mRNAs available for translation. More recently, an ever-expanding diversity of new mechanisms by which nuclear RNA decay factors finely tune the expression of protein-encoding genes have been uncovered. Here, we review the current understanding of how mammalian cells shape their protein-encoding potential by regulating the decay of pre-mRNAs and mRNAs in the nucleus.

Toll like receptor 4 (TLR4) gene polymorphism and its association with somatic cell score and milk production traits in Indian dromedary camels

Our study aimed to explore the genetic variation in the Toll-like receptor 4 (TLR4) gene and establish its association with somatic cell score (SCS) and milk production traits in four Indian camel breeds namely Bikaneri, Kachchhi, Jaisalmeri and Mewari. TLR4 gene fragment of 573 bp spanning 5' UTR, exon-1 and partial intron-1 region was amplified and genotyped using the PCR-sequence based typing method. Only one SNP located at position C472T was identified. Genotyping revealed two alleles (C and T) and three genotypes: CC, CT and TT. The genotype frequencies for CC, CT and TT were 0.116, 0.326 and 0.558 and allele frequencies for C and T alleles were 0.279 and 0.721, respectively. Association study inferred that the effect of genotype on SCS, lactation yield (LY) and peak yield (PY) was non-significant however heterozygote (CT) genotypes recorded lower SCS and higher LY and PY. It can be concluded that the TLR4 gene possesses limited genetic variation, depicting polymorphism at a single locus in Indian camel breeds with a predominance of the TT genotype. The association study indicated that heterozygote animals possess better udder health and production performance, the statistical significance of which needs to be established using a large data set.

RNA Pol II-dependent transcription efficiency fine-tunes A-to-I editing levels

A-to-I RNA editing is a widespread epitranscriptomic phenomenon leading to the conversion of adenosines to inosines, which are primarily interpreted as guanosines by cellular machines. Consequently, A-to-I editing can alter splicing or lead to recoding of transcripts. As misregulation of editing can cause a variety of human diseases, A-to-I editing requires tight regulation of the extent of deamination, particularly in protein-coding regions. The bulk of A-to-I editing occurs cotranscriptionally. Thus, we studied A-to-I editing regulation in the context of transcription and pre-mRNA processing. We show that stimulation of transcription impacts editing levels. Activation of the transcription factor MYC leads to an up-regulation of A-to-I editing, particularly in transcripts that are suppressed upon MYC activation. Moreover, low pre-mRNA synthesis rates and low pre-mRNA expression levels support high levels of editing. We also show that editing levels greatly differ between nascent pre-mRNA and mRNA in a cellular system, as well as in mouse tissues. Editing levels can increase or decrease from pre-mRNA to mRNA and can vary across editing targets and across tissues, showing that pre-mRNA processing is an important layer of editing regulation. Several lines of evidence suggest that the differences emerge during pre-mRNA splicing. Moreover, actinomycin D treatment of primary neuronal cells and editing level analysis suggests that regulation of editing levels also depends on transcription.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

The eukaryotic translation initiation factor eIF4E unexpectedly acts in splicing thereby coupling mRNA processing with translation: eIF4E induces widescale splicing reprogramming providing system-wide connectivity between splicing, nuclear mRNA export and translation

Recent findings position the eukaryotic translation initiation factor eIF4E as a novel modulator of mRNA splicing, a process that impacts the form and function of resultant proteins. eIF4E physically interacts with the spliceosome and with some intron-containing transcripts implying a direct role in some splicing events. Moreover, eIF4E drives the production of key components of the splicing machinery underpinning larger scale impacts on splicing. These drive eIF4E-dependent reprogramming of the splicing signature. This work completes a series of studies demonstrating eIF4E acts in all the major mRNA maturation steps whereby eIF4E drives production of the RNA processing machinery and escorts some transcripts through various maturation steps. In this way, eIF4E couples the mRNA processing-export-translation axis linking nuclear mRNA processing to cytoplasmic translation. eIF4E elevation is linked to worse outcomes in acute myeloid leukemia patients where these activities are dysregulated. Understanding these effects provides new insight into post-transcriptional control and eIF4E-driven cancers.

The haplotype-resolved telomere-to-telomere carnation (Dianthus caryophyllus) genome reveals the correlation between genome architecture and gene expression

Carnation (Dianthus caryophyllus) is one of the most valuable commercial flowers, due to its richness of color and form, and its excellent storage and vase life. The diverse demands of the market require faster breeding in carnations. A full understanding of carnations is therefore required to guide the direction of breeding. Hence, we assembled the haplotype-resolved gap-free carnation genome of the variety 'Baltico', which is the most common white standard variety worldwide. Based on high-depth HiFi, ultra-long nanopore, and Hi-C sequencing data, we assembled the telomere-to-telomere (T2T) genome to be 564 479 117 and 568 266 215 bp for the two haplotypes Hap1 and Hap2, respectively. This T2T genome exhibited great improvement in genome assembly and annotation results compared with the former version. The improvements were seen when different approaches to evaluation were used. Our T2T genome first informs the analysis of the telomere and centromere region, enabling us to speculate about specific centromere characteristics that cannot be identified by high-order repeats in carnations. We analyzed allele-specific expression in three tissues and the relationship between genome architecture and gene expression in the haplotypes. This demonstrated that the length of the genes, coding sequences, and introns, the exon numbers and the transposable element insertions correlate with gene expression ratios and levels. The insertions of transposable elements repress expression in gene regulatory networks in carnation. This gap-free finished T2T carnation genome provides a valuable resource to illustrate the genome characteristics and for functional genomics analysis in further studies and molecular breeding.

Genetic variants of CSN1S1, CSN2, CSN3, and BLG genes and their association with dairy production traits in Sahiwal cattle and Nili-Ravi buffaloes

Milk protein genes are associated with milk yield and composition in dairy animals. The present study aimed to identify milk protein genes (CSN1S1, CSN2, CSN3, and BLG) genetic variants and their association with milk yield in Sahiwal cattle and Nili-Ravi buffaloes. One hundred animals from each species were selected to collect blood samples and milk production records. Primers were designed for these milk protein genes for PCR amplification. Sequencing of resultant PCR products revealed a higher number of SNPs (13 vs. 7, 5 vs. 1, and 6 vs. 2) in Sahiwal as compared to Nili-Ravi animals in CSN1S1, CSN2, and CSN3 genes, respectively. However, a single SNP was observed in BLG gene of both species. Association analysis revealed that one SNP in BLG gene of Nili-Ravi was associated (p < 0.05) with 305-day milk yield. Two SNPs at CSN1S1 gene in Sahiwal were associated with dry-period. Similarly, one SNP at CSN1S1 and two SNPs at CSN3 gene showed significant association (p < 0.05) with average calving-interval in Sahiwal while two SNPs in CSN1S1 gene were associated (p < 0.05) with this trait in Nili-Ravi. These SNPs could be helpful as candidate variants for marker-assisted selection in cattle and buffaloes for improvement of lactation performance.

A dual role of RBM42 in modulating splicing and translation of CDKN1A/p21 during DNA damage response

p53-mediated cell cycle arrest during DNA damage is dependent on the induction of p21 protein, encoded by the CDKN1A gene. p21 inhibits cyclin-dependent kinases required for cell cycle progression to guarantee accurate repair of DNA lesions. Hence, fine-tuning of p21 levels is crucial to preserve genomic stability. Currently, the multilayered regulation of p21 levels during DNA damage is not fully understood. Herein, we identify the human RNA binding motif protein 42 (RBM42) as a regulator of p21 levels during DNA damage. Genome-wide transcriptome and interactome analysis reveals that RBM42 alters the expression of p53-regulated genes during DNA damage. Specifically, we demonstrate that RBM42 facilitates CDKN1A splicing by counteracting the splicing inhibitory effect of RBM4 protein. Unexpectedly, we also show that RBM42, underpins translation of various splicing targets, including CDKN1A. Concordantly, transcriptome-wide mapping of RBM42-RNA interactions using eCLIP further substantiates the dual function of RBM42 in regulating splicing and translation of its target genes, including CDKN1A. Collectively, our data show that RBM42 couples splicing and translation machineries to fine-tune gene expression during DNA damage response.

The polyA tail facilitates splicing of last introns with weak 3' splice sites via PABPN1

The polyA tail of mRNAs is important for many aspects of RNA metabolism. However, whether and how it regulates pre-mRNA splicing is still unknown. Here, we report that the polyA tail acts as a splicing enhancer for the last intron via the nuclear polyA binding protein PABPN1 in HeLa cells. PABPN1-depletion induces the retention of a group of introns with a weaker 3' splice site, and they show a strong 3'-end bias and mainly locate in nuclear speckles. The polyA tail is essential for PABPN1-enhanced last intron splicing and functions in a length-dependent manner. Tethering PABPN1 to nonpolyadenylated transcripts also promotes splicing, suggesting a direct role for PABPN1 in splicing regulation. Using TurboID-MS, we construct the PABPN1 interactome, including many spliceosomal and RNA-binding proteins. Specifically, PABPN1 can recruit RBM26&27 to promote splicing by interacting with the coiled-coil and RRM domain of RBM27. PABPN1-regulated terminal intron splicing is conserved in mice. Together, our study establishes a novel mode of post-transcriptional splicing regulation via the polyA tail and PABPN1.

Hepatotoxicity evaluation and possible mechanisms of decabrominated diphenyl ethers (BDE-209) in broilers: Oxidative stress, inflammatory, and transcriptomics

Decabrominated diphenyl ether (BDE-209), a persistent organic pollutant, is linked to a great number of health problems, the most severe of which impact the liver due to its role in the elimination and degradation of exogenous harmful substances. Though the hepatotoxicity of BDE-209 has been observed, its underlying mechanism is yet unknown. The purpose of this study is to thoroughly investigate the hepatotoxicity of BDE-209 and its molecular processes in broilers by subjecting 120 male broilers to varied concentrations of BDE-209 for 42 days. We observed that the bioaccumulation of BDE-209 in the liver in a dose-dependent manner, and that BDE-209 exposure can raise the concentrations of ALT, AST, and GGT, accompanied by hepatocyte fatty degeneration and inflammatory foci. In the hepatic homogenates, oxidative stress was evidenced by elevated levels of MDA and ROS and decreased activies of SOD and CAT. Additionally, pro-inflammatory cytokines including IL-1, IL-1β, TNF-α, IL-8 levels were increased, whereas anti-inflammatory cytokine IL-4 level was declined. Furthermore, RNA sequencing revealed that genes involved in inflammation were considerably dysregulated, and real-time PCR verified the expressed alterations of numerous genes related to the MAPK and WNT signaling pathways. The protein concentrations of NF-κB, β-catenin, and WNT5A, and the phosphorylation levels of JNK and ERK were all dramatically enhanced. The current study indicates that BDE-209 exposure can cause hepatotoxicity in broilers via bioaccumulation and oxidative stress, which then activates the MAPK and WNT signaling pathways, subsequently generating inflammation and hepatic injury.

Causal ALS genes impact the MHC class II antigen presentation pathway

Mutations in RNA/DNA-binding proteins cause amyotrophic lateral sclerosis (ALS), but the underlying disease mechanisms remain unclear. Here, we report that a set of ALS-associated proteins, namely FUS, EWSR1, TAF15, and MATR3, impact the expression of genes encoding the major histocompatibility complex II (MHC II) antigen presentation pathway. Both subunits of the MHC II heterodimer, HLA-DR, are down-regulated in ALS gene knockouts/knockdown in HeLa and human microglial cells, due to loss of the MHC II transcription factor CIITA. Importantly, hematopoietic progenitor cells (HPCs) derived from human embryonic stem cells bearing the FUSR495X mutation and HPCs derived from C9ORF72 ALS patient induced pluripotent stem cells also exhibit disrupted MHC II expression. Given that HPCs give rise to numerous immune cells, our data raise the possibility that loss of the MHC II pathway results in global failure of the immune system to protect motor neurons from damage that leads to ALS.

A Case-Control Study of the Relationship Between Genetic Polymorphism and Cretinism in Xinjiang

Background

Cretinism is a subtype of congenital hypothyroidism, an endocrine disorder resulting from inadequate thyroid hormone production or receptor deficiency. Genetic abnormalities play a major role in the development of thyroid dysfunction.

Methods

We recruited 183 participants with cretinism and 119 healthy participants from the Xinjiang Uyghur Autonomous Region and randomly selected 29 tag single nucleotide polymorphisms (tSNPs) in TSHB, PAX8, TPO, NKX2-5, and TSHR in all participants. We compared genotype and allele frequencies between cases and controls utilizing the chi-squared test, logistic regression analysis, and haplotype analysis.

Results

Using the chi-squared test, a single SNP was found to be associated with cretinism (recessive model: rs3754363, OR = 0.46, 95% CI = 0.27-0.80, P = 0.00519; genotype model: P = 0.01677). We stratified neurological, myxedematous, and mixed type and determined that another SNP was associated with a higher risk when comparing myxedematous type to the neurological type (rs2277923).

Conclusion

rs3754363 has a statistically significant protective effect on people with cretinism, while rs2277923 may play a greater role in promoting the development of neurocretinism.

Monomerization of TDP-43 is a key determinant for inducing TDP-43 pathology in amyotrophic lateral sclerosis

The cytoplasmic aggregation of TAR DNA binding protein-43 (TDP-43), also known as TDP-43 pathology, is the pathological hallmark of amyotrophic lateral sclerosis (ALS). However, the mechanism underlying TDP-43 cytoplasmic mislocalization and subsequent aggregation remains unclear. Here, we show that TDP-43 dimerization/multimerization is impaired in the postmortem brains and spinal cords of patients with sporadic ALS and that N-terminal dimerization-deficient TDP-43 consists of pathological inclusion bodies in ALS motor neurons. Expression of N-terminal dimerization-deficient mutant TDP-43 in Neuro2a cells and induced pluripotent stem cell-derived motor neurons recapitulates TDP-43 pathology, such as Nxf1-dependent cytoplasmic mislocalization and aggregate formation, which induces seeding effects. Furthermore, TDP-DiLuc, a bimolecular luminescence complementation reporter assay, could detect decreased N-terminal dimerization of TDP-43 before TDP-43 pathological changes caused by the transcription inhibition linked to aberrant RNA metabolism in ALS. These findings identified TDP-43 monomerization as a critical determinant inducing TDP-43 pathology in ALS.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Tumor educated platelet: the novel BioSource for cancer detection

Platelets, involved in the whole process of tumorigenesis and development, constantly absorb and enrich tumor-specific substances in the circulation during their life span, thus called "Tumor Educated Platelets" (TEPs). The alterations of platelet mRNA profiles have been identified as tumor markers due to the regulatory mechanism of post-transcriptional splicing. Small nuclear RNAs (SnRNAs), the important spliceosome components in platelets, dominate platelet RNA splicing and regulate the splicing intensity of pre-mRNA. Endogenous variation at the snRNA levels leads to widespread differences in alternative splicing, thereby driving the development and progression of neoplastic diseases. This review systematically expounds the bidirectional tumor-platelets interactions, especially the tumor induced alternative splicing in TEP, and further explores whether molecules related to alternative splicing such as snRNAs can serve as novel biomarkers for cancer diagnostics.

Mechanisms of RNA export and nuclear retention

With the identification of huge amount of noncoding RNAs in recent years, the concept of RNA localization has extended from traditional mRNA export to RNA export of mRNA and ncRNA as well as nuclear retention of ncRNA. This review aims to summarize the recent findings from studies on the mechanisms of export of different RNAs and nuclear retention of some lncRNAs in higher eukaryotes, with a focus on splicing-dependent TREX recruitment for the export of spliced mRNA and the sequence-dependent mechanism of mRNA export in the absence of splicing. In addition, evidence to support the involvement of m⁶ A modification in RNA export with the coordination between the methylase complex and TREX complex as well as sequence-dependent nuclear retention of lncRNA is recapitulated. Finally, a model of sequence-dependent RNA localization is proposed along with the many questions that remain to be answered. This article is categorized under: RNA Export and Localization > RNA Localization RNA Export and Localization > Nuclear Export/Import.

An integrative analysis of enhancer of yellow 2 homolog (ENY2) as a molecular biomarker in pan-cancer

ENY2 (Enhancer of yellow 2 transcription factor) is a transcription nuclear protein and primarily participates in the course of mRNA export and histone deubiquitination to influence gene expression. Current studies have shown that the expression of ENY2 is significantly upregulated in multiple cancers. However, the exact association between ENY2 and pan-cancers has not been fully established. Here, we comprehensively analyzed ENY2 from the online public database and The Cancer Genome Atlas (TCGA) database, including gene expression level in pan-cancer, comparison of ENY2 expression in different molecular and immune subtypes of pan-cancer, targeted protein, biological functions, molecular signatures, diagnostic and prognostic value in pan-cancer. Moreover, we focused on head and neck squamous cell carcinoma (HNSC) and explored ENY2 from the perspective of the correlations with clinical characteristics, prognosis, co-expression genes, differentially expressed genes (DEGs) and immune Infiltration. Our findings showed that the expression of ENY2 differed enormously not only in most cancer types but also in different molecular and immune subtypes of cancers. High accuracy in predicting cancers and notable correlations with prognosis of certain cancers suggested that ENY2 might be a potential diagnostic and prognostic biomarker of cancers. In addition, ENY2 was identified to be significantly correlated with clinical stage, gender, histologic grade and lymphovascular invasion in HNSC. Overexpression of ENY2 could lead to a worse overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) in HNSC, especially in different clinical subgroups of HNSC. Taken together, ENY2 showed strong correlation with the diagnosis and prognosis of pan-cancer, and was an independent prognostic risk factor of HNSC, which may serve as a potential target for cancer management.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

RNA modification-mediated translational control in immune cells

RNA modifications play a vital role in multiple pathways of mRNA metabolism, and translational regulation is essential for immune cells to promptly respond to stimuli and adapt to the microenvironment. N6-methyladenosine (m6A) methylation, which is the most abundant mRNA modification in eukaryotes, primarily functions in the regulation of RNA splicing and degradation. However, the role of m6Amethylation in translational control and its underlying mechanism remain controversial. The role of m6A methylation in translation regulation in immune cells has received relatively limited attention. In this review, we aim to provide a comprehensive summary of current studies on the translational regulation of m6A modifications and recent advances in understanding the translational control regulated by RNA modifications during the immune response. Furthermore, we envision the possible pathways through which m6A modifications may be involved in the regulation of immune cell function via translational control.

Roles of minor spliceosome in intron recognition and the convergence with the better understood major spliceosome

Catalyzed by spliceosomes in the nucleus, RNA splicing removes intronic sequences from precursor RNAs in eukaryotes to generate mature RNA, which also significantly increases proteome complexity and fine-tunes gene expression. Most metazoans have two coexisting spliceosomes; the major spliceosome, which removes >99.5% of introns, and the minor spliceosome, which removes far fewer introns (only 770 at present have been predicted in the human genome). Both spliceosomes are large and dynamic machineries, each consisting of five small nuclear RNAs (snRNAs) and more than 100 proteins. However, the dynamic assembly, catalysis, and protein composition of the minor spliceosome are still poorly understood. With different splicing signals, minor introns are rare and usually distributed alone and flanked by major introns in genes, raising questions of how they are recognized by the minor spliceosome and how their processing deals with the splicing of neighboring major introns. Due to large numbers of introns and close similarities between the two machinery, cooperative, and competitive recognition by the two spliceosomes has been investigated. Functionally, many minor-intron-containing genes are evolutionarily conserved and essential. Mutations in the minor spliceosome exhibit a variety of developmental defects in plants and animals and are linked to numerous human diseases. Here, we review recent progress in the understanding of minor splicing, compare currently known components of the two spliceosomes, survey minor introns in a wide range of organisms, discuss cooperation and competition of the two spliceosomes in splicing of minor-intron-containing genes, and contributions of minor splicing mutations in development and diseases. This article is categorized under: RNA Processing > Processing of Small RNAs RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Three-Dimensional Printing of an Apigenin-Loaded Mucoadhesive Film for Tailored Therapy to Oral Leukoplakia and the Chemopreventive Effect on a Rat Model of Oral Carcinogenesis

Oral leukoplakia, which presents as white lesions in the oral cavity, including on the tongue, is precancerous in nature. Conservative treatment is preferable, since surgical removal can markedly reduce the patient’s quality of life. In the present study, we focused on the flavonoid apigenin as a potential compound for preventing carcinogenesis, and an apigenin-loaded mucoadhesive oral film was prepared using a three-dimensional (3D) bioprinter (semi-solid extrusion-type 3D printer). Apigenin-loaded printer inks are composed of pharmaceutical excipients (HPMC, CARBOPOL, and Poloxamer), water, and ethanol to dissolve apigenin, and the appropriate viscosity of printer ink after adjusting the ratios allowed for the successful 3D printing of the film. After drying the 3D-printed object, the resulting film was characterized. The chemopreventive effect of the apigenin-loaded film was evaluated using an experimental rat model that had been exposed to 4-nitroquinoline 1-oxide (4NQO) to induce oral carcinogenesis. Treatment with the apigenin-loaded film showed a remarkable chemopreventive effect based on an analysis of the specimen by immunohistostaining. These results suggest that the apigenin-loaded mucoadhesive film may help prevent carcinogenesis. This successful preparation of apigenin-loaded films by a 3D printer provides useful information for automatically fabricating other tailored films (with individual doses and shapes) for patients with oral leukoplakia in a future clinical setting.

Finding Correlations Between mRNA and Protein Levels in Leishmania Development: Is There a Discrepancy?

Multiple genes and proteins have been identified as differentially expressed in the stages of the Leishmania life cycle. The differentiation processes are implicated in specific transcriptional and proteomic adjustments driven by gene expression regulation mechanisms. Leishmania parasites lack gene-specific transcriptional control, and gene expression regulation mostly depends on posttranscriptional mechanisms. Due to the lack of transcriptional regulation, criticism regarding the relevance of transcript quantification as a possible and efficient prediction of protein levels is recurrent in studies that use transcriptomic information. The advent of high-throughput technologies has improved the analysis of genomes, transcriptomes and proteomes for different organisms under several conditions. Nevertheless, defining the correlation between transcriptional and proteomic profiles requires arduous and expensive work and remains a challenge in Leishmania. In this review, we analyze transcriptomic and proteomic data for several Leishmania species in two different stages of the parasite life cycle: metacyclogenesis and amastigogenesis (amastigote differentiation). We found a correlation between mRNA and protein levels of 60.9% and 69.8% for metacyclogenesis and amastigogenesis, respectively; showing that majority mRNA and protein levels increase or decrease concomitantly. Among the analyzed genes that did not present correlation indicate that transcriptomic data should be carefully interpreted as protein expression. We also discuss possible explanations and mechanisms involved for this lack of correlation.

From Isles of Königsberg to Islets of Langerhans: Examining the Function of the Endocrine Pancreas Through Network Science

Islets of Langerhans are multicellular microorgans located in the pancreas that play a central role in whole-body energy homeostasis. Through secretion of insulin and other hormones they regulate postprandial storage and interprandial usage of energy-rich nutrients. In these clusters of hormone-secreting endocrine cells, intricate cell-cell communication is essential for proper function. Electrical coupling between the insulin-secreting beta cells through gap junctions composed of connexin36 is particularly important, as it provides the required, most important, basis for coordinated responses of the beta cell population. The increasing evidence that gap-junctional communication and its modulation are vital to well-regulated secretion of insulin has stimulated immense interest in how subpopulations of heterogeneous beta cells are functionally arranged throughout the islets and how they mediate intercellular signals. In the last decade, several novel techniques have been proposed to assess cooperation between cells in islets, including the prosperous combination of multicellular imaging and network science. In the present contribution, we review recent advances related to the application of complex network approaches to uncover the functional connectivity patterns among cells within the islets. We first provide an accessible introduction to the basic principles of network theory, enumerating the measures characterizing the intercellular interactions and quantifying the functional integration and segregation of a multicellular system. Then we describe methodological approaches to construct functional beta cell networks, point out possible pitfalls, and specify the functional implications of beta cell network examinations. We continue by highlighting the recent findings obtained through advanced multicellular imaging techniques supported by network-based analyses, giving special emphasis to the current developments in both mouse and human islets, as well as outlining challenges offered by the multilayer network formalism in exploring the collective activity of islet cell populations. Finally, we emphasize that the combination of these imaging techniques and network-based analyses does not only represent an innovative concept that can be used to describe and interpret the physiology of islets, but also provides fertile ground for delineating normal from pathological function and for quantifying the changes in islet communication networks associated with the development of diabetes mellitus.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Nuclear mRNA Export and Aging

The relationship between transcription and aging is one that has been studied intensively and experimentally with diverse attempts. However, the impact of the nuclear mRNA export on the aging process following its transcription is still poorly understood, although the nuclear events after transcription are coupled closely with the transcription pathway because the essential factors required for mRNA transport, namely TREX, TREX-2, and nuclear pore complex (NPC), physically and functionally interact with various transcription factors, including the activator/repressor and pre-mRNA processing factors. Dysregulation of the mediating factors for mRNA export from the nucleus generally leads to the aberrant accumulation of nuclear mRNA and further impairment in the vegetative growth and normal lifespan and the pathogenesis of neurodegenerative diseases. The optimal stoichiometry and density of NPC are destroyed during the process of cellular aging, and their damage triggers a defect of function in the nuclear permeability barrier. This review describes recent findings regarding the role of the nuclear mRNA export in cellular aging and age-related neurodegenerative disorders.

Genome-Wide Characterization, Evolutionary Analysis of ARF Gene Family, and the Role of SaARF4 in Cd Accumulation of Sedum alfredii Hance

Auxin response factors (ARFs) play important roles in plant development and environmental adaption. However, the function of ARFs in cadmium (Cd) accumulation are still unknown. Here, 23 SaARFs were detected in the genome of hyperaccumulating ecotype of Sedum alfredii Hance (HE), and they were not evenly distributed on the chromosomes. Their protein domains remained highly conservative. SaARFs in the phylogenetic tree can be divided into three groups. Genes in the group Ⅰ contained three introns at most. However, over ten introns were found in other two groups. Collinearity relationships were exhibited among ten SaARFs. The reasons for generating SaARFs may be segmental duplication and rearrangements. Collinearity analysis among different species revealed that more collinear genes of SaARFs can be found in the species with close relationships of HE. A total of eight elements in SaARFs promoters were related with abiotic stress. The qRT-PCR results indicated that four SaARFs can respond to Cd stress. Moreover, that there may be functional redundancy among six SaARFs. The adaptive selection and functional divergence analysis indicated that SaARF4 may undergo positive selection pressure and an adaptive-evolution process. Overexpressing SaARF4 effectively declined Cd accumulation. Eleven single nucleotide polymorphism (SNP) sites relevant to Cd accumulation can be detected in SaARF4. Among them, only one SNP site can alter the sequence of the SaARF4 protein, but the SaARF4 mutant of this site did not cause a significant difference in cadmium content, compared with wild-type plants. SaARFs may be involved in Cd-stress responses, and SaARF4 may be applied for decreasing Cd accumulation of plants.

Comprehensive Evolutionary Analysis of CPP Genes in Brassica napus L. and Its Two Diploid Progenitors Revealing the Potential Molecular Basis of Allopolyploid Adaptive Advantage Under Salt Stress

Allopolyploids exist widely in nature and have strong environmental adaptability. The typical allopolyploid Brassica napus L. is a widely cultivated crop, but whether it is superior to its diploid progenitors in abiotic stress resistance and the key genes that may be involved are not fully understood. Cystein-rich polycomb-like protein (CPP) genes encode critical transcription factors involved in the response of abiotic stress, including salt stress. To explore the potential molecular basis of allopolyploid adaptation to salt stress, we comprehensively analyzed the characteristics and salt stress response of the CPP genes in B. napus and its two diploid progenitors in this study. We found some molecular basis that might be associated with the adaptability of B. napus, including the expansion of the CPP gene family, the acquisition of introns by some BnCPPs, and abundant cis-acting elements upstream of BnCPPs. We found two duplication modes (whole genome duplication and transposed duplication) might be the main reasons for the expansion of CPP gene family in B. napus during allopolyploidization. CPP gene expression levels and several physiological indexes were changed in B. napus and its diploid progenitors after salt stress, suggesting that CPP genes might play important roles in the response of salt stress. We found that some BnCPPs might undergo new functionalization or subfunctionalization, and some BnCPPs also show biased expression, which might contribute to the adaptation of B. napus under saline environment. Compared with diploid progenitors, B. napus showed stronger physiological responses, and BnCPP gene expression also showed higher changes after salt stress, indicating that the allopolyploid B. napus had an adaptive advantage under salt stress. This study could provide evidence for the adaptability of polyploid and provide important clues for the study of the molecular mechanism of salt stress resistance in B. napus.

Non-Darwinian Molecular Biology

With the discovery of the double helical structure of DNA, a shift occurred in how biologists investigated questions surrounding cellular processes, such as protein synthesis. Instead of viewing biological activity through the lens of chemical reactions, this new field used biological information to gain a new profound view of how biological systems work. Molecular biologists asked new types of questions that would have been inconceivable to the older generation of researchers, such as how cellular machineries convert inherited biological information into functional molecules like proteins. This new focus on biological information also gave molecular biologists a way to link their findings to concepts developed by genetics and the modern synthesis. However, by the late 1960s this all changed. Elevated rates of mutation, unsustainable genetic loads, and high levels of variation in populations, challenged Darwinian evolution, a central tenant of the modern synthesis, where adaptation was the main driver of evolutionary change. Building on these findings, Motoo Kimura advanced the neutral theory of molecular evolution, which advocates that selection in multicellular eukaryotes is weak and that most genomic changes are neutral and due to random drift. This was further elaborated by Jack King and Thomas Jukes, in their paper “Non-Darwinian Evolution”, where they pointed out that the observed changes seen in proteins and the types of polymorphisms observed in populations only become understandable when we take into account biochemistry and Kimura’s new theory. Fifty years later, most molecular biologists remain unaware of these fundamental advances. Their adaptionist viewpoint fails to explain data collected from new powerful technologies which can detect exceedingly rare biochemical events. For example, high throughput sequencing routinely detects RNA transcripts being produced from almost the entire genome yet are present less than one copy per thousand cells and appear to lack any function. Molecular biologists must now reincorporate ideas from classical biochemistry and absorb modern concepts from molecular evolution, to craft a new lens through which they can evaluate the functionality of transcriptional units, and make sense of our messy, intricate, and complicated genome.

SUMOylation modulates the function of DDX19 in mRNA export

Nuclear export of mRNAs is a critical regulatory step in eukaryotic gene expression. The mRNA transcript undergoes extensive processing, and is loaded with a set of RNA-binding proteins (RBPs) to form export-competent messenger ribonucleoprotein particles (mRNPs) in the nucleus. During the transit of mRNPs through the nuclear pore complex (NPC), the DEAD-box ATPase - DDX19 - remodels mRNPs at the cytoplasmic side of the NPC, by removing a subset of RNA-binding proteins to terminate mRNP export. This requires the RNA-dependent ATPase activity of DDX19 and its dynamic interactions with Gle1 and Nup214. However, the regulatory mechanisms underlying these interactions are unclear. We find that DDX19 gets covalently attached with a small ubiquitin-like modifier (SUMO) at lysine 26, which enhances its interaction with Gle1. Furthermore, a SUMOylation-defective mutant of human DDX19B, K26R, failed to provide a complete rescue of the mRNA export defect caused by DDX19 depletion. Collectively, our results suggest that SUMOylation fine-tunes the function of DDX19 in mRNA export by regulating its interaction with Gle1. This study identifies SUMOylation of DDX19 as a modulatory mechanism during the mRNA export process.

Methyl-lysine readers PHF20 and PHF20L1 define two distinct gene expression-regulating NSL complexes

The methyl-lysine readers plant homeodomain finger protein 20 (PHF20) and its homolog PHF20-like protein 1 (PHF20L1) are known components of the nonspecific lethal (NSL) complex that regulates gene expression through its histone acetyltransferase activity. In the current model, both PHF homologs coexist in the same NSL complex, although this was not formally tested; nor have the functions of PHF20 and PHF20L1 regarding NSL complex integrity and transcriptional regulation been investigated. Here, we perform an in-depth biochemical and functional characterization of PHF20 and PHF20L1 in the context of the NSL complex. Using mass spectrometry, genome-wide chromatin analysis, and protein-domain mapping, we identify the existence of two distinct NSL complexes that exclusively contain either PHF20 or PHF20L1. We show that the C-terminal domains of PHF20 and PHF20L1 are essential for complex formation with NSL, and the Tudor 2 domains are required for chromatin binding. The genome-wide chromatin landscape of PHF20–PHF20L1 shows that these proteins bind mostly to the same genomic regions, at promoters of highly expressed/housekeeping genes. Yet, deletion of PHF20 and PHF20L1 does not abrogate gene expression or impact the recruitment of the NSL complex to those target gene promoters, suggesting the existence of an alternative mechanism that compensates for the transcription of genes whose sustained expression is important for critical cellular functions. This work shifts the current paradigm and lays the foundation for studies on the differential roles of PHF20 and PHF20L1 in regulating NSL complex activity in physiological and diseases states.

Interrogating the transcriptome with metabolically incorporated ribonucleosides

RNA is a central player in biological processes, but there remain major gaps in our understanding of transcriptomic processes and the underlying biochemical mechanisms regulating RNA in cells. A powerful strategy to facilitate molecular analysis of cellular RNA is the metabolic incorporation of chemical probes. In this review, we discuss current approaches for RNA metabolic labeling with modified ribonucleosides and their integration with Next-Generation Sequencing, mass spectrometry-based proteomics, and fluorescence microscopy in order to interrogate RNA behavior in its native context.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

Dynamic self-reinforcement of gene expression determines acquisition of cellular mechanical memory

Mechanotransduction describes activation of gene expression by changes in the cell's physical microenvironment. Recent experiments show that mechanotransduction can lead to long-term "mechanical memory," in which cells cultured on stiff substrates for sufficient time (priming phase) maintain altered phenotype after switching to soft substrates (dissipation phase) as compared to unprimed controls. The timescale of memory acquisition and retention is orders of magnitude larger than the timescale of mechanosensitive cellular signaling, and memory retention time changes continuously with priming time. We develop a model that captures these features by accounting for positive reinforcement in mechanical signaling. The sensitivity of reinforcement represents the dynamic transcriptional state of the cell composed of protein lifetimes and three-dimensional chromatin organization. Our model provides a single framework connecting microenvironment mechanical history to cellular outcomes ranging from no memory to terminal differentiation. Predicting cellular memory of environmental changes can help engineer cellular dynamics through changes in culture environments.

To export, or not to export: how nuclear export factor variants resolve Piwi's dilemma

Piwi-interacting RNAs (piRNAs) defend animal gonads by guiding PIWI-clade Argonaute proteins to silence transposons. The nuclear Piwi/piRNA complexes confer transcriptional repression of transposons, which is accompanied with heterochromatin formation at target loci. On the other hand, piRNA clusters, genomic loci that transcribe piRNA precursors composed of transposon fragments, are often recognized by piRNAs to define their heterochromatic identity. Therefore, Piwi/piRNA complexes must resolve this conundrum of silencing transposons while allowing the expression of piRNA precursors, at least in Drosophila germlines. This review is focused on recent advances how the piRNA pathway deals with this genetic conflict.

The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape

Viral infections are known to hijack the transcription and translation of the host cell. However, the extent to which viral proteins coordinate these perturbations remains unclear. Here we used a model system, the human T-cell leukemia virus type 1 (HTLV-1), and systematically analyzed the transcriptome and interactome of key effectors oncoviral proteins Tax and HBZ. We showed that Tax and HBZ target distinct but also common transcription factors. Unexpectedly, we also uncovered a large set of interactions with RNA-binding proteins, including the U2 auxiliary factor large subunit (U2AF2), a key cellular regulator of pre-mRNA splicing. We discovered that Tax and HBZ perturb the splicing landscape by altering cassette exons in opposing manners, with Tax inducing exon inclusion while HBZ induces exon exclusion. Among Tax- and HBZ-dependent splicing changes, we identify events that are also altered in Adult T cell leukemia/lymphoma (ATLL) samples from two independent patient cohorts, and in well-known cancer census genes. Our interactome mapping approach, applicable to other viral oncogenes, has identified spliceosome perturbation as a novel mechanism coordinated by Tax and HBZ to reprogram the transcriptome.

Tax and HBZ are two viral regulatory proteins encoded by the human T-cell leukemia virus type 1 (HTLV-1) via sense and antisense transcripts, respectively. Both proteins are known to drive oncogenic processes that culminate in a T-cell neoplasm, known as Adult T cell leukemia/lymphoma (ATLL). We measured the effects of Tax and HBZ on host gene expression pathway by analyzing the interactome with cellular transcriptional and post-transcriptional regulators, and the transcriptome and mRNA splicing of cell lines expressing either Tax or HBZ. We compared our results with data obtained from independent cohorts of Japanese and Afro-Caribbean patients, and identified common splicing changes that might represent clinically useful biomarkers for ATLL. Finally, we provide evidence that the viral protein Tax can reprogram initial steps of the T-cell transcriptome diversification by hijacking the U2AF complex, a key cellular regulator of pre-mRNA splicing.

Intron-retained radish (Raphanus sativus L.) RsMYB1 transcripts found in colored-taproot lines enhance anthocyanin accumulation in transgenic Arabidopsis plants

KEY MESSAGE

Overexpression of the naturally occurring intron-retained (IR) forms of radish RsMYB1 and RsTT8 transcripts in Arabidopsis causes a substantial increase in anthocyanin accumulation. The production of anthocyanins in plants is tightly controlled by the MYB-bHLH-WD40 (MBW) complex. In this study, analysis of four radish (Raphanus sativus L.) inbred lines with different colored taproots revealed that regulatory genes of anthocyanin biosynthesis, RsMYB1 and RsTT8, produce three transcripts, one completely spliced and two intron retention (IR1 and IR2) forms. Transcripts RsMYB1-IR1 and RsMYB1-IR2 retained the 1st (380 nt) and 2nd (149 nt) introns, respectively; RsTT8-IR1 retained the 4th intron (113 nt); RsTT8-IR2 retained both the 3rd (128 nt) and 4th introns. Levels of most IR forms were substantially low in radish samples, but the RsTT8-IR2 level was higher than RsTT8 in red skin/red flesh (RsRf) root. Since all IR forms contained a stop codon within the intron, they were predicted to encode truncated proteins with defective interaction domains, resulting in the inability to form the MBW complex in vivo. However, tobacco leaves transiently co-expressing RsMYB1-IRs and RsTT8-IRs showed substantially higher anthocyanin accumulation than those co-expressing their spliced forms. Consistently, co-expression of constructs encoding truncated proteins with spliced or IR forms of their interaction partner in tobacco leaves did not result in anthocyanin accumulation. Compared with RsMYB1, the overexpression of RsMYB1-IRs in Arabidopsis pap1 mutant increased anthocyanin accumulation by > sevenfold and upregulated the expression of Arabidopsis flavonoid biosynthesis genes including AtTT8. Our results suggest that the stable co-expression of RsMYB1-IRs in fruit trees and vegetable crops could be used to increase their anthocyanin contents.

Expression of a human cDNA in moss results in spliced mRNAs and fragmentary protein isoforms

Production of biopharmaceuticals relies on the expression of mammalian cDNAs in host organisms. Here we show that the expression of a human cDNA in the moss Physcomitrium patens generates the expected full-length and four additional transcripts due to unexpected splicing. This mRNA splicing results in non-functional protein isoforms, cellular misallocation of the proteins and low product yields. We integrated these results together with the results of our analysis of all 32,926 protein-encoding Physcomitrella genes and their 87,533 annotated transcripts in a web application, physCO, for automatized optimization. A thus optimized cDNA results in about twelve times more protein, which correctly localizes to the ER. An analysis of codon preferences of different production hosts suggests that similar effects occur also in non-plant hosts. We anticipate that the use of our methodology will prevent so far undetected mRNA heterosplicing resulting in maximized functional protein amounts for basic biology and biotechnology.

Deciphering the association of intronic single nucleotide polymorphisms of crystallin gene family with congenital cataract

Purpose

Introns play an important role in gene regulation and expression. Single nucleotide polymorphisms (SNPs) in introns have the potential to cause disease and alter the genotype-phenotype association. Hence, this study aimed to decipher the association of SNPs in the introns of the crystallin gene in congenital cataracts.

Methods

SNPs in the introns of crystallin gene family - CRYAA (rs3788059), CRYAB (rs2070894), CRYBA4 (rs2071861), and CRYBB2 (rs5752083, rs5996863) - were genotyped in 248 participants consisting of 141 congenital cataracts and 107 healthy controls by allele-specific oligonucleotide polymerase chain reaction method. Around 10% of samples for each SNPs were sequenced to confirm the genotypes. The allele, genotype, and haplotype frequency were evaluated by the SHEsis online tool.

Results

Using dominant model, the "A" allele of rs3788059 was found to have an increased risk toward congenital cataract development whereas the "G" allele was found to be protective (AA + AG vs. GG; odds ratio [OR] 95% confidence interval [CI] = 3.73 [1.71, 8.15], P = 0.0009). The "A" allele of both rs2070894 (AA + AG vs. GG; OR [95% CI] = 0.49 [0.29, 0.84], P = 0.012) and rs5752083 (AA + AC vs. CC; OR [95% CI] = 0.25 [0.08, 0.76], P = 0.016) were suggested to have a protective role by the dominant model. The A-C-T haplotype (rs2071861, rs5752083, and rs5996863) was found to be a significant risk factor for the development of congenital cataract.

Conclusion

Intronic SNPs in crystallin genes may play a role in the predisposition toward congenital cataract. However, the present findings need to be replicated in a large cohort with more number of samples.

VEGFA's distal enhancer regulates its alternative splicing in CML

Enhancer demethylation in leukemia has been shown to lead to overexpression of genes which promote cancer characteristics. The vascular endothelial growth factor A (VEGFA) enhancer, located 157 Kb downstream of its promoter, is demethylated in chronic myeloid leukemia (CML). VEGFA has several alternative splicing isoforms with different roles in cancer progression. Since transcription and splicing are coupled, we wondered whether VEGFA enhancer activity can also regulate the gene's alternative splicing to contribute to the pathology of CML. Our results show that mutating the VEGFA +157 enhancer promotes exclusion of exons 6a and 7 and activating the enhancer by tethering a chromatin activator has the opposite effect. In line with these results, CML patients present with high expression of +157 eRNA and inclusion of VEGFA exons 6a and 7. In addition, our results show that the positive regulator of RNAPII transcription elongation, CCNT2, binds VEGFA's promoter and enhancer, and its silencing promotes exclusion of exons 6a and 7 as it slows down RNAPII elongation rate. Thus our results suggest that VEGFA's +157 enhancer regulates its alternative splicing by increasing RNAPII elongation rate via CCNT2. Our work demonstrates for the first time a connection between an endogenous enhancer and alternative splicing regulation of its target gene.

The Role of RNA Modifications and RNA-modifying Proteins in Cancer Therapy and Drug Resistance

The advent of new genome-wide sequencing technologies has uncovered abnormal RNA modifications and RNA editing in a variety of human cancers. The discovery of reversible RNA N6-methyladenosine (RNA: m⁶A) by fat mass and obesity-associated protein (FTO) demethylase has led to exponential publications on the pathophysiological functions of m⁶A and its corresponding RNA modifying proteins (RMPs) in the past decade. Some excellent reviews have summarized the recent progress in this field. Compared to the extent of research into RNA: m⁶A and DNA 5-methylcytosine (DNA: m⁵C), much less is known about other RNA modifications and their associated RMPs, such as the role of RNA: m⁵C and its RNA cytosine methyltransferases (RCMTs) in cancer therapy and drug resistance. In this review, we will summarize the recent progress surrounding the function, intramolecular distribution and subcellular localization of several major RNA modifications, including 5' cap N7-methylguanosine (m7G) and 2'-O-methylation (Nm), m⁶A, m⁵C, A-to-I editing, and the associated RMPs. We will then discuss dysregulation of those RNA modifications and RMPs in cancer and their role in cancer therapy and drug resistance.

Co-transcriptional splicing efficiencies differ within genes and between cell types

Pre-mRNA splicing is carried out by the spliceosome and involves splice site recognition, removal of introns, and ligation of exons. Components of the spliceosome have been shown to interact with the elongating RNA polymerase II (RNAPII) which is thought to allow splicing to occur concurrently with transcription. However, little is known about the regulation and efficiency of co-transcriptional splicing in human cells. In this study, we used Bru-seq and BruChase-seq to determine the co-transcriptional splicing efficiencies of 17,000 introns expressed across 6 human cell lines. We found that less than half of all introns across these 6 cell lines were co-transcriptionally spliced. Splicing efficiencies for individual introns showed variations across cell lines, suggesting that splicing may be regulated in a cell-type specific manner. Moreover, the splicing efficiency of introns varied within genes. The efficiency of co-transcriptional splicing did not correlate with gene length, intron position, splice site strengths, or the intron/neighboring exons GC content. However, we identified binding signals from multiple RNA binding proteins (RBPs) that correlated with splicing efficiency, including core spliceosomal machinery components-such as SF3B4, U2AF1 and U2AF2 showing higher binding signals in poorly spliced introns. In addition, multiple RBPs, such as BUD13, PUM1 and SND1, showed preferential binding in exons that flank introns with high splicing efficiencies. The nascent RNA splicing patterns presented here across multiple cell types add to our understanding of the complexity in RNA splicing, wherein RNA-binding proteins may play important roles in determining splicing outcomes in a cell type- and intron-specific manner.