Folding free energies of 5'-UTRs impact post-transcriptional regulation on a genomic scale in yeast

Posttranscriptional regulation of the T-box gene midline via the 3'UTR in Drosophila is complex and cell- and tissue-dependent

The T-box (Tbx) proteins have a 180-230 amino acid DNA-binding domain, first reported in the Brachyury (T) protein. They are highly conserved among metazoans. They regulate a multitude of cellular functions in development and disease. Here, we report posttranscriptional and translational regulation of midline (mid), a Tbx member in Drosophila. We found that the 3'UTR of mid has mRNA degradation elements and AT-rich sequences. In Schneider S2 cells, mid-mRNA could be detected only when the transgene was without the 3'UTR. Similarly, the 3'UTR linked to the Renilla luciferase reporter significantly reduced the activity of the Luciferase, whereas deleting only the degradation elements from the 3'UTR resulted in reduced activity, but not as much. Overexpression of mid in MP2, an embryonic neuroblast, showed no significant difference in the levels of mid-mRNA between the 2 transgenes, with and without the 3'UTR, indicating the absence of posttranscriptional regulation of mid in MP2. Moreover, while elevated mid-RNA was detected in MP2 in nearly all hemisegments, only a fifth of those hemisegments had elevated levels of the protein. Overexpression of the 2 transgenes resulted in MP2-lineage defects at about the same frequency. These results indicate a translational/posttranslational regulation of mid in MP2. The regulation of ectopically expressed mid in the wing imaginal disc was complex. In the wing disc, where mid is not expressed, the ectopic expression of the transgene lacking the 3'UTR had a higher level of mid-RNA and the protein had a stronger phenotypic effect. These results indicate that the 3'UTR can subject mid-mRNA to degradation in a cell- and tissue-specific manner. We further report a balancer-mediated transgenerational modifier effect on the expression and gain of function effects of the 2 transgenes.

Determination of codon pattern and evolutionary forces acting on genes linked to inflammatory bowel disease

Inflammatory bowel disease (IBD) is an inflammatory disorder of the gastrointestinal tract. The present study attempted to understand the codon usage preferences in genes associated with IBD progression. Compositional analysis, codon usage bias (CUB), Relative synonymous codon usage (RSCU), RNA structure, and expression analysis were performed to obtain a comprehensive picture of codon usage in IBD genes. Compositional analysis of 62 IBD-associated genes revealed that G and T are the most and least abundant nucleotides, respectively. ApG, CpA, and TpG dinucleotides were overrepresented or randomly used, while ApC, CpG, GpT, and TpA dinucleotides were either underrepresented or randomly used in genes related to IBD. The codons influencing the codon usage the most in IBD genes were CGC and AGG. A comparison of codon usage between IBD, and pancreatitis (non-IBD inflammatory disease) indicated that only codon CTG codon usage was significantly different between IBD and pancreatitis. At the same time, there were codons ATA, ACA, CGT, CAA, GTA, CCT, ATT, GCT, CGG, TTG, and CAG for whom codon usage was significantly different for IBD and housekeeping gene sets. The results suggest similar codon usage in at least two inflammatory disorders, IBD and pancreatitis. The analysis helps understand the codon biology, factors affecting gene expression of IBD-associated genes, and the evolution of these genes. The study helps reveal the molecular patterns associated with IBD.

Codon-optimization in gene therapy: promises, prospects and challenges

Codon optimization has evolved to enhance protein expression efficiency by exploiting the genetic code's redundancy, allowing for multiple codon options for a single amino acid. Initially observed in E. coli, optimal codon usage correlates with high gene expression, which has propelled applications expanding from basic research to biopharmaceuticals and vaccine development. The method is especially valuable for adjusting immune responses in gene therapies and has the potenial to create tissue-specific therapies. However, challenges persist, such as the risk of unintended effects on protein function and the complexity of evaluating optimization effectiveness. Despite these issues, codon optimization is crucial in advancing gene therapeutics. This study provides a comprehensive review of the current metrics for codon-optimization, and its practical usage in research and clinical applications, in the context of gene therapy.

Compositional Features and Codon Usage Pattern of Genes Associated with Parkinson's Disease

Codon usage bias (CUB) is the phenomenon of non-uniform usage of synonymous codons in which some codons are more used than others and it helps in understanding the molecular organization of genome. Bioinformatic approach was used to analyze the protein-coding sequences of genes associated with Parkinson's disease (PD) to explore compositional features and codon usage pattern as no details work was reported yet. The average improved effective number of codons (Nc) and Nc prime were 42.74 and 44.26 respectively, indicated that CUB was low in these genes. In most of the genes, the overall GC content was almost 50% and GC content at the 1st codon position was the highest while GC content at the 2nd codon position was lowest. Relative synonymous codon usage (RSCU) analysis elucidated over-represented (p > 1.6) and under-represented codons (p < 0.6). The GTG (Val) is the only codon over-represented in all genes. Over-represented codons except (GTG) were A or T ending while under-represented codons (except ACT) were G or C ending. The codons namely TTA (Leu), CTA (Leu), ATC (Ile), ATA (Ile), AGT (Ser), AAC (Asn), TGT (Cys), TGC (Cys), CGC (Arg), AGA (Arg), and AGG (Arg) were absent in SNCA1 to SNCA8 genes. The codon TCG (Ser) was absent in all genes except UCHL1 and PINK1. Correspondence analysis (COA) revealed that the pattern of codon usage differs among genes associated with PD. Neutrality plot analysis indicated some of the points are diagonal distribution suggested that mutation pressure influenced the CUB in genes associated with PD.

Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A proxitome-RNA-capture approach reveals that processing bodies repress coregulated hub genes

Cellular condensates are usually ribonucleoprotein assemblies with liquid- or solid-like properties. Because these subcellular structures lack a delineating membrane, determining their compositions is difficult. Here we describe a proximity-biotinylation approach for capturing the RNAs of the condensates known as processing bodies (PBs) in Arabidopsis (Arabidopsis thaliana). By combining this approach with RNA detection, in silico, and high-resolution imaging approaches, we studied PBs under normal conditions and heat stress. PBs showed a much more dynamic RNA composition than the total transcriptome. RNAs involved in cell wall development and regeneration, plant hormonal signaling, secondary metabolism/defense, and RNA metabolism were enriched in PBs. RNA-binding proteins and the liquidity of PBs modulated RNA recruitment, while RNAs were frequently recruited together with their encoded proteins. In PBs, RNAs follow distinct fates: in small liquid-like PBs, RNAs get degraded while in more solid-like larger ones, they are stored. PB properties can be regulated by the actin-polymerizing SCAR (suppressor of the cyclic AMP)-WAVE (WASP family verprolin homologous) complex. SCAR/WAVE modulates the shuttling of RNAs between PBs and the translational machinery, thereby adjusting ethylene signaling. In summary, we provide an approach to identify RNAs in condensates that allowed us to reveal a mechanism for regulating RNA fate.

