Structural and functional features of eukaryotic mRNA untranslated regions

Non-Viral RNA Therapies for Non-Small Cell Lung Cancer and Their Corresponding Clinical Trials

Ribonucleic acid (RNA)-based therapies represent a promising class of drugs for the treatment of non-small cell lung cancer (NSCLC) due to their ability to modulate gene expression. Therapies leveraging small interfering RNA (siRNA), messenger RNA (mRNA), microRNA (miRNA), and antisense oligonucleotides (ASOs) offer various advantages over conventional treatments, including the ability to target specific genetic mutations and the potential for personalized medicine approaches. However, the clinical translation of these therapeutics for the treatment of NSCLC faces challenges in delivery due to their immunogenicity, negative charge, and large size, which can be mitigated with delivery platforms. In this review, we provide a description of the pathophysiology of NSCLC and an overview of RNA-based therapeutics, specifically highlighting their potential application in the treatment of NSCLC. We discuss relevant classes of RNA and their therapeutic potential for NSCLC. We then discuss challenges in delivery and non-viral delivery strategies such as lipid- and polymer-based nanoparticles that have been developed to address these issues in preclinical models. Furthermore, we provide a summary table of clinical trials that leverage RNA therapies for NSCLC [which includes their National Clinical Trial (NCT) numbers] to highlight the current progress in NSCLC. We also discuss how these NSCLC therapies can be integrated with existing treatment modalities to enhance their efficacy and improve patient outcomes. Overall, we aim to highlight non-viral strategies that tackle RNA delivery challenges while showcasing RNA's potential as a next-generation therapy for NSCLC treatment.

Multiplexed assays of human disease-relevant mutations reveal UTR dinucleotide composition as a major determinant of RNA stability

Untranslated regions (UTRs) contain crucial regulatory elements for RNA stability, translation and localization, so their integrity is indispensable for gene expression. Approximately 3.7% of genetic variants associated with diseases occur in UTRs, yet a comprehensive understanding of UTR variant functions remains limited due to inefficient experimental and computational assessment methods. To systematically evaluate the effects of UTR variants on RNA stability, we established a massively parallel reporter assay on 6555 UTR variants reported in human disease databases. We examined the RNA degradation patterns mediated by the UTR library in two cell lines, and then applied LASSO regression to model the influential regulators of RNA stability. We found that UA dinucleotides and UA-rich motifs are the most prominent destabilizing element. Gain of UA dinucleotide outlined mutant UTRs with reduced stability. Studies on endogenous transcripts indicate that high UA-dinucleotide ratios in UTRs promote RNA degradation. Conversely, elevated GC content and protein binding on UA dinucleotides protect high-UA RNA from degradation. Further analysis reveals polarized roles of UA-dinucleotide-binding proteins in RNA protection and degradation. Furthermore, the UA-dinucleotide ratio of both UTRs is a common characteristic of genes in innate immune response pathways, implying a coordinated stability regulation through UTRs at the transcriptomic level. We also demonstrate that stability-altering UTRs are associated with changes in biobank-based health indices, underscoring the importance of precise UTR regulation for wellness. Our study highlights the importance of RNA stability regulation through UTR primary sequences, paving the way for further exploration of their implications in gene networks and precision medicine.

UTR-Insight: integrating deep learning for efficient 5' UTR discovery and design

The 5' UTR is critical for mRNA stability and translation efficiency in therapeutics. We developed UTR-Insight, a model integrating a pretrained language model with a CNN-Transformer architecture, explaining 89.1% of the mean ribosome load (MRL) variation in random 5' UTRs and 82.8% in endogenous 5' UTRs, surpassing existing models. Using UTR-Insight, we performed high-throughput in silico screening of hundreds of thousands of endogenous 5' UTRs from primates, mice, and viruses. The screened sequences increased protein expression by up to 319% compared to the human α-globin 5' UTR, and UTR-Insight-designed sequences achieved even greater expression levels than high-performing endogenous 5' UTRs.

From Sequence to System: Enhancing IVT mRNA Vaccine Effectiveness through Cutting-Edge Technologies

The COVID-19 pandemic has spotlighted the potential of in vitro transcribed (IVT) mRNA vaccines with their demonstrated efficacy, safety, cost-effectiveness, and rapid manufacturing. Numerous IVT mRNA vaccines are now under clinical trials for a range of targets, including infectious diseases, cancers, and genetic disorders. Despite their promise, IVT mRNA vaccines face hurdles such as limited expression levels, nonspecific targeting beyond the liver, rapid degradation, and unintended immune activation. Overcoming these challenges is crucial to harnessing the full therapeutic potential of IVT mRNA vaccines for global health advancement. This review provides a comprehensive overview of the latest research progress and optimization strategies for IVT mRNA molecules and delivery systems, including the application of artificial intelligence (AI) models and deep learning techniques for IVT mRNA structure optimization and mRNA delivery formulation design. We also discuss recent development of the delivery platforms, such as lipid nanoparticles (LNPs), polymers, and exosomes, which aim to address challenges related to IVT mRNA protection, cellular uptake, and targeted delivery. Lastly, we offer insights into future directions for improving IVT mRNA vaccines, with the hope to spur further progress in IVT mRNA vaccine research and development.

First Report of SNPs Detection in TMEM154 Gene in Sheep from Poland and Their Association with SRLV Infection Status

Small ruminant lentiviruses (SRLVs) infect sheep, causing a multiorganic disease called maedi-visna or ovine progressive pneumonia, which significantly affects the production and welfare of sheep, generating serious economic losses. Although not all infected animals develop fully symptomatic disease, they constantly spread the virus in the flock. Since the infection is incurable and no vaccine is available, another approach is necessary to control SRLV infections. In recent years, an alternative for culling infected animals has become the approach based on identifying genetic markers for selecting SRLV-resistant individuals. Recent reports revealed several candidates, including gene encoding transmembrane protein 154 (TMEM154). Several single nucleotide polymorphisms (SNPs) are found within this gene in sheep of different breeds, but only some can be considered as resistant markers. This study aimed to investigate the presence of single polymorphic sites in TMEM154 gene in sheep of selected Polish flocks and assess their association with the infection and proviral load in the context of susceptibility to SRLV infection. In total 107 sheep, representing three breeds, were screened for their SRLV infection status by serological and PCR testing. All these animals were also genotyped to characterize the presence of SNPs in TMEM154 gene and estimate their potential of being the SRLV-resistance marker. The frequency of identified alleles differed among breeds. Moreover, the positive association between TMEM154 genotype and SRLV status was found for E35K polymorphism and two polymorphic sites in 5'UTR in one of analyzed breed. However, when the relationship between SNPs and SRLV proviral load was analyzed, five had a strong association, considering the whole population of tested sheep. Presented data, for the first time, identified the presence of SNPs in TMEM154 gene in sheep housed in Polish flocks and suggested that selecting SRLV-resistant animals based on this analysis might be possible, but further validation in a larger group of sheep is required.

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

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

Advances in non-viral mRNA delivery to the spleen

Developing safe and effective delivery strategies for localizing messenger RNA (mRNA) payloads to the spleen is an important goal in the field of genetic medicine. Accomplishing this goal is challenging due to the instability, size, and charge of mRNA payloads. Here, we provide an analysis of non-viral delivery technologies that have been developed to deliver mRNA payloads to the spleen. Specifically, our review begins by outlining the unique anatomy and potential targets for mRNA delivery within the spleen. Next, we describe approaches in mRNA sequence engineering that can be used to improve mRNA delivery to the spleen. Then, we describe advances in non-viral carrier systems that can package and deliver mRNA payloads to the spleen, highlighting key advances in the literature in lipid nanoparticle (LNP) and polymer nanoparticle (PNP) technology platforms. Finally, we provide commentary and outlook on how splenic mRNA delivery may afford next-generation treatments for autoimmune disorders and cancers. In undertaking this approach, our goal with this review is to both establish a fundamental understanding of drug delivery challenges associated with localizing mRNA payloads to the spleen, while also broadly highlighting the potential to use these genetic medicines to treat disease.

Upstream open reading frames: new players in the landscape of cancer gene regulation

