Search

Alu-elements comprise a large part of the human genome and some insertions have been shown to cause diseases. Here, we illuminate the protective role of an Alu-element in the 3'UTR of the human Factor 9 gene and its ability to ameliorate a poly(A) site mutation in a hemophilia B patient, preventing him from developing a severe disease. Using a minigene, we examined the disease-causing mutation and the modifying effect of the transposon in cellulo. Further, we simulated evolutionary scenarios regarding alternative polyadenylation before and after Alu insertion. A sequence analysis revealed that Old World monkeys displayed a highly conserved polyadenylation sites in this Alu-element, whereas New World monkeys lacked this motif, indicating a selective pressure. We conclude that this transposon has inserted shortly before the separation of Old and New World monkeys and thus also serves as a molecular landmark in primate evolution., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Kopp et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The translational repressor Nanos (Nos) regulates a single target, maternal hunchback (hb) mRNA, to govern abdominal segmentation in the early Drosophila embryo. Nos is recruited to sites in the 3' UTR of hb mRNA in collaboration with the sequence-specific RNA-binding protein Pumilio (Pum); on its own, Nos has no binding specificity. Nos is expressed at other stages of development, but very few mRNA targets that might mediate its action at these stages have been described. Nor has it been clear whether Nos is targeted to other mRNAs in concert with Pum or via other mechanisms. In this report, we identify mRNAs targeted by Nos via 2 approaches. First, we identify mRNAs depleted upon expression of a chimera bearing Nos fused to the nonsense mediated decay (NMD) factor Upf1. We find that, in addition to hb, Upf1-Nos depletes approximately 2,600 mRNAs from the maternal transcriptome in early embryos. Virtually all of these appear to be targeted in a canonical, hb-like manner in concert with Pum. In a second, more conventional approach, we identify mRNAs that are stabilized during the maternal zygotic transition (MZT) in embryos from nos- females. Most (86%) of the 1,185 mRNAs regulated by Nos are also targeted by Upf1-Nos, validating use of the chimera. Previous work has shown that 60% of the maternal transcriptome is degraded in early embryos. We find that maternal mRNAs targeted by Upf1-Nos are hypoadenylated and inefficiently translated at the ovary-embryo transition; they are subsequently degraded in the early embryo, accounting for 59% of all destabilized maternal mRNAs. We suggest that the late ovarian burst of Nos represses a large fraction of the maternal transcriptome, priming it for later degradation by other factors in the embryo., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Marhabaie et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

A recently reported Schizophrenia-associated genetic variant in the 3'UTR of the human furin gene, a homolog of C. elegans kpc-1, highlights an important role of the furin 3'UTR in neuronal development. We isolate three kpc-1 mutants that display abnormal dendrite arborization in PVD neurons and defective male mating behaviors. We show that the kpc-1 3'UTR participates in dendrite branching and self-avoidance. The kpc-1 3'UTR facilitates mRNA localization to branching points and contact points between sibling dendrites and promotes translation efficiency. A predicted secondary structural motif in the kpc-1 3'UTR is required for dendrite self-avoidance. Animals with over-expression of DMA-1, a PVD dendrite receptor, exhibit similar dendrite branching and self-avoidance defects that are suppressed with kpc-1 over-expression. Our results support a model in which KPC-1 proteins are synthesized at branching points and contact points to locally down-regulate DMA-1 receptors to promote dendrite branching and self-avoidance of a mechanosensory neuron important for male courtship., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Shih et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The eukaryotic mRNA life cycle includes transcription, nuclear mRNA export and degradation. To quantify all these processes simultaneously, we perform thiol-linked alkylation after metabolic labeling of RNA with 4-thiouridine (4sU), followed by sequencing of RNA (SLAM-seq) in the nuclear and cytosolic compartments of human cancer cells. We develop a model that reliably quantifies mRNA-specific synthesis, nuclear export, and nuclear and cytosolic degradation rates on a genome-wide scale. We find that nuclear degradation of polyadenylated mRNA is negligible and nuclear mRNA export is slow, while cytosolic mRNA degradation is comparatively fast. Consequently, an mRNA molecule generally spends most of its life in the nucleus. We also observe large differences in the nuclear export rates of different 3'UTR transcript isoforms. Furthermore, we identify genes whose expression is abruptly induced upon metabolic labeling. These transcripts are exported substantially faster than average mRNAs, suggesting the existence of alternative export pathways. Our results highlight nuclear mRNA export as a limiting factor in mRNA metabolism and gene regulation., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Müller et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Research has identified the large conductance voltage- and calcium-activated potassium channel (BK) as a key regulator of neuronal excitability genetically associated to behavioral alcohol tolerance. Sensitivity to ethanol at the molecular level is characterized by acute potentiation of channel activity. BK isoforms show variations in alcohol sensitivity and are differentially distributed on the plasma membrane surface in response to prolonged exposure. MicroRNA (MiRNA) targeting of alcohol-sensitive isoforms coupled with active internalization of BK channels in response to ethanol are believed to be key in establishing homeostatic adaptations that produce persistent changes within the plasma membrane of neurons. In fact, microRNA 9 (miR-9) upregulated expression is a key event in persistent alcohol tolerance mediating acute EtOH desensitization of BK channels. The exact nature of these interactions remains a current topic of discussion. To further study the effects of miR-9 on the expression and distribution of BK channel isoforms we designed an experimental model by transfecting human BK channel isoforms ZERO heterologous constructs in human embryonic kidney cells 293 (HEK293) cells respectively expressing 2.1 (miR-9 responsive), 2.2 (unresponsive) and control (no sequence) 3'untranslated region (3'UTR) miRNA recognition sites. We used imaging techniques to characterize the stably transfected monoclonal cell lines, and electrophysiology to validate channel activity. Finally, we used immunocytochemistry to validate isoform responsiveness to miR-9. Our findings suggest the cell lines were successfully transfected to express either the 2.1 or 2.2 version of ZERO. Patch clamp recordings confirm that these channels retain their functionality and immunohistochemistry shows differential responses to miR-9, making these cells viable for use in future alcohol dependence studies., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Cordero Padilla et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Alternative polyadenylation (APA) is a widespread mechanism of gene regulation that generates mRNA isoforms with alternative 3' untranslated regions (3' UTRs). Our previous study has revealed the global 3' UTR shortening of host mRNAs through APA upon viral infection. However, how the dynamic changes in the APA landscape occur upon viral infection remains largely unknown. Here we further found that, the reduced protein abundance of CPSF6, one of the core 3' processing factors, promotes the usage of proximal poly(A) sites (pPASs) of many immune related genes in macrophages and fibroblasts upon viral infection. Shortening of the 3' UTR of these transcripts may improve their mRNA stability and translation efficiency, leading to the promotion of type I IFN (IFN-I) signalling-based antiviral immune responses. In addition, dysregulated expression of CPSF6 is also observed in many immune related physiological and pathological conditions, especially in various infections and cancers. Thus, the global APA dynamics of immune genes regulated by CPSF6, can fine-tune the antiviral response as well as the responses to other cellular stresses to maintain the tissue homeostasis, which may represent a novel regulatory mechanism for antiviral immunity., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Ge et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Mammalian orthoreovirus (MRV) is a prototypic member of the Spinareoviridae family and has ten double-stranded RNA segments. One copy of each segment must be faithfully packaged into the mature virion, and prior literature suggests that nucleotides (nts) at the terminal ends of each gene likely facilitate their packaging. However, little is known about the precise packaging sequences required or how the packaging process is coordinated. Using a novel approach, we have determined that 200 nts at each terminus, inclusive of untranslated regions (UTR) and parts of the open reading frame (ORF), are sufficient for packaging S gene segments (S1-S4) individually and together into replicating virus. Further, we mapped the minimal sequences required for packaging the S1 gene segment into a replicating virus to 25 5' nts and 50 3' nts. The S1 UTRs, while not sufficient, were necessary for efficient packaging, as mutations of the 5' or 3' UTRs led to a complete loss of virus recovery. Using a second novel assay, we determined that 50 5' nts and 50 3' nts of S1 are sufficient to package a non-viral gene segment into MRV. The 5' and 3' termini of the S1 gene are predicted to form a panhandle structure and specific mutations within the stem of the predicted panhandle region led to a significant decrease in viral recovery. Additionally, mutation of six nts that are conserved across the three major serotypes of MRV that are predicted to form an unpaired loop in the S1 3' UTR, led to a complete loss of viral recovery. Overall, our data provide strong experimental proof that MRV packaging signals lie at the terminal ends of the S gene segments and offer support that the sequence requirements for efficient packaging of the S1 segment include a predicted panhandle structure and specific sequences within an unpaired loop in the 3' UTR., Competing Interests: The authors have declared no competing interests exist., (Copyright: © 2024 Dhar et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Localization of oskar mRNA to the posterior of the Drosophila oocyte is essential for abdominal patterning and germline development. oskar localization is a multi-step process involving temporally and mechanistically distinct transport modes. Numerous cis-acting elements and trans-acting factors have been identified that mediate earlier motor-dependent transport steps leading to an initial accumulation of oskar at the posterior. Little is known, however, about the requirements for the later localization phase, which depends on cytoplasmic flows and results in the accumulation of large oskar ribonucleoprotein granules, called founder granules, by the end of oogenesis. Using super-resolution microscopy, we show that founder granules are agglomerates of smaller oskar transport particles. In contrast to the earlier kinesin-dependent oskar transport, late-phase localization depends on the sequence as well as on the structure of the spliced oskar localization element (SOLE), but not on the adjacent exon junction complex deposition. Late-phase localization also requires the oskar 3' untranslated region (3' UTR), which targets oskar to founder granules. Together, our results show that 3' UTR-mediated targeting together with SOLE-dependent agglomeration leads to accumulation of oskar in large founder granules at the posterior of the oocyte during late stages of oogenesis. In light of previous work showing that oskar transport particles are solid-like condensates, our findings indicate that founder granules form by a process distinct from that of well-characterized ribonucleoprotein granules like germ granules, P bodies, and stress granules. Additionally, they illustrate how an individual mRNA can be adapted to exploit different localization mechanisms depending on the cellular context., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Eichler et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Introduction: Altered circadian rhythms underlie manifestation of several cardiovascular disorders, however a little is known about the mediating biomolecules. Multiple transcriptional-translational feedback loops control circadian-clockwork wherein; micro RNAs (miRNAs) are known to manifest post transcriptional regulation. This study assesses miR34a-5p as a mediating biomolecule., Method: 8-10-week-old male C57BL/6J mice (n = 6/group) were subjected to photoperiodic manipulation induced chronodisruption and thoracic aortae were examined for miRNA, gene (qPCR) and protein (Immunoblot) expression studies. Histomorphological changes were assessed for pro-atherogenic manifestations (fibrillar arrangement, collagen/elastin ratio, intima-media thickening). Computational studies for miRNA-mRNA target prediction were done using TargetScan and miRDB. Correlative in vitro studies were done in serum synchronized HUVEC cells. Time point based studies were done at five time points (ZT 0, 6, 12, 18, 24) in 24h., Results: Chronodisruption induced hypomethylation in the promoter region of miR34a-5p, in the thoracic aortae, culminating in elevated miRNA titers. In a software-based detection of circadian-clock-associated targets of miR34a-5p, Clock and Sirt1 genes were identified. Moreover, miR34a-5p exhibited antagonist circadian oscillations to that of its target genes CLOCK and SIRT1 in endothelial cells. Luciferase reporter gene assay further showed that miR34a-5p interacts with the 3'UTR of the Clock gene to lower its expression, disturbing the operation of positive arm of circadian clock system. Elevated miR34a-5p and impeded SIRT1 expression in a chronodisruptive aortae exhibited pro-atherogenic changes observed in form of gene expression, increased collagen/elastin ratio, fibrillar derangement and intimal-media thickening., Conclusion: The study reports for the first time chronodisruption mediated miR34a-5p elevation, its circadian expression and interaction with the 3'UTR of Clock gene to impede its expression. Moreover, elevated miR34a-5p and lowered SIRT1 expression in the chronodisruptive aortae lead off cause-consequence relationship of chronodisruption mediated proatherogenic changes., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Vyas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Ets-related gene (ERG) is overexpressed as a fusion protein in prostate cancer. During metastasis, the pathological role of ERG is associated with cell proliferation, invasion, and angiogenesis. Here, we hypothesized that miRNAs regulate ERG expression through its 3'UTR. Several bioinformatics tools were used to identify miRNAs and their binding sites on 3'UTR of ERG. The selected miRNAs expression was analyzed in prostate cancer samples by qPCR. The miRNAs overexpression was induced in prostate cancer cells (VCaP) to analyze ERG expression. Reporter gene assay was performed to evaluate the ERG activity in response to selected miRNAs. The expression of ERG downstream target genes was also investigated through qPCR after miRNAs overexpression. To observe the effects of selected miRNAs on cell proliferation and migration, scratch assay was performed to calculate the cell migration rate. miR-4482 and miR-3912 were selected from bioinformatics databases. miR-4482 and -3912 expression were decreased in prostate cancer samples, as compared to controls (p<0.05 and p<0.001), respectively. Overexpression of miR-4482 and miR-3912 significantly reduced ERG mRNA (p<0.001 and p<0.01), respectively) and protein (p<0.01) in prostate cancer cells. The transcriptional activity of ERG was significantly reduced (p<0.01) in response to miR-4482 and-3912. ERG angiogenic targets and cell migration rate was also reduced significantly (p<0.001) after miR-4482 and -3912 over-expression. This study indicates that miR-4482 and -3912 can suppress the ERG expression and its target genes, thereby, halt prostate cancer progression. These miRNAs may be employed as a potential therapeutic target for the miRNA-based therapy against prostate cancer., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Mumtaz et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

