Search

Previous research suggests that genetic variants in the interleukin-6 (IL-6) gene contribute to sport-related traits and athletic performance. We aimed to identify sequence variants in the IL-6 gene region comprising the 3' untranslated region (UTR) in the Turkish triathletes and sedentary individuals and assessed their putative roles in tendency to athletic phenotype. Sequence variants were identified in the Turkish triathletes (n = 47) and sedentary individu-als (n = 46) by Sanger sequencing. Allele/genotype frequencies and linkage disequilibrium (LD) patterns were calculated by the Haploview program. The functional significance of the detected variants was analyzed using in silico prediction tools. Four single nucleotide vari-ants (rs13306435, rs747302620, rs2069849, rs13306436) were detected in saliva samples of the participants by sequencing the target region. Notably, rs13306436-3'UTR/IL-6 was only seen in the triathletes, while the exonic rs747302620 was observed in only sedentary group. Also, rs13306436G>A causes loss/gain sites for binding multiple miRNAs that may be associated with athletic performance. Our findings indicate that the 3'UTR/IL-6 may have functional relevance in determining sports talent. Future comprehensive studies focusing on the IL-6 gene in athletes may pave the way for not only determining the athletic status of the individuals but also have implications for translational medicine.

IntroductionTobacco vein mottling virus (TVMV) is a member of the family Potyviridae. The 3’ untranslated region (3’UTR) of viral genomic RNA has been reported to significantly impact viral infection. Nevertheless, the role of the TVMV 3’UTR during viral infection remains unknown.MethodsHere, a 3’UTR-GFP expression vector was transiently expressed in Nicotiana benthamiana, in which the 3’UTR of TVMV was introduced upstream of the green fluorescent protein (GFP) gene. Transcriptome sequencing was performed to analyze the genes associated with plant resistance. The effect of the TVMV 3’UTR on GFP expression was studied using an Agrobacterium-mediated transient expression assay, revealing that the TVMV 3’UTR significantly inhibited GFP expression. Transcriptome analysis of differentially expressed genes in 3’UTR-GFP in N. benthamiana was performed to elucidate the why the TVMV 3’UTR inhibited GFP expression.ResultsEighty genes related to plant disease resistance were differentially expressed, including 29 upregulated and 51 downregulated genes. Significantly upregulated genes included those encoding the calcium-binding protein CML24, leucine-rich repeat receptor-like tyrosine-protein kinase, and respiratory burst oxidase homolog protein E. The significantly downregulated genes included calcium-binding protein 7, ethylene-responsive transcription factor 10, endoglucanase 5, and receptor-like protein kinase.DiscussionThese findings indicate that the 3’UTR of TVMV may inhibit the expression of GFP gene by inducing the expression of plant resistance genes. This study provides a theoretical basis for further research on the function and mechanism of the TVMV 3’UTR.

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

The striatal D1 dopamine receptor (D1R) and A2a adenosine receptor (A2aR) signaling pathways play important roles in drug-related behaviors. These receptors activate the Golf protein comprised of a specific combination of αolfβ2γ7 subunits. During assembly, the γ7 subunit sets the cellular level of the Golf protein. In turn, the amount of Golf protein determines the collective output from both D1R and A2aR signaling pathways. This study shows the Gng7 gene encodes multiple γ7 transcripts differing only in their non-coding regions. In striatum, Transcript 1 is the predominant isoform. Preferentially expressed in the neuropil, Transcript 1 is localized in dendrites where it undergoes post-transcriptional regulation mediated by regulatory elements in its 3′ untranslated region that contribute to translational suppression of the γ7 protein. Earlier studies on gene-targeted mice demonstrated loss of γ7 protein disrupts assembly of the Golf protein. In the current study, morphological analysis reveals the loss of the Golf protein is associated with altered dendritic morphology of medium spiny neurons. Finally, behavioral analysis of conditional knockout mice with cell-specific deletion of the γ7 protein in distinct populations of medium spiny neurons reveals differential roles of the Golf protein in mediating behavioral responses to cocaine. Altogether, these findings provide a better understanding of the regulation of γ7 protein expression, its impact on Golf function, and point to a new potential target and mechanisms for treating addiction and related disorders.

Abstract BRAFV600E comes as two main splicing variants. The well-studied ref isoform and the recently discovered X1 isoform are co-expressed in cancer cells and differ in terms of 3’UTR length and sequence, as well as C-term protein sequence. Here, we use a melanoma model in zebrafish to study the role played by each isoform in larval pigmentation, nevi formation, and their progression into melanoma tumours. We show that both BRAFV600E-ref and BRAFV600E-X1 proteins promote larval pigmentation and nevi formation, while melanoma-free survival curves performed in adult fish indicate that BRAFV600E-ref protein is a much stronger melanoma driver that BRAFV600E-X1 protein. Crucially, we also show that the presence of the 3’UTR suppresses the effect of ref protein. Our data highlight the necessity to undertake a systematic study of BRAFV600E isoforms, in order to uncover the full spectrum of their kinase-(in)dependent and coding-(in)dependent functions, hence to develop more informed strategies for therapeutic targeting.

The biotechnological landscape has witnessed significant growth in biological therapeutics particularly in the field of recombinant protein production. Here we investigate the function of 3′UTR cis-regulatory elements in increasing mRNA and protein levels in different biological therapeutics and model systems, spanning from monoclonal antibodies to mRNA vaccines. We explore the regulatory function of iPLUS - a universal sequence capable of consistently augmenting recombinant protein levels. By incorporating iPLUS in a vector to express a monoclonal antibody used in immunotherapy, in a mammalian cell line used by the industry (ExpiCHO), trastuzumab production increases by 2-fold. As yeast Pichia pastoris is widely used in the manufacture of industrial enzymes and pharmaceuticals, we then used iPLUS in tandem (3x) and iPLUSv2 (a variant of iPLUS) to provide proof-of-concept data that it increases the production of a reporter protein more than 100-fold. As iPLUS functions by also increasing mRNA levels, we hypothesize that these sequences could be used as an asset in the mRNA vaccine industry. In fact, by including iPLUSv2 downstream of Spike we were able to double its production. Moreover, the same effect was observed when we introduced iPLUSv2 downstream of MAGEC2, a tumor-specific antigen tested for cancer mRNA vaccines. Taken together, our study provides data (TLR4) showing that iPLUS may be used as a valuable asset in a variety of systems used by the biotech and biopharmaceutical industry. Our results underscore the critical role of non-coding sequences in controlling gene expression, offering a promising avenue to accelerate, enhance, and cost-effectively optimize biopharmaceutical production processes.

MicroRNAs (miRNAs) are small non-coding RNAs that participate as powerful genetic regulators. MiRNAs can interfere with cellular processes by interacting with a broad spectrum of target genes under physiological and pathological states, including cancer development and progression. Major histocompatibility complex major histocompatibility complex class I-related chain A (MICA) belongs to a family of proteins that bind the natural-killer group 2, member D (NKG2D) receptor on Natural Killer cells and other cytotoxic lymphocytes. MICA plays a crucial role in the host’s innate immune response to several disease settings, including cancer. MICA harbors various single nucleotide polymorphisms (SNPs) located in its 3′-untranslated region (3′UTR), a characteristic that increases the complexity of MICA regulation, favoring its post-transcriptional modulation by miRNAs under physiological and pathological conditions. Here, we conducted an in-depth analysis of MICA 3′UTR sequences according to each MICA allele described to date using NCBI database. We also systematically evaluated interactions between miRNAs and their putative targets on MICA 3′UTR containing SNPs using in silico analysis. Our in silico results showed that MICA SNPs rs9266829, rs 1880, and rs9266825, located in the target sequence of miRNAs hsa-miR-106a-5p, hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, hsa-miR-93, hsa-miR-1207.5p, and hsa-miR-711 could modify the binding free energy between −8.62 and −18.14 kcal/mol, which may affect the regulation of MICA expression. We believe that our results may provide a starting point for further exploration of miRNA regulatory effects depending on MICA allelic variability; they may also be a guide to conduct miRNA in silico analysis for other highly polymorphic genes.

Abstract The body muscle is an important tissue used in organisms for proper viability and locomotion. Although this tissue is generally well studied and characterized, and many pathways have been elucidated throughout the years, we still lack a comprehensive understanding of its transcriptome and how it controls muscle development and function. Here, we have updated a nuclear FACS sorting-based methodology to isolate and sequence a high-quality muscle transcriptome from Caenorhabditis elegans mixed-stage animals. We have identified 2848 muscle-specific protein-coding genes, including 78 transcription factors and 206 protein-coding genes containing an RNA binding domain. We studied their interaction network, performed a detailed promoter analysis, and identified novel muscle-specific cis-acting elements. We have also identified 16 high-quality muscle-specific miRNAs, studied their function in vivo using fluorochrome-based analyses, and developed a high-quality C. elegans miRNA interactome incorporating other muscle-specific datasets produced by our lab and others. Our study expands our understanding of how muscle tissue functions in C. elegans and in turn provides results that can in the future be applied to humans to study muscular-related diseases.