The translation of RNA by ribosomes represents a central biological process and one of the most dysregulated processes in cancer. While translation is traditionally thought to occur exclusively in the protein-coding regions of messenger RNAs (mRNAs), recent transcriptome-wide approaches have shown abundant ribosome activity across diverse stretches of RNA transcripts. The most common type of this kind of ribosome activity occurs in gene leader sequences, also known as 5' untranslated regions (UTRs) of the mRNA, that precede the main coding sequence. Translation of these upstream open reading frames (uORFs) is now known to occur in upwards of 25% of all protein-coding genes. With diverse functions from RNA regulation to microprotein generation, uORFs are rapidly igniting a new arena of cancer biology, where they are linked to cancer genetics, cancer signaling, and tumor-immune interactions. This review focuses on the contributions of uORFs and their associated 5'UTR sequences to cancer biology.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Differences in 5'untranslated regions highlight the importance of translational regulation of dosage sensitive genes

BACKGROUND

Untranslated regions (UTRs) are important mediators of post-transcriptional regulation. The length of UTRs and the composition of regulatory elements within them are known to vary substantially across genes, but little is known about the reasons for this variation in humans. Here, we set out to determine whether this variation, specifically in 5'UTRs, correlates with gene dosage sensitivity.

RESULTS

We investigate 5'UTR length, the number of alternative transcription start sites, the potential for alternative splicing, the number and type of upstream open reading frames (uORFs) and the propensity of 5'UTRs to form secondary structures. We explore how these elements vary by gene tolerance to loss-of-function (LoF; using the LOEUF metric), and in genes where changes in dosage are known to cause disease. We show that LOEUF correlates with 5'UTR length and complexity. Genes that are most intolerant to LoF have longer 5'UTRs, greater TSS diversity, and more upstream regulatory elements than their LoF tolerant counterparts. We show that these differences are evident in disease gene-sets, but not in recessive developmental disorder genes where LoF of a single allele is tolerated.

CONCLUSIONS

Our results confirm the importance of post-transcriptional regulation through 5'UTRs in tight regulation of mRNA and protein levels, particularly for genes where changes in dosage are deleterious and lead to disease. Finally, to support gene-based investigation we release a web-based browser tool, VuTR, that supports exploration of the composition of individual 5'UTRs and the impact of genetic variation within them.

A plant virus manipulates the long-winged morph of insect vectors

Wing dimorphism of insect vectors is a determining factor for viral long-distance dispersal and large-area epidemics. Although plant viruses affect the wing plasticity of insect vectors, the potential underlying molecular mechanisms have seldom been investigated. Here, we found that a planthopper-vectored rice virus, rice stripe virus (RSV), specifically induces a long-winged morph in male insects. The analysis of field populations demonstrated that the long-winged ratios of male insects are closely associated with RSV infection regardless of viral titers. A planthopper-specific and testis-highly expressed gene, Encounter, was fortuitously found to play a key role in the RSV-induced long-winged morph. Encounter resembles malate dehydrogenase in the sequence, but it does not have corresponding enzymatic activity. Encounter is upregulated to affect male wing dimorphism at early larval stages. Encounter is closely connected with the insulin/insulin-like growth factor signaling pathway as a downstream factor of Akt, of which the transcriptional level is activated in response to RSV infection, resulting in the elevated expression of Encounter. In addition, an RSV-derived small interfering RNA directly targets Encounter to enhance its expression. Our study reveals an unreported mechanism underlying the direct regulation by a plant virus of wing dimorphism in its insect vectors, providing the potential way for interrupting viral dispersal.

Exploring the influence of DNA methylation and single nucleotide polymorphisms of the Myostatin gene on growth traits in the hybrid grouper (Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male))

Investigations into the correlation between growth characteristics and DNA methylation levels, along with genetic variations, can provide fundamental insights to enhance growth performance in groupers. The Myostatin (mstn) gene plays a vital role in regulating skeletal muscle development and growth. This study scrutinized the DNA methylation levels of the mstn gene across hybrid groupers (E. fuscoguttatus (♀) × E. polyphekadion (♂)) and their parental species, to evaluate its impact on growth attributes in grouper fish. The nucleotide sequence of the mstn gene was directly sequenced in the hybrid grouper, exhibiting different growth performance to identify the single nucleotide polymorphisms (SNPs) of the mstn gene and explore their correlation with growth characteristics. The findings revealed no significant differences in global DNA methylation levels within muscle tissue among the hybrid grouper and parents. However, significant differences in DNA methylation sites were discovered between the hybrid grouper and E. polyphekadion at sites 824 and 1521 (located at exon 2 and intron 2, respectively), and between E. fuscoguttatus and E. polyphekadion at site 1521. These variations could potentially influence the mRNA expression of the mstn gene. The study also identified that SNP g.1003 T > C in exon 2 of the mstn gene was significantly associated with various growth traits including body weight, total length, body length, head length, caudal peduncle height, and body height (p < 0.01). Specimens with the TT genotype at site 1003 demonstrated superior growth performance compared to those with the TC genotype. Furthermore, microstructural analyses of muscle tissue showed that the average area and diameter of muscle fibers in TT genotype individuals were significantly greater than those in TC genotype individuals. Therefore, this research provides robust evidence linking the DNA methylation level and polymorphisms of the mstn gene with growth traits, which could be beneficial for grouper breeding programs.

Effects of Combinations of Untranslated-Region Sequences on Translation of mRNA

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.

Translational fidelity screens in mammalian cells reveal eIF3 and eIF4G2 as regulators of start codon selectivity

The translation initiation machinery and the ribosome orchestrate a highly dynamic scanning process to distinguish proper start codons from surrounding nucleotide sequences. Here, we performed genome-wide CRISPRi screens in human K562 cells to systematically identify modulators of the frequency of translation initiation at near-cognate start codons. We observed that depletion of any eIF3 core subunit promoted near-cognate start codon usage, though sensitivity thresholds of each subunit to sgRNA-mediated depletion varied considerably. Double sgRNA depletion experiments suggested that enhanced near-cognate usage in eIF3D depleted cells required canonical eIF4E cap-binding and was not driven by eIF2A or eIF2D-dependent leucine tRNA initiation. We further characterized the effects of eIF3D depletion and found that the N-terminus of eIF3D was strictly required for accurate start codon selection, whereas disruption of the cap-binding properties of eIF3D had no effect. Lastly, depletion of eIF3D activated TNFα signaling via NF-κB and the interferon gamma response. Similar transcriptional profiles were observed upon knockdown of eIF1A and eIF4G2, which also promoted near-cognate start codon usage, suggesting that enhanced near-cognate usage could potentially contribute to NF-κB activation. Our study thus provides new avenues to study the mechanisms and consequences of alternative start codon usage.

A highly contiguous genome assembly reveals sources of genomic novelty in the symbiotic fungus Rhizophagus irregularis

The root systems of most plant species are aided by the soil-foraging capacities of symbiotic arbuscular mycorrhizal (AM) fungi of the Glomeromycotina subphylum. Despite recent advances in our knowledge of the ecology and molecular biology of this mutualistic symbiosis, our understanding of the AM fungi genome biology is just emerging. Presented here is a close to T2T genome assembly of the model AM fungus Rhizophagus irregularis DAOM197198, achieved through Nanopore long-read DNA sequencing and Hi-C data. This haploid genome assembly of R. irregularis, alongside short- and long-read RNA-Sequencing data, was used to produce a comprehensive annotation catalog of gene models, repetitive elements, small RNA loci, and DNA cytosine methylome. A phylostratigraphic gene age inference framework revealed that the birth of genes associated with nutrient transporter activity and transmembrane ion transport systems predates the emergence of Glomeromycotina. While nutrient cycling in AM fungi relies on genes that existed in ancestor lineages, a burst of Glomeromycotina-restricted genetic innovation is also detected. Analysis of the chromosomal distribution of genetic and epigenetic features highlights evolutionarily young genomic regions that produce abundant small RNAs, suggesting active RNA-based monitoring of genetic sequences surrounding recently evolved genes. This chromosome-scale view of the genome of an AM fungus genome reveals previously unexplored sources of genomic novelty in an organism evolving under an obligate symbiotic life cycle.

The dark side of mRNA translation and the translation machinery in glioblastoma