RNA viral genomes compact information into functional RNA structures. Here, using chikungunya virus as a model, we investigated the structural requirements of conserved RNA elements in the 3' untranslated region (3'UTR) for viral replication in mosquito and mammalian cells. Using structural predictions and co-variation analysis, we identified a highly stable and conserved Y-shaped structure (SLY) at the end of the 3'UTR that is duplicated in the Asian lineage. Functional studies with mutant viruses showed that the SLY has host-specific functions during viral replication and evolution. The SLY positively modulates viral replication in mosquito cells but has the opposite effect in mammalian cells. Additional structural/functional analyses showed that maintaining the Y-shaped fold and specific nucleotides in the loop are critical for full SLY functionality and optimal viral replication in mosquito cells. Experimental adaptation of viruses with duplicated SLYs to mammalian cells resulted in the generation of heterogeneous viral populations comprising variants with diverse 3'UTRs, contrasting with the homogeneous populations from viruses without SLY copies. Altogether, our findings constitute the first evidence of an RNA secondary structure in the 3'UTR of chikungunya virus genome that plays host-dependent functions., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Bardossy et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Aim: In a retrospective and exploratory case-control study, we examined the genetic association of two common polymorphisms in the 3' untranslated region (UTR) of DICER1 (rs3742330) and DROSHA (rs10719) genes in primary open-angle glaucoma (POAG) and primary angle-closure glaucoma (PACG), and its related clinical phenotypes in a Saudi cohort., Methods: DNA genotyping was performed using TaqMan real-time PCR assays in 500 participants, including 152 POAG, 102 PACG, and 246 non-glaucomatous controls. Statistical analyses were performed to examine the association(s)., Results: Allele and genotype frequency of rs3742330 and rs10719 did not vary significantly in POAG and PACG compared to controls. No significant deviation was observed from Hardy-Weinberg Equilibrium (p > 0.05). Gender stratification revealed no significant allelic/genotype association with glaucoma types. Also, these polymorphisms showed no significant genotype effect on clinical markers such as intraocular pressure, cup/disc ratio, and the number of antiglaucoma medications. Logistic regression showed no effect of age, sex, rs3742330, and rs10719 genotypes on the risk of disease outcome. We also examined a combined allelic effect of rs3742330 (A>G) and rs10719 (A>G). However, none of the allelic combinations significantly affected POAG and PACG., Conclusions: The 3' UTR polymorphisms rs3742330 and rs10719 of DICER1 and DROSHA genes are not associated with POAG and PACG or its related glaucoma indices in this Middle-Eastern cohort of Saudi Arab ethnicity. However, there is a need to validate the results on a broader population and other ethnicities., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Kondkar et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