Abstract Background Folates are crucial for the biosynthesis of nucleotides and amino acids, essential for cell proliferation and development. Folate deficiency induces DNA damage, developmental defects, and tumorigenicity. The obligatory enzyme folylpolyglutamate synthetase (FPGS) mediates intracellular folate retention via cytosolic and mitochondrial folate polyglutamylation. Our previous paper demonstrated the association of the cytosolic FPGS (cFPGS) with the cytoskeleton and various cell protrusion proteins. Based on these recent findings, the aim of the current study was to investigate the potential role of cFPGS at cell protrusions. Results Here we uncovered a central role for two G-quadruplex (GQ) motifs in the 3′UTR of FPGS mediating the localization of cFPGS mRNA and protein at cell protrusions. Using the MBSV6-loop reporter system and fluorescence microscopy, we demonstrate that following folate deprivation, cFPGS mRNA is retained in the endoplasmic reticulum, whereas upon 15 min of folate repletion, this mRNA is rapidly translocated to cell protrusions in a 3′UTR- and actin-dependent manner. The actin dependency of this folate-induced mRNA translocation is shown by treatment with Latrunculin B and inhibitors of the Ras homolog family member A (RhoA) pathway. Upon folate repletion, the FPGS 3′UTR GQs induce an amoeboid/mesenchymal hybrid cell phenotype during migration and invasion through a collagen gel matrix. Targeted disruption of the 3′UTR GQ motifs by introducing point mutations or masking them by antisense oligonucleotides abrogated cell protrusion targeting of cFPGS mRNA. Conclusions Collectively, the GQ motifs within the 3′UTR of FPGS regulate its transcript and protein localization at cell protrusions in response to a folate cue, inducing cancer cell invasive phenotype. These novel findings suggest that the 3′UTR GQ motifs of FPGS constitute an attractive druggable target aimed at inhibition of cancer invasion and metastasis.

Decompression sickness (DCS) with neurological disorders includes an inappropriate inflammatory response which degenerates slowly, even after the disappearance of the bubbles. There is high inter-individual variability in terms of the occurrence of DCS that could have been mastered by the selection and then the breeding of DCS-resistant rats. We hypothesized the selection of single-nucleotide polymorphisms (SNPs) linked to autoimmunity operated upon a generation of a DCS-resistant strain of rats. We used the candidate gene approach and targeted SNPs linked to the signaling cascade that directly regulates inflammation of innate immunity transiting by the Toll-like receptors. Twenty candidate SNPs were investigated in 36 standard rats and 33 DCS-resistant rats. For the first time, we identify a diplotype (i.e., with matched haplotypes)—when coinherited—that strengthens protection against DCS, which is not strictly homozygous and suggests that a certain tolerance may be considered. We deduced an ideal haplotype of six variants from it (MyD88_50-T, _49-A, _97-C coupled to NFKB_85-T, _69-T, _45-T) linked to the resistant phenotype. Four among the six identified variants are located in pre- and/or post-transcriptional areas regulating MyD88 or NFKB1 expression. Because of missense mutations, the other two variants induce a structural change in the NFKB1 protein complex including one damage alteration according to the Missense3D algorithm. In addition to the MyD88/NFKB1 haplotype providing rats with a strong resistance to DCS, this also highlights the importance that the immune response, here linked to the genetic heritage, can have in the development of DCS and offer a new perspective for therapeutic strategies.

IntroductionMacrophages are essential cells of the immune system that alter their inflammatory profile depending on their microenvironment. Alternative polyadenylation in the 3’UTR (3’UTR-APA) and intronic polyadenylation (IPA) are mechanisms that modulate gene expression, particularly in cancer and activated immune cells. Yet, how polarization and colorectal cancer (CRC) cells affect 3’UTR-APA and IPA in primary human macrophages was unclear.MethodsIn this study, we isolated primary human monocytes from healthy donors, differentiated and polarized them into a pro-inflammatory state and performed indirect co-cultures with CRC cells. ChrRNA-Seq and 3’RNA-Seq was performed to quantify gene expression and characterize new 3’UTR-APA and IPA mRNA isoforms.ResultsOur results show that polarization of human macrophages from naïve to a pro-inflammatory state causes a marked increase of proximal polyA site selection in the 3’UTR and IPA events in genes relevant to macrophage functions. Additionally, we found a negative correlation between differential gene expression and IPA during pro-inflammatory polarization of primary human macrophages. As macrophages are abundant immune cells in the CRC microenvironment that either promote or abrogate cancer progression, we investigated how indirect exposure to CRC cells affects macrophage gene expression and 3’UTR-APA and IPA events. Co-culture with CRC cells alters the inflammatory phenotype of macrophages, increases the expression of pro-tumoral genes and induces 3’UTR-APA alterations. Notably, some of these gene expression differences were also found in tumor-associated macrophages of CRC patients, indicating that they are physiologically relevant. Upon macrophage pro-inflammatory polarization, SRSF12 is the pre-mRNA processing gene that is most upregulated. After SRSF12 knockdown in M1 macrophages there is a global downregulation of gene expression, in particular in genes involved in gene expression regulation and in immune responses.DiscussionOur results reveal new 3’UTR-APA and IPA mRNA isoforms produced during pro-inflammatory polarization of primary human macrophages and CRC co-culture that may be used in the future as diagnostic or therapeutic tools. Furthermore, our results highlight a function for SRSF12 in pro-inflammatory macrophages, key cells in the tumor response.

mRNA vaccines have been shown to be effective in combating the COVID-19 pandemic. The amount of research on the use of mRNAs as preventive and therapeutic modalities has undergone explosive growth in the last few years. Nonetheless, the issue of the stability of mRNA molecules and their translation efficiency remains incompletely resolved. These characteristics of mRNA directly affect the expression level of a desired protein. Regulatory elements of RNA—5′ and 3′ untranslated regions (UTRs)—are responsible for translation efficiency. An optimal combination of the regulatory sequences allows mRNA to significantly increase the target protein’s expression. We assessed the translation efficiency of mRNA encoding of firefly luciferase with various 5′ and 3′UTRs in vitro on cell lines DC2.4 and THP1. We found that mRNAs containing 5′UTR sequences from eukaryotic genes HBB, HSPA1A, Rabb, or H4C2, or from the adenoviral leader sequence TPL, resulted in higher levels of luciferase bioluminescence 4 h after transfection of DC2.4 cells as compared with 5′UTR sequences used in vaccines mRNA-1273 and BNT162b2 from Moderna and BioNTech. mRNA containing TPL as the 5′UTR also showed higher efficiency (as compared with the 5′UTR from Moderna) at generating a T-cell response in mice immunized with mRNA vaccines encoding a multiepitope antigen. By contrast, no effects of various 5′UTRs and 3′UTRs were detectable in THP1 cells, suggesting that the observed effects are cell type specific. Further analyses enabled us to identify potential cell type-specific RNA-binding proteins that differ in landing sites within mRNAs with various 5′UTRs and 3′UTRs. Taken together, our data indicate high translation efficiency of TPL as a 5′UTR, according to experiments on DC2.4 cells and C57BL/6 mice.

The 3' untranslated regions (UTRs) modulate gene expression levels by regulating mRNA stability and translation. We previously showed that the replacement of the negative regulatory elements from the 3'UTR of glial cell line-derived neurotrophic factor (GDNF) resulted in increased endogenous GDNF expression while retaining its normal spatiotemporal expression pattern. Here, we have developed a methodology for the generation of in vivo hyper- and hypomorphic alleles via 3'UTR targeting using the CRISPR/Cas9 system. We demonstrate that CRISPR/Cas9-mediated excision of a long inhibitory sequence from Gdnf native 3'UTR in mouse zygotes increases the levels of endogenous GDNF with similar phenotypic alterations in embryonic kidney development as we described in GDNF constitutive and conditional hypermorphic mice. Furthermore, we show that CRISPR/Cas9-mediated targeting of 3’UTRs in vivo allows the modulation of the expression levels of two other morphogens, Gdf11 and Bdnf. Together, our work demonstrates the power of in vivo 3’UTR editing using the CRISPR/Cas9 system to create hyper- and hypomorphic alleles, suggesting wide applicability in studies on gene function and potentially, in gene therapy.