Tailor made: the art of therapeutic mRNA design

mRNA medicine is a new and rapidly developing field in which the delivery of genetic information in the form of mRNA is used to direct therapeutic protein production in humans. This approach, which allows for the quick and efficient identification and optimization of drug candidates for both large populations and individual patients, has the potential to revolutionize the way we prevent and treat disease. A key feature of mRNA medicines is their high degree of designability, although the design choices involved are complex. Maximizing the production of therapeutic proteins from mRNA medicines requires a thorough understanding of how nucleotide sequence, nucleotide modification and RNA structure interplay to affect translational efficiency and mRNA stability. In this Review, we describe the principles that underlie the physical stability and biological activity of mRNA and emphasize their relevance to the myriad considerations that factor into therapeutic mRNA design.

The regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

Highlights

In vivo MPRA systematically interrogates the regulatory potential of endogenous 5' UTRs The 5' UTR alone is sufficient to regulate the dynamics of ribosome recruitment during early embryogenesis The MPRA identifies 5' UTR cis -regulatory motifs for translation initiation control 5' UTR length, upstream AUGs and motif grammar contribute to the differential regulatory capability of 5' UTR switching isoforms.

Dynamic Translational Landscape Revealed by Genome-Wide Ribosome Profiling under Drought and Heat Stress in Potato

The yield and quality of potatoes, an important staple crop, are seriously threatened by high temperature and drought stress. In order to deal with this adverse environment, plants have evolved a series of response mechanisms. However, the molecular mechanism of potato's response to environmental changes at the translational level is still unclear. In this study, we performed transcriptome- and ribosome-profiling assays with potato seedlings growing under normal, drought, and high-temperature conditions to reveal the dynamic translational landscapes for the first time. The translational efficiency was significantly affected by drought and heat stress in potato. A relatively high correlation (0.88 and 0.82 for drought and heat stress, respectively) of the fold changes of gene expression was observed between the transcriptional level and translational level globally based on the ribosome-profiling and RNA-seq data. However, only 41.58% and 27.69% of the different expressed genes were shared by transcription and translation in drought and heat stress, respectively, suggesting that the transcription or translation process can be changed independently. In total, the translational efficiency of 151 (83 and 68 for drought and heat, respectively) genes was significantly changed. In addition, sequence features, including GC content, sequence length, and normalized minimal free energy, significantly affected the translational efficiencies of genes. In addition, 28,490 upstream open reading frames (uORFs) were detected on 6463 genes, with an average of 4.4 uORFs per gene and a median length of 100 bp. These uORFs significantly affected the translational efficiency of downstream major open reading frames (mORFs). These results provide new information and directions for analyzing the molecular regulatory network of potato seedlings in response to drought and heat stress.

Rapid 40S scanning and its regulation by mRNA structure during eukaryotic translation initiation

How the eukaryotic 43S preinitiation complex scans along the 5' untranslated region (5' UTR) of a capped mRNA to locate the correct start codon remains elusive. Here, we directly track yeast 43S-mRNA binding, scanning, and 60S subunit joining by real-time single-molecule fluorescence spectroscopy. 43S engagement with mRNA occurs through a slow, ATP-dependent process driven by multiple initiation factors including the helicase eIF4A. Once engaged, 43S scanning occurs rapidly and directionally at ∼100 nucleotides per second, independent of multiple cycles of ATP hydrolysis by RNA helicases post ribosomal loading. Scanning ribosomes can proceed through RNA secondary structures, but 5' UTR hairpin sequences near start codons drive scanning ribosomes at start codons backward in the 5' direction, requiring rescanning to arrive once more at a start codon. Direct observation of scanning ribosomes provides a mechanistic framework for translational regulation by 5' UTR structures and upstream near-cognate start codons.

A 5′ UTR Mutation Contributes to Down-Regulation of Bbs7 in the Berlin Fat Mouse

The Bardet–Biedl Syndrome 7 (Bbs7) gene was identified as the most likely candidate gene causing juvenile obesity in the Berlin Fat Mouse Inbred (BFMI) line. Bbs7 expression is significantly lower in the brain, adipose tissue, and liver of BFMI mice compared to lean C57BL/6NCrl (B6N) mice. A DNA sequence comparison between BFMI and B6N revealed 16 sequence variants in the Bbs7 promoter region. Here, we tested if these mutations contribute to the observed differential expression of Bbs7. In a cell-based dual-luciferase assay, we compared the effects of the BFMI and the B6N haplotypes of different regions of the Bbs7 promotor on the reporter gene expression. A single-nucleotide polymorphism (SNP) was identified causing a significant reduction in the reporter gene expression. This SNP (rs29947545) is located in the 5′ UTR of Bbs7 at Chr3:36.613.350. The SNP is not unique to BFMI mice but also occurs in several other mouse strains, where the BFMI allele is not associated with lower Bbs7 transcript amounts. Thus, we suggest a compensatory mutation in the other mouse strains that keeps Bbs7 expression at the normal level. This compensatory mechanism is missing in BFMI mice and the cell lines tested.

Analysis of the Compositional Features and Codon Usage Pattern of Genes Involved in Human Autophagy

Autophagy plays an intricate role in paradigmatic human pathologies such as cancer, and neurodegenerative, cardiovascular, and autoimmune disorders. Autophagy regulation is performed by a set of autophagy-related (ATG) genes, first recognized in yeast genome and subsequently identified in other species, including humans. Several other genes have been identified to be involved in the process of autophagy either directly or indirectly. Studying the codon usage bias (CUB) of genes is crucial for understanding their genome biology and molecular evolution. Here, we examined the usage pattern of nucleotide and synonymous codons and the influence of evolutionary forces in genes involved in human autophagy. The coding sequences (CDS) of the protein coding human autophagy genes were retrieved from the NCBI nucleotide database and analyzed using various web tools and software to understand their nucleotide composition and codon usage pattern. The effective number of codons (ENC) in all genes involved in human autophagy ranges between 33.26 and 54.6 with a mean value of 45.05, indicating an overall low CUB. The nucleotide composition analysis of the autophagy genes revealed that the genes were marginally rich in GC content that significantly influenced the codon usage pattern. The relative synonymous codon usage (RSCU) revealed 3 over-represented and 10 under-represented codons. Both natural selection and mutational pressure were the key forces influencing the codon usage pattern of the genes involved in human autophagy.