SNAP-25 protein is a key protein of the SNARE complex that is involved in synaptic vesicles fusion with plasma membranes and neurotransmitter release, playing a fundamental role in neural plasticity. Recently the concentration of three specific miRNAs-miR-27b-3p, miR-181a-5p and miR-23a-3p -was found to be associated with a specific SNAP-25 polymorphism (rs363050). in silico analysis showed that all the three miRNAs target SNAP-25, but the effect of the interaction between these miRNAs and the 3'UTR of SNAP-25 mRNA is currently unknown. For this reason, we verified in vitro whether miR-27b-3p, miR-181a-5p and miR-23a-3p modulate SNAP-25 gene and protein expression. Initial experiments using miRNAs-co-transfected Vero cells and SNAP-25 3'UTR luciferase reporter plasmids showed that miR-181a-5p (p≤0.01) and miR-23a-3p (p<0.05), but not miR-27b-3p, modulate the luciferase signal, indicating that these two miRNAs bind the SNAP-25 3'UTR. Results obtained using human oligodendroglial cell line (MO3.13) transfected with miR-181a-5p or miR-27b-3p confirmed that miR-181a-5p and miR-23a-3p regulate SNAP-25 gene and protein expression. Interestingly, the two miRNAs modulate in an opposite way SNAP-25, as miR-181a-5p significantly increases (p<0.0005), whereas miR-23a-3p decreases (p<0.0005) its expression. These results for the first time describe the ability of miR-181a-5p and miR-23a-3p to modulate SNAP-25 expression, suggesting their possible use as biomarkers or as therapeutical targets for diseases in which SNAP-25 expression is altered., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Agostini et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

T cells are key players in our defence against infections and malignancies. When T cells differentiate or become activated, they undergo substantial alterations in gene expression. Even though RNA expression levels are now well documented throughout different stages of T cells, it is not well understood how mRNA expression translates into the protein landscape. By combining paired RNA sequencing and mass spectrometry data of primary human CD8+ T cells, we report that mRNA expression is a poor proxy for the overall protein output, irrespective of the differentiation or activation status. Yet, gene class stratification revealed a function-specific correlation of mRNA with protein expression. This gene class-specific expression pattern associated with differences in gene characteristics such as sequence conservation and untranslated region (UTR) lengths. In addition, the presence of AU-rich elements in the 3'UTR associated with alterations in mRNA and protein abundance T cell activation dependent, gene class-specific manner. In conclusion, our study highlights the role of gene characteristics as a determinant for gene expression in T cells., Competing Interests: The authors have declared that no competing interests exist.

Parkinson's disease is the second most common age-related, neurodegenerative disease. A small collection of genes has been linked to Parkinson's disease including LRRK2, SAT1, and SNCA, the latter of which encodes the protein alpha-synuclein that aggregates in Lewy bodies as a hallmark of the disease. Overexpression of even wild-type versions of these genes can lead to pathogenesis, yet the regulatory mechanisms that control protein production of the genes are not fully understood. Pumilio proteins belong to the highly conserved PUF family of eukaryotic RNA-binding proteins that post-transcriptionally regulate gene expression through binding conserved motifs in the 3' untranslated region (UTR) of mRNA targets known as PUF Recognition Elements (PREs). The 3'UTRs of LRRK2, SNCA and SAT1 each contain multiple putative PREs. Knockdown (KD) of the two human Pumilio homologs (Pumilio 1 and Pumilio 2) in a neurodegenerative model cell line, SH-SY5Y, resulted in increased SNCA and LRRK2 mRNA, as well as alpha-synuclein levels, suggesting these genes are normally repressed by the Pumilio proteins. Some studies have indicated a relationship between Pumilio and microRNA activities on the same target, especially when their binding sites are close together. LRRK2, SNCA, and SAT1 each contain several putative microRNA-binding sites within the 3'UTR, some of which reside near PREs. Small RNA-seq and microRNA qPCR assays were performed in both wild type and Pumilio KD SH-SY5Y cells to analyze global and differential microRNA expression. One thousand four hundred and four microRNAs were detected across wild type and Pumilio KD cells. Twenty-one microRNAs were differentially expressed between treatments, six of which were previously established to be altered in Parkinson's disease patient samples or research models. Expression of ten miRs predicted to target LRRK2 and SNCA was verified by RT-qPCR. Collectively, our results demonstrate that Pumilios and microRNAs play a multi-faceted role in regulating Parkinson's disease-associated genes., Competing Interests: The authors have declared that no competing interests exist.

Adenosine bases of RNA can be transiently modified by the deposition of a methyl-group to form N6-methyladenosine (m6A). This adenosine-methylation is an ancient process and the enzymes involved are evolutionary highly conserved. A genetic screen designed to identify suppressors of late flowering transgenic Arabidopsis plants overexpressing the miP1a microProtein yielded a new allele of the FIONA1 (FIO1) m6A-methyltransferase. To characterize the early flowering phenotype of fio1 mutant plants we employed an integrative approach of mRNA-seq, Nanopore direct RNA-sequencing and meRIP-seq to identify differentially expressed transcripts as well as differentially methylated RNAs. We provide evidence that FIO1 is the elusive methyltransferase responsible for the 3'-end methylation of the FLOWERING LOCUS C (FLC) transcript. Furthermore, our genetic and biochemical data suggest that 3'-methylation stabilizes FLC mRNAs and non-methylated FLC is a target for rapid degradation., Competing Interests: The authors have declared that no competing interests exist.

Dengue viruses (DENV) are expanding global pathogens that are transmitted through the bite of mosquitoes, mostly Aedes aegypti. As RNA viruses, DENV rely on RNA-binding proteins (RBPs) to complete their life cycle. Alternatively, RBPs can act as restriction factors that prevent DENV multiplication. While the importance of RBPs is well-supported in humans, there is a dearth of information about their influence on DENV transmission by mosquitoes. Such knowledge could be harnessed to design novel, effective interventions against DENV. Here, we successfully adapted RNA-affinity chromatography coupled with mass spectrometry-a technique initially developed in mammalian cells-to identify RBPs in Ae. aegypti cells. We identified fourteen RBPs interacting with DENV serotype 2 3'UTR, which is involved in the viral multiplication and produces subgenomic flaviviral RNA (sfRNA). We validated the RNA affinity results for two RBPs by confirming that AePur binds the 3'UTR, whereas AeStaufen interacts with both 3'UTR and sfRNA. Using in vivo functional evaluation, we determined that RBPs like AeRan, AeExoRNase, and AeRNase have pro-viral functions, whereas AeGTPase, AeAtu, and AePur have anti-viral functions in mosquitoes. Furthermore, we showed that human and mosquito Pur homologs have a shared affinity to DENV2 RNA, although the anti-viral effect is specific to the mosquito protein. Importantly, we revealed that AeStaufen mediates a reduction of gRNA and sfRNA copies in several mosquito tissues, including the salivary glands and that AeStaufen-mediated sfRNA reduction diminishes the concentration of transmission-enhancing sfRNA in saliva, thereby revealing AeStaufen's role in DENV transmission. By characterizing the first RBPs that associate with DENV2 3'UTR in mosquitoes, our study unravels new pro- and anti-viral targets for the design of novel therapeutic interventions as well as provides foundation for studying the role of RBPs in virus-vector interactions., Competing Interests: The authors have declared that no competing interests exist.

Selenoproteins contain the 21st amino acid, selenocysteine (Sec), which is incorporated at select UGA codons when a specialized hairpin sequence, the Sec insertion sequence (SECIS) element, is present in the 3' UTR. Aside from the SECIS, selenoprotein mRNA 3' UTRs are not conserved between different selenoproteins within a species. In contrast, the 3'-UTR of a given selenoprotein is often conserved across species, which supports the hypothesis that cis-acting elements in the 3'-UTR other than the SECIS exert post-transcriptional control on selenoprotein expression. In order to determine the function of one such SECIS context, we chose to focus on the plasma selenoprotein, SELENOP, which is required to maintain selenium homeostasis as a selenium transport protein that contains 10 Sec residues. It is unique in that its mRNA contains two SECIS elements in the context of a highly conserved 843-nucleotide 3' UTR. Here we have used RNA affinity chromatography and identified PTBP1 as the major RNA binding protein that specifically interacts with the sequence between the two SECIS elements. We then used CRISPR/Cas9 genome editing to delete two regions surrounding the first SECIS element. We found that these sequences are involved in regulating SELENOP mRNA and protein levels, which are inversely altered as a function of selenium concentrations., Competing Interests: The authors have declared that no competing interests exist.

Feed efficiency (FE) is an important trait in the broiler industry due to its direct correlation to efficient muscle growth instead of fat deposition. The present study characterized and compared gene expression profiles in abdominal fat from broiler chickens of different FE levels to enhance the understanding of FE biology. Specifically, traditional whole-transcript RNA-sequencing (RNA-seq) and 3' UTR-sequencing (3' UTR-seq) were applied to 22 and 61 samples, respectively. Overall, these two sequencing techniques shared a high correlation (0.76) between normalized counts, although 3' UTR-seq showed a higher variance in sequencing and mapping performance statistics across samples and a lower rate of uniquely mapped reads. A higher percentage of 3' UTR-seq reads mapped to introns suggested the frequent presence of cleavage sites in introns, thus warranting future research to study its regulatory function. Differential expression analysis identified 1198 differentially expressed genes (DEGs) between high FE (HFE) and intermediate FE (IFE) chickens with False Discovery Rate < 0.05 and fold change > 1.2. The processes that were significantly enriched by the DEGs included extracellular matrix remodeling and mechanisms impacting gene expression at the transcriptional and translational levels. Gene ontology enrichment analysis suggested that the divergence in fat deposition and FE in broiler chickens could be associated with peroxisome and lipid metabolism possibly regulated by G0/G1 switch gene 2 (G0S2)., Competing Interests: The authors have declared that no competing interests exist.