Among the different types of cancer affecting the central nervous system (CNS), glioblastoma (GB) is classified by the World Health Organization (WHO) as the most common and aggressive CNS cancer in adults. GB incidence is more frequent among persons aged 45–55 years old. GB treatments are based on tumor resection, radiation, and chemotherapies. The current development of novel molecular biomarkers (MB) has led to a more accurate prediction of GB progression. Moreover, clinical, epidemiological, and experimental studies have established genetic variants consistently associated with the risk of suffering GB. However, despite the advances in these fields, the survival expectancy of GB patients is still shorter than 2 years. Thus, fundamental processes inducing tumor onset and progression remain to be elucidated. In recent years, mRNA translation has been in the spotlight, as its dysregulation is emerging as a key cause of GB. In particular, the initiation phase of translation is most involved in this process. Among the crucial events, the machinery performing this phase undergoes a reconfiguration under the hypoxic conditions in the tumor microenvironment. In addition, ribosomal proteins (RPs) have been reported to play translation-independent roles in GB development. This review focuses on the research elucidating the tight relationship between translation initiation, the translation machinery, and GB. We also summarize the state-of-the-art drugs targeting the translation machinery to improve patients’ survival. Overall, the recent advances in this field are shedding new light on the dark side of translation in GB.

Identification of TSGA10 and GGNBP2 splicing variants in 5' untranslated region with distinct expression profiles in brain tumor samples

INTRODUCTION

Brain tumors (BTs) are perceived as one of the most common malignancies among children. The specific regulation of each gene can play a critical role in cancer progression. The present study aimed to determine the transcripts of the TSGA10 and GGNBP2 genes, considering the alternative 5'UTR region, and investigating the expression of these different transcripts in BTs.

MATERIAL AND METHODS

Public data on brain tumor microarray datasets in GEO were analyzed with R software to evaluate the expression levels of TSGA10 and GGNBP2 genes (the Pheatmap package in R was also used to plot DEGs in a heat map). In addition, to validate our in-silico data analysis, RT-PCR was performed to determine the splicing variants of TSGA10 and GGNBP2 genes in testis and brain tumor samples. The expression levels of splice variants of these genes were analyzed in 30 brain tumor samples and two testicular tissue samples as a positive control.

RESULTS

In silico results show that the differential expression levels of TSGA10 and GGNBP2 were significant in the GEO datasets of BTs compared to normal samples (with adjusted p-value<0.05 and log fold change > 1). This study's experimental results showed that the TSGA10 gene produces four different transcripts with two distinct promoter regions and splicing exon 4. The relative mRNA expression of transcripts without exon 4 was higher than transcripts with exon 4 in BT samples (p-value<001). In GGNBP2, exon 2 in the 5'UTR region and exon 6 in the coding sequence were spliced. The expression analysis results showed that the relative mRNA expression of transcript variants without exon 2 was higher than other transcript variants with exon 2 in BT samples (p-value<001).

CONCLUSION

The decreased expression levels of transcripts with longer 5'UTR in BT samples than in testicular or low-grade brain tumor samples may decrease their translation efficiency. Therefore, decreased amounts of TSGA10 and GGNBP2 as potential tumor suppressor proteins, especially in high-grade brain tumors, may cause cancer development by angiogenesis and metastasis.

New insights into genome annotation in Podospora anserina through re-exploiting multiple RNA-seq data

BACKGROUND

Publicly available RNA-seq datasets are often underused although being helpful to improve functional annotation of eukaryotic genomes. This is especially true for filamentous fungi genomes which structure differs from most well annotated yeast genomes. Podospora anserina is a filamentous fungal model, which genome has been sequenced and annotated in 2008. Still, the current annotation lacks information about cis-regulatory elements, including promoters, transcription starting sites and terminators, which are instrumental to integrate epigenomic features into global gene regulation strategies.

RESULTS

Here we took advantage of 37 RNA-seq experiments that were obtained in contrasted developmental and physiological conditions, to complete the functional annotation of P. anserina genome. Out of the 10,800 previously annotated genes, 5'UTR and 3'UTR were defined for 7554, among which, 3328 showed differential transcriptional signal starts and/or transcriptional end sites. In addition, alternative splicing events were detected for 2350 genes, mostly due alternative 3'splice sites and 1732 novel transcriptionally active regions (nTARs) in unannotated regions were identified.

CONCLUSIONS

Our study provides a comprehensive genome-wide functional annotation of P. anserina genome, including chromatin features, cis-acting elements such as UTRs, alternative splicing events and transcription of non-coding regions. These new findings will likely improve our understanding of gene regulation strategies in compact genomes, such as those of filamentous fungi. Characterization of alternative transcripts and nTARs paves the way to the discovery of putative new genes, alternative peptides or regulatory non-coding RNAs.

SFPQ promotes RAS-mutant cancer cell growth by modulating 5'-UTR mediated translational control of CK1α

Oncogenic mutations in the RAS family of small GTPases are commonly found in human cancers and they promote tumorigenesis by altering gene expression networks. We previously demonstrated that Casein Kinase 1α (CK1α), a member of the CK1 family of serine/threonine kinases, is post-transcriptionally upregulated by oncogenic RAS signaling. Here, we report that the CK1α mRNA contains an exceptionally long 5'-untranslated region (UTR) harbouring several translational control elements, implicating its involvement in translational regulation. We demonstrate that the CK1α 5'-UTR functions as an IRES element in HCT-116 colon cancer cells to promote cap-independent translation. Using tobramycin-affinity RNA-pulldown assays coupled with identification via mass spectrometry, we identified several CK1α 5'-UTR-binding proteins, including SFPQ. We show that RNA interference targeting SFPQ reduced CK1α protein abundance and partially blocked RAS-mutant colon cancer cell growth. Importantly, transcript and protein levels of SFPQ and other CK1α 5'-UTR-associated RNA-binding proteins (RBPs) are found to be elevated in early stages of RAS-mutant cancers, including colorectal and lung adenocarcinoma. Taken together, our study uncovers a previously unappreciated role of RBPs in promoting RAS-mutant cancer cell growth and their potential to serve as promising biomarkers as well as tractable therapeutic targets in cancers driven by oncogenic RAS.

Identification, Characterization and Comparison of the Genome-Scale UTR Introns from Six Citrus Species

Ever since their discovery, introns within the coding sequence (CDS) of transcripts have been paid great attention. However, the introns located in the untranslated regions (UTRs) are often ignored. Here, we identified, characterized and compared the UTR introns (UIs) from six citrus species. Results showed that the average intron number of UTRs is greatly lower than that of CDSs. Among all six citrus species, the number and density of 5′UTR introns (5UIs) are higher than those of 3′UTR introns (3UIs). The UI densities varied greatly among different citrus species. There are 11 and 9 types of splice site (SS) pairs for the UIs of C. sinensis and C. medica, respectively. However, the UIs of the other four citrus species all own only three kinds of SS pairs. The ‘GT-AG’, accounting for more than 95% of both 5UIs and 3UIs SS pairs for all the six species, is the most popular type. Moreover, 81 5UIs and 26 3UIs were identified as common UIs among the six citrus species, and the transcripts containing these common UIs were mostly involved in gene expression or gene expression regulation. Our study revealed that the UIs’ length, abundance, density and SS pair types varied among different citrus species and that many UI-containing genes play important roles in gene expression regulation. Our findings have great implications for future citrus UI function research.

Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle

Fertility is a key driver of economic profitability in cattle production. A number of studies have identified genes associated with fertility using genome wide association studies and differential gene expression analysis; however, the genes themselves are poorly characterized in cattle. Here, we selected 13 genes from the literature which have previously been shown to have strong evidence for an association with fertility in Brahman cattle (Bos taurus indicus) or closely related breeds. We examine the expression variation of the 13 genes that are associated with cattle fertility using RNA-seq, CAGE-seq, and ISO-seq data from 11 different tissue samples from an adult Brahman cow and a Brahman fetus. Tissues examined include blood, liver, lung, kidney, muscle, spleen, ovary, and uterus from the cow and liver and lung from the fetus. The analysis revealed several novel isoforms, including seven from SERPINA7. The use of three expression characterization methodologies (5′ cap selected ISO-seq, CAGE-seq, and RNA-seq) allowed the identification of isoforms that varied in their length of 5′ and 3′ untranslated regions, variation otherwise undetectable (collapsed as degraded RNA) in generic isoform identification pipelines. The combinations of different sequencing technologies allowed us to overcome the limitations of relatively low sequence depth in the ISO-seq data. The lower sequence depth of the ISO-seq data was also reflected in the lack of observed expression of some genes that were observed in the CAGE-seq and RNA-seq data from the same tissue. We identified allele specific expression that was tissue-specific in AR, IGF1, SOX9, STAT3, and TAF9B. Finally, we characterized an exon of TAF9B as partially nested within the neighboring gene phosphoglycerate kinase 1. As this study only examined two animals, even more transcriptional variation may be present in a genetically diverse population. This analysis reveals the large amount of transcriptional variation within mammalian fertility genes and illuminates the fact that the transcriptional landscape cannot be fully characterized using a single technology alone.