Abstract Background Genome-wide association studies (GWAS) have uncovered thousands of genetic variants that are associated with complex human traits and diseases. miRNAs are single-stranded non-coding RNAs. In particular, genetic variants located in the 3’UTR region of mRNAs may play an important role in gene regulation through their interaction with miRNAs. Existing studies have not been thoroughly conducted to elucidate 3’UTR variants discovered through GWAS. The goal of this study is to analyze patterns of GWAS functional variants located in 3’UTRs about their relevance in the network between hosting genes and targeting miRNAs, and elucidate the association between the genes harboring these variants and genetic traits. Methods We employed MIGWAS, ANNOVAR, MEME, and DAVID software packages to annotate the variants obtained from GWAS for 31 traits and elucidate the association between their harboring genes and their related traits. We identified variants that occurred in the motif regions that may be functionally important in affecting miRNA binding. We also conducted pathway analysis and functional annotation on miRNA targeted genes harboring 3’UTR variants for a trait with the highest percentage of 3’UTR variants occurring. Results The Child Obesity trait has the highest percentage of 3’UTR variants (75%). Of the 16 genes related to the Child Obesity trait, 5 genes (ETV7, GMEB1, NFIX, ZNF566, ZBTB40) had a significant association with the term DNA-Binding (p

Dita Maria Virginia,1 Christine Patramurti,2 Fenty,1,3 Christianus Heru Setiawan,1 Jeffry Julianus,2 Phebe Hendra,1 Nicholas Adi Perdana Susanto3 1Department of Pharmacology and Clinical Pharmacy, Faculty of Pharmacy, Sanata Dharma University, Yogyakarta, Indonesia; 2Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Sanata Dharma University, Yogyakarta, Indonesia; 3Bethesda Lempuyangwangi Hospital, Yogyakarta, IndonesiaCorrespondence: Dita Maria Virginia, Universitas Sanata Dharma, Mrican, Tromol Pos 29, Yogyakarta, 55281, Indonesia, Tel +62 274 513301, Fax +62 274 562383, Email virginia@usd.ac.idPurpose: This study aimed to explore the association of rs857148 A>C as 3ʹUTR variants with blood pressure, HbA1c profile, and lipid profiles as cardiometabolic parameters among patients with T2DM receiving metformin.Patients and Methods: This cross-sectional analytic research was conducted with 114 consecutively selected patients with T2DM. Polymerase chain reaction-restriction fragment length polymorphism was conducted to determine rs857148. A total of 108 patients fulfilled inclusion and exclusion criteria.Results: Genotype distribution agreed with the Hardy Weinberg Equation for Equilibrium (p> 0.05) but wildtype allele was found as the minor allele. Subjects with CC genotype and C allele had enhanced HbA1c levels (OR=7.12; 95% CI=1.05– 48.26; p=0.04; OR=1.66; 95% CI=1.06– 2.60; p=0.03, respectively). It was confirmed by dominant model whereas subjects with AA tended to have reduced HbA1c compared to AC+CC genotype (OR=0.15; 95% CI=0.02– 0.97; p=0.047). AC genotype had significant correlation to total cholesterol (OR=1.05; 95% CI=1.01– 1.10; p=0.03) compared to AA genotype.Conclusion: We conclude that polymorphism of rs87148, specifically CC genotype and C allele, has a significant association with HbA1c and total cholesterol after considering oral hypoglycemia agent dose, age, gender, and combination therapy, compared to AA genotype. Future studies that involve a larger sample population and more rigorous selection criteria are required.Keywords: PRKAA2, rs857148, 3ʹUTR, cardiometabolic, type 2 diabetes mellitus

Inducing senescence in cancer cells is an emerging strategy for cancer therapy. The dysregulation and mutation of genes encoding cyclin‐dependent kinases (CDKs) have been implicated in various human cancers. However, whether CDK can induce cancer cell senescence remains poorly understood. We observed that CDK16 expression was high in multiple cancer types, including lung cancer, whereas various replicative senescence models displayed low CDK16 expression. CDK16 knockdown caused senescence‐associated phenotypes in lung cancer cell lines. Interestingly, the CDK16 3′ UTR was shortened in cancer and lengthened in senescence models, which was regulated by alternative polyadenylation (APA). The longer 3′UTR [using the distal polyA (pA) site] generated less protein than the shorter one (using the proximal pA site). Since microRNAs (miRNAs) usually bind to the 3′UTR of target genes to suppress their expression, we investigated whether miRNAs targeting the region between the shortened and longer 3′UTR are responsible for the reduced expression. We found that miR‐485‐5p targeted the 3′UTR between the distal and proximal pA site and caused senescence‐associated phenotypes by reducing protein production from the longer CDK16 transcript. Of note, CDK16 knockdown led to a reduced expression of MYC proto‐oncogene, bHLH transcription factor (MYC) and CD274 molecule (PD‐L1), which in turn enhanced the tumor‐suppressive effects of senescent cancer cells. The present study discovered that CDK16, whose expression is under the regulation of APA and miR‐485‐5p, is a potential target for prosenescence therapy for lung cancer.

mRNA-based therapeutics have been found to be a promising treatment strategy in immunotherapy, gene therapy, and cancer treatments. Effectiveness of mRNA therapeutics depends on the level and duration of a desired protein’s expression, which is determined by various cis- and trans-regulatory elements of the mRNA. Sequences of 5′ and 3′ untranslated regions (UTRs) are responsible for translational efficiency and stability of mRNA. An optimal combination of the regulatory sequences allows researchers to significantly increase the target protein’s expression. Using both literature data and previously obtained experimental data, we chose six sequences of 5′UTRs (adenoviral tripartite leader [TPL], HBB, rabbit β-globin [Rabb], H4C2, Moderna, and Neo2) and five sequences of 3′UTRs (mtRNR-EMCV, mtRNR-AES, mtRNR-mtRNR, BioNTech, and Moderna). By combining them, we constructed 30 in vitro transcribed RNAs encoding firefly luciferase with various combinations of 5′- and 3′UTRs, and the resultant bioluminescence was assessed in the DC2.4 cell line at 4, 8, 24, and 72 h after transfection. The cellular data enabled us to identify the best seven combinations of 5′- and 3′UTRs, whose translational efficiency was then assessed in BALB/c mice. Two combinations of 5′- and 3′UTRs (5′Rabb-3′mtRNR-EMCV and 5′TPL-3′Biontech) led to the most pronounced increase in the luciferase amount in the in vivo experiment in mice. Subsequent analysis of the stability of the mRNA indicated that the increase in luciferase expression is explained primarily by the efficiency of translation, not by the number of RNA molecules. Altogether, these findings suggest that 5′UTR-and-3′UTR combinations 5′Rabb-3′mtRNR- EMCV and 5′TPL-3′Biontech lead to high expression of target proteins and may be considered for use in preventive and therapeutic modalities based on mRNA.

BackgroundHuman leukocyte antigen G (HLA-G) is an immune checkpoint molecule with relevance in several cancers. The aim of this study was to evaluate the potential role of soluble HLA-G (sHLA-G), its genetic polymorphisms and its haplotype structure in the susceptibility and prognosis of primary cervical cancer in a Chinese Han population.MethodsWe investigated sHLA-G plasma levels and 3’ untranslated region (3’UTR) polymorphisms through ELISA and direct DNA sequencing, respectively, in cervical cancer patients (120 cases) and healthy control women (96 cases). The data were analyzed for associations using PowerMarker, Haploview, and GraphPad Prism.ResultsIn this study, 8 polymorphic sites, 16 haplotypes and 23 diplotypes in the HLA-G 3’UTR were identified in our study population. We observed that each pair of 8 polymorphic sites exhibited linkage disequilibrium. The heterozygote CT genotype at position +3422 (rs17875408) was more common in cervical cancer patients than in healthy women (OR=5.285, P0.7). Among patients, mean sHLA-G levels increased with increasing FIGO stages but were not related to the overall survival time.ConclusionsThe results of the present study enhance our understanding of how HLA-G 3’UTR polymorphisms can influence the peripheral sHLA-G plasma level and play a key role in cervical carcinogenesis. This study further confirmed that sHLA-G may represent a novel plasma biomarker for the prognosis and potential therapeutic target of cervical cancer.

Glutamine synthetase (GS) is the key enzyme of nitrogen assimilation induced under nitrogen limiting conditions. The carbon skeleton of glutamate and glutamine, 2-oxoglutarate, is supplied from the TCA cycle, but how this metabolic flow is controlled in response to nitrogen availability remains unknown. We show that the expression of the E1o component of 2-oxoglutarate dehydrogenase, SucA, is repressed under nitrogen limitation in Salmonella enterica and Escherichia coli. The repression is exerted at the post-transcriptional level by an Hfq-dependent sRNA GlnZ generated from the 3′UTR of the GS-encoding glnA mRNA. Enterobacterial GlnZ variants contain a conserved seed sequence and primarily regulate sucA through base-pairing far upstream of the translation initiation region. During growth on glutamine as the nitrogen source, the glnA 3′UTR deletion mutants expressed SucA at higher levels than the S. enterica and E. coli wild-type strains, respectively. In E. coli, the transcriptional regulator Nac also participates in the repression of sucA. Lastly, this study clarifies that the release of GlnZ from the glnA mRNA by RNase E is essential for the post-transcriptional regulation of sucA. Thus, the mRNA coordinates the two independent functions to balance the supply and demand of the fundamental metabolites.

ABSTRACT Tembusu virus (TMUV), an avian mosquito-borne flavivirus, was first identified from Culex tritaeniorhynchus in 1955. To validate the effects of the 3′-untranslated region (3′UTR) in viral host-specific adaptation, we generated a set of chimeric viruses using CQW1 (duck strain) and MM 1775 (mosquito strain) as backbones with heterogeneous 3′UTRs. Compared with rMM 1775, rMM-CQ3′UTR (recombinant MM 1775 virus carrying the 3′UTR of CQW1) exhibited enhanced proliferation in vitro, with peak titers increasing by 5-fold in duck embryonic fibroblast (DEF) cells or 12-fold in baby hamster kidney (BHK-21) cells; however, the neurovirulence of rMM-CQ3′UTR was attenuated in 14-day-old Kunming mice via intracranial injection, with slower weight loss, lower mortality, and reduced viral loads. In contrast, rCQ-MM3′UTR showed similar growth kinetics in vitro and neurovirulence in mice compared with those of rCQW1. Then, the Stem-loop I (SLI) structure, which showed the highest variation within the 3′UTR between CQW1 and MM 1775, was further chosen for making chimeric viruses. The peak titers of rMM-CQ3′UTRSLI displayed a 15- or 4-fold increase in vitro, and the neurovirulence in mice was attenuated, compared with that of rMM 1775; rCQ-MM3′UTRSLI displayed comparable multiplication ability in vitro but was significantly attenuated in mice, in contrast with rCQW1. In conclusion, we demonstrated that the TMUV SLI structure of the 3′UTR was responsible for viral host-specific adaptation of the mosquito-derived strain in DEF and BHK-21 cells and regulated viral pathogenicity in 14-day-old mice, providing a new understanding of the functions of TMUV 3′UTR in viral host switching and the pathogenicity changes in mice. IMPORTANCE Mosquito-borne flaviviruses (MBFVs) constitute a large number of mosquito-transmitted viruses. The 3′-untranslated region (3′UTR) of MBFV has been suggested to be relevant to viral host-specific adaptation. However, the evolutionary strategies for host-specific fitness among MBFV are different, and the virulence-related structures within the 3′UTR are largely unknown. Here, using Tembusu virus (TMUV), an avian MBFV as models, we observed that the duck-derived SLI of the 3′UTR significantly enhanced the proliferation ability of mosquito-derived TMUV in baby hamster kidney (BHK-21) and duck embryonic fibroblast (DEF) cells, suggesting that the SLI structure was crucial for viral host-specific adaptation of mosquito-derived TMUVs in mammalian and avian cells. In addition, all SLI mutant viruses exhibited reduced viral pathogenicity in mice, indicating that SLI structure was a key factor for the pathogenicity in mice. This study provides a new insight into the functions of the MBFV 3′UTR in viral host switching and pathogenicity changes in mice.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has becoming globally public health threat. Recently studies were focus on SARS-CoV-2 RNA to design vaccine and drugs. It was demonstrated that virus RNA could play as sponge to host noncoding RNAs to regulate cellular processes. Bioinformatic research predicted a series of motif on SARS-CoV-2 genome where are targets of human miRNAs. In this study, we used dual-luciferase reporter assays to validate the interaction between 3’UTR of SARS-CoV-2 S (S-3’UTR) gene and bioinformatic predicted targeting miRNAs. The growth of 293T cells and HUVECs with overexpressed S-3’UTR was determined, while miRNAs and IL6, TNF-α levels were checked in this condition. Then, miR-296 and miR-602 mimic were introduced into 293T cells and HUVECs with overexpressed S-3’UTR, respectively, to reveal the underlying regulation mechanism. In results, we screened 19 miRNAs targeting the S-3’UTR, including miR-296 and miR-602. In 293T cell, S-3’UTR could inhibit 293T cell growth through down-regulation of miR-296. By reducing miR-602, S-3’UTR could induce HUVECs cell proliferation, alter the cell cycle, reduce apoptosis, and enhanced IL6 and TNF-αlevel. In conclusion, SARS-CoV-2 RNA could play as sponge of host miRNA to disturb cell growth and cytokine signaling. It suggests an important clue for designing COVID-19 drug and vaccine.

Breast cancer is the most common malignancy in women and has a high incidence rate and mortality. Abnormal regulation of gene expression plays an important role in breast cancer occurrence and development. RNA-binding proteins (RBPs) are one kind of the key regulators for gene expression. By interacting with RNA, RBPs are widely involved in RNA cutting, transport, editing, intracellular localization, and translation regulation. RBPs are important during breast cancer occurrence and progression by engaging in many aspects, like proliferation, migration, invasion, and stemness. Therefore, comprehensively understanding the role of RBPs in breast cancer progression can facilitate early diagnosis, timely treatment, and long-term survival and quality of life of breast cancer patients.

FRAXopathies are caused by the expansion of the CGG repeat in the 5′UTR of the FMR1 gene, which encodes the protein responsible for the synthesis of FMRP. This mutation leads to dramatic changes in FMRP expression at both the mRNA and protein levels. Evidence is emerging that changes in FMR1 mRNA expression can lead to the dysregulation of the miRNAs that target its 3′UTR. In the present work, B-lymphocyte cell lines obtained from patients with FRAXopathies were used, and a wide variety of FMR1 gene activities were observed, allowing the identification of the relationships between FMR1 dysregulation and miRNA activity. We studied the expression levels of eight miRNAs that target the FMR1 gene. To prove the interaction of the studied miRNAs with FMR1, a plasmid was constructed that possesses three primary structures: the miRNA gene, with expression driven by an inducible promoter; a constitutively expressed FusionRed reporter; and an eGFP reporter followed by the 3′UTR of the FMR1 gene. We evaluated changes in miRNA expression in response to alterations in FMR1 gene activity in a model cell line as well as interactions with some miRNAs with the FMR1 3′UTR.

Abstract Occurring in over 60% of human genes, alternative polyadenylation (APA) results in numerous transcripts with differing 3’ends, thus greatly expanding the diversity of mRNAs and of proteins derived from a single gene. As a key molecular mechanism, APA is involved in various gene regulation steps including mRNA maturation, mRNA stability, cellular RNA decay, and protein diversification. APA is frequently dysregulated in cancers leading to changes in oncogenes and tumor suppressor gene expressions. Recent studies have revealed various APA regulatory mechanisms that promote the development and progression of a number of human diseases, including cancer. Here, we provide an overview of four types of APA and their impacts on gene regulation. We focus particularly on the interaction of APA with microRNAs, RNA binding proteins and other related factors, the core pre-mRNA 3’end processing complex, and 3’UTR length change. We also describe next-generation sequencing methods and computational tools for use in poly(A) signal detection and APA repositories and databases. Finally, we summarize the current understanding of APA in cancer and provide our vision for future APA related research.

Definitive diagnosis to sudden cardiac death (SCD) is often challenging since the postmortem examination on SCD victims could hardly demonstrate an adequate cause of death. It is therefore important to uncover the inherited risk component to SCD. Signal transducer and activators of transcription 5 A (STAT5A) is a member of the STAT family and a transcription factor that is activated by many cell ligands and associated with various cardiovascular processes. In this study, we performed a systematic variant screening on the STAT5A to filter potential functional genetic variations. Based on the screening results, an insertion/deletion polymorphism (rs3833144) in 3’UTR of STAT5A was selected as the candidate variant. A total of 159 SCD cases and 668 SCD matched healthy controls was enrolled to perform a case-control study and evaluate the association between rs3833144 and SCD susceptibility in Chinese populations. Logistic regression analysis showed that the deletion allele of rs3833144 had significantly increased the SCD risk (odds ratio (OR) = 1.54; 95% confidence interval (CI) = 1.18–2.01; P = 0.000955). Further genotype-expression eQTL analysis showed that samples with deletion allele appeared to lower expression of STAT5A, and in silico prediction suggested the local 3 D structure changes of STAT5A mRNA caused by the variant. On the other hand, the bioinformatic analysis presented that promoters of RARA and PTGES3L-AARSD1 could interact with rs3833144, and eQTL analysis showed the higher expression of both genes in samples with deletion allele. Dual-luciferase activity assays also suggested the significant regulatory role of rs3833144 in gene transcription. Our current data thus suggested a possible involvement of rs3833144 to SCD predisposition in Chinese populations and rs3833144 with potential function roles may become a candidate marker for SCD diagnosis and prevention.

Mannose-binding lectin (MBL) and lectin complement pathway have become targets of increasing clinical interest. Many aspects of MBL have been recently explored, including the structural properties that allow it to distinguish self from non-self/altered-self structures. Experimental evidences have declared the additional 5′- and 3′-variants that in amalgamation with well-known secretor polymorphisms change MBL function and concentration. Moreover, the current review highlights the differential behavior of MBL on exposure with extra/intracellular pathogens and in autoimmune diseases, stressing the fact that “high MBL levels can increase diseases susceptibility,” a paradox that needs justification. Attributable to these discrepancies, no absolute level of MBL deficiency could be defined so far and thus must be interpreted for specific diseases through case–control population-specific designs. Overall, it is evident that further research is needed about MBL and the lectin pathway of complement. Particularly, the transformative role of MBL over evolution is of interest and its role with regard to pathogenesis of different diseases and potential therapeutic targets within the respective pathways should be further explored. Apart from this, it is necessary to adopt an extensive locus-wide methodology to apprehend the clinical significance of MBL2 polymorphisms in a variety of infectious diseases by the future studies.

Abstract Background Understanding gene regulation is important but difficult. Elucidating tissue-specific gene regulation mechanism is even more challenging and requires gene co-expression network assembled from protein–protein interaction, transcription factor and gene binding, and post-transcriptional regulation (e.g., miRNA targeting) information. The miRNA binding affinity could therefore be changed by SNP(s) located at the 3′ untranslated regions (3′UTR) of the target messenger RNA (mRNA) which miRNA(s) interacts with. Genome-wide association study (GWAS) has reported significant numbers of loci hosting SNPs associated with many traits. The goal of this study is to pinpoint GWAS functional variants located in 3′UTRs and elucidate if the genes harboring these variants along with their targeting miRNAs are associated with genetic traits relevant to certain tissues. Methods By applying MIGWAS, CoCoNet, ANNOVAR, and DAVID bioinformatics software and utilizing the gene expression database (e.g. GTEx data) to study GWAS summary statistics for 43 traits from 28 GWAS studies, we have identified a list of miRNAs and targeted genes harboring 3′UTR variants, which could contribute to trait-relevant tissue over miRNA-target gene network. Results Our result demonstrated that strong association between traits and tissues exists, and in particular, the Primary Biliary Cirrhosis (PBC) trait has the most significant p-value for all 180 tissues among all 43 traits used for this study. We reported SNPs located in 3′UTR regions of genes (SFMBT2, ZC3HAV1, and UGT3A1) targeted by miRNAs for PBC trait and its tissue association network. After employing Gene Ontology (GO) analysis for PBC trait, we have also identified a very important miRNA targeted gene over miRNA-target gene network, PFKL, which encodes the liver subunit of an enzyme. Conclusions The non-coding variants identified from GWAS studies are casually assumed to be not critical to translated protein product. However, 3′ untranslated regions (3′UTRs) of genes harbor variants can often change the binding affinity of targeting miRNAs playing important roles in protein translation degree. Our study has shown that GWAS variants could play important roles on miRNA-target gene networks by contributing the association between traits and tissues. Our analysis expands our knowledge on trait-relevant tissue network and paves way for future human disease studies.

RNA stability plays an important role in gene expression. Here, using 3′ end sequencing of newly made and pre-existing poly(A)+ RNAs, we compare transcript stability in multiple human cell lines, including HEK293T, HepG2, and SH-SY5Y. We show that while mRNA stability is generally conserved across the cell lines, specific transcripts having a high GC content and possibly more stable secondary RNA structures are relatively more stable in SH-SY5Y cells compared to the other 2 cell lines. These features also differentiate stability levels of alternative polyadenylation (APA) 3′UTR isoforms in a cell type-specific manner. Using differentiation of a neural stem cell line as a model, we show that mRNA stability difference could contribute to gene expression changes in neurogenesis and confirm the neuronal identity of SH-SY5Y cells at both gene expression and APA levels. In addition, compared to transcripts using 3′-most exon cleavage/polyadenylation sites (PASs), those using intronic PASs are generally less stable, especially when the PAS-containing intron is large and has a strong 5′ splice site, suggesting that intronic polyadenylation mostly plays a negative role in gene expression. Interestingly, the differential mRNA stability among APA isoforms appears to buffer PAS choice in these cell lines. Moreover, we found that several other poly(A)+ RNA species, including promoter-associated long noncoding RNAs and transcripts encoded by the mitochondrial genome, are more stable in SH-SY5Y cells than the other 2 cell lines, further highlighting distinct RNA metabolism in neuronal cells. Together, our results indicate that distinct RNA stability control in neuronal cells may contribute to the gene expression and APA programs that define their cell identity.

Yong-ping Liu, Xue Wu, Jing-hua Meng, Jun Yao, Bao-jie Wang School of Forensic Medicine, China Medical University, Shenyang 110122, People’s Republic of ChinaCorrespondence: Bao-jie WangSchool of Forensic Medicine, China Medical University, No. 77 Puhe Road, Shenbei New District, Shenyang 110122, People’s Republic of ChinaTel +8618842462979Email wangbaojie777@163.comPurpose: Abnormal expression of the NR1 subunit of the N-methyl-d-aspartate (NMDA) receptor may potentially increase the susceptibility to neuropsychiatric diseases. The purpose of this study was to investigate the functional sequence of the 3ʹUTR of the human GRIN1 gene, which encodes the GluN1 receptor to determine the effect on the expression of GluN1 receptor.Methods: We transferred seven recombinant pmirGLO recombinant vectors containing the 3ʹUTR truncated fragment of the GRIN1 gene into HEK-293, SK-N-SH, and U87 cell lines and compared the relative fluorescence intensity of adjacent length fragments. The TargetScan database was used to predict miRNAs. Then, miRNA mimics/inhibitors were co-transfected into the three cell lines with the 3ʹUTR of GRIN1 (pmirGLO - GRIN1), to investigate their influence on GRIN1 gene expression.Results: Compared with the pmirGLo-Basic vector, the relative fluorescence intensity of the complete GRIN1 gene 3ʹUTR recombinant sequence − 27 bp – +1284 bp (the next base of the stop codon is +1) was significantly decreased in all three cell lines. The relative fluorescence intensities were significantly different between − 27 bp – +294 bp and − 27 bp – +497 bp regions, and between − 27 bp – +708 bp and − 27 bp – +907 bp regions. According to the prediction of the TargetScan database and analysis, miR-212-5p, miR-324-3p and miR-326 may bind to +295 bp – +497 bp, while miR-491-5p may bind to +798 bp – +907 bp. After co-transfection of miRNA mimic/inhibitor or mimic/inhibitor NC with a recombinant vector in the 3ʹUTR region of GRIN1 gene, we found that has-miR-491-5p inhibited GRIN1 expression significantly in all three cell lines, while has-miR-326 inhibitor upregulated GRIN1 expression in HEK-293 and U87 cells.Conclusion: miR-491-5p may bind to the 3ʹUTR of the GRIN1 gene (+799 bp – +805 bp, the next base of the stop codon is +1) and down-regulate gene expression in HEK-293, SK-N-SH, and U87 cell lines, which implicates a potential role of miR-491-5p in central nervous system diseases.Keywords: human GRIN1 gene, 3ʹUTR, miRNA, GluN1 receptor, Parkinson’s disease

Abstract MicroRNA are small non-coding RNAs that post-transcriptionally regulate the expression levels of messenger RNAs. MicroRNA regulation activity depends on the recognition of binding sites located on mRNA molecules. ComiR is a web tool realized to predict the targets of a set of microRNAs, starting from their expression profile. ComiR was trained with the information regarding binding sites in the 3’utr region, by using a reliable dataset containing the targets of endogenously expressed microRNA in D. melanogaster S2 cells. This dataset was obtained by comparing the results from two different experimental approaches, i.e., inhibition, and immunoprecipitation of the AGO1 protein--a component of the microRNA induced silencing complex. In this work, we tested whether including coding region binding sites in ComiR algorithm improves the performance of the tool in predicting microRNA targets. We focused the analysis on the D. melanogaster species and updated the ComiR underlying database with the currently available releases of mRNA and microRNA sequences. As a result, we find that ComiR algorithm trained with the information related to the coding regions is more efficient in predicting the microRNA targets, with respect to the algorithm trained with 3’utr information. On the other hand, we show that 3’utr based predictions can be seen as complementary to the coding region based predictions, which suggests that both predictions, from 3’utr and coding regions, should be considered in comprehensive analysis. Furthermore, we observed that the lists of targets obtained by analyzing data from one experimental approach only, that is, inhibition or immunoprecipitation of AGO1, are not reliable enough to test the performance of our microRNA target prediction algorithm. Further analysis will be conducted to investigate the effectiveness of the tool with data from other species, provided that validated datasets, as obtained from the comparison of RISC proteins inhibition and immunoprecipitation experiments, will be available for the same samples. Finally, we propose to upgrade the existing ComiR web-tool by including the coding region based trained model, available together with the 3’utr based one.

Abstract Background 3′-Untranslated regions (3′UTRs) play crucial roles in mRNA metabolism, such as by controlling mRNA stability, translation efficiency, and localization. Intriguingly, in some genes the 3′UTR is longer than their coding regions, pointing to additional, unknown functions. Here, we describe a protein-coding function of 3′UTRs upon frameshift-inducing alternative splicing in more than 10% of human and mouse protein-coding genes. Results 3′UTR-encoded amino acid sequences show an enrichment of PxxP motifs and lead to interactome rewiring. Furthermore, an elevated proline content increases protein disorder and reduces protein stability, thus allowing splicing-controlled regulation of protein half-life. This could also act as a surveillance mechanism for erroneous skipping of penultimate exons resulting in transcripts that escape nonsense mediated decay. The impact of frameshift-inducing alternative splicing on disease development is emphasized by a retinitis pigmentosa-causing mutation leading to translation of a 3′UTR-encoded, proline-rich, destabilized frameshift-protein with altered protein-protein interactions. Conclusions We describe a widespread, evolutionarily conserved mechanism that enriches the mammalian proteome, controls protein expression and protein-protein interactions, and has important implications for the discovery of novel, potentially disease-relevant protein variants.

Abstract Reduced biotinidase activity is associated with a spectrum of deficiency ranging from total deficiency to heterozygous levels, a finding that is not always explained by the pathogenic variants observed in the BTD gene. The investigation of miRNAs, regulatory elements and variants in the 3’UTR region may present relevance in understanding the genotype-phenotype association. The aims of the study were to characterize the regulatory elements of the 3’UTR of the BTD gene and identify variants and miRNAs which may explain the discrepancies observed between genotype and biochemical phenotype. We evaluated 92 individuals with reduced biotinidase activity (level of heterozygotes = 33, borderline = 35, partial DB = 20 or total DB= 4) with previously determined BTD genotype. The 3’UTR of the BTD gene was Sanger sequenced. In silico analysis was performed to identify miRNAs and regulatory elements. No variants were found in the 3’UTR. We found 97 possible miRNAs associated with the BTD gene, 49 predicted miRNAs involved in the alanine, biotin, citrate and pyruvate metabolic pathways and 5 genes involved in biotin metabolism. Six AU-rich elements were found. Our data suggest variants in the 3'UTR of BTD do not explain the genotype-phenotype discrepancies found in Brazilian individuals with reduced biotinidase.

Duck hepatitis A virus 1 (DHAV-1) is a highly contagious etiological agent that causes acute hepatitis in young ducklings. MicroRNAs (miRNAs) play important regulatory roles in response to pathogens. However, the interplay between DHAV-1 infection and miRNAs remains ambiguous. We characterized and compared miRNA and mRNA expression profiles in duck embryo fibroblasts cells (DEFs) infected with DHAV-1. In total, 36 and 96 differentially expressed (DE) miRNAs, and 4110 and 2595 DE mRNAs, were identified at 12 and 24 h after infection. In particular, 126 and 275 miRNA–mRNA pairs with a negative correlation were chosen to construct an interaction network. Subsequently, we identified the functional annotation of DE mRNAs and target genes of DE miRNAs enriched in diverse Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, which may be important for virus resistance, cell proliferation, and metabolism. Moreover, upregulated miR-222a could negatively regulate DHAV-1 replication in DEFs and downregulate integrin subunit beta 3 (ITGB3) expression by targeting the 3′ untranslated region (3′UTR), indicating that miR-222a may modulate DHAV-1 replication via interaction with ITGB3. In conclusion, the results reveal changes of mRNAs and miRNAs during DHAV-1 infection and suggest miR-222a as an antiviral factor against DHAV-1.

Rapid responses to changes within subcellular compartments of highly polarized cells, such as neuron axons, depend on local translation and post-transcriptional regulation. The mechanism by which microRNAs (miRNAs) regulate this process is not fully understood. Here, using live cell imaging and RNA sequencing analysis, we demonstrated how miRNAs can differentially control hundreds of transcripts at the subcellular level. We demonstrated that the seed match length of the miRNA target-sequence regulates both mRNA stability and protein translation rates. While longer seed matches have an increased inhibitory effect, transcriptome analysis did not reveal differences in seed match length between axonal and somata mRNAs of motor neurons. However, mRNA variants with longer 3′UTR are enriched in axons and contain multiple repeats of specific miRNA target sequences. Finally, we demonstrated that the long 3′UTR mRNA variant of the motor protein Kif5b is enriched explicitly in motor neuron axons and contains multiple sequence repeats for binding miR-129-5p. This subsequently results in the differential post-transcriptional regulation of kif5b and its synthesis in axons. Thus, we suggest that the number of miRNA binding sites at the 3′UTR of the mRNA, rather than the miRNA seed match length, regulates the axonal transcriptome.

Multiple tumor exacerbations and treatment procedures, such as extracellular matrix remodeling, metabolic reprogramming, immunological evasion, and resistance to chemotherapy and radiotherapy, are influenced by intratumoral hypoxia. It is becoming increasingly clear how hypoxia interacts with the extracellular matrix and how this affects the growth of cancer. We analyzed the published sequencing results of hypoxia-stressed mouse kidney tumor cells and found that the expression of miR-29b was significantly downregulated. There are several sites that are complementary to the miR-29 seed sequence in the 3’ non-coding regions (3’UTRs) of various extracellular matrix-related genes, including collagen IV. We analyzed the sequences of the 3’UTRs of different subunits of collagen IV in different species and constructed the corresponding phylogenetic trees. We found that the 3’UTRs of Col4a1 and Col4a4 may have been subjected to particular evolutionary pressures. By cloning the 3’UTRs of collagen IV subunits into the psiCHECKTM-2 vector, we found that seven of the eight sites in the Col4a3–Col4a6 gene complementary to miR-29 were significantly repressed by miR-29a, b (except for the 7774–7781 of Col4a3 gene). The inhibitory efficiency of miR-29a, b on these seven sites was between 27% and 57%. The research on the regulation of miR-29 and extracellular matrix by hypoxia can provide a theoretical basis for tumor and fibrosis research and treatment.

Ficolin-2 is regarded as an important innate immunity factor endowed with both lectin (carbohydrate recognition) qualities and ability to induce complement activation. The aim of this study was to investigate the association of the FCN2 3’-untranslated region (3’UTR) polymorphisms with ficolin-2 expression and perinatal complications in preterm neonates. The sequencing analysis allowed us to identify six 3’UTR polymorphisms with minor allele frequency (MAF) >1%: rs4521835, rs73664188, rs11103564, rs11103565, rs6537958 and rs6537959. Except for rs4521835, all adhered to Hardy-Weinberg expectations. Moreover, rs6537958 and rs6537959 were shown to be in perfect linkage disequilibrium (LD) with nine other genetic polymorphisms: rs7040372, rs7046516, rs747422, rs7847431, rs6537957, rs6537960, rs6537962, rs11462298 and rs7860507 together stretched on a distance of 1242 bp and very high LD with rs11103565. The 3’UTR region was shown to bind nuclear extract proteins. The polymorphisms at rs4521835 and rs73664188 were found to influence serum ficolin-2 concentration significantly. All polymorphisms identified create (together with exon 8 polymorphism, rs7851696) two haplotype blocks. Among 49 diplotypes (D1-D49) created from rs7851696 (G>T), rs4521835 (T>G), rs73664188 (T>C), rs11103564 (T>C), rs11103565 (G>A) and rs6537959 (T>A), twenty two occurred with frequency >1%. Two diplotypes: D13 (GTTTGT/GGTCGT) and D10 (GTTTGT/GGTCGA), were significantly more frequent among preterm neonates with early onset of infection and pneumonia, compared with newborns with no infectious complications (OR 2.69 and 2.81, respectively; both p

Abstract Background Lung cancer is the second most common cancer with an extremely poor overall survival rate. Post-transcriptional regulation of gene expression play many important roles in human cancer, and one of the potential mechanisms underlying this is alternative mRNA maturation at its 3′ untranslated regions (3′-UTRs). Methods Cancer tissues and paired adjacent normal lung tissues from 26 patients diagnosed with non-small cell lung cancer (NSCLC) were analyzed by in vitro transcription-sequencing alternative polyadenylation sites (IVT-SAPAS). 41,773,101 reads in average were obtained from each paired sample. A potential regulation of Cleavage Stimulation Factor Subunit 2 (CSTF2) on 3′UTR length of genes was tested in H460 cells. Results 1439 (10.26%) genes showed up-regulated expression and 1364 (9.72%) genes showed down-regulated expression in lung cancer tissue versus normal lung tissue, and shorten 3′UTR in cancer tissue was detected in cancer tissues collected from 96.2% (25/26) patients, indicating lung cancer tend to have shortened 3′UTRs of these identified genes. KEGG analysis showed 1855 genes with shorten 3′UTR were enriched in mTOR signaling, ubiquitin mediated proteolysis and RNA degradation. Knocking down CSTF2 expression in H460 cells results in 3′UTR elongation of genes that was identified to be with shortened length in cancer tissues. Conclusion Alternative polyadenylation (APA) site-switching of 3′UTRs is prevalent in NSCLC, and CSTF2 may serve as an oncogene regulates the 3′UTR length of cancer related genes in NSCLC.

We have previously reported the presence of m6A in the AMV (Alfamovirus, Bromoviridae) genome. Interestingly, two of these putative m6A-sites are in hairpin (hp) structures in the 3’UTR of the viral RNA3. One site (2012AAACU2016) is in the loop of hpB, within the coat protein binding site 1 (CPB1), while the other (1900UGACC1904) is in the lower stem of hpE, a loop previously associated with AMV negative-strand RNA synthesis. In this work, we have performed in vivo experiments to assess the role of these two regions, containing the putative m6A-sites in the AMV cycle, by introducing compensatory point mutations to interfere with or abolish the m6A-tag of these sites. Our results suggest that the loop of hpB could be involved in viral replication/accumulation. Meanwhile, in the 1900UGACC1904 motif of the hpE, the maintenance of the adenosine residue and the lower stem hpE structure are necessary for in vivo plus-strand accumulation. These results extend our understanding of the requirements for hpE in the AMV infection cycle, indicating that both the residue identity and the base-pairing capacity in this structure are essential for viral accumulation.

RNA molecules are increasingly being identified as facilitating or impeding the interaction of proteins and nucleic acids, serving as so-called scaffolds or decoys. Long non-coding RNAs have been commonly implicated in such roles, particularly in the regulation of nuclear processes including chromosome topology, regulation of chromatin state and gene transcription, and assembly of nuclear biomolecular condensates such as paraspeckles. Recently, an increased awareness of cytoplasmic RNA scaffolds and decoys has begun to emerge, including the identification of non-coding regions of mRNAs that can also function in a scaffold-like manner to regulate interactions of nascently translated proteins. Collectively, cytoplasmic RNA scaffolds and decoys are now implicated in processes such as mRNA translation, decay, protein localization, protein degradation and assembly of cytoplasmic biomolecular condensates such as P-bodies. Here, we review examples of RNA scaffolds and decoys in both the nucleus and cytoplasm, illustrating common themes, the suitability of RNA to such roles, and future challenges in identifying and better understanding RNA scaffolding and decoy functions.

Bacterial messenger RNAs (mRNAs) are composed of 5′ and 3′ untranslated regions (UTRs) that flank the coding sequences (CDSs). In eukaryotes, 3′UTRs play key roles in post-transcriptional regulatory mechanisms. Shortening or deregulation of these regions is associated with diseases such as cancer and metabolic disorders. Comparatively, little is known about the functions of 3′UTRs in bacteria. Over the past few years, 3′UTRs have emerged as important players in the regulation of relevant bacterial processes such as virulence, iron metabolism, and biofilm formation. This MiniReview is an update for the different 3′UTR-mediated mechanisms that regulate gene expression in bacteria. Some of these include 3′UTRs that interact with the 5′UTR of the same transcript to modulate translation, 3′UTRs that are targeted by specific ribonucleases, RNA-binding proteins and small RNAs (sRNAs), and 3′UTRs that act as reservoirs of trans-acting sRNAs, among others. In addition, recent findings regarding a differential evolution of bacterial 3′UTRs and its impact in the species-specific expression of orthologous genes are also discussed.

Background: Tourette Syndrome (TS) is a neurodevelopmental disorder that presents with motor and vocal tics early in childhood. The aim of this study was to investigate genetic variants in the 3' untranslated region (3'UTR) of TS candidate genes with a putative link to microRNA (miRNA) mediated regulation or gene expression.Methods: We used an in silico approach to identify 32 variants in the 3'UTR of 18 candidate genes putatively changing the binding site for miRNAs. In a sample composed of TS cases and controls (n = 290), as well as TS family trios (n = 148), we performed transmission disequilibrium test (TDT) and meta-analysis.Results: We found positive association of rs3750486 in the LIM homeobox 6 (LHX6) gene (p = 0.021) and rs7795011 in the inner mitochondrial membrane peptidase subunit 2 (IMMP2L) gene (p = 0.029) with TS in our meta-analysis. The TDT showed an over-transmission of the A allele of rs1042201 in the arylacetamide deacetylase (AADAC) gene in TS patients (p = 0.029).Conclusion: This preliminary study provides further support for the involvement of LHX6, IMMP2L, and AADAC genes, as well as epigenetic mechanisms, such as altered miRNA mediated gene expression regulation in the etiology of TS.

BackgroundColorectal cancer (CRC) is one of the leading causes of cancer death worldwide. microRNAs (miRNAs) repress gene expression by binding to complementary sequences in the 3' untranslated region (3'UTR) of target mRNAs. Alternative polyadenylation (APA) are relevant to the variability of the 3'UTR of mRNA. However, the posttranscriptional dysregulation of miRNAs and APA in CRC are poorly understood.MethodIn this study, we conducted small RNA sequencing to identify differentially expressed miRNAs (DERs) and their target genes. Function analysis on DER-target genes can explain the regulation roles of miRNAs in CRC. The mutual regulation of miRNAs and APA was analyzed by combining miRNA data to 3'UTR alteration using 3' termini of polyadenylated RNAs sequencing (3T-seq) technique, and this was validated using TCGA gene expression data.ResultsOur results showed 64 significant differentially expressed miRNAs (DERs) in CRC patients. Their target genes were related to cell adhesion and transcription regulation and were prevailingly involved in the CRC-related pathway. Integrative analysis of the miRNA and APA profile revealed 16 DERs were correlated with 12 polyadenylation factors, and six of them were significantly differently expressed in CRC. We also found four DERs that lost binding sites due to APA and showed a positive correlation between the miRNA and gene expression.ConclusionOur study found that miRNAs regulated APA by modulating key polyadenylation factors, and several miRNAs lost their suppression on mRNA due to APA. Associating this with gene expression may provide some important clues for a deeper study of posttranscriptional cellular regulation and biomarker research in CRC. Our data provided the first evidence that the interaction between miRNAs and APA associated with gene expression could serve as biomarkers for CRC, suggesting that hsa-miR-133a-3p and MLEC, hsa-miR-145-5p and SET, hsa-miR-1-3p and PPIA, and hsa-miR-378d and YY1 might be novel and potential biomarkers in improving the diagnosis of CRC.