Gene activity is controlled by multiple molecular mechanisms, for instance through transcription factors or by microRNAs (miRNAs), among others. Established bioinformatics tools for the prediction of miRNA target genes face the challenge of ensuring accuracy, due to high false positive rates. Further, these tools present poor overlap. However, we demonstrated that it is possible to filter good predictions of miRNA targets from the bulk of all predictions by using information from the gene regulatory network. Here, we take advantage of this strategy that selects a large subset of predicted microRNA binding sites as more likely to possess less false-positives because of their over-representation in RE1 silencing transcription factor (REST)-regulated genes from the background of TargetScanHuman 7.2 predictions to identify useful features for the prediction of microRNA targets. These enriched miRNA families would have silencing activity for neural transcripts overlapping the repressive activity on neural genes of REST. We analyze properties of associated microRNA binding sites and contrast the outcome to the background. We found that the selected subset presents significant differences respect to the background: (i) lower GC-content in the vicinity of the predicted miRNA binding site, (ii) more target genes with multiple identical microRNA binding sites and (iii) a higher density of predicted microRNA binding sites close to the 3' terminal end of the 3'-UTR. These results suggest that network selection of miRNA-mRNA pairs could provide useful features to improve microRNA target prediction., Competing Interests: The authors have declared that no competing interests exist.

RNA stability control contributes to the proper expression of gene products. Messenger RNAs (mRNAs) in eukaryotic cells possess a 5' cap structure and the 3' poly(A) tail which are important for mRNA stability and efficient translation. The Ccr4-Not complex is a major cytoplasmic deadenylase and functions in mRNA degradation. The CLB1-6 genes in Saccharomyces cerevisiae encode B-type cyclins which are involved in the cell cycle progression together with the cyclin-dependent kinase Cdc28. The CLB genes consist of CLB1/2, CLB3/4, and CLB5/6 whose gene products accumulate at the G2-M, S-G2, and late G1 phase, respectively. These Clb protein levels are thought to be mainly regulated by the transcriptional control and the protein stability control. Here we investigated regulation of CLB1-6 expression by Ccr4. Our results show that all CLB1-6 mRNA levels were significantly increased in the ccr4Δ mutant compared to those in wild-type cells. Clb1, Clb4, and Clb6 protein levels were slightly increased in the ccr4Δ mutant, but the Clb2, Clb3, and Clb5 protein levels were similar to those in wild-type cells. Since both CLB6 mRNA and Clb6 protein levels were most significantly increased in the ccr4Δ mutant, we further analyzed the cis-elements for the Ccr4-mediated regulation within CLB6 mRNA. We found that there were destabilizing sequences in both coding sequence and 3' untranslated region (3' UTR). The destabilizing sequences in the coding region were found to be both within and outside the sequences corresponding the cyclin domain. The CLB6 3' UTR was sufficient for mRNA destabilization and decrease of the reporter GFP gene and this destabilization involved Ccr4. Our results suggest that CLB6 expression is regulated by Ccr4 through the coding sequence and 3' UTR of CLB6 mRNA., Competing Interests: The authors have declared that no competing interests exist.

Within the immune system, microRNAs (miRNAs) exert key regulatory functions. However, what are the mRNA targets regulated by miRNAs and how miRNAs are transcriptionally regulated themselves remain for the most part unknown. We found that in primary human memory T helper lymphocytes, miR-150 was the most abundantly expressed miRNA, and its expression decreased drastically upon activation, suggesting regulatory roles. Constitutive MIR150 gene expression required the RFX family of transcription factors, and its activation-induced down-regulation was linked to their reduced expression. By performing miRNA pull-down and sequencing experiments, we identified PDGFA-associated protein 1 (PDAP1) as one main target of miR-150 in human T lymphocytes. PDAP1 acted as an RNA-binding protein (RBP), and its CRISPR/Cas-9-mediated deletion revealed that it prominently contributed to the regulation of T-cell proliferation. Overall, using an integrated approach involving quantitative analysis, unbiased genomics, and genome editing, we identified RFX factors, miR-150, and the PDAP1 RBP as the components of a regulatory axis that restrains proliferation of primary human T lymphocytes., Competing Interests: The authors have declared that no competing interests exist.

Translation of messenger RNAs (mRNAs) with premature termination codons produces truncated proteins with potentially deleterious effects. This is prevented by nonsense-mediated mRNA decay (NMD) of these mRNAs. NMD is triggered by ribosomes terminating upstream of a splice site marked by an exon-junction complex (EJC), but also acts on many mRNAs lacking a splice junction after their termination codon. We developed a genome-wide CRISPR flow cytometry screen to identify regulators of mRNAs with premature termination codons in K562 cells. This screen recovered essentially all core NMD factors and suggested a role for EJC factors in degradation of PTCs without downstream splicing. Among the strongest hits were the translational repressors GIGYF2 and EIF4E2. GIGYF2 and EIF4E2 mediate translational repression but not mRNA decay of a subset of NMD targets and interact with NMD factors genetically and physically. Our results suggest a model wherein recognition of a stop codon as premature can lead to its translational repression through GIGYF2 and EIF4E2., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: RG is a shareholder and scientific advisory board member of Moderna Inc. and a consultant for FL63, Celgene, Vertex, Arrakis and Initial Therapeutics.

Papaya ringspot virus biotype-P is a detrimental pathogen of economically important papaya and cucurbits worldwide. The mutation prone feature of this virus perhaps accounts for its geographical dissemination. In this study, investigations of the atypical PRSV-P strain was conducted based on phylogenetic, recombination and genetic differentiation analyses considering of it's likely spread across India and Bangladesh. Full length genomic sequences of 38 PRSV isolates and 35 CP gene sequences were subjected to recombination analysis. A total of 61 recombination events were detected in aligned complete PRSV genome sequences. 3 events were detected in complete genome of PRSV strain PK whereas one was in its CP gene sequence. The PRSV-PK appeared to be recombinant of a major parent from Bangladesh. However, the genetic differentiation based on full length genomic sequences revealed less frequent gene flow between virus PRSV-PK and the population from America, India, Colombia, other Asian Countries and Australia. Whereas, frequent gene flow exists between Pakistan and Bangladesh virus populations. These results provided evidence correlating geographical position and genetic distances. We speculate that the genetic variations and evolutionary dynamics of this virus may challenge the resistance developed in papaya against PRSV and give rise to virus lineage because of its atypical emergence where geographic spread is already occurring., Competing Interests: The authors have declared that no competing interests exist.

RNA regulation is essential to successful reproduction. Messenger RNAs delivered from parent to progeny govern early embryonic development. RNA-binding proteins (RBPs) are the key effectors of this process, regulating the translation and stability of parental transcripts to control cell fate specification events prior to zygotic gene activation. The KH-domain RBP MEX-3 is conserved from nematode to human. It was first discovered in Caenorhabditis elegans, where it is essential for anterior cell fate and embryo viability. Here, we show that loss of the endogenous mex-3 3´UTR disrupts its germline expression pattern. An allelic series of 3´UTR deletion variants identify repressing regions of the UTR and demonstrate that repression is not precisely coupled to reproductive success. We also show that several RBPs regulate mex-3 mRNA through its 3´UTR to define its unique germline spatiotemporal expression pattern. Additionally, we find that both poly(A) tail length control and the translation initiation factor IFE-3 contribute to its expression pattern. Together, our results establish the importance of the mex-3 3´UTR to reproductive health and its expression in the germline. Our results suggest that additional mechanisms control MEX-3 function when 3´UTR regulation is compromised., Competing Interests: The authors have declared that no competing interests exist.