Direct epitranscriptomic regulation of mammalian translation initiation through N4-acetylcytidine

mRNA function is influenced by modifications that modulate canonical nucleobase behavior. We show that a single modification mediates distinct impacts on mRNA translation in a position-dependent manner. Although cytidine acetylation (ac4C) within protein-coding sequences stimulates translation, ac4C within 5' UTRs impacts protein synthesis at the level of initiation. 5' UTR acetylation promotes initiation at upstream sequences, competitively inhibiting annotated start codons. Acetylation further directly impedes initiation at optimal AUG contexts: ac4C within AUG-flanking Kozak sequences reduced initiation in base-resolved transcriptome-wide HeLa results and in vitro utilizing substrates with site-specific ac4C incorporation. Cryo-EM of mammalian 80S initiation complexes revealed that ac4C in the -1 position adjacent to an AUG start codon disrupts an interaction between C and hypermodified t6A at nucleotide 37 of the initiator tRNA. These findings demonstrate the impact of RNA modifications on nucleobase function at a molecular level and introduce mRNA acetylation as a factor regulating translation in a location-specific manner.

Variation in upstream open reading frames contributes to allelic diversity in maize protein abundance

SignificanceProteins are the machinery which execute essential cellular functions. However, measuring their abundance within an organism can be difficult and resource-intensive. Cells use a variety of mechanisms to control protein synthesis from mRNA, including short open reading frames (uORFs) that lie upstream of the main coding sequence. Ribosomes can preferentially translate uORFs instead of the main coding sequence, leading to reduced translation of the main protein. In this study, we show that uORF sequence variation between individuals can lead to different rates of protein translation and thus variable protein abundances. We also demonstrate that natural variation in uORFs occurs frequently and can be linked to whole-plant phenotypes, indicating that uORF sequence variation likely contributes to plant adaptation.

Increased expression of the fluorescent reporter protein ymNeonGreen in Saccharomyces cerevisiae by reducing RNA secondary structure near the start codon

•

Stable secondary RNA structure 3′ of the start codon inhibits ymNeonGreen expression.

•

Removing secondary RNA structure increased expression in S. cerevisiae and E. coli.

•

Expression was increased 2-fold in S. cerevisiae.

•

Expression was increased 3.8-fold in E. coli.

•

Increased expression in S. cerevisiae was promoter independent.

Expression of a new fluorescent reporter protein called mNeonGreen, that is not based on the jellyfish green fluorescent protein (GFP) sequence, shows increased brightness and folding speed compared to enhanced GFP. However, in vivo brightness of mNeonGreen and its yeast-optimized variant ymNeonGreen in S. cerevisiae is lower than expected, limiting the use of this high quantum yield, fast-folding reporter in budding yeast. This study shows that secondary RNA structure near the start codon in the ymNeonGreen ORF inhibits expression in S. cerevisiae. Removing secondary structure, without altering the ymNeonGreen protein sequence, led to a 2 and 4-fold increase in fluorescence when expressed in S. cerevisiae and E. coli, respectively. In S. cerevisiae, increased fluorescence was seen with strong and weak promoters and led to higher transcript levels suggesting greater transcript stability and improved expression in the absence of stable secondary RNA structure near the start codon.

Rational design and construction of multi-copy biomanufacturing islands in mammalian cells

Cell line development is a critical step in the establishment of a biopharmaceutical manufacturing process. Current protocols rely on random transgene integration and amplification. Due to considerable variability in transgene integration profiles, this workflow results in laborious screening campaigns before stable producers can be identified. Alternative approaches for transgene dosage increase and integration are therefore highly desirable. In this study, we present a novel strategy for the rapid design, construction, and genomic integration of engineered multiple-copy gene constructs consisting of up to 10 gene expression cassettes. Key to this strategy is the diversification, at the sequence level, of the individual gene cassettes without altering their protein products. We show a computational workflow for coding and regulatory sequence diversification and optimization followed by experimental assembly of up to nine gene copies and a sentinel reporter on a contiguous scaffold. Transient transfections in CHO cells indicates that protein expression increases with the gene copy number on the scaffold. Further, we stably integrate these cassettes into a pre-validated genomic locus. Altogether, our findings point to the feasibility of engineering a fully mapped multi-copy recombinant protein 'production island' in a mammalian cell line with greatly reduced screening effort, improved stability, and predictable product titers.

Predicting sequence and structural features of effective piRNA target binding sites

Piwi-interacting RNA (piRNA) targets are usually identified through base pairing between the piRNA seed region and complementary bases on the target mRNAs, which often results in false predictions. Crosslinking immunoprecipitation (CLIP) study emerges as a promising method that enables accurate identification of PIWI-clade-based targets containing RNA-binding sites. In the present study, we have analyzed the piRNA-target CLIP-seq datasets to uncover the additional characteristic features of piRNA targets. We studied important sequence and structural features using IP⁺ and IP^- set targets that might enhance the accuracy of target site predictions. Analysis has revealed substantial enrichment of AU in target sites as well as in and around the 30 nts upstream and downstream of target sites in IP⁺ set relative to IP^- set that might be contributing to lowering the minimal folding energy of target sites of IP⁺ + set that might be easing the base pairing between piRNA and their targets. We have also found a lower MFE threshold (en) and higher miRanda score for piRNA targets. Interestingly, we have found that majority of the target sites are residing within 3'UTR, suggesting 3'UTR as a preferential target site like that of miRNA targets. Thus, we hypothesize that our findings on additional key features of piRNA target sites might be valuable in identifying the potential targets of piRNA accurately, which will aid in decrypting their functional importance. This article is protected by copyright. All rights reserved.

Combinations of slow-translating codon clusters can increase mRNA half-life in Saccharomyces cerevisiae