MDM4 polymorphisms associated with the risk but not the prognosis of esophageal cancer in Cixian high-incidence region from northern China

AIM

The mouse double minute 4 (MDM4) may contribute to tumorgenesis by inhibiting p53 tumor suppressor activity. This study was designed to investigate whether MDM4 polymorphisms could affect susceptibility to esophageal squamous cell carcinoma (ESCC) and the survival of ESCC patients in a population from Cixian high-incidence region of northern China, which has not been explored.

METHODS

MDM4 rs1380576 and rs4245739 were genotyped by polymerase chain reaction-ligase detection reaction (PCR-LDR) in 568 ESCC patients and 578 controls.

RESULTS

Compared to rs1380576 C/C genotype, C/G genotype was associated with decreased risk of ESCC (odds ratio [OR] = 0.761, 95% confidence interval [CI] = 0.595-0.973). Compared to rs4245739 A/A genotype, A/C or C/C genotype was related to increased susceptibility to ESCC (OR = 1.551, 95% CI = 1.001-2.402). Individuals with GC haplotype had significantly higher risk of ESCC than those with CA or GA haplotype (OR = 1.598, 95% CI = 1.048-2.438). Neither rs1380576 nor rs4245739 influenced the survival of ESCC patients.

CONCLUSION

rs1380576 and rs4245739 may be used to predict susceptibility to ESCC for population in Cixian high-incidence region.

Identification of functional tumor necrosis factor-alpha promoter variants associated with Helicobacter pylori infection in the Sudanese population: Computational approach

BACKGROUND

Helicobacter pylori (H. pylori) is a ubiquitous bacterium that affects nearly half of the world's population with a high morbidity and mortality rate. Polymorphisms within the tumor necrosis factor-alpha (TNF-A) promoter region are considered a possible genetic basis for this disease.

AIM

To functionally characterize the genetic variations in the TNF-A 5'-region (-584 to +107) of Sudanese patients infected with H. pylori using in silico tools.

METHODS

An observational study was carried out in major public and private hospitals in Khartoum state. A total of 122 gastric biopsies were taken from patients who had been referred for endoscopy. Genomic DNA was extracted. Genotyping of the TNF-A-1030 polymorphism was performed using PCR with confronting two-pair primer to investigate its association with the susceptibility to H. pylori infection in the Sudanese population. Furthermore, Sanger sequencing was applied to detect single nucleotide polymorphisms in the 5'-region (-584 to +107) of TNF-A in H. pylori-infected patients. Bioinformatics analyses were used to predict whether these mutations would alter transcription factor binding sites or composite regulatory elements in this region. A comparative profiling analysis was conducted in 11 species using the ECR browser and multiple-sequence local alignment and visualization search engine to investigate the possible conservation. Also, a multivariate logistic regression model was constructed to estimate odds ratios and their 95% confidence intervals for the association between TNF-A-1030, sociodemographic characteristics and H. pylori infection. Differences were statistically significant if P < 0.05. Statistical analyses were performed using Stata version 11 software.

RESULTS

A total of seven single nucleotide polymorphisms were observed in the TNF-A 5'-region of Sudanese patients infected with H. pylori. Only one of them (T > A, -76) was located at the in silico-predicted promoter region (-146 to +10), and it was predicted to alter transcription factor binding sites and composite regulatory elements. A novel mutation (A > T, +27) was detected in the 5' untranslated region, and it could affect the post-transcriptional regulatory pathways. Genotyping of TNF-A-1030 showed a lack of significant association between -1030T and susceptibility to H. pylori and gastric cancer in the studied population (P = 0.1756) and (P = 0.8116), respectively. However, a significant association was detected between T/C genotype and H. pylori infection (39.34% vs 19.67%, odds ratio = 2.69, 95% confidence interval: 1.17-6.17, P = 0.020). Mammalian conservation was observed for the (-146 to +10) region in chimpanzee (99.4%), rhesus monkey (95.6%), cow (91.8%), domesticated dog (89.3%), mouse (84.3%), rat (82.4%) and opossum (78%).

CONCLUSION

Computational analysis was a valuable method for understanding TNF-A gene expression patterns and guiding further in vitro and in vivo experimental validation.

Genome-Scale Computational Identification and Characterization of UTR Introns in Atalantia buxifolia

Accumulated evidence has shown that CDS introns (CIs) play important roles in regulating gene expression. However, research on UTR introns (UIs) is limited. In this study, UIs (including 5′UTR and 3′UTR introns (5UIs and 3UIs)) were identified from the Atalantia buxifolia genome. The length and nucleotide distribution characteristics of both 5UIs and 3UIs and the distributions of cis-acting elements and transcription factor binding sites (TFBSs) in 5UIs were investigated. Moreover, PageMan enrichment analysis was applied to show the possible roles of transcripts containing UIs (UI-Ts). In total, 1077 5UIs and 866 3UIs were identified from 897 5UI-Ts and 670 3UI-Ts, respectively. Among them, 765 (85.28%) 5UI-Ts and 527 (78.66%) 3UI-Ts contained only one UI, and 94 (6.38%) UI-Ts contained both 5UI and 3UI. The UI density was lower than that of CDS introns, but their mean and median intron sizes were ~2 times those of the CDS introns. The A. buxifolia 5UIs were rich in gene-expression-enhancement-related elements and contained many TFBSs for BBR-BPC, MIKC_MADS, AP2 and Dof TFs, indicating that 5UIs play a role in regulating or enhancing the expression of downstream genes. Enrichment analysis revealed that UI-Ts involved in ‘not assigned’ and ‘RNA’ pathways were significantly enriched. Noteworthily, 119 (85.61%) of the 3UI-Ts were genes encoding pentatricopeptide (PPR) repeat-containing proteins. These results will be helpful for the future study of the regulatory roles of UIs in A. buxifolia.

The Landscape of Alternative Splicing Regulating Potassium Use Efficiency in Nicotiana tabacum

Alternative splicing (AS) occurs extensively in eukaryotes as an essential mechanism for regulating transcriptome complexity and diversity, but the AS landscape regulating potassium (K) use efficiency in plants is unclear. In this study, we performed high-throughput transcriptome sequencing of roots and shoots from allopolyploid Nicotiana tabacum under K⁺ deficiency. Preliminary physiological analysis showed that root system architecture was dramatically changed due to potassium deficiency and that IAA content was significantly reduced in root and shoot. AS analysis showed that a total of 28,179 genes exhibited 54,457 AS events, and 1,510 and 1,732 differentially alternatively spliced (DAS) events were identified in shoots and roots under low K⁺ stress. Nevertheless, only 120 DAS events occurred in both shoots and roots, implying that most DAS events were tissue-specific. Both in shoot and the root, the proportion of DAS genes in differentially expressed (DE) genes equaled that in non-DE genes, which indicated that AS might play a unique regulatory role in response to low potassium. Gene ontology analysis further indicated that transcription regulation and AS modulation worked independently in response to low K⁺ stress in tobacco, as their target biological processes were different. Totally 45 DAS transcription factors (TFs) were found, which were involved in 18 TF families. Five Auxin response factor (ARF) TFs were significantly DAS in root, suggesting that response to auxin was probably subject to AS regulation in the tobacco root. Our study shows that AS variation occurs extensively and has a particular regulatory mechanism under K⁺ deficiency in tobacco. The study also links changes in root system architecture with the changes in AS of ARF TFs, which implied the functional significance of these AS events for root growth and architecture.

Determination of the Affinity of Eukaryotic DDX3 RNA Helicase to the Characteristic Elements of mRNA Secondary Structure

DDX3 RNA helicase is involved in many processes of RNA metabolism in eukaryotic cells. Many studies of DDX3 have shown that it is also involved in the translation initiation process, both cap-dependent and IRES-dependent. However, the specificity of the secondary structure of mRNA 5'-UTRs, which require DDX3 RNA helicase for effective translation, has not yet been determined. We performed a bioinformatic analysis of the 5'-UTR secondary structures in the pool of DDX3-dependent mRNAs in silico and predicted that the length of 5'-UTRs for such mRNAs is less than the average for the genome and that there are also characteristic hairpin structures in the region of the first 50 nucleotides from the 5'-end of the mRNA.

The Evolution of G-quadruplex Structure in mRNA Untranslated Region