Adenosine deaminase acting on RNA 1 (ADAR1), an enzyme responsible for adenosine-to-inosine RNA editing, is composed of two isoforms: nuclear p110 and cytoplasmic p150. Deletion of Adar1 or Adar1 p150 genes in mice results in embryonic lethality with overexpression of interferon-stimulating genes (ISGs), caused by the aberrant recognition of unedited endogenous transcripts by melanoma differentiation-associated protein 5 (MDA5). However, among numerous RNA editing sites, how many RNA sites require editing, especially by ADAR1 p150, to avoid MDA5 activation and whether ADAR1 p110 contributes to this function remains elusive. In particular, ADAR1 p110 is abundant in the mouse brain where a subtle amount of ADAR1 p150 is expressed, whereas ADAR1 mutations cause Aicardi-Goutières syndrome, in which the brain is one of the most affected organs accompanied by the elevated expression of ISGs. Therefore, understanding RNA editing-mediated prevention of MDA5 activation in the brain is especially important. Here, we established Adar1 p110-specific knockout mice, in which the upregulated expression of ISGs was not observed. This result suggests that ADAR1 p150-mediated RNA editing is enough to suppress MDA5 activation. Therefore, we further created Adar1 p110/Adar2 double knockout mice to identify ADAR1 p150-mediated editing sites. This analysis demonstrated that although the elevated expression of ISGs was not observed, only less than 2% of editing sites were preserved in the brains of Adar1 p110/Adar2 double knockout mice. Of note, we found that some sites were highly edited, which was comparable to those found in wild-type mice, indicating the presence of ADAR1 p150-specific sites. These data suggest that RNA editing at a very limited sites, which is mediated by a subtle amount of ADAR1 p150, is sufficient to prevents MDA5 activation, at least in the mouse brain., Competing Interests: The authors have declared that no competing interests exist.

Tomato ringspot virus (ToRSV, genus Nepovirus, family Secoviridae, order Picornavirales) is a bipartite positive-strand RNA virus, with each RNA encoding one large polyprotein. ToRSV RNAs are linked to a 5'-viral genome-linked protein (VPg) and have a 3' polyA tail, suggesting a non-canonical cap-independent translation initiation mechanism. The 3' untranslated regions (UTRs) of RNA1 and RNA2 are unusually long (~1.5 kb) and share several large stretches of sequence identities. Several putative in-frame start codons are present in the 5' regions of the viral RNAs, which are also highly conserved between the two RNAs. Using reporter transcripts containing the 5' region and 3' UTR of the RNA2 of ToRSV Rasp1 isolate (ToRSV-Rasp1) and in vitro wheat germ extract translation assays, we provide evidence that translation initiates exclusively at the first AUG, in spite of a poor codon context. We also show that both the 5' region and 3' UTR of RNA2 are required for efficient cap-independent translation of these transcripts. We identify translation-enhancing elements in the 5' proximal coding region of the RNA2 polyprotein and in the RNA2 3' UTR. Cap-dependent translation of control reporter transcripts was inhibited when RNAs consisting of the RNA2 3' UTR were supplied in trans. Taken together, our results suggest the presence of a CITE in the ToRSV-Rasp1 RNA2 3' UTR that recruits one or several translation factors and facilitates efficient cap-independent translation together with the 5' region of the RNA. Non-overlapping deletion mutagenesis delineated the putative CITE to a 200 nts segment (nts 773-972) of the 1547 nt long 3' UTR. We conclude that the general mechanism of ToRSV RNA2 translation initiation is similar to that previously reported for the RNAs of blackcurrant reversion virus, another nepovirus. However, the position, sequence and predicted structures of the translation-enhancing elements differed between the two viruses., Competing Interests: The authors have declared that no competing interests exist.

Once loaded onto Argonaute proteins, microRNAs form a silencing complex called miRISC that targets mostly the 3'UTR of mRNAs to silence their translation. How microRNAs are transported to and from their target mRNA remains poorly characterized. While some reports linked intracellular trafficking to microRNA activity, it is still unclear how these pathways coordinate for proper microRNA-mediated gene silencing and turnover. Through a forward genetic screen using Caenorhabditis elegans, we identified the RabGAP tbc-11 as an important factor for the microRNA pathway. We show that TBC-11 acts mainly through the small GTPase RAB-6 and that its regulation is required for microRNA function. The absence of functional TBC-11 increases the pool of microRNA-unloaded Argonaute ALG-1 that is likely associated to endomembranes. Furthermore, in this condition, this pool of Argonaute accumulates in a perinuclear region and forms a high molecular weight complex. Altogether, our data suggest that the alteration of TBC-11 generates a fraction of ALG-1 that cannot bind to target mRNAs, leading to defective gene repression. Our results establish the importance of intracellular trafficking for microRNA function and demonstrate the involvement of a small GTPase and its GAP in proper Argonaute localization in vivo., Competing Interests: The authors have declared that no competing interests exist.

ELAV/Hu factors are conserved RNA binding proteins (RBPs) that play diverse roles in mRNA processing and regulation. The founding member, Drosophila Elav, was recognized as a vital neural factor 35 years ago. Nevertheless, little was known about its impacts on the transcriptome, and potential functional overlap with its paralogs. Building on our recent findings that neural-specific lengthened 3' UTR isoforms are co-determined by ELAV/Hu factors, we address their impacts on splicing. While only a few splicing targets of Drosophila are known, ectopic expression of each of the three family members (Elav, Fne and Rbp9) alters hundreds of cassette exon and alternative last exon (ALE) splicing choices. Reciprocally, double mutants of elav/fne, but not elav alone, exhibit opposite effects on both classes of regulated mRNA processing events in larval CNS. While manipulation of Drosophila ELAV/Hu RBPs induces both exon skipping and inclusion, characteristic ELAV/Hu motifs are enriched only within introns flanking exons that are suppressed by ELAV/Hu factors. Moreover, the roles of ELAV/Hu factors in global promotion of distal ALE splicing are mechanistically linked to terminal 3' UTR extensions in neurons, since both processes involve bypass of proximal polyadenylation signals linked to ELAV/Hu motifs downstream of cleavage sites. We corroborate the direct action of Elav in diverse modes of mRNA processing using RRM-dependent Elav-CLIP data from S2 cells. Finally, we provide evidence for conservation in mammalian neurons, which undergo broad programs of distal ALE and APA lengthening, linked to ELAV/Hu motifs downstream of regulated polyadenylation sites. Overall, ELAV/Hu RBPs orchestrate multiple broad programs of neuronal mRNA processing and isoform diversification in Drosophila and mammalian neurons., Competing Interests: The authors have declared that no competing interests exist.

Faithful replication of the entire genome requires replication forks to copy large contiguous tracts of DNA, and sites of persistent replication fork stalling present a major threat to genome stability. Understanding the distribution of sites at which replication forks stall, and the ensuing fork processing events, requires genome-wide methods that profile replication fork position and the formation of recombinogenic DNA ends. Here, we describe Transferase-Activated End Ligation sequencing (TrAEL-seq), a method that captures single-stranded DNA 3' ends genome-wide and with base pair resolution. TrAEL-seq labels both DNA breaks and replication forks, providing genome-wide maps of replication fork progression and fork stalling sites in yeast and mammalian cells. Replication maps are similar to those obtained by Okazaki fragment sequencing; however, TrAEL-seq is performed on asynchronous populations of wild-type cells without incorporation of labels, cell sorting, or biochemical purification of replication intermediates, rendering TrAEL-seq far simpler and more widely applicable than existing replication fork direction profiling methods. The specificity of TrAEL-seq for DNA 3' ends also allows accurate detection of double-strand break sites after the initiation of DNA end resection, which we demonstrate by genome-wide mapping of meiotic double-strand break hotspots in a dmc1Δ mutant that is competent for end resection but not strand invasion. Overall, TrAEL-seq provides a flexible and robust methodology with high sensitivity and resolution for studying DNA replication and repair, which will be of significant use in determining mechanisms of genome instability., Competing Interests: The authors have declared that no competing interests exist.

Mutations in RanBP2 (also known as Nup358), one of the main components of the cytoplasmic filaments of the nuclear pore complex, contribute to the overproduction of acute necrotizing encephalopathy (ANE1)-associated cytokines. Here we report that RanBP2 represses the translation of the interleukin 6 (IL6) mRNA, which encodes a cytokine that is aberrantly up-regulated in ANE1. Our data indicates that soon after its production, the IL6 messenger ribonucleoprotein (mRNP) recruits Argonautes bound to let-7 microRNA. After this mRNP is exported to the cytosol, RanBP2 sumoylates mRNP-associated Argonautes, thereby stabilizing them and enforcing mRNA silencing. Collectively, these results support a model whereby RanBP2 promotes an mRNP remodelling event that is critical for the miRNA-mediated suppression of clinically relevant mRNAs, such as IL6., Competing Interests: The authors have declared that no competing interests exist.