The presence of a single cluster of nonoptimal codons was found to decrease a transcript's half-life through the interaction of the ribosome-associated quality control machinery with stalled ribosomes in Saccharomyces cerevisiae The impact of multiple nonoptimal codon clusters on a transcript's half-life, however, is unknown. Using a kinetic model, we predict that inserting a second nonoptimal cluster near the 5' end can lead to synergistic effects that increase a messenger RNA's (mRNA's) half-life in S. cerevisiae Specifically, the 5' end cluster suppresses the formation of ribosome queues, reducing the interaction of ribosome-associated quality control factors with stalled ribosomes. We experimentally validate this prediction by introducing two nonoptimal clusters into three different genes and find that their mRNA half-life increases up to fourfold. The model also predicts that in the presence of two clusters, the cluster closest to the 5' end is the primary determinant of mRNA half-life. These results suggest the "translational ramp," in which nonoptimal codons are located near the start codon and increase translational efficiency, may have the additional biological benefit of allowing downstream slow-codon clusters to be present without decreasing mRNA half-life. These results indicate that codon usage bias plays a more nuanced role in controlling cellular protein levels than previously thought.

A trans locus causes a ribosomopathy in hypertrophic hearts that affects mRNA translation in a protein length-dependent fashion

BACKGROUND

Little is known about the impact of trans-acting genetic variation on the rates with which proteins are synthesized by ribosomes. Here, we investigate the influence of such distant genetic loci on the efficiency of mRNA translation and define their contribution to the development of complex disease phenotypes within a panel of rat recombinant inbred lines.

RESULTS

We identify several tissue-specific master regulatory hotspots that each control the translation rates of multiple proteins. One of these loci is restricted to hypertrophic hearts, where it drives a translatome-wide and protein length-dependent change in translational efficiency, altering the stoichiometric translation rates of sarcomere proteins. Mechanistic dissection of this locus across multiple congenic lines points to a translation machinery defect, characterized by marked differences in polysome profiles and misregulation of the small nucleolar RNA SNORA48. Strikingly, from yeast to humans, we observe reproducible protein length-dependent shifts in translational efficiency as a conserved hallmark of translation machinery mutants, including those that cause ribosomopathies. Depending on the factor mutated, a pre-existing negative correlation between protein length and translation rates could either be enhanced or reduced, which we propose to result from mRNA-specific imbalances in canonical translation initiation and reinitiation rates.

CONCLUSIONS

We show that distant genetic control of mRNA translation is abundant in mammalian tissues, exemplified by a single genomic locus that triggers a translation-driven molecular mechanism. Our work illustrates the complexity through which genetic variation can drive phenotypic variability between individuals and thereby contribute to complex disease.

Analysis of codon usage of Horseshoe Bat Hepatitis B virus and its host

In the present analysis, codon usage strategies and base distribution of Horseshoe bat hepatitis B virus (HBHBV) were analyzed and compared with its host Rhinolophus sinicus, as no work was yet reported. The magnitude of synonymous codon usage bias (CUB) in the virus and its host was low with higher proportion of the base C. Notably, 21 more frequently used codons, 19 less frequently used codons and 3 underrepresented codons (TCG, ACG and GCG) were found to be similar in both virus and its host coding sequences. Neutrality plot analysis reported greater role of natural selection in HBHBV (67.84%) and R. sinicus (76.90%) over mutation pressure. Base skewness and protein properties also influenced the CUB of genes. Further, codon usage analysis depicted, HBHBV and R. sinicus had many similarities in codon usage patterns that might reflect viral adaptation to its host.

In vivo structure of the Ty1 retrotransposon RNA genome

Long terminal repeat (LTR)-retrotransposons constitute a significant part of eukaryotic genomes and influence their function and evolution. Like other RNA viruses, LTR-retrotransposons efficiently utilize their RNA genome to interact with host cell machinery during replication. Here, we provide the first genome-wide RNA secondary structure model for a LTR-retrotransposon in living cells. Using SHAPE probing, we explore the secondary structure of the yeast Ty1 retrotransposon RNA genome in its native in vivo state and under defined in vitro conditions. Comparative analyses reveal the strong impact of the cellular environment on folding of Ty1 RNA. In vivo, Ty1 genome RNA is significantly less structured and more dynamic but retains specific well-structured regions harboring functional cis-acting sequences. Ribosomes participate in the unfolding and remodeling of Ty1 RNA, and inhibition of translation initiation stabilizes Ty1 RNA structure. Together, our findings support the dual role of Ty1 genomic RNA as a template for protein synthesis and reverse transcription. This study also contributes to understanding how a complex multifunctional RNA genome folds in vivo, and strengthens the need for studying RNA structure in its natural cellular context.

A comparative genomics-based study of positive strand RNA viruses emphasizing on SARS-CoV-2 utilizing dinucleotide signature, codon usage and codon context analyses

The novel corona virus disease or COVID-19 caused by a positive strand RNA virus (PRV) called SARS-CoV-2 is plaguing the entire planet as we conduct this study. In this study a multifaceted analysis was carried out employing dinucleotide signature, codon usage and codon context to compare and unravel the genomic as well as genic characteristics of the SARS-CoV-2 isolates and how they compare to other PRVs which represents some of the most pathogenic human viruses. The main emphasis of this study was to comprehend the codon biology of the SARS-CoV-2 in the backdrop of the other PRVs like Poliovirus, Japanese encephalitis virus, Hepatitis C virus, Norovirus, Rubella virus, Semliki Forest virus, Zika virus, Dengue virus, Human rhinoviruses and the Betacoronaviruses since codon usage pattern along with the nucleotide composition prevalent within the viral genome helps to understand the biology and evolution of viruses. Our results suggest discrete genomic dinucleotide signature within the PRVs. Some of the genes from the different SARS-CoV-2 isolates were also found to demonstrate heterogeneity in terms of their dinucleotide signature. The SARS-CoV-2 isolates also demonstrated a codon context trend characteristically dissimilar to the other PRVs. The findings of this study are expected to contribute to the developing global knowledge base in countering COVID-19.

STARRPeaker: uniform processing and accurate identification of STARR-seq active regions

STARR-seq technology has employed progressively more complex genomic libraries and increased sequencing depths. An issue with the increased complexity and depth is that the coverage in STARR-seq experiments is non-uniform, overdispersed, and often confounded by sequencing biases, such as GC content. Furthermore, STARR-seq readout is confounded by RNA secondary structure and thermodynamic stability. To address these potential confounders, we developed a negative binomial regression framework for uniformly processing STARR-seq data, called STARRPeaker. Moreover, to aid our effort, we generated whole-genome STARR-seq data from the HepG2 and K562 human cell lines and applied STARRPeaker to comprehensively and unbiasedly call enhancers in them.

Polymorphisms of the ACSL1 Gene Influence Milk Production Traits and Somatic Cell Score in Chinese Holstein Cows

Simple Summary

Milk production traits of cows are important economic indicators of the livestock industry. Many dairy farms strive to improve the quality of their milk. Long-chain acyl-CoA synthetase 1 (ACSL1) is a gene related to lipid metabolism. It is widely found in various organisms and can affect fat content and protein content in milk. Single nucleotide polymorphisms (SNP) refers to the polymorphism of DNA sequence caused by a single nucleotide variation at the gene level, which plays a vital function in the genetic study of milk production traits in dairy cows. Our study identified six SNPs of the ACSL1 gene in Chinese Holstein cows, which were related to milk yield, milk fat content, milk protein content and somatic cell score (SCS) to some extent. In summary, the pleiotropic effects of bovine ACSL1 for milk production traits were found in this paper, which will provide a reference for Chinese Holstein cow breeding selection and high economic benefits.

Abstract

Improving the quality of milk is a challenge for zootechnicians and dairy farms across the globe. Long-chain acyl-CoA synthetase 1 (ACSL1) is a significant member of the long-chain acyl-CoA synthetase gene family. It is widely found in various organisms and influences the lactation performance of cows, including fat percentage, milk protein percentage etc. Our study was aimed to investigate the genetic effects of single nucleotide polymorphisms (SNPs) in ACSL1 on milk production traits. Twenty Chinese Holstein cows were randomly selected to extract DNA from their blood samples for PCR amplification and sequencing to identify SNPs of the bovine ACSL1 gene, and six SNPs (5’UTR-g.20523C>G, g.35446C>T, g.35651G>A, g.35827C>T, g.35941G>A and g.51472C>T) were discovered. Then, Holstein cow genotyping (n = 992) was performed by Sequenom MassARRAY based on former SNP information. Associations between SNPs and milk production traits and somatic cell score (SCS) were analyzed by the least-squares method. The results showed that SNP g.35827C>T was in high linkage disequilibrium with g.35941G>A. Significant associations were found between SNPs and test-day milk yield (TDMY), fat content (FC), protein content (PC) and SCS (p < 0.05). Among these SNPs, SNP 5’UTR-g.20523C>G showed an extremely significant effect on PC and SCS (p < 0.01). The SNP g.35446C>T showed a statistically significant effect on FC, PC, and SCS (p < 0.01), and also TDMY (p < 0.05). The SNP g.35651G>A had a statistically significant effect on PC (p < 0.01). The SNP g.35827C>T showed a highly significant effect on TDMY, FC, and SCS (p < 0.01) and significantly influenced PC (p < 0.05). Lastly, SNP g.51472C>T was significantly associated with TDMY, FC, and SCS (p < 0.05). In summary, the pleiotropic effects of bovine ACSL1 for milk production traits were found in this paper, but further investigation will be required on the intrinsic correlation to provide a theoretical basis for the research on molecular genetics of milk quality traits of Holstein cows.

Improving the Expression of Human Granulocyte Colony Stimulating Factor in Escherichia coli by Reducing the GC-content and Increasing mRNA Folding Free Energy at 5'-Terminal End

Purpose: Strategy for improving the production of biopharmaceutical protein continues to develop due to increasing market demand. Human granulocyte colony stimulating factor (hG-CSF) is one of biopharmaceutical proteins that has many applications, and easily produced in Escherichia coli expression system. Previous studies reported that codon usage, rare codon, mRNA folding and GC-content at 5’-terminal end were crucial for protein production in E. coli. In the present study, the effect of reducing the GC-content and increasing the mRNA folding free energy at the 5’-terminal end on the expression level of hG-CSF proteins was investigated.

Methods: Synonymous codon substitutions were performed to generate mutant variants of open reading frame (ORF) with lower GC-content at 5’-terminal ends. Oligoanalyzer tool was used to calculate the GC content of eight codons sequence after ATG. Whereas, mRNA folding free energy was predicted using KineFold and RNAfold tools. The template DNA was amplified using three variant forward primers and one same reverse primer. Those DNA fragments were individually cloned into pJexpress414 expression vector and were confirmed using restriction and DNA sequencing analyses. The confirmed constructs were transformed into E. coli NiCo21(DE3) host cells and the recombinant protein was expressed using IPTG-induction. Total protein obtained were characterized using SDS-PAGE, Western blot and ImageJ software analyses.

Results: The result showed that the mutant variant with lower GC-content and higher mRNA folding free energy near the translation initiation region (TIR) could produce a higher amount of hG-CSF proteins compared to the original gene sequence.

Conclusion: This study emphasized the important role of the nucleotide composition immediately downstream the start codon to achieve high-yield protein product on heterologous expression in E. coli.

On the Way to Understanding the Interplay between the RNA Structure and Functions in Cells: A Genome-Wide Perspective

RNAs adopt specific structures in order to perform their biological activities. The structure of RNA is an important layer of gene expression regulation, and can impact a plethora of cellular processes, starting with transcription, RNA processing, and translation, and ending with RNA turnover. The development of high-throughput technologies has enabled a deeper insight into the sophisticated interplay between the structure of the cellular transcriptome and the living cells environment. In this review, we present the current view on the RNA structure in vivo resulting from the most recent transcriptome-wide studies in different organisms, including mammalians, yeast, plants, and bacteria. We focus on the relationship between the mRNA structure and translation, mRNA stability and degradation, protein binding, and RNA posttranscriptional modifications.

Ribosome profiling reveals a functional role for autophagy in mRNA translational control

Autophagy promotes protein degradation, and therefore has been proposed to maintain amino acid pools to sustain protein synthesis during metabolic stress. To date, how autophagy influences the protein synthesis landscape in mammalian cells remains unclear. Here, we utilize ribosome profiling to delineate the effects of genetic ablation of the autophagy regulator, ATG12, on translational control. In mammalian cells, genetic loss of autophagy does not impact global rates of cap dependent translation, even under starvation conditions. Instead, autophagy supports the translation of a subset of mRNAs enriched for cell cycle control and DNA damage repair. In particular, we demonstrate that autophagy enables the translation of the DNA damage repair protein BRCA2, which is functionally required to attenuate DNA damage and promote cell survival in response to PARP inhibition. Overall, our findings illuminate that autophagy impacts protein translation and shapes the protein landscape.

Codon usage pattern and its influencing factors in different genomes of hepadnaviruses

Codon usage bias (CUB) arises from the preference for a codon over codons for the same amino acid. The major factors contributing to CUB are evolutionary forces, compositional properties, gene expression, and protein properties. The present analysis was performed to investigate the compositional properties and the extent of CUB across the genomes of members of the family Hepadnaviridae, as previously no work using bioinformatic tools has been reported. The viral genes were found to be AT rich with low CUB. Analysis of relative synonymous codon usage (RSCU) was used to identify overrepresented and underrepresented codons for each amino acid. Correlation analysis of overall nucleotide composition and its composition at the third codon position suggested that mutation pressure might influence the CUB. A highly significant correlation was observed between GC12 and GC3 (r = 0.910, p < 0.01), indicating that directional mutation affected all three codon positions across the genome. Translational selection (P2) and mutational responsive index (MRI) values of genes suggested that mutation plays a more important role than translational selection in members of the family Hepadnaviridae.

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.

Modular 5'-UTR hexamers for context-independent tuning of protein expression in eukaryotes

Functional characterization of regulatory DNA elements in broad genetic contexts is a prerequisite for forward engineering of biological systems. Translation initiation site (TIS) sequences are attractive to use for regulating gene activity and metabolic pathway fluxes because the genetic changes are minimal. However, limited knowledge is available on tuning gene outputs by varying TISs in different genetic and environmental contexts. Here, we created TIS hexamer libraries in baker’s yeast Saccharomyces cerevisiae directly 5′ end of a reporter gene in various promoter contexts and measured gene activity distributions for each library. Next, selected TIS sequences, resulted in almost 10-fold changes in reporter outputs, were experimentally characterized in various environmental and genetic contexts in both yeast and mammalian cells. From our analyses, we observed strong linear correlations (R² = 0.75–0.98) between all pairwise combinations of TIS order and gene activity. Finally, our analysis enabled the identification of a TIS with almost 50% stronger output than a commonly used TIS for protein expression in mammalian cells, and selected TISs were also used to tune gene activities in yeast at a metabolic branch point in order to prototype fitness and carotenoid production landscapes. Taken together, the characterized TISs support reliable context-independent forward engineering of translation initiation in eukaryotes.

Dietary restriction induces posttranscriptional regulation of longevity genes

Dietary restriction (DR) increases life span through adaptive changes in gene expression. To understand more about these changes, we analyzed the transcriptome and translatome of Caenorhabditis elegans subjected to DR. Transcription of muscle regulatory and structural genes increased, whereas increased expression of amino acid metabolism and neuropeptide signaling genes was controlled at the level of translation. Evaluation of posttranscriptional regulation identified putative roles for RNA-binding proteins, RNA editing, miRNA, alternative splicing, and nonsense-mediated decay in response to nutrient limitation. Using RNA interference, we discovered several differentially expressed genes that regulate life span. We also found a compensatory role for translational regulation, which offsets dampened expression of a large subset of transcriptionally down-regulated genes. Furthermore, 3' UTR editing and intron retention increase under DR and correlate with diminished translation, whereas trans-spliced genes are refractory to reduced translation efficiency compared with messages with the native 5' UTR. Finally, we find that smg-6 and smg-7, which are genes governing selection and turnover of nonsense-mediated decay targets, are required for increased life span under DR.

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.

The landscape of transcription initiation across latent and lytic KSHV genomes

Precise promoter annotation is required for understanding the mechanistic basis of transcription initiation. In the context of complex genomes, such as herpesviruses where there is extensive genic overlap, identification of transcription start sites (TSSs) is particularly problematic and cannot be comprehensively accessed by standard RNA sequencing approaches. Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic gammaherpesvirus and the etiological agent of Kaposi's sarcoma and the B cell lymphoma primary effusion lymphoma (PEL). Here, we leverage RNA annotation and mapping of promoters for analysis of gene expression (RAMPAGE) and define KSHV TSSs transcriptome-wide and at nucleotide resolution in two widely used models of KSHV infection, namely iSLK.219 cells and the PEL cell line TREx-BCBL1-RTA. By mapping TSSs over a 96 h time course of reactivation we confirm 48 of 50 previously identified TSSs. Moreover, we identify over 100 novel transcription start site clusters (TSCs) in each cell line. Our analyses identified cell-type specific differences in TSC positions as well as promoter strength, and defined motifs within viral core promoters. Collectively, by defining TSSs at high resolution we have greatly expanded the transcriptional landscape of the KSHV genome and identified transcriptional control mechanisms at play during KSHV lytic reactivation.

Neomycin-dependent hammerhead ribozymes for the direct control of gene expression in Saccharomyces cerevisiae

Hammerhead ribozyme-based RNA switches have been proven to be powerful tools for conditional gene regulation in various organisms. We present neomycin-dependent hammerhead ribozymes (HHR) that influence gene expression in a ligand- and dose-dependent manner in S. cerevisiae. We utilized a novel design of fusing the aptamer domain to the HHR enabling for the first time the identification of genetic ON- and OFF-switches within the same library. For this purpose a neomycin aptamer was fused to stem 1 of a type 3 hammerhead ribozyme via an addressable three-way junction that shows high flexibility at the connection site. An in vivo screening approach identified sequences that allow to induce or repress gene expression 2- to 3-fold in response to neomycin addition. The ribozyme switches operate at neomycin concentrations that show no toxic effect on cell growth and pose powerful genetic tools to study and modulate cellular function in yeast.

Design of RNA hairpin modules that predictably tune translation in yeast

Modular parts for tuning translation are prevalent in prokaryotic synthetic biology but lacking for eukaryotic synthetic biology. Working in Saccharomyces cerevisiae yeast, we here describe how hairpin RNA structures inserted into the 5′ untranslated region (5′UTR) of mRNAs can be used to tune expression levels by 100-fold by inhibiting translation. We determine the relationship between the calculated free energy of folding in the 5′UTR and in vivo protein abundance, and show that this enables rational design of hairpin libraries that give predicted expression outputs. Our approach is modular, working with different promoters and protein coding sequences, and outperforms promoter mutation as a way to predictably generate a library where a protein is induced to express at a range of different levels. With this new tool, computational RNA sequence design can be used to predictably fine-tune protein production for genes expressed in yeast.

Virus-Induced Changes in mRNA Secondary Structure Uncover cis-Regulatory Elements that Directly Control Gene Expression

mRNAs carry two layers of information, the genetic code and the information that dictates their post-transcriptional fate. The latter function relies on a complex interplay between cis-elements and trans-regulators, and unbiased identification of these elements is still challenging. To identify cis-elements that control gene expression, we use dimethyl sulfate (DMS) mutational profiling with sequencing and map changes in mRNA secondary structure following viral infection. Our dynamic structural data reveal a major role for ribosomes in unwinding secondary structures, which is further supported by the relationship we uncover between structure and translation efficiency. Moreover, our analysis revealed dozens of regions in viral and cellular mRNAs that exhibit changes in secondary structure. In-depth analysis of these regions reveals cis-elements in 3' UTRs that regulate mRNA stability and elements within coding sequences that control translation. Overall, our study demonstrates how mapping dynamic changes in mRNA structure allows unbiased identification of functional regulatory elements.

Bombyx mori cyclin-dependent kinase inhibitor is involved in regulation of the silkworm cell cycle

Cyclin-dependent kinase inhibitors (CKIs) are negative regulators of the cell cycle. They can bind to cyclin-dependent kinase (CDK)-cyclin complexes and inhibit CDK activities. We identified a single homologous gene of the CDK interacting protein/kinase inhibitory protein (Cip/Kip) family, BmCKI, in the silkworm, Bombyx mori. The gene transcribes two splice variants: a 654-bp-long BmCKI-L (the longer splice variant) encoding a protein with 217 amino acids and a 579-bp-long BmCKI-S (the shorter splice variant) encoding a protein with 192 amino acids. BmCKI-L and BmCKI-S contain the Cip/Kip family conserved cyclin-binding domain and the CDK-binding domain. They are localized in the nucleus and have an unconventional bipartite nuclear localization signal at amino acid residues 181-210. Overexpression of BmCKI-L or BmCKI-S affected cell cycle progression; the cell cycle was arrested in the first gap phase of cell cycle (G1). RNA interference of BmCKI-L or BmCKI-S led to cells accumulating in the second gap phase and the mitotic phase of cell cycle (G2/M). Both BmCKI-L and BmCKI-S are involved in cell cycle regulation and probably have similar effects. The transgenic silkworm with BmCKI-L overexpression (BmCKI-L-OE), exhibited embryonic lethal, larva developmental retardation and lethal phenotypes. These results suggest that BmCKI-L might regulate the growth and development of silkworm. These findings clarify the function of CKIs and increase our understanding of cell cycle regulation in the silkworm.

Toward Predictable 5'UTRs in Saccharomyces cerevisiae: Development of a yUTR Calculator

Fine-tuning biosynthetic pathways is crucial for the development of economic feasible microbial cell factories. Therefore, the use of computational models able to predictably design regulatory sequences for pathway engineering proves to be a valuable tool, especially for modifying genes at the translational level. In this study we developed a computational approach for the de novo design of 5'-untranslated regions (5'UTRs) in Saccharomyces cerevisiae with a predictive outcome on translation initiation rate. On the basis of existing data, a partial least-squares (PLS) regression model was trained and showed good performance on predicting protein abundances of an independent test set. This model was further used for the construction of a "yUTR calculator" that can design 5'UTR sequences with a diverse range of desired translation efficiencies. The predictive power of our yUTR calculator was confirmed in vivo by different representative case studies. As such, these results show the great potential of data driven approaches for reliable pathway engineering in S. cerevisiae.