Abstract Background Whole transcriptome RNA variant analyses have shown that adenosine deaminases acting on RNA (ADAR) enzymes modify a large proportion of cellular RNAs, contributing to transcriptome diversity and cancer evolution. Despite the advances in the understanding of ADAR function in breast cancer, ADAR RNA editing functional consequences are not fully addressed. Results We characterized A to G(I) mRNA editing in 81 breast cell lines, showing increased editing at 3′UTR and exonic regions in breast cancer cells compared to immortalized non-malignant cell lines. In addition, tumors from the BRCA TCGA cohort show a 24% increase in editing over normal breast samples when looking at 571 well-characterized UTRs targeted by ADAR1. Basal-like subtype breast cancer patients with high level of ADAR1 mRNA expression shows a worse clinical outcome and increased editing in their 3′UTRs. Interestingly, editing was particularly increased in the 3′UTRs of ATM, GINS4 and POLH transcripts in tumors, which correlated with their mRNA expression. We confirmed the role of ADAR1 in this regulation using a shRNA in a breast cancer cell line (ZR-75-1). Conclusions Altogether, these results revealed a significant association between the mRNA editing in genes related to cancer-relevant pathways and clinical outcomes, suggesting an important role of ADAR1 expression and function in breast cancer.

Abstract Background Circular RNAs (circRNAs) are a new member of noncoding RNAs (ncRNAs) that have recently been described as key regulators of gene expression. Our previous study had identified the negative correlation between circHIPK3 and bladder cancer grade, invasion, as well as lymph node metastasis. However, the roles of circRNAs in cellular proliferation in bladder cancer remain largely unknown. Methods We had analyzed circRNA high-throughout sequencing from human tissues and determined bladder cancer related circRNA-3 (BCRC-3, GenBank: KU921434.1) as a new candidate circRNA derived from PSMD1 gene. The expression levels of circRNAs, mRNAs and miRNAs in human tissues and cells were detected by quantitative real-time PCR (qRT-PCR). The effects of BCRC-3 on cancer cells were explored by transfecting with plasmids in vitro and in vivo. RNA pull down assay, luciferase reporter assay and fluorescence in situ hybridization were applied to verify the interaction between BCRC-3 and microRNAs. Anticancer effects of methyl jasmonate (MJ) were measured by flow cytometry assay, western blot and qRT-PCR. Results BCRC-3 was lowly expressed in bladder cancer tissues and cell lines. Proliferation of BC cells was suppressed by ectopic expression of BCRC-3 in vitro and in vivo. Mechanistically, overexpression of BCRC-3 induced the expression of cyclin-dependent kinase inhibitor 1B (p27). Importantly, BCRC-3 could directly interact with miR-182-5p, and subsequently act as a miRNA sponge to promote the miR-182-5p-targeted 3’UTR activity of p27. Furthermore, MJ significantly increased the expression of BCRC-3, resulting in an obvious up-regulation of p27. Conclusions BCRC-3 functions as a tumor inhibitor to suppress BC cell proliferation through miR-182-5p/p27 axis, which would be a novel target for BC therapy.