Background: Brain metastasis (BM) is a major cause of morbidity and mortality in breast cancer (BC) and its molecular mechanism remains poorly understood. Transmigration of metastatic cells through the brain endothelium is an essential step in BM. Metalloproteinase-1 (MMP-1) overexpression plays a key role in promoting trans-endothelial migration by degrading the inter-endothelial junctions and disrupting the endothelial integrity. However, little is known about the molecular mechanisms that induce MMP-1 in metastatic cells granting them a brain invasive phenotype. MiR-202-3p is downregulated in brain metastases compared to primary breast tumors and directly targets MMP-1. Here, we unraveled a critical role of miR-202-3p loss in MMP-1 upregulation promoting transmigration of metastatic cells through the brain endothelium., Methods: A variant of the MDA-MB-231 human BC cell line (MDA-MB-231-BrM2) selected for its propensity to form brain metastases was found to express high levels of MMP-1 and low levels of miR-202-3p compared to the parental cells. Using a gain-and-loss of function approach, we modulated levels of miR-202-3p and examined the resultant effect on MMP-1 expression. Effect of miR-202-3p modulation on integrity of the brain endothelium and the transmigrative ability of BC cells were also examined., Results: Loss of miR-202-3p in breast cancer cells enhanced their transmigration through the brain endothelium by upregulating MMP-1 and disrupting the inter-endothelial junctions (claudin-5, ZO-1 and ß-catenin). Restoring miR-202-3p exerted a metastasis-suppressive effect and preserved the endothelial barrier integrity., Conclusions: Our study identified a critical regulatory role of miR-202-3p in brain metastasis and shed light on miR-202-3p/MMP-1 axis as a novel prognostic and therapeutic target that can be exploited to predict and prevent brain metastasis in breast cancer patients., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Background: The world is going through the critical phase of COVID-19 pandemic, caused by human coronavirus, SARS-CoV-2. Worldwide concerted effort to identify viral genomic changes across different sub-types has identified several strong changes in the coding region. However, there have not been many studies focusing on the variations in the 5' and 3' untranslated regions and their consequences. Considering the possible importance of these regions in host mediated regulation of viral RNA genome, we wanted to explore the phenomenon., Methods: To have an idea of the global changes in 5' and 3'-UTR sequences, we downloaded 8595 complete and high-coverage SARS-CoV-2 genome sequence information from human host in FASTA format from Global Initiative on Sharing All Influenza Data (GISAID) from 15 different geographical regions. Next, we aligned them using Clustal Omega software and investigated the UTR variants. We also looked at the putative host RNA binding protein (RBP) and microRNA binding sites in these regions by 'RBPmap' and 'RNA22 v2' respectively. Expression status of selected RBPs and microRNAs were checked in lungs tissue., Results: We identified 28 unique variants in SARS-CoV-2 UTR region based on a minimum variant percentage cut-off of 0.5. Along with 241C>T change the important 5'-UTR change identified was 187A>G, while 29734G>C, 29742G>A/T and 29774C>T were the most familiar variants of 3'UTR among most of the continents. Furthermore, we found that despite the variations in the UTR regions, binding of host RBP to them remains mostly unaltered, which further influenced the functioning of specific miRNAs., Conclusion: Our results, shows for the first time in SARS-Cov-2 infection, a possible cross-talk between host RBPs-miRNAs and viral UTR variants, which ultimately could explain the mechanism of escaping host RNA decay machinery by the virus. The knowledge might be helpful in developing anti-viral compounds in future., Competing Interests: The authors have declared that no competing interests exist.

Based on microRNA (miR) microarray analysis, we previously found that miR22-5p expression is decreased in the mid-luteal endometrium of women with minimal/mild endometriosis. Bioinformatics analysis predicted that miR22-5p targets ten-eleven translocation (TET2) 3'-untranslated region. This study aimed to determine the regulation and roles of miR22-5p in the pathogenesis of minimal/mild endometriosis-associated infertility. MiR22-5p and TET2 expression in the mid-luteal endometrium from women with or without minimal/mild endometriosis was analyzed. After transfection with miR22-5p mimics or inhibitor, TET2 expression was analyzed by quantitative reverse transcription (RT-q) PCR, western blotting and immunohistochemistry. 5-Hydroxymethylcytosine was determined by immunofluorescence and dot blotting. Expression and promoter methylation of estrogen receptor 2 (ESR2) was measured by RT-qPCR and western blotting, and by bisulfite sequencing, respectively. We first established that miR22-5p expression decreased and TET2 expression increased in minimal/mild endometriosis during implantation window. TET2 was found to be a direct target of miR22-5p. MiR22-5p regulated the expression of ESR2, but did not directly affect methylation of its promoter region (-197/+359). Our results suggest that an imbalance in miR22-5p expression in the mid-luteal endometrium may be involved in minimal/mild endometriosis-associated infertility., Competing Interests: The authors have declared that no competing interests exist.

Cellular stress response mechanisms normally function to enhance survival and allow for cells to return to homeostasis following an adverse event. Cancer cells often co-opt these same mechanisms as a means to evade apoptosis and mitigate a state of constant cellular stress. Activating transcription factor 5 (ATF5) is upregulated under diverse stress conditions and is overexpressed in a variety of cancers. It was demonstrated ATF5 is a survival factor in transformed, but not normal cells. However, the regulation of ATF5 is not fully understood. The purpose of the present study was to investigate miRNA regulation at the 3' untranslated region (UTR) of ATF5, with the goal of demonstrating a reversal of the upregulation of ATF5 induced under diverse cellular stress in cancer cells. A multifactorial approach using in silico analysis was employed to identify miRNAs 433-3p, 520b-3p, and 129-5p as potential regulators of ATF5, based on their predicted binding sites over the span of the ATF5 3' UTR. Luciferase reporter assay data validated all three miRNA candidates by demonstrating direct binding to the target ATF5 3'. However, functional studies revealed miR-520b-3p as the sole candidate able to reverse the upregulation of ATF5 protein under diverse cellular stress. Additionally, miR-520b-3p levels were inversely related to ATF5 mRNA under endoplasmic reticulum stress and amino acid deprivation. This is the first evidence that regulation at the 3' UTR is involved in modulating ATF5 levels under cellular stress and suggests an important role for miRNA-520b-3p in the regulation of ATF5., Competing Interests: The authors have declared that no competing interests exist.

Poly(A) tails at the 3' end of eukaryotic messenger RNAs control mRNA stability and translation efficiency. Facilitated by various NGS methods, alternative polyadenylation sites determining the 3'-UTR length of gene transcripts have been extensively studied. However, poly(A) lengths demonstrating dynamic and developmental regulation remain largely unexplored. The recently developed NGS-based methods for genome-wide poly(A) profiling have promoted the study of genom-wide poly(A) dynamics. Here we present a straight forward NGS-method for poly(A) profiling, which applies a direct 3'-end adaptor ligation and the template switching for 5'-end adaptor ligation for cDNA library construction. Poly(A) lengths are directly calculated from base call data using a self-developed pipeline pA-finder. The libraries were directly sequenced from the 3'-UTR regions into the followed poly(A) tails, firstly on NextSeq 500 to produce single-end 300-nt reads, demonstrating the method feasibility and that optimization of the fragmented RNA size for cDNA library construction could detecting longer poly (A) tails. We next applied Poly(A)-seq cDNA libraries containing 40-nt and 120-nt poly(A) tail spike-in RNAs on HiSeq X-ten and NovaSeq 6000 to obtain 150-nt and 250-nt pair-end reads. The sequencing profiles of the spike-in RNAs demonstrated both high accuracy and high quality score in reading poly(A) tails. The poly(A) signal bleeding into the 3' adaptor sequence and a sharp decreased quality score at the junction were observed, allowing the modification of pA-finder to remove homopolymeric signal bleeding. We hope that wide applications of Poly(A)-seq help facilitate the study of the development- and disease-related poly(A) dynamics and regulation, and of the recent emerging mixed tailing regulation., Competing Interests: This study was supported by ABLife Inc.. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

MicroRNAs (miRNAs) are important regulators of gene expression, and their expression is associated with many physiological conditions. Here, we investigated potential associations between expression levels of miRNAs in human placenta and the onset of spontaneous term birth. Using RNA sequencing, we identified 54 miRNAs differentially expressed during spontaneous term labor compared to elective term births. Expression levels of 23 miRNAs were upregulated, whereas 31 were downregulated at least 1.5-fold. The upregulated miRNA miR-371a-5p putatively targets CPPED1, expression of which decreases during spontaneous birth. We used a luciferase reporter-based assay to test whether a miR-371a-5p mimic affected translation when it bound to the 3' untranslated region of CPPED1. In this setting, the miR-371a-5p mimic resulted in lower luciferase activity, which suggests that miR-371a-5p regulates levels of CPPED1. In conclusion, inversely correlated levels of miR-371a-5p and CPPED1 suggest a role for both in spontaneous delivery., Competing Interests: The authors declare no competing interests.