Gene length as a regulator for ribosome recruitment and protein synthesis: theoretical insights

Protein synthesis rates are determined, at the translational level, by properties of the transcript’s sequence. The efficiency of an mRNA can be tuned by varying the ribosome binding sites controlling the recruitment of the ribosomes, or the codon usage establishing the speed of protein elongation. In this work we propose transcript length as a further key determinant of translation efficiency. Based on a physical model that considers the kinetics of ribosomes advancing on the mRNA and diffusing in its surrounding, as well as mRNA circularisation and ribosome drop-off, we explain how the transcript length may play a central role in establishing ribosome recruitment and the overall translation rate of an mRNA. According to our results, the proximity of the 3′ end to the ribosomal recruitment site of the mRNA could induce a feedback in the translation process that would favour the recycling of ribosomes. We also demonstrate how this process may be involved in shaping the experimental ribosome density-gene length dependence. Finally, we argue that cells could exploit this mechanism to adjust and balance the usage of its ribosomal resources.

Translation initiation events on structured eukaryotic mRNAs generate gene expression noise

Gene expression stochasticity plays a major role in biology, creating non-genetic cellular individuality and influencing multiple processes, including differentiation and stress responses. We have addressed the lack of knowledge about posttranscriptional contributions to noise by determining cell-to-cell variations in the abundance of mRNA and reporter protein in yeast. Two types of structural element, a stem–loop and a poly(G) motif, not only inhibit translation initiation when inserted into an mRNA 5΄ untranslated region, but also generate noise. The noise-enhancing effect of the stem–loop structure also remains operational when combined with an upstream open reading frame. This has broad significance, since these elements are known to modulate the expression of a diversity of eukaryotic genes. Our findings suggest a mechanism for posttranscriptional noise generation that will contribute to understanding of the generally poor correlation between protein-level stochasticity and transcriptional bursting. We propose that posttranscriptional stochasticity can be linked to cycles of folding/unfolding of a stem–loop structure, or to interconversion between higher-order structural conformations of a G-rich motif, and have created a correspondingly configured computational model that generates fits to the experimental data. Stochastic events occurring during the ribosomal scanning process can therefore feature alongside transcriptional bursting as a source of noise.