The pluripotent transcription factor NANOG is essential for maintaining embryonic stem cells and driving tumorigenesis. We previously showed that PKC activity is involved in the regulation of NANOG expression. To explore the possible involvement of microRNAs in regulating the expression of key pluripotency factors, we performed a genome-wide analysis of microRNA expression in the embryonal carcinoma cell line NT2/D1 in the presence of the PKC activator, PMA. We found that MIR630 was significantly upregulated in PMA-treated cells. Experimentally, we showed that transfection of MIR630 mimic into embryonal carcinoma cell lines directly targeted the 3′UTR of OCT4, SOX2, and NANOG and markedly suppressed their expression. RNAhybrid and RNA22 algorithms were used to predict miRNA target sites in the NANOG 3′UTR, four possible target sites of MIR630 were identified. To examine the functional interaction between MIR630 and NANOG mRNA, the predicted MIR630 target sites in the NANOG 3′UTR were deleted and the activity of the reporters were compared. After targeted mutation of the predicted MIR630 target sites, the MIR630 mimic inhibited NANOG significantly less than the wild-type reporters. It is worth noting that mutation of a single putative binding site in the 3′UTR of NANOG did not completely abolish MIR630-mediated suppression, suggesting that MIR630 in the NANOG 3′UTR may have multiple binding sites and act together to maximally repress NANOG expression. Interestingly, MIR630 mimics significantly downregulated NANOG gene transcription. Exogenous expression of OCT4, SOX2, and NANOG lacking the 3′UTR almost completely rescued the reduced transcriptional activity of MIR630. MIR630 mediated the expression of differentiation markers in NT2/D1 cells, suggesting that MIR630 leads to the differentiation of NT2/D1 cell. Our findings show that MIR630 represses NANOG through transcriptional and post-transcriptional regulation, suggesting a direct link between core pluripotency factors and MIR630.