Many post-transcriptional mechanisms operate via mRNA 3'UTRs to regulate protein expression, and such controls are crucial for development. We show that homozygous mutations in two zebrafish exon junction complex (EJC) core genes rbm8a and magoh leads to muscle disorganization, neural cell death, and motor neuron outgrowth defects, as well as dysregulation of mRNAs subjected to nonsense-mediated mRNA decay (NMD) due to translation termination ≥ 50 nts upstream of the last exon-exon junction. Intriguingly, we find that EJC-dependent NMD also regulates a subset of transcripts that contain 3'UTR introns (3'UI) < 50 nts downstream of a stop codon. Some transcripts containing such stop codon-proximal 3'UI are also NMD-sensitive in cultured human cells and mouse embryonic stem cells. We identify 167 genes that contain a conserved proximal 3'UI in zebrafish, mouse and humans. foxo3b is one such proximal 3'UI-containing gene that is upregulated in zebrafish EJC mutant embryos, at both mRNA and protein levels, and loss of foxo3b function in EJC mutant embryos significantly rescues motor axon growth defects. These data are consistent with EJC-dependent NMD regulating foxo3b mRNA to control protein expression during zebrafish development. Our work shows that the EJC is critical for normal zebrafish development and suggests that proximal 3'UIs may serve gene regulatory function in vertebrates., Competing Interests: The authors have declared that no competing interests exist.

MicroRNAs (miRNAs) control the abundance of the majority of the vertebrate transcriptome. The recognition sequences, or target sites, for bilaterian miRNAs are found predominantly in the 3' untranslated regions (3'UTRs) of mRNAs, and are amongst the most highly conserved motifs within 3'UTRs. However, little is known regarding the evolutionary pressures that lead to loss and gain of such target sites. Here, we quantify the selective pressures that act upon miRNA target sites. Notably, selective pressure extends beyond deeply conserved binding sites to those that have undergone recent substitutions. Our approach reveals that even amongst ancient animal miRNAs, which exert the strongest selective pressures on 3'UTR sequences, there are striking differences in patterns of target site evolution between miRNAs. Considering only ancient animal miRNAs, we find three distinct miRNA groups, each exhibiting characteristic rates of target site gain and loss during mammalian evolution. The first group both loses and gains sites rarely. The second group shows selection only against site loss, with site gains occurring at a neutral rate, whereas the third loses and gains sites at neutral or above expected rates. Furthermore, mutations that alter the strength of existing target sites are disfavored. Applying our approach to individual transcripts reveals variation in the distribution of selective pressure across the transcriptome and between miRNAs, ranging from strong selection acting on a small subset of targets of some miRNAs, to weak selection on many targets for other miRNAs. miR-20 and miR-30, and many other miRNAs, exhibit broad, deeply conserved targeting, while several other comparably ancient miRNAs show a lack of selective constraint, and a small number, including mir-146, exhibit evidence of rapidly evolving target sites. Our approach adds valuable perspective on the evolution of miRNAs and their targets, and can also be applied to characterize other 3'UTR regulatory motifs., Competing Interests: The authors have declared that no competing interests exist.

There is a reciprocal relationship between the circadian and the reward systems. Polymorphisms in several circadian rhythm-related (clock) genes were associated with drug addiction. This study aims to search for associations between 895 variants in 39 circadian rhythm-related genes and opioid addiction (OUD). Genotyping was performed with the Smokescreen® array. Ancestry was verified by principal/MDS component analysis and the sample was limited to European Americans (EA) (OUD; n = 435, controls; n = 138). Nominally significant associations (p < 0.01) were detected for several variants in genes encoding vasoactive intestinal peptide receptor 2 (VIPR2), period circadian regulator 2 (PER2), casein kinase 1 epsilon (CSNK1E), and activator of transcription and developmental regulator (AUTS2), but no signal survived correction for multiple testing. There was intriguing association signal for the untranslated region (3' UTR) variant rs885863 in VIPR2, (p = .0065; OR = 0.51; 95% CI 0.31-0.51). The result was corroborated in an independent EA OUD sample (n = 398, p = 0.0036; for the combined samples). Notably, this SNP is an expression quantitative trait locus (cis-eQTL) for VIPR2 and a long intergenic non-coding RNA, lincRNA 689, in a tissue-specific manner, based on the Genotype-Tissue Expression (GTEx) project. Vasoactive intestinal peptide (VIP) is an important peptide of light-activated suprachiasmatic nucleus cells. It regulates diverse physiological processes including circadian rhythms, learning and memory, and stress response. This is the first report of an association of a VIPR2 variant and OUD. Additionally, analysis of combinations of single nucleotide polymorphisms (SNPs) genotypes revealed an association of PER2 SNP rs80136044, and SNP rs4128839, located 41.6 kb downstream of neuropeptide Y receptor type 1 gene, NPY1R (p = 3.4 × 10-6, OR = 11.4, 95% CI 2.7-48.2). The study provides preliminary insight into the relationship between genetic variants in circadian rhythm genes and long non-coding RNA (lncRNAs) in their vicinity, and opioid addiction., Competing Interests: The authors have declared that no competing interests exist.

Eastern equine encephalitis virus (EEEV), a mosquito-borne RNA virus, is one of the most acutely virulent viruses endemic to the Americas, causing between 30% and 70% mortality in symptomatic human cases. A major factor in the virulence of EEEV is the presence of four binding sites for the hematopoietic cell-specific microRNA, miR-142-3p, in the 3' untranslated region (3' UTR) of the virus. Three of the sites are "canonical" with all 7 seed sequence residues complimentary to miR-142-3p while one is "non-canonical" and has a seed sequence mismatch. Interaction of the EEEV genome with miR-142-3p limits virus replication in myeloid cells and suppresses the systemic innate immune response, greatly exacerbating EEEV neurovirulence. The presence of the miRNA binding sequences is also required for efficient EEEV replication in mosquitoes and, therefore, essential for transmission of the virus. In the current studies, we have examined the role of each binding site by point mutagenesis of the seed sequences in all combinations of sites followed by infection of mammalian myeloid cells, mosquito cells and mice. The resulting data indicate that both canonical and non-canonical sites contribute to cell infection and animal virulence, however, surprisingly, all sites are rapidly deleted from EEEV genomes shortly after infection of myeloid cells or mice. Finally, we show that the virulence of a related encephalitis virus, western equine encephalitis virus, is also dependent upon miR-142-3p binding sites., Competing Interests: The authors have declared that no competing interests exist

microRNAs (miRNAs) are potent regulators of gene expression that function in diverse developmental and physiological processes. Argonaute proteins loaded with miRNAs form the miRNA Induced Silencing Complexes (miRISCs) that repress gene expression at the post-transcriptional level. miRISCs target genes through partial sequence complementarity between the miRNA and the target mRNA's 3' UTR. In addition to being targeted by miRNAs, these mRNAs are also extensively regulated by RNA-binding proteins (RBPs) through RNA processing, transport, stability, and translation regulation. While the degree to which RBPs and miRISCs interact to regulate gene expression is likely extensive, we have only begun to unravel the mechanisms of this functional cooperation. An RNAi-based screen of putative ALG-1 Argonaute interactors has identified a role for a conserved RNA binding protein, HRPK-1, in modulating miRNA activity during C. elegans development. Here, we report the physical and genetic interaction between HRPK-1 and ALG-1/miRNAs. Specifically, we report the genetic and molecular characterizations of hrpk-1 and its role in C. elegans development and miRNA-mediated target repression. We show that loss of hrpk-1 causes numerous developmental defects and enhances the mutant phenotypes associated with reduction of miRNA activity, including those of lsy-6, mir-35-family, and let-7-family miRNAs. In addition to hrpk-1 genetic interaction with these miRNA families, hrpk-1 is required for efficient regulation of lsy-6 target cog-1. We report that hrpk-1 plays a role in processing of some but not all miRNAs and is not required for ALG-1/AIN-1 miRISC assembly. We suggest that HRPK-1 may functionally interact with miRNAs by both affecting miRNA processing and by enhancing miRNA/miRISC gene regulatory activity and present models for its activity., Competing Interests: The authors have declared that no competing interests exist.

3' downstream inhibitory stem plays a crucial role in locking rpoS mRNA 5' untranslated region in a self-inhibitory state. Here, we used optical tweezers to study the unfolding/refolding of rpoS inhibitory stem in the absence and presence of Mg2+. We found adding Mg2+ decreased the free energy of the RNA junction without re-arranging its secondary structure, through confirming that this RNA formed a canonical RNA three-way junction. We suspected increased free energy might change the relative orientation of different stems of rpoS and confirmed this by small angle X-ray scattering. Such changed conformation may improve Hfq-bridged annealing between sRNA and rpoS RNA inhibitory stem. We established a convenient route to analyze the changes of RNA conformation and folding dynamics by combining optical tweezers with X-ray scattering methods. This route can be easily applied in the studies of other RNA structure and ligand-RNA., Competing Interests: The authors have declared that no competing interests exist.