The RNA G-quadruplex (rG4) is a kind of non-canonical high-order secondary structure with important biological functions and is enriched in untranslated regions (UTRs) of protein-coding genes. However, how rG4 structures evolve is largely unknown. Here, we systematically investigated the evolution of RNA sequences around UTR rG4 structures in 5 eukaryotic organisms. We found universal selection on UTR sequences, which facilitated rG4 formation in all the organisms that we analyzed. While G-rich sequences were preferred in the rG4 structural region, C-rich sequences were selectively not preferred. The selective pressure acting on rG4 structures in the UTRs of genes with higher G content was significantly smaller. Furthermore, we found that rG4 structures experienced smaller evolutionary selection near the translation initiation region in the 5' UTR, near the polyadenylation signals in the 3' UTR, and in regions flanking the miRNA targets in the 3' UTR. These results suggest universal selection for rG4 formation in the UTRs of eukaryotic genomes and the selection may be related to the biological functions of rG4s.

Sequence Analysis and Ontogenetic Expression Patterns of Cone Opsin Genes in the Bluefin Killifish (Lucania goodei)

Sensory systems allow for the transfer of environmental stimuli into internal cues that can alter physiology and behavior. Many studies of visual systems focus on opsins to compare spectral sensitivity among individuals, populations, and species living in different lighting environments. This requires an understanding of the cone opsins, which can be numerous. The bluefin killifish is a good model for studying the interaction between environments and visual systems as they are found in both clear springs and tannin-stained swamps. We conducted a genome-wide screening and demonstrated that the bluefin killifish has 9 cone opsins: 1 SWS1 (354 nm), 2 SWS2 (SWS2B: 359 nm, SWS2A: 448 nm), 2 RH2 (RH2-2: 476 nm, RH2-1: 537 nm), and 4 LWS (LWS-1: 569 nm, LWS-2: 524 nm, LWS-3: 569 nm, LWS-R: 560 or 569 nm). These 9 cone opsins were located on 4 scaffolds. One scaffold contained the 2 SWS2 and 3 of the 4 LWS opsins in the same syntenic order as found in other cyprinodontoid fishes. We also compared opsin expression in larval and adult killifish under clear water conditions, which mimic springs. Two of the newly discovered opsins (LWS-2 and LWS-3) were expressed at low levels (<0.2%). Whether these opsins make meaningful contributions to visual perception in other contexts (i.e., swamp conditions) is unclear. In contrast, there was an ontogenetic change from using LWS-R to LWS-1 opsin. Bluefin killifish adults may be slightly more sensitive to longer wavelengths, which might be related to sexual selection and/or foraging preferences.

Dissection of a rice OsMac1 mRNA 5' UTR to uncover regulatory elements that are responsible for its efficient translation

The untranslated regions (UTRs) of mRNAs are involved in many posttranscriptional regulatory pathways. The rice OsMac1 mRNA has three splicing variants of the 5' UTR (UTRa, UTRb, and UTRc), which include a CU-rich region and three upstream open reading frames (uORFs). UTRc contains an additional 38-nt sequence, termed sp38, which acts as a strong translational enhancer of the downstream ORF; reporter analysis revealed translational efficiencies >15-fold higher with UTRc than with the other splice variants. Mutation analysis of UTRc demonstrated that an optimal sequence length of sp38, rather than its nucleotide sequence is essential for UTRc to promote efficient translation. In addition, the 5' 100 nucleotides of CU-rich region contribute to UTRc translational enhancement. Strikingly, three uORFs did not reveal their inhibitory potential within the full-length leader, whereas deletion of the 5' leader fragment preceding the leader region with uORFs nearly abolished translation. Computational prediction of UTRc structural motifs revealed stem-loop structures, termed SL1-SL4, and two regions, A and B, involved in putative intramolecular interactions. Our data suggest that SL4 binding to Region-A and base pairing between Region-B and the UTRc 3'end are critically required for translational enhancement. Since UTRc is not capable of internal initiation, we presume that the three-dimensional leader structures can allow translation of the leader downstream ORF, likely allowing the bypass of uORFs.

Vps34 and TOR Kinases Coordinate HAC1 mRNA Translation in the Presence or Absence of Ire1-Dependent Splicing

In the budding yeast Saccharomyces cerevisiae, an mRNA, called HAC1, exists in a translationally repressed form in the cytoplasm. Under conditions of cellular stress, such as when unfolded proteins accumulate inside the endoplasmic reticulum (ER), an RNase Ire1 removes an intervening sequence (intron) from the HAC1 mRNA by nonconventional cytosolic splicing. Removal of the intron results in translational derepression of HAC1 mRNA and production of a transcription factor that activates expression of many enzymes and chaperones to increase the protein-folding capacity of the cell. Here, we show that Ire1-mediated RNA cleavage requires Watson-Crick base pairs in two RNA hairpins, which are located at the HAC1 mRNA exon-intron junctions. Then, we show that the translational derepression of HAC1 mRNA can occur independent of cytosolic splicing. These results are obtained from HAC1 variants that translated an active Hac1 protein from the unspliced mRNA. Additionally, we show that the phosphatidylinositol-3-kinase Vps34 and the nutrient-sensing kinases TOR and GCN2 are key regulators of HAC1 mRNA translation and consequently the ER stress responses. Collectively, our data suggest that the cytosolic splicing and the translational derepression of HAC1 mRNA are coordinated by unique and parallel networks of signaling pathways.

Translational Control of Serrate Expression in Drosophila Cells

BACKGROUND/AIM

The DSL proteins, Serrate and Delta, which act as Notch receptor ligands, mediate signalling between adjacent cells, when a ligand-expressing cell binds to Notch on an adjacent receiving cell. Notch is ubiquitously expressed and DSL protein mis-expression can have devastating developmental consequences. Although transcriptional regulation of Delta and Serrate has been amply documented, we examined whether they are also regulated at the level of translation.

MATERIALS AND METHODS

We generated a series of deletions to investigate the initiation codon usage for Serrate using Drosophila S2 cells.

RESULTS

Serrate mRNA contains three putative ATG initiation codons spanning a 60-codon region upstream of its signal peptide; we found that each one can act as an initiation codon, however, with a different translational efficiency.

CONCLUSION

Serrate expression is strictly regulated at the translational level.

Characterization of Rheumatoid Arthritis Risk-Associated SNPs and Identification of Novel Therapeutic Sites Using an In-Silico Approach

Simple Summary

Rheumatoid arthritis (RA) is a complex disease resulting from multiple genetic and environmental pathogenic factors. The genetic factors include single-nucleotide polymorphisms (SNPs), which have been reported to be associated with RA, but their specific role in the pathogenesis of RA remains unexplained. This study explains the potential role of RA risk-associated SNPs in its pathogenesis in order to provide a basis for understanding the genetic complexity of RA. Several roles of these SNPs are described in this study, and may also aid in the design of a therapeutic strategy for RA. Furthermore, novel potential therapeutic sites have also been researched, resulting in the identification of three novel therapeutic targets. The therapeutic strategies for the treatment of RA include inflammatory pathway-targeting drugs, which alleviate inflammation in joints. There is always a need for novel therapeutic targets that can play a role in alleviating inflammation in autoimmune diseases including RA. Therefore, these novel therapeutic sites are very important, and further experimental studies are required.

Abstract

Single-nucleotide polymorphisms (SNPs) are reported to be associated with many diseases, including autoimmune diseases. In rheumatoid arthritis (RA), about 152 SNPs are reported to account for ~15% of its heritability. These SNPs may result in the alteration of gene expression and may also affect the stability of mRNA, resulting in diseased protein. Therefore, in order to predict the underlying mechanism of these SNPs and identify novel therapeutic sites for the treatment of RA, several bioinformatics tools were used. The damaging effect of 23 non-synonymous SNPs on proteins using different tools suggested four SNPs, including rs2476601 in PTPN22, rs5029941 and rs2230926 in TNFAIP3, and rs34536443 in TYK2, to be the most damaging. In total, 42 of 76 RA-associated intronic SNPs were predicted to create or abolish potential splice sites. Moreover, the analysis of 11 RA-associated UTR SNPs indicated that only one SNP, rs1128334, located in 3′UTR of ETS1, caused functional pattern changes in BRD-BOX. For the identification of novel therapeutics sites to treat RA, extensive gene–gene interaction network interactive pathways were established, with the identification of 13 potential target sites for the development of RA drugs, including three novel target genes. The anticipated effect of these findings on RA pathogenesis may be further validated in both in vivo and in vitro studies.