Abstract Natural antisense transcripts (NAT) and alternative polyadenylation (APA) of messenger RNA (mRNA) are important contributors of transcriptome complexity, each playing a critical role in multiple biological processes. However, whether they have crosstalk and function collaboratively is unclear. We discovered that APA enriched in human sense-antisense (S-AS) gene pairs, and finally focused on RNASEH2C-KAT5 S-AS pair for further study. In cis but not in trans over-expression of the antisense KAT5 gene promoted the usage of distal polyA (pA) site in sense gene RNASEH2C, which generated longer 3′ untranslated region (3′UTR) and produced less protein, accompanying with slowed cell growth. Mechanistically, elevated Pol II occupancy coupled with SRSF3 could explain the higher usage of distal pA site. Finally, NAT-mediated downregulation of sense gene’s protein level in RNASEH2C-KAT5 pair was specific for human rather than mouse, which lacks the distal pA site of RNASEH2C. We provided the first evidence to support that certain gene affected phenotype may not by the protein of its own, but by affecting the expression of its overlapped gene through APA, implying an unexpected view for understanding the link between genotype and phenotype.

Abstract Background Arabidopsis has 5 paralogs of the S-adenosylmethionine decarboxylase (SAMDC) gene. Neither their specific role in development nor the role of positive/purifying selection in genetic divergence of this gene family is known. While some data are available on organ-specific expression of AtSAMDC1, AtSAMDC2, AtSAMDC3 and AtSAMDC4, not much is known about their promoters including AtSAMDC5, which is believed to be non-functional. Results (1) Phylogenetic analysis of the five AtSAMDC genes shows similar divergence pattern for promoters and coding sequences (CDSs), whereas, genetic divergence of 5’UTRs and 3’UTRs was independent of the promoters and CDSs; (2) while AtSAMDC1 and AtSAMDC4 promoters exhibit high activity (constitutive in the former), promoter activities of AtSAMDC2, AtSAMDC3 and AtSAMDC5 are moderate to low in seedlings (depending upon translational or transcriptional fusions), and are localized mainly in the vascular tissues and reproductive organs in mature plants; (3) based on promoter activity, it appears that AtSAMDC5 is both transcriptionally and translationally active, but based on it’s coding sequence it seems to produce a non-functional protein; (4) though 5’-UTR based regulation of AtSAMDC expression through upstream open reading frames (uORFs) in the 5’UTR is well known, no such uORFs are present in AtSAMDC4 and AtSAMDC5; (5) the promoter regions of all five AtSAMDC genes contain common stress-responsive elements and hormone-responsive elements; (6) at the organ level, the activity of AtSAMDC enzyme does not correlate with the expression of specific AtSAMDC genes or with the contents of spermidine and spermine. Conclusions Differential roles of positive/purifying selection were observed in genetic divergence of the AtSAMDC gene family. All tissues express one or more AtSAMDC gene with significant redundancy, and concurrently, there is cell/tissue-specificity of gene expression, particularly in mature organs. This study provides valuable information about AtSAMDC promoters, which could be useful in future manipulation of crop plants for nutritive purposes, stress tolerance or bioenergy needs. The AtSAMDC1 core promoter might serve the need of a strong constitutive promoter, and its high expression in the gametophytic cells could be exploited, where strong male/female gametophyte-specific expression is desired; e.g. in transgenic modification of crop varieties.