Errors throughout gene expression are likely deleterious, hence genomes are under selection to ameliorate their consequences. Additional stop codons (ASCs) are in-frame nonsense 'codons' downstream of the primary stop which may be read by translational machinery should the primary stop have been accidentally read through. Prior evidence in several eukaryotes suggests that ASCs are selected to prevent potentially-deleterious consequences of read-through. We extend this evidence showing that enrichment of ASCs is common but not universal for single cell eukaryotes. By contrast, there is limited evidence as to whether the same is true in other taxa. Here, we provide the first systematic test of the hypothesis that ASCs act as a fail-safe mechanism in eubacteria, a group with high read-through rates. Contra to the predictions of the hypothesis we find: there is paucity, not enrichment, of ASCs downstream; substitutions that degrade stops are more frequent in-frame than out-of-frame in 3' sequence; highly expressed genes are no more likely to have ASCs than lowly expressed genes; usage of the leakiest primary stop (TGA) in highly expressed genes does not predict ASC enrichment even compared to usage of non-leaky stops (TAA) in lowly expressed genes, beyond downstream codon +1. Any effect at the codon immediately proximal to the primary stop can be accounted for by a preference for a T/U residue immediately following the stop, although if anything, TT- and TC- starting codons are preferred. We conclude that there is no compelling evidence for ASC selection in eubacteria. This presents an unusual case in which the same error could be solved by the same mechanism in eukaryotes and prokaryotes but is not. We discuss two possible explanations: that, owing to the absence of nonsense mediated decay, bacteria may solve read-through via gene truncation and in eukaryotes certain prion states cause raised read-through rates., Competing Interests: The authors have declared that no competing interests exist.

Growth hormone (GH) is an important hormone released by the pituitary gland that plays a key role in the growth and development of organisms. In our study, TargetScan analysis and the dual luciferase reporter assays were used to predict and screen for miRNAs that might act on the rat Gh1 gene, and we identified miR-543-5p. Then, the GH3 cell line and the primary rat pituitary cells were transfected with miRNA mimic, inhibitor, and siRNA. We detected the Gh1 gene expression and the GH secretion by real-time PCR and ELISAs, respectively, to verify the regulatory effect of miR-543-5p on GH secretion. The results showed that miR-543-5p can inhibit Gh1 mRNA expression and reduce GH secretion. MiR-543-5p inhibitor upregulated Gh1 mRNA expression and increased GH secretion compared with the negative control. In summary, miR-543-5p downregulates Gh1 expression, resulting in a decrease in GH synthesis and secretion, which demonstrates the important role of miRNAs in regulating GH and animal growth and development., Competing Interests: The authors have declared that no competing interests exist.

Tauopathies are neurodegenerative diseases that affect millions of people worldwide including those with Alzheimer's disease. While many efforts have focused on understanding the role of tau protein in neurodegeneration, there has been little done to systematically analyze and study the structures within tau's encoding RNA and their connection to disease pathology. Knowledge of RNA structure can provide insights into disease mechanisms and how to affect protein production for therapeutic benefit. Using computational methods based on thermodynamic stability and evolutionary conservation, we identified structures throughout the tau pre-mRNA, especially at exon-intron junctions and within the 5' and 3' untranslated regions (UTRs). In particular, structures were identified at twenty exon-intron junctions. The 5' UTR contains one structured region, which lies within a known internal ribosome entry site. The 3' UTR contains eight structured regions, including one that contains a polyadenylation signal. A series of functional experiments were carried out to assess the effects of mutations associated with mis-regulation of alternative splicing of exon 10 and to identify regions of the 3' UTR that contain cis-regulatory elements. These studies defined novel structural regions within the mRNA that affect stability and pre-mRNA splicing and may lead to new therapeutic targets for treating tau-associated diseases., Competing Interests: The authors have declared that no competing interests exist.

Background: Great progress has been made in the understanding of inflammatory processes in psoriasis. However, clarifying the role of genetic variability in processes regulating inflammation, including post-transcriptional regulation by microRNA (miRNA), remains a challenge., Objectives: We therefore investigated single nucleotide polymorphisms (SNPs) with a predicted change in the miRNA/mRNA interaction of genes involved in the psoriasis inflammatory processes., Methods: Studied SNPs rs2910164 C/G-miR-146a, rs4597342 T/C-ITGAM, rs1368439 G/T-IL12B, rs1468488 C/T-IL17RA were selected using a bioinformatics analysis of psoriasis inflammation-associated genes. These SNPs were then genotyped using a large cohort of women with psoriasis (n = 241) and healthy controls (n = 516)., Results: No significant association with psoriasis was observed for rs2910164, rs1368439, and rs1468488 genotypes. However, the major allele T of rs4597342 -ITGAM was associated with approximately 28% higher risk for psoriasis in comparison to the patients with the C allele (OR = 1.28, 95% CI 1.01-1.61, p = 0.037). In case of genotypes, the effect of the T allele indicates the dominant model of disease penetrance as the CT and TT genotypes increase the chance of psoriasis up to 42% in comparison to CC homozygotes of rs4597342 (OR = 1.42, 95% CI = 1.05-1.94, p = 0.025)., Conclusion: SNP rs4597342 in 3'UTR of ITGAM influencing miR-21 binding may be considered a risk factor for psoriasis development. Upregulated miR-21 in psoriasis is likely to inhibit CD11b production in the case of the rs4597342 T allele which may lead to Mac-1 dysfunction, resulting in an aberrant function of innate immune cells and leading to the production of cytokines involved in psoriasis pathogenesis., Competing Interests: The authors have declared that no competing interests exist.

Non-small cell lung cancer (NSCLC) is the most common cancer worldwide. MicroRNAs have been shown to be correlated with biological processes of various tumors. In this study, we observed that the expression of miR-582-5p was lower in NSCLC tissues than that in para-carcinoma tissues. Ectopic expression of miR-582-5p significantly inhibited NCI-H358 cell proliferation and invasion. Knockdown of miR-582-5p showed the opposite results, with cell growth rate and the invasive capacity of PC-9 cells enhanced. Furthermore, we elucidated that NOTCH1 is a target of miR-582-5p and there is an inverse correlation between miR-582-5p and NOTCH1 expression in NSCLC tissues. Overexpression of NOTCH1 in miR-582-5p-overexpressing NCI-H358 cells could partially reverse the inhibition of cell proliferation and invasion by miR-582-5p. Thus, our research demonstrated that miR-582-5p suppresses NSCLC cell lines' growth and invasion via targeting oncoprotein NOTCH1 and restoration of miR-582-5p might be feasible therapeutic strategy for NSCLC., Competing Interests: The authors have declared that no competing interests exist.

Alternative polyadenylation (APA) in 3' untranslated regions (3' UTR) plays an important role in regulating transcript abundance, localization, and interaction with microRNAs. Length-variation of 3'UTRs by APA contributes to efficient proliferation of cancer cells. In this study, we investigated APA in single cancer cells and tumor microenvironment cells to understand the physiological implication of APA in different cell types. We analyzed APA patterns and the expression level of genes from the 515 single-cell RNA sequencing (scRNA-seq) dataset from 11 breast cancer patients. Although the overall 3'UTR length of individual genes was distributed equally in tumor and non-tumor cells, we found a differential pattern of polyadenylation in gene sets between tumor and non-tumor cells. In addition, we found a differential pattern of APA across tumor types using scRNA-seq data from 3 glioblastoma patients and 1 renal cell carcinoma patients. In detail, 1,176 gene sets and 53 genes showed the distinct pattern of 3'UTR shortening and over-expression as signatures for five cell types including B lymphocytes, T lymphocytes, myeloid cells, stromal cells, and breast cancer cells. Functional categories of gene sets for cellular proliferation demonstrated concordant regulation of APA and gene expression specific to cell types. The expression of APA genes in breast cancer was significantly correlated with the clinical outcome of earlier stage breast cancer patients. We identified cell type-specific APA in single cells, which allows the identification of cell types based on 3'UTR length variation in combination with gene expression. Specifically, an immune-specific APA signature in breast cancer could be utilized as a prognostic marker of early stage breast cancer., Competing Interests: W.P. is employed by a commercial company, GENINUS Inc. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome., Competing Interests: The authors have declared that no competing interests exist.