Association of Vitamin D Metabolism Gene Polymorphisms with Autoimmunity: Evidence in Population Genetic Studies

A high prevalence of vitamin D (calcidiol) serum deficiency has been described in several autoimmune diseases, including multiple sclerosis (MS), rheumatoid arthritis (AR), and systemic lupus erythematosus (SLE). Vitamin D is a potent immunonutrient that through its main metabolite calcitriol, regulates the immunomodulation of macrophages, dendritic cells, T and B lymphocytes, which express the vitamin D receptor (VDR), and they produce and respond to calcitriol. Genetic association studies have shown that up to 65% of vitamin D serum variance may be explained due to genetic background. The 90% of genetic variability takes place in the form of single nucleotide polymorphisms (SNPs), and SNPs in genes related to vitamin D metabolism have been linked to influence the calcidiol serum levels, such as in the vitamin D binding protein (VDBP; rs2282679 GC), 25-hydroxylase (rs10751657 CYP2R1), 1α-hydroxylase (rs10877012, CYP27B1) and the vitamin D receptor (FokI (rs2228570), BsmI (rs1544410), ApaI (rs7975232), and TaqI (rs731236) VDR). Therefore, the aim of this comprehensive literature review was to discuss the current findings of functional SNPs in GC, CYP2R1, CYP27B1, and VDR associated to genetic risk, and the most common clinical features of MS, RA, and SLE.

Association of resistin (rs3745367) and urotensin II (rs228648 and rs2890565) gene polymorphisms with risk of type 2 diabetes mellitus in Indian population

Insulin resistance may become the most powerful predictor of future development of type 2 diabetes mellitus (T2DM) and a therapeutic target for the treatment of the same. Both Resistin, an adipose derived peptide hormone and Urotensin II a potent vasoconstrictor, are reported to be involved in the development of insulin resistance and T2DM but the results remain contradictory. Therefore, investigations were carried out to study the association of T2DM and single nucleotide polymorphism (SNP) in Resistin (RETN) gene at rs3745367 (+ 299 G > A) and Urotensin II (UTS2) gene at rs228648 (+ 143 G > A) and rs2890565 (+ 3836 C > T) in a North Indian population. Method: The present case-control study, conducted from August 2017 to July 2020, involved 168 T2DM patients and 102 healthy controls. SNPs rs3745367, rs228648 and rs2890565 were amplified from genomic DNA in the studied samples by polymerase chain reaction (PCR) using specific primers. The amplified products were genotyped by restriction fragment length polymorphism (RFLP) using particular restriction endonucleases. Clinical parameters viz. glycosylated haemoglobin (HbA1c), fasting blood glucose (FBG), high density lipoprotein cholesterol (HDL-C), triglycerides (TG), total cholesterol (CHL) and fasting insulin were determined by enzymatic methods. Result and conclusion: A statistically significant association between T2DM and RETN gene at SNP rs3745367 (p = 0.001) and UTS2 gene at SNP rs2890565 (p = 0.001) was observed. In RETN gene SNP rs3745367, insulin and homeostasis model assessment of insulin resistance (HOMA-IR) were found to be higher in GA + AA combined genotype than in GG genotype for T2DM subjects. Regression analysis revealed that SNP rs2890565 and HOMA-IR were independently associated with the risk of development of T2DM when three SNPs were taken as independent variable adjusted for clinical variables. Among four haplotypes, A/T was found associated with increased risk of T2DM as determined for rs228648 and rs2890565 of UTS2 gene. It can be concluded from these results that polymorphism at rs3745367 of RETN gene and at rs2890565 of UTS2 gene are associated with risk of T2DM in North Indian population.

5' Untranslated Region Elements Show High Abundance and Great Variability in Homologous ABCA Subfamily Genes

The 12 members of the ABCA subfamily in humans are known for their ability to transport cholesterol and its derivatives, vitamins, and xenobiotics across biomembranes. Several ABCA genes are causatively linked to inborn diseases, and the role in cancer progression and metastasis is studied intensively. The regulation of translation initiation is implicated as the major mechanism in the processes of post-transcriptional modifications determining final protein levels. In the current bioinformatics study, we mapped the features of the 5' untranslated regions (5'UTR) known to have the potential to regulate translation, such as the length of 5'UTRs, upstream ATG codons, upstream open-reading frames, introns, RNA G-quadruplex-forming sequences, stem loops, and Kozak consensus motifs, in the DNA sequences of all members of the subfamily. Subsequently, the conservation of the features, correlations among them, ribosome profiling data as well as protein levels in normal human tissues were examined. The 5'UTRs of ABCA genes contain above-average numbers of upstream ATGs, open-reading frames and introns, as well as conserved ones, and these elements probably play important biological roles in this subfamily, unlike RG4s. Although we found significant correlations among the features, we did not find any correlation between the numbers of 5'UTR features and protein tissue distribution and expression scores. We showed the existence of single nucleotide variants in relation to the 5'UTR features experimentally in a cohort of 105 breast cancer patients. 5'UTR features presumably prepare a complex playground, in which the other elements such as RNA binding proteins and non-coding RNAs play the major role in the fine-tuning of protein expression.

How Many Messenger RNAs Can Be Translated by the START Mechanism?

Translation initiation is a key step in the protein synthesis stage of the gene expression pathway of all living cells. In this important process, ribosomes have to accurately find the AUG start codon in order to ensure the integrity of the proteome. “Structure Assisted RNA Translation”, or “START”, has been proposed to use stable secondary structures located in the coding sequence to augment start site selection by steric hindrance of the progression of pre-initiation complex on messenger RNA. This implies that such structures have to be located downstream and at on optimal distance from the AUG start codon (i.e., downstream nucleotide +16). In order to assess the importance of the START mechanism in the overall mRNA translation process, we developed a bioinformatic tool to screen coding sequences for such stable structures in a 50 nucleotide-long window spanning the nucleotides from +16 to +65. We screened eight bacterial genomes and six eukaryotic genomes. We found stable structures in 0.6–2.5% of eukaryotic coding sequences. Among these, approximately half of them were structures predicted to form G-quadruplex structures. In humans, we selected 747 structures. In bacteria, the coding sequences from Gram-positive bacteria contained 2.6–4.2% stable structures, whereas the structures were less abundant in Gram-negative bacteria (0.2–2.7%). In contrast to eukaryotes, putative G-quadruplex structures are very rare in the coding sequence of bacteria. Altogether, our study reveals that the START mechanism seems to be an ancient strategy to facilitate the start codon recognition that is used in different kingdoms of life.

Selection Shapes Synonymous Stop Codon Use in Mammals

Phylogenetic models of the evolution of protein-coding sequences can provide insights into the selection pressures that have shaped them. In the application of these models synonymous nucleotide substitutions, which do not alter the encoded amino acid, are often assumed to have limited functional consequences and used as a proxy for the neutral rate of evolution. The ratio of nonsynonymous to synonymous substitution rates is then used to categorize the selective regime that applies to the protein (e.g., purifying selection, neutral evolution, diversifying selection). Here, we extend the Muse and Gaut model of codon evolution to explore the extent of purifying selection acting on substitutions between synonymous stop codons. Using a large collection of coding sequence alignments, we estimate that a high proportion (approximately 57%) of mammalian genes are affected by selection acting on stop codon preference. This proportion varies substantially by codon, with UGA stop codons far more likely to be conserved. Genes with evidence of selection acting on synonymous stop codons have distinctive characteristics, compared to unconserved genes with the same stop codon, including longer [Formula: see text] untranslated regions (UTRs) and shorter mRNA half-life. The coding regions of these genes are also much more likely to be under strong purifying selection pressure. Our results suggest that the preference for UGA stop codons found in many multicellular eukaryotes is selective rather than mutational in origin.

Optimization of 5' Untranslated Region of Modified mRNA for Use in Cardiac or Hepatic Ischemic Injury