Deep learning of the regulatory grammar of yeast 5' untranslated regions from 500,000 random sequences

Our ability to predict protein expression from DNA sequence alone remains poor, reflecting our limited understanding of cis-regulatory grammar and hampering the design of engineered genes for synthetic biology applications. Here, we generate a model that predicts the protein expression of the 5′ untranslated region (UTR) of mRNAs in the yeast Saccharomyces cerevisiae. We constructed a library of half a million 50-nucleotide-long random 5′ UTRs and assayed their activity in a massively parallel growth selection experiment. The resulting data allow us to quantify the impact on protein expression of Kozak sequence composition, upstream open reading frames (uORFs), and secondary structure. We trained a convolutional neural network (CNN) on the random library and showed that it performs well at predicting the protein expression of both a held-out set of the random 5′ UTRs as well as native S. cerevisiae 5′ UTRs. The model additionally was used to computationally evolve highly active 5′ UTRs. We confirmed experimentally that the great majority of the evolved sequences led to higher protein expression rates than the starting sequences, demonstrating the predictive power of this model.

Melting temperature highlights functionally important RNA structure and sequence elements in yeast mRNA coding regions

Secondary structure elements in the coding regions of mRNAs play an important role in gene expression and regulation, but distinguishing functional from non-functional structures remains challenging. Here we investigate the dependence of sequence–structure relationships in the coding regions on temperature based on the recent PARTE data by Wan et al. Our main finding is that the regions with high and low thermostability (high T_m and low T_m regions) are under evolutionary pressure to preserve RNA secondary structure and primary sequence, respectively. Sequences of low T_m regions display a higher degree of evolutionary conservation compared to high T_m regions. Low T_m regions are under strong synonymous constraint, while high T_m regions are not. These findings imply that high T_m regions contain thermo-stable functionally important RNA structures, which impose relaxed evolutionary constraint on sequence as long as the base-pairing patterns remain intact. By contrast, low thermostability regions contain single-stranded functionally important conserved RNA sequence elements accessible for binding by other molecules. We also find that theoretically predicted structures of paralogous mRNA pairs become more similar with growing temperature, while experimentally measured structures tend to diverge, which implies that the melting pathways of RNA structures cannot be fully captured by current computational approaches.