Mesencephalic astrocyte derived neurotrophic factor (MANF) and cerebral dopamine neurotrophic factor (CDNF) are novel evolutionary conserved trophic factors, which exhibit cytoprotective activity via negative regulation of unfolded protein response (UPR) and inflammation. Despite multiple reports demonstrating detrimental effect of MANF/CDNF downregulation, little is known about the control of their expression. miRNAs—small non-coding RNAs—are important regulators of gene expression. Their dysregulation was demonstrated in multiple pathological processes and their ability to modulate levels of other neurotrophic factors, glial cell line-derived neurotrophic factor (GDNF) and brain-derived neurotrophic factor (BDNF), was previously reported. Here, for the first time we demonstrated direct regulation of MANF and CDNF by miRNAs. Using bioinformatic tools, reporter assay and analysis of endogenous MANF and CDNF, we identified that miR-144 controls MANF expression, and miR-134 and miR-141 downregulate CDNF levels. We also demonstrated that this effect is human-specific and is executed via predicted binding sites of corresponding miRNAs. Finally, we found that miR-382 suppressed hCDNF expression indirectly. In conclusion, we demonstrate for the first time direct regulation of MANF and CDNF expression by specific miRNAs, despite the fact their binding sites are not strongly evolutionary conserved. Furthermore, we demonstrate a functional effect of miR-144 mediated regulation of MANF on ER stress response markers. These findings emphasize that (1) prediction of miRNA targets based on evolutionary conservation may miss biologically meaningful regulatory pairs; and (2) interpretation of miRNA regulatory effects in animal models should be cautiously validated.

Complete 3′UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3′UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3′UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3′UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3′UTR motifs.

FUS is an RNA-binding protein associated with frontotemporal lobar degeneration (FTLD) and amyotrophic lateral sclerosis (ALS). Previous reports have demonstrated intrinsic roles of FUS in synaptic function. However, the mechanism underlying FUS’s regulation of synaptic morphology has remained unclear. We found that reduced mature spines after FUS depletion were associated with the internalization of PSD-95 within the dendritic shaft. Mass spectrometry of PSD-95-interacting proteins identified SynGAP, whose expression decreased after FUS depletion. Moreover, FUS and the ELAV-like proteins ELAVL4 and ELAVL1 control SynGAP mRNA stability in a 3′UTR length-dependent manner, resulting in the stable expression of the alternatively spliced SynGAP isoform α2. Finally, abnormal spine maturation and FTLD-like behavioral deficits in FUS-knockout mice were ameliorated by SynGAP α2. Our findings establish an important link between FUS and ELAVL proteins for mRNA stability control and indicate that this mechanism is crucial for the maintenance of synaptic morphology and cognitive function.