Modified mRNA (modRNA) is a gene-delivery platform for transiently introducing a single gene or several genes of interest to different cell types and tissues. modRNA is considered to be a safe vector for gene transfer, as it negligibly activates the innate immune system and does not compromise the genome integrity. The use of modRNA in basic and translational science is rising, due to the clinical potential of modRNA. We are currently using modRNA to induce cardiac regeneration post-ischemic injury. Major obstacles in using modRNA for cardiac ischemic disease include the need for the direct and single administration of modRNA to the heart and the inefficient translation of modRNA due to its short half-life. Modulation of the 5' untranslated region (5' UTR) to enhance translation efficiency in ischemic cardiac disease has great value, as it can reduce the amount of modRNA needed per delivery and will achieve higher and longer protein production post-single delivery. Here, we identified that 5' UTR, from the fatty acid metabolism gene carboxylesterase 1D (Ces1d), enhanced the translation of firefly luciferase (Luc) modRNA by 2-fold in the heart post-myocardial infarction (MI). Moreover, we identified, in the Ces1d, a specific RNA element (element D) that is responsible for the improvement of modRNA translation and leads to a 2.5-fold translation increment over Luc modRNA carrying artificial 5' UTR, post-MI. Importantly, we were able to show that 5' UTR Ces1d also enhances modRNA translation in the liver, but not in the kidney, post-ischemic injury, indicating that Ces1d 5' UTR and element D may play a wider role in translation of protein under an ischemic condition.

Glucagon-Dependent Suppression of mTORC1 is Associated with Upregulation of Hepatic FGF21 mRNA Translation

Fibroblast growth factor 21 (FGF21) is a peptide hormone that acts to enhance insulin sensitivity and reverse many of the metabolic defects associated with consumption of a high-fat diet. Recent studies show that the liver is the primary source of FGF21 in the blood, and that hepatic FGF21 expression is upregulated by glucagon. Interestingly, glucagon acts to upregulate FGF21 production by primary cultures of rat hepatocytes and H4IIE and HepG2 hepatocarcinoma cells independent of changes in FGF21 mRNA abundance, suggesting that FGF21 protein expression is regulated post-transcriptionally. Based on these observations, the goal of the present study was to assess whether or not FGF21 mRNA is translationally regulated. The results show that FGF21 mRNA translation and secretion of the hormone are significantly upregulated in H4IIE cells exposed to 25 nM glucagon, independent of changes in FGF21 mRNA abundance. Furthermore, the glucagon-induced upregulation of FGF21 mRNA translation is associated with suppressed activity of the mechanistic target of rapamycin in complex 1 (mTORC1). Similarly, the results show that rapamycin-induced suppression of mTORC1 leads to upregulation of FGF21 mRNA translation with no change in FGF21 mRNA abundance. In contrast, activation of mTORC1 by refreshing the culture medium leads to downregulation of FGF21 mRNA translation. Notably, re-feeding fasted rats also leads to downregulation of FGF21 mRNA translation concomitantly with activation of mTORC1 in the liver. Overall, the findings support a model in which glucagon acts to upregulate FGF21 production by hepatocytes through suppression of mTORC1 and subsequent upregulation of FGF21 mRNA translation.

Genome-Wide Identification and Characterization of UTR-Introns of Citrus sinensis

Introns exist not only in coding sequences (CDSs) but also in untranslated regions (UTRs) of a gene. Recent studies in animals and model plants such as Arabidopsis have revealed that the UTR-introns (UIs) are widely presented in most genomes and involved in regulation of gene expression or RNA stability. In the present study, we identified introns at both 5'UTRs (5UIs) and 3'UTRs (3UIs) of sweet orange genes, investigated their size and nucleotide distribution characteristics, and explored the distribution of cis-elements in the UI sequences. Functional category of genes with predicted UIs were further analyzed using GO, KEGG, and PageMan enrichment. In addition, the organ-dependent splicing and abundance of selected UI-containing genes in root, leaf, and stem were experimentally determined. Totally, we identified 825 UI- and 570 3UI-containing transcripts, corresponding to 617 and 469 genes, respectively. Among them, 74 genes contain both 5UI and 3UI. Nucleotide distribution analysis showed that 5UI distribution is biased at both ends of 5'UTR whiles 3UI distribution is biased close to the start site of 3'UTR. Cis- elements analysis revealed that 5UI and 3UI sequences were rich of promoter-enhancing related elements, indicating that they might function in regulating the expression through them. Function enrichment analysis revealed that genes containing 5UI are significantly enriched in the RNA transport pathway. While, genes containing 3UI are significantly enriched in splicesome. Notably, many pentatricopeptide repeat-containing protein genes and the disease resistance genes were identified to be 3UI-containing. RT-PCR result confirmed the existence of UIs in the eight selected gene transcripts whereas alternative splicing events were found in some of them. Meanwhile, qRT-PCR result showed that UIs were differentially expressed among organs, and significant correlation was found between some genes and their UIs, for example: The expression of VPS28 and its 3UI was significantly negative correlated. This is the first report about the UIs in sweet orange from genome-wide level, which could provide evidence for further understanding of the role of UIs in gene expression regulation.

Identification of transcription factors that bind to the 5'-UTR of the barley PHO2 gene

In barley and other higher plants, phosphate homeostasis is maintained by a regulatory network involving the PHO2 (PHOSPHATE2) encoding ubiquitin-conjugating (UBC) E2 enzyme, the PHR1 (PHOSPHATE STARVATION RESPONSE 1) transcription factor (TF), IPS1 (INDUCED BYPHOSPHATESTARVATION1) RNA, and miR399. During phosphate ion (Pi) deprivation, PHR1 positively regulates MIR399 expression, after transcription and processing mature miR399 guides the Ago protein to the 5'-UTR of PHO2 transcripts. Non-coding IPS1 RNA is highly expressed during Pi starvation, and the sequestration of miR399 molecules protects PHO2 mRNA from complete degradation. Here, we reveal new cis- and trans-regulatory elements that are crucial for efficient PHO2 gene expression in barley. We found that the 5'-UTR of PHO2 contains two PHR1 binding sites (P1BSs) and one Pi-responsive PHO element. Using a yeast one-hybrid (Y1H) assay, we identified two candidate proteins that might mediate this transcriptional regulation: a barley PHR1 ortholog and a TF containing an uncharacterized MYB domain. Additional results classified this new potential TF as belonging to the APL (ALTERED PHLOEM DEVELOPMENT) protein family, and we observed its nuclear localization in barley protoplasts. Pi starvation induced the accumulation of barley APL transcripts in both the shoots and roots. Interestingly, the deletion of the P1BS motif from the first intron of the barley 5'-UTR led to a significant increase in the transcription of a downstream β-glucuronidase (GUS) reporter gene in tobacco leaves. Our work extends the current knowledge about putative cis- and trans-regulatory elements that may affect the expression of the barley PHO2 gene.

5' untranslated regions: the next regulatory sequence in yeast synthetic biology

When developing industrial biotechnology processes, Saccharomyces cerevisiae (baker's yeast or brewer's yeast) is a popular choice as a microbial host. Many tools have been developed in the fields of synthetic biology and metabolic engineering to introduce heterologous pathways and tune their expression in yeast. Such tools mainly focus on controlling transcription, whereas post-transcriptional regulation is often overlooked. Herein we discuss regulatory elements found in the 5' untranslated region (UTR) and their influence on protein synthesis. We provide not only an overall picture, but also a set of design rules on how to engineer a 5' UTR. The reader is also referred to currently available models that allow gene expression to be tuned predictably using different 5' UTRs.

Prediction of Deleterious Non-synonymous SNPs of Human STK11 Gene by Combining Algorithms, Molecular Docking, and Molecular Dynamics Simulation

Serine-threonine kinase11 (STK11) is a tumor suppressor gene which plays a key role in regulating cell growth and apoptosis. It is widely known as a multitasking kinase and engaged in cell polarity, cell cycle arrest, chromatin remodeling, energy metabolism, and Wnt signaling. The substitutions of single amino acids in highly conserved regions of the STK11 protein are associated with Peutz-Jeghers syndrome (PJS), which is an autosomal dominant inherited disorder. The abnormal function of the STK11 protein is still not well understood. In this study, we classified disease susceptible single nucleotide polymorphisms (SNPs) in STK11 by using different computational algorithms. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking and molecular dynamics analysis. Our results show that W239R and W308C variants are likely to be highly deleterious mutations found in the catalytic kinase domain, which may destabilize structure and disrupt the activation of the STK11 protein as well as reduce its catalytic efficiency. The W239R mutant is likely to have a greater impact on destabilizing the protein structure compared to the W308C mutant. In conclusion, these mutants can help to further realize the large pool of disease susceptibilities linked with catalytic kinase domain activation of STK11 and assist to develop an effective drug for associated diseases.

PTBP1 mRNA isoforms and regulation of their translation