Antisense oligonucleotides targeting translation inhibitory elements in 5' UTRs can selectively increase protein levels

A variety of diseases are caused by deficiencies in amounts or activity of key proteins. An approach that increases the amount of a specific protein might be of therapeutic benefit. We reasoned that translation could be specifically enhanced using trans-acting agents that counter the function of negative regulatory elements present in the 5' UTRs of some mRNAs. We recently showed that translation can be enhanced by antisense oligonucleotides (ASOs) that target upstream open reading frames. Here we report the amount of a protein can also be selectively increased using ASOs designed to hybridize to other translation inhibitory elements in 5' UTRs. Levels of human RNASEH1, LDLR, and ACP1 and of mouse ACP1 and ARF1 were increased up to 2.7-fold in different cell types and species upon treatment with chemically modified ASOs targeting 5' UTR inhibitory regions in the mRNAs encoding these proteins. The activities of ASOs in enhancing translation were sequence and position dependent and required helicase activity. The ASOs appear to improve the recruitment of translation initiation factors to the target mRNA. Importantly, ASOs targeting ACP1 mRNA significantly increased the level of ACP1 protein in mice, suggesting that this approach has therapeutic and research potentials.

α+-Thalassemia Caused by an 811 bp Deletion in Individuals from Nanning, Guangxi: A Report of Two Cases