Polypyrimidine tract-binding proteins (PTBPs) are RNA binding proteins that regulate a number of posttranscriptional events. Human PTBP1 transits between the nucleus and cytoplasm and is thought to regulate RNA processes in both. However, information about PTBP1 mRNA isoforms and regulation of PTPB1 expression remains incomplete. Here we mapped the major PTBP1 mRNA isoforms in HEK293T cells and identified alternative 5' and 3' untranslated regions (5'-UTRs, 3'-UTRs), as well as alternative splicing patterns in the protein coding region. We also assessed how the observed PTBP1 mRNA isoforms contribute to PTBP1 expression in different phases of the cell cycle. Previously, PTBP1 mRNAs were shown to crosslink to eukaryotic translation initiation factor 3 (eIF3). We find that eIF3 binds differently to each PTBP1 mRNA isoform in a cell cycle dependent manner. We also observe a strong correlation between eIF3 binding to PTBP1 mRNAs and repression of PTBP1 levels during the S phase of the cell cycle. Our results provide evidence of translational regulation of PTBP1 protein levels during the cell cycle, which may affect downstream regulation of alternative splicing and translation mediated by PTBP1 protein isoforms.

An intron-derived motif strongly increases gene expression from transcribed sequences through a splicing independent mechanism in Arabidopsis thaliana

Certain introns significantly increase mRNA accumulation by a poorly understood mechanism. These introns have no effect when located upstream, or more than ~1 Kb downstream, of the start of transcription. We tested the ability of a formerly non-stimulating intron containing 11 copies of the sequence TTNGATYTG, which is over-represented in promoter-proximal introns in Arabidopsis thaliana, to affect expression from various positions. The activity profile of this intron at different locations was similar to that of a natural intron from the UBQ10 gene, suggesting that the motif increases mRNA accumulation by the same mechanism. A series of introns with different numbers of this motif revealed that the effect on expression is linearly dependent on motif copy number up to at least 20, with each copy adding another 1.5-fold increase in mRNA accumulation. Furthermore, 6 copies of the motif stimulated mRNA accumulation to a similar degree from within an intron or when introduced into the 5'-UTR and coding sequences of an intronless construct, demonstrating that splicing is not required for this sequence to boost expression. The ability of this motif to substantially elevate expression from several hundred nucleotides downstream of the transcription start site reveals a novel type of eukaryotic gene regulation.

Post-transcriptional Regulation by 3' UTRs Can Be Masked by Regulatory Elements in 5' UTRs

In mRNA sequences, 3' UTRs are thought to contain most elements that specifically regulate localization, turnover, and translation. Although high-throughput experiments indicate that many RNA-binding proteins (RBPs) also bind 5' UTRs, much less is known about specific post-transcriptional control exerted by 5' UTRs. GLD-1 is a conserved RBP and a translational repressor with essential roles in Caenorhabditis elegans germ cell development. Previously, we showed that GLD-1 binds highly conserved sites in both 3' and 5' UTRs. Here, by targeted single-copy insertion of transgenes, we systematically tested in vivo functionality of 5' and 3' UTR binding sites individually and in combination. Our data show that sites in 5' UTRs mediate specific and strong translational repression, independent of exact position. Intriguingly, we found that the functionality of 3' UTR sites can be masked by 5' UTR sites and vice versa. We conclude that it is important to study both UTRs simultaneously.

Extensions, Extra Factors, and Extreme Complexity: Ribosomal Structures Provide Insights into Eukaryotic Translation

Although the basic aspects of protein synthesis are preserved in all kingdoms of life, there are many important structural and functional differences between bacterial and the more complex eukaryotic ribosomes. High-resolution cryo-electron microscopy (cryo-EM) and X-ray crystallography structures of eukaryotic ribosomes have revealed the complex architectures of eukaryotic ribosomes and species-specific variations in protein and ribosomal RNA (rRNA) extensions. They also enabled structural studies of a range of eukaryotic ribosomal complexes involved in translation initiation, elongation, and termination, revealing unique mechanistic features of the eukaryotic translation process, especially with respect to the identification and recognition of translation start and stop codons on messenger RNAs (mRNAs). Most recently, structural biology has provided insights into the eukaryotic ribosomal biogenesis pathway by visualizing several of its complex intermediates. This review highlights the past decade's structural work on eukaryotic ribosomes and its implications on our understanding of eukaryotic translation.

Functionally Significant Features in the 5' Untranslated Region of the ABCA1 Gene and Their Comparison in Vertebrates

Single nucleotide polymorphisms located in 5′ untranslated regions (5′UTRs) can regulate gene expression and have clinical impact. Recognition of functionally significant sequences within 5′UTRs is crucial in next-generation sequencing applications. Furthermore, information about the behavior of 5′UTRs during gene evolution is scarce. Using the example of the ATP-binding cassette transporter A1 (ABCA1) gene (Tangier disease), we describe our algorithm for functionally significant sequence finding. 5′UTR features (upstream start and stop codons, open reading frames (ORFs), GC content, motifs, and secondary structures) were studied using freely available bioinformatics tools in 55 vertebrate orthologous genes obtained from Ensembl and UCSC. The most conserved sequences were suggested as hot spots. Exon and intron enhancers and silencers (sc35, ighg2 cgamma2, ctnt, gh-1, and fibronectin eda exon), transcription factors (TFIIA, TATA, NFAT1, NFAT4, and HOXA13), some of them cancer related, and microRNA (hsa-miR-4474-3p) were localized to these regions. An upstream ORF, overlapping with the main ORF in primates and possibly coding for a small bioactive peptide, was also detected. Moreover, we showed several features of 5′UTRs, such as GC content variation, hairpin structure conservation or 5′UTR segmentation, which are interesting from a phylogenetic point of view and can stimulate further evolutionary oriented research.

Evolutionary and functional implications of 3' untranslated region length of mRNAs by comprehensive investigation among four taxonomically diverse metazoan species

Background

In eukaryotic organisms, it has been well acknowledged that 3′ untranslated regions (3′ UTRs) of mRNA are actively involved in post-transcriptional regulations of gene expression. Although both shortening and lengthening of 3′ UTRs of specific candidate genes were explicitly documented to have functional consequences, landscape of 3′ UTR lengths in relation to evolutionary dynamics and biological meanings remains to be elucidated when large-scale data become available.

Objectives

The primary objective of this study was to revealed different inter- and intra-species patterns on length distribution of 3′ UTRs in comparison with 5′ UTRs and coding regions.

Methods

In the present study, we investigated 3′ UTR lengths in a highly curated set of 57,135 mRNA sequences among four well-studied and taxonomically diverse metazoan species (fruit fly, zebrafish, mouse and human).

Results

The average length ratio of 3′–5′ UTRs considerably increased from fruit fly (twofold) to human (fivefold). Moreover, genes can be characterized by the obviously different patterns of evolutionary change on 3′ UTR lengths. By utilizing the Gene Ontology annotations, genes with differential lengths of 3′ UTRs were suggested to have the divergent functional implications. In human, we further revealed that ubiquitously transcribed genes had higher median lengths of 3′ UTRs than the genes of tissue-restricted expressions.

Conclusion

We conducted a comprehensive analysis and provided an overview regarding 3′ UTRs biology of mRNAs in animals, whereas the mechanistic explanations of 3′ UTRs length variation in correlation to regulation of gene expression still remain to be further studied.

Electronic supplementary material

The online version of this article (10.1007/s13258-019-00808-8) contains supplementary material, which is available to authorized users.

Characterization of mammalian Lipocalin UTRs in silico: Predictions for their role in post-transcriptional regulation

The Lipocalin family is a group of homologous proteins characterized by its big array of functional capabilities. As extracellular proteins, they can bind small hydrophobic ligands through a well-conserved β-barrel folding. Lipocalins evolutionary history sprawls across many different taxa and shows great divergence even within chordates. This variability is also found in their heterogeneous tissue expression pattern. Although a handful of promoter regions have been previously described, studies on UTR regulatory roles in Lipocalin gene expression are scarce. Here we report a comprehensive bioinformatic analysis showing that complex post-transcriptional regulation exists in Lipocalin genes, as suggested by the presence of alternative UTRs with substantial sequence conservation in mammals, alongside a high diversity of transcription start sites and alternative promoters. Strong selective pressure could have operated upon Lipocalins UTRs, leading to an enrichment in particular sequence motifs that limit the choice of secondary structures. Mapping these regulatory features to the expression pattern of early and late diverging Lipocalins suggests that UTRs represent an additional phylogenetic signal, which may help to uncover how functional pleiotropy originated within the Lipocalin family.