A novel 811 bp deletion was first identified in two families of Nanning City, Guangxi Zhuang Autonomous Region, People's Republic of China (PRC). The molecular basis of this anomaly is a deletion from NG_000006.1: g.32945_33755, and is 20 bp upstream of the translation initiation codon of HBA2. From analyses of the blood indices of the two probands, the 811 bp deletion is an α⁺-thalassemia (α⁺-thal). This is the first report of this deletional thalassemia anywhere in the world.

Deleting the first disulphide bond in an arenicin derivative enhances its expression in Pichia pastoris

The marine antimicrobial peptide NZ17074, a variant of arenicin-3 from Arenicola marina that has broad antimicrobial activity and high bioavailability, can be designed to treat bacterial and fungal diseases. To reduce the toxicity of NZ17074, N6 was designed by replacing a cysteine in positions 3 and 20 with alanine, fused to the C-terminus of the small ubiquitin-like modifier tag (SUMO), and expressed in yeast. SUMO-N6 yielded as much as 921 mg l^-1 at 72 h after induction in a fermentor and increased 1·8-fold over SUMO-NZ17074. After cleavage with 30% formic acid and purification by a Sephadex G-25 column, 9·7 mg of the recombinant peptide N6 (rN6) was obtained from one-litre fermentation broth, increasing 1·4-fold over NZ17074. Compared to NZ17074, rN6 displayed almost identical antimicrobial activity with a minimal inhibitory concentration of 0·5, 0·25-0·5, 4, 0·25-16 and 16 μg ml^-1 against Escherichia, Salmonella, Pseudomonas, Staphylococcus and Streptococcus strains. Our results indicate that the first disulphide bond, Cys3-Cys20, in NZ17074 is not necessary for antimicrobial activity and that its deletion might reduce toxicity to host cells. These findings may help design new antimicrobial peptides harbouring fewer disulphide bridges and may have more potent activity.

SIGNIFICANCE AND IMPACT OF THE STUDY

Disulphide bond formation is an important step in the protein expression and can also influence protein secretion. A deletion of the first disulphide bond in NZ17074 increased the secreted level of target protein, and its antimicrobial activity was almost unaffected by the deletion of the first disulphide bond. The first disulphide bond in NZ17074 is favourable for correctly forming another disulphide bond during expression but not necessary for its activity. This may help design and produce a novel class of antimicrobial peptides harbouring fewer disulphide bridges to save the cost.

The Evolutionary Basis of Translational Accuracy in Plants

Mistranslation errors compromise fitness by wasting resources on nonfunctional proteins. In order to reduce the cost of mistranslations, natural selection chooses the most accurately translated codons at sites that are particularly important for protein structure and function. We investigated the determinants underlying selection for translational accuracy in several species of plants belonging to three clades: Brassicaceae, Fabidae, and Poaceae. Although signatures of translational selection were found in genes from a wide range of species, the underlying factors varied in nature and intensity. Indeed, the degree of synonymous codon bias at evolutionarily conserved sites varied among plant clades while remaining uniform within each clade. This is unlikely to solely reflect the diversity of tRNA pools because there is little correlation between synonymous codon bias and tRNA abundance, so other factors must affect codon choice and translational accuracy in plant genes. Accordingly, synonymous codon choice at a given site was affected not only by the selection pressure at that site, but also its participation in protein domains or mRNA secondary structures. Although these effects were detected in all the species we analyzed, their impact on translation accuracy was distinct in evolutionarily distant plant clades. The domain effect was found to enhance translational accuracy in dicot and monocot genes with a high GC content, but to oppose the selection of more accurate codons in monocot genes with a low GC content.