Relationships among stop codon usage bias, its context, isochores, and gene expression level in various eukaryotes

Codon usage bias and phylogenetic analysis of chloroplast genome in 36 gracilariaceae species

Gracilariaceae is a group of marine large red algae and main source of agar with important economic and ecological value. The codon usage patterns of chloroplast genomes in 36 species from Graciliaceae show that GC range from 0.284 to 0.335, the average GC3 range from 0.135 to 0.243 and the value of ENC range from 35.098 to 42.327, which indicates these genomes are rich in AT and prefer to use codons ending with AT in these species. Nc plot, PR2 plot, neutrality plot analyses and correlation analysis indicate that these biases may be caused by multiple factors, such as natural selection and mutation pressure, but prolonged natural selection is the main driving force influencing codon usage preference. The cluster analysis and phylogenetic analysis show that the differentiation relationship of them is different and indicate that codons with weak or unbiased preferences may also play an irreplaceable role in these species' evolution. In addition, we identified 26 common high-frequency codons and 8-18 optimal codons all ending in A/U in these 36 species. Our results will not only contribute to carrying out transgenic work in Gracilariaceae species to maximize the protein yield in the future, but also lay a theoretical foundation for further exploring systematic classification of them.

Condensed, microtubule-coating thin organelles for orthogonal translation in mammalian cells

Membraneless organelles are capable of selectively performing complex tasks in living cells despite dynamically exchanging with their surroundings. This is an exquisite example how self-organization of proteins and RNAs can lead to more complex functionalities in living systems. Importantly, the absence of a membrane boundary can enable easier access to larger macromolecular complexes that can be challenging to be transported across a membrane. We previously formed orthogonally translating designer membraneless organelles by combining phase separation with kinesin motor proteins to highly enrich engineered translational factors in large organelles. We also showed that even submicron thick designer organelles can be formed, by mounting them onto membranes, which, presumable assisted by 2D condensation, leads to thin film-like condensates. In this study we show that orthogonal translation can also be built with fiber-like appearing organelles. Here, the microtubule-end binding protein EB1 was used to form fiber-like OT organelles along the microtubule cytoskeleton that perform highly selective and efficient orthogonal translation. We also show an improved simplified design of OT organelles. Together this extends OT technology and demonstrates that the microtubule cytoskeleton is a powerful platform for advanced synthetic organelle engineering.

Codon usage bias

Codon usage bias is the preferential or non-random use of synonymous codons, a ubiquitous phenomenon observed in bacteria, plants and animals. Different species have consistent and characteristic codon biases. Codon bias varies not only with species, family or group within kingdom, but also between the genes within an organism. Codon usage bias has evolved through mutation, natural selection, and genetic drift in various organisms. Genome composition, GC content, expression level and length of genes, position and context of codons in the genes, recombination rates, mRNA folding, and tRNA abundance and interactions are some factors influencing codon bias. The factors shaping codon bias may also be involved in evolution of the universal genetic code. Codon-usage bias is critical factor determining gene expression and cellular function by influencing diverse processes such as RNA processing, protein translation and protein folding. Codon usage bias reflects the origin, mutation patterns and evolution of the species or genes. Investigations of codon bias patterns in genomes can reveal phylogenetic relationships between organisms, horizontal gene transfers, molecular evolution of genes and identify selective forces that drive their evolution. Most important application of codon bias analysis is in the design of transgenes, to increase gene expression levels through codon optimization, for development of transgenic crops. The review gives an overview of deviations of genetic code, factors influencing codon usage or bias, codon usage bias of nuclear and organellar genes, computational methods to determine codon usage and the significance as well as applications of codon usage analysis in biological research, with emphasis on plants.

Identification of permissive amber suppression sites for efficient non-canonical amino acid incorporation in mammalian cells

The genetic code of mammalian cells can be expanded to allow the incorporation of non-canonical amino acids (ncAAs) by suppressing in-frame amber stop codons (UAG) with an orthogonal pyrrolysyl-tRNA synthetase (PylRS)/tRNAPylCUA (PylT) pair. However, the feasibility of this approach is substantially hampered by unpredictable variations in incorporation efficiencies at different stop codon positions within target proteins. Here, we apply a proteomics-based approach to quantify ncAA incorporation rates at hundreds of endogenous amber stop codons in mammalian cells. With these data, we compute iPASS (Identification of Permissive Amber Sites for Suppression; available at www.bultmannlab.eu/tools/iPASS), a linear regression model to predict relative ncAA incorporation efficiencies depending on the surrounding sequence context. To verify iPASS, we develop a dual-fluorescence reporter for high-throughput flow-cytometry analysis that reproducibly yields context-specific ncAA incorporation efficiencies. We show that nucleotides up- and downstream of UAG synergistically influence ncAA incorporation efficiency independent of cell line and ncAA identity. Additionally, we demonstrate iPASS-guided optimization of ncAA incorporation rates by synonymous exchange of codons flanking the amber stop codon. This combination of in silico analysis followed by validation in living mammalian cells substantially simplifies identification as well as adaptation of sites within a target protein to confer high ncAA incorporation rates.

ANAP: A versatile, fluorescent probe of ion channel gating and regulation

Fluorescence spectroscopy and microscopy are non-destructive methods that provide real-time measurements of ion channel structural dynamics. As such, they constitute a direct path linking the high-resolution structural models from X-ray crystallography and cryo-electron microscopy with the high-resolution functional data from ionic current measurements. The utility of fluorescence as a reporter of channel structure is limited by the palette of available fluorophores. Thiol-reactive fluorophores are small and bright, but are restricted in terms of the positions on a protein that can be labeled and present significant issues with background incorporation. Genetically encoded fluorescent protein tags are specific to a protein of interest, but are very large and usually only used to label the free N- and C-termini of proteins. L-3-(6-acetylnaphthalen-2-ylamino)-2-aminopropionic acid (ANAP) is a fluorescent amino acid that can be specifically incorporated into virtually any site on a protein of interest using amber stop-codon suppression. Due to its environmental sensitivity and potential as a donor in fluorescence resonance energy transfer experiments, it has been adopted by numerous investigators to study voltage, ligand, and temperature-dependent activation of a host of ion channels. Simultaneous measurements of ionic currents and ANAP fluorescence yield exceptional mechanistic insights into channel function. In this chapter, I will summarize the current literature regarding ANAP and ion channels and discuss the practical aspects of using ANAP, including potential pitfalls and confounds.

Disrupting upstream translation in mRNAs is associated with human disease

Ribosome-profiling has uncovered pervasive translation in non-canonical open reading frames, however the biological significance of this phenomenon remains unclear. Using genetic variation from 71,702 human genomes, we assess patterns of selection in translated upstream open reading frames (uORFs) in 5’UTRs. We show that uORF variants introducing new stop codons, or strengthening existing stop codons, are under strong negative selection comparable to protein-coding missense variants. Using these variants, we map and validate gene-disease associations in two independent biobanks containing exome sequencing from 10,900 and 32,268 individuals, respectively, and elucidate their impact on protein expression in human cells. Our results suggest translation disrupting mechanisms relating uORF variation to reduced protein expression, and demonstrate that translation at uORFs is genetically constrained in 50% of human genes.

The significance of translated upstream open reading frames is not well known. Here, the authors investigate genetic variants in these regions, finding that they are under high evolutionary constraint and may contribute to disease.

Sensing through Non-Sensing Ocular Ion Channels

Ion channels are membrane-spanning integral proteins expressed in multiple organs, including the eye. In the eye, ion channels are involved in various physiological processes, like signal transmission and visual processing. A wide range of mutations have been reported in the corresponding genes and their interacting subunit coding genes, which contribute significantly to an array of blindness, termed ocular channelopathies. These mutations result in either a loss- or gain-of channel functions affecting the structure, assembly, trafficking, and localization of channel proteins. A dominant-negative effect is caused in a few channels formed by the assembly of several subunits that exist as homo- or heteromeric proteins. Here, we review the role of different mutations in switching a “sensing” ion channel to “non-sensing,” leading to ocular channelopathies like Leber’s congenital amaurosis 16 (LCA16), cone dystrophy, congenital stationary night blindness (CSNB), achromatopsia, bestrophinopathies, retinitis pigmentosa, etc. We also discuss the various in vitro and in vivo disease models available to investigate the impact of mutations on channel properties, to dissect the disease mechanism, and understand the pathophysiology. Innovating the potential pharmacological and therapeutic approaches and their efficient delivery to the eye for reversing a “non-sensing” channel to “sensing” would be life-changing.

Enabling large-scale genome editing at repetitive elements by reducing DNA nicking

To extend the frontier of genome editing and enable editing of repetitive elements of mammalian genomes, we made use of a set of dead-Cas9 base editor (dBE) variants that allow editing at tens of thousands of loci per cell by overcoming the cell death associated with DNA double-strand breaks and single-strand breaks. We used a set of gRNAs targeting repetitive elements-ranging in target copy number from about 32 to 161 000 per cell. dBEs enabled survival after large-scale base editing, allowing targeted mutations at up to ∼13 200 and ∼12 200 loci in 293T and human induced pluripotent stem cells (hiPSCs), respectively, three orders of magnitude greater than previously recorded. These dBEs can overcome current on-target mutation and toxicity barriers that prevent cell survival after large-scale genome engineering.

Computational discovery and modeling of novel gene expression rules encoded in the mRNA

The transcript is populated with numerous overlapping codes that regulate all steps of gene expression. Deciphering these codes is very challenging due to the large number of variables involved, the non-modular nature of the codes, biases and limitations in current experimental approaches, our limited knowledge in gene expression regulation across the tree of life, and other factors. In recent years, it has been shown that computational modeling and algorithms can significantly accelerate the discovery of novel gene expression codes. Here, we briefly summarize the latest developments and different approaches in the field.

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.

Increased gene translation stringency in mammalian cells by nonsense suppression at multiple permissive sites with a single noncanonical amino acid

The considerable potential of engineered cells compels the development of strategies for the stringent control of gene expression. A promising approach is the introduction of a premature stop codon (PTC) into a selected gene that is expressed only in the presence of noncanonical amino acids through nonsense suppression. Here, different strategies of amber PTC readthrough in mammalian cells were tested. The use of a tRNA synthetase together with a TAG codon-specific tRNA achieved PTC readthrough depending on the addition of a noncanonical amino acid (4-benzoyl-L-phenylalanine; Bpa). While single TAG codon incorporation exhibited detectable expression of the reporter protein even in the absence of Bpa, the use of a double PTC enabled virtually leakage-free functional gene translation. The introduction of an additional 5'-PTC, therefore, represents a generally applicable strategy to increase stringency in gene translation.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

The novel aminoglycoside, ELX-02, permits CTNSW138X translational read-through and restores lysosomal cystine efflux in cystinosis

Background

Cystinosis is a rare disorder caused by recessive mutations of the CTNS gene. Current therapy decreases cystine accumulation, thus slowing organ deterioration without reversing renal Fanconi syndrome or preventing eventual need for a kidney transplant.15-20% of cystinosis patients harbour at least one nonsense mutation in CTNS, leading to premature end of translation of the transcript. Aminoglycosides have been shown to permit translational read-through but have high toxicity level, especially in the kidney and inner ear. ELX-02, a modified aminoglycoside, retains it read-through ability without the toxicity.

Methods and findings

We ascertained the toxicity of ELX-02 in cells and in mice as well as the effect of ELX-02 on translational read-through of nonsense mutations in cystinotic mice and human cells. ELX-02 was not toxic in vitro or in vivo, and permitted read-through of nonsense mutations in cystinotic mice and human cells.

Conclusions

ELX-02 has translational read-through activity and produces a functional CTNS protein, as evidenced by reduced cystine accumulation. This reduction is comparable to cysteamine treatment. ELX-02 accumulates in the kidney but neither cytotoxicity nor nephrotoxicity was observed.

Evaluation and optimisation of unnatural amino acid incorporation and bioorthogonal bioconjugation for site-specific fluorescent labelling of proteins expressed in mammalian cells

Many biophysical techniques that are available to study the structure, function and dynamics of cellular constituents require modification of the target molecules. Site-specific labelling of a protein is of particular interest for fluorescence-based single-molecule measurements including single-molecule FRET or super-resolution microscopy. The labelling procedure should be highly specific but minimally invasive to preserve sensitive biomolecules. The modern molecular engineering toolkit provides elegant solutions to achieve the site-specific modification of a protein of interest often necessitating the incorporation of an unnatural amino acid to introduce a unique reactive moiety. The Amber suppression strategy allows the site-specific incorporation of unnatural amino acids into a protein of interest. Recently, this approach has been transferred to the mammalian expression system. Here, we demonstrate how the combination of unnatural amino acid incorporation paired with current bioorthogonal labelling strategies allow the site-specific engineering of fluorescent dyes into proteins produced in the cellular environment of a human cell. We describe in detail which parameters are important to ensure efficient incorporation of unnatural amino acids into a target protein in human expression systems. We furthermore outline purification and bioorthogonal labelling strategies that allow fast protein preparation and labelling of the modified protein. This way, the complete eukaryotic proteome becomes available for single-molecule fluorescence assays.

Forces acting on codon bias in malaria parasites

Malaria parasite genomes have a range of codon biases, with Plasmodium falciparum one of the most AT-biased genomes known. We examined the make up of synonymous coding sites and stop codons in the core genomes of representative malaria parasites, showing first that local DNA context influences codon bias similarly across P. falciparum, P. vivax and P. berghei, with suppression of CpG dinucleotides and enhancement of CpC dinucleotides, both within and aross codons. Intense asexual phase gene expression in P. falciparum and P. berghei is associated with increased A3:G3 bias but reduced T3:C3 bias at 2-fold sites, consistent with adaptation of codons to tRNA pools and avoidance of wobble tRNA interactions that potentially slow down translation. In highly expressed genes, the A3:G3 ratio can exceed 30-fold while the T3:C3 ratio can be less than 1, according to the encoded amino acid and subsequent base. Lysine codons (AAA/G) show distinctive behaviour with substantially reduced A3:G3 bias in highly expressed genes, perhaps because of selection against frameshifting when the AAA codon is followed by another adenine. Intense expression is also associated with a strong bias towards TAA stop codons (found in 94% and 89% of highly expressed P. falciparum and P. berghei genes respectively) and a proportional rise in the TAAA stop ‘tetranucleotide’. The presence of these expression-linked effects in the relatively AT-rich malaria parasite species adds weight to the suggestion that AT-richness in the Plasmodium genus might be a fitness adaptation. Potential explanations for the relative lack of codon bias in P. vivax include the distinct features of its lifecycle and its effective population size over evolutionary time.

Humanin Prevents Age-Related Cognitive Decline in Mice and is Associated with Improved Cognitive Age in Humans

Advanced age is associated with a decline in cognitive function, likely caused by a combination of modifiable and non-modifiable factors such as genetics and lifestyle choices. Mounting evidence suggests that humanin and other mitochondrial derived peptides play a role in several age-related conditions including neurodegenerative disease. Here we demonstrate that humanin administration has neuroprotective effects in vitro in human cell culture models and is sufficient to improve cognition in vivo in aged mice. Furthermore, in a human cohort, using mitochondrial GWAS, we identified a specific SNP (rs2854128) in the humanin-coding region of the mitochondrial genome that is associated with a decrease in circulating humanin levels. In a large, independent cohort, consisting of a nationally-representative sample of older adults, we find that this SNP is associated with accelerated cognitive aging, supporting the concept that humanin is an important factor in cognitive aging.

Site-Specific Labelling of Native Mammalian Proteins for Single-Molecule FRET Measurements

Human cells are complex entities in which molecular recognition and selection are critical for cellular processes often driven by structural changes and dynamic interactions. Biomolecules appear in different chemical states, and modifications, such as phosphorylation, affect their function. Hence, using proteins in their chemically native state in biochemical and biophysical assays is essential. Single-molecule FRET measurements allow exploration of the structure, function and dynamics of biomolecules but cannot be fully exploited for the human proteome, as a method for the site-specific coupling of organic dyes into native, non-recombinant mammalian proteins is lacking. We address this issue showing the site-specific engineering of fluorescent dyes into human proteins on the basis of bioorthogonal reactions. We show the applicability of the method to study functional and post-translationally modified proteins on the single-molecule level, among them the hitherto inaccessible human Argonaute 2.

Eukaryotic translational termination efficiency is influenced by the 3' nucleotides within the ribosomal mRNA channel

When a stop codon is at the 80S ribosomal A site, there are six nucleotides (+4 to +9) downstream that are inferred to be occupying the mRNA channel. We examined the influence of these downstream nucleotides on translation termination success or failure in mammalian cells at the three stop codons. The expected hierarchy in the intrinsic fidelity of the stop codons (UAA>UAG>>UGA) was observed, with highly influential effects on termination readthrough mediated by nucleotides at position +4 and position +8. A more complex influence was observed from the nucleotides at positions +5 and +6. The weakest termination contexts were most affected by increases or decreases in the concentration of the decoding release factor (eRF1), indicating that eRF1 binding to these signals was rate-limiting. When termination efficiency was significantly reduced by cognate suppressor tRNAs, the observed influence of downstream nucleotides was maintained. There was a positive correlation between experimentally measured signal strength and frequency of the signal in eukaryotic genomes, particularly in Saccharomyces cerevisiae and Drosophila melanogaster. We propose that termination efficiency is not only influenced by interrogation of the stop signal directly by the release factor, but also by downstream ribosomal interactions with the mRNA nucleotides in the entry channel.

Running the Stop Sign: Readthrough of a Premature UAG Termination Signal in the Translation of a Zebrafish (Danio rerio) Taurine Biosynthetic Enzyme

The UAG termination codon is generally recognized as the least efficient and least frequently used of the three universal stop codons. This is substantiated by numerous studies in an array of organisms. We present here evidence of a translational readthrough of a mutant nonsense UAG codon in the transcript from the cysteine sulfinic acid decarboxylase (csad) gene (ENSDARG00000026348) in zebrafish. The csad gene encodes the terminal enzyme in the taurine biosynthetic pathway. Taurine is a critical amino acid for all animals, playing several essential roles throughout the body, including modulation of the immune system. The sa9430 zebrafish strain (ZDB-ALT-130411-5055) has a point mutation leading to a premature stop codon (UAG) 20 amino acids 5' of the normal stop codon, UGA. Data from immunoblotting, enzyme activity assays, and mass spectrometry provide evidence that the mutant is making a CSAD protein identical to that of the wild-type (XP_009295318.1) in terms of size, activity, and amino acid sequence. UAG readthrough has been described in several species, but this is the first presentation of a case in fish. Also presented are the first data substantiating the ability of a fish CSAD to utilize cysteic acid, an alternative to the standard substrate cysteine sulfinic acid, to produce taurine.

Click-MS: Tagless Protein Enrichment Using Bioorthogonal Chemistry for Quantitative Proteomics

Epitope-tagging is an effective tool to facilitate protein enrichment from crude cell extracts. Traditionally, N- or C-terminal fused tags are employed, which, however, can perturb protein function. Unnatural amino acids (UAAs) harboring small reactive handles can be site-specifically incorporated into proteins, thus serving as a potential alternative for conventional protein tags. Here, we introduce Click-MS, which combines the power of site-specific UAA incorporation, bioorthogonal chemistry, and quantitative mass spectrometry-based proteomics to specifically enrich a single protein of interest from crude mammalian cell extracts. By genetic encoding of p-azido-l-phenylalanine, the protein of interest can be selectively captured using copper-free click chemistry. We use Click-MS to enrich proteins that function in different cellular compartments, and we identify protein-protein interactions, showing the great potential of Click-MS for interaction proteomics workflows.

Selective forces and mutational biases drive stop codon usage in the human genome: a comparison with sense codon usage

BACKGROUND

The three stop codons UAA, UAG, and UGA signal the termination of mRNA translation. As a result of a mechanism that is not adequately understood, they are normally used with unequal frequencies.

RESULTS

In this work, we showed that selective forces and mutational biases drive stop codon usage in the human genome. We found that, in respect to sense codons, stop codon usage was affected by stronger selective forces but was less influenced by neutral mutational biases. UGA is the most frequent termination codon in human genome. However, UAA was the preferred stop codon in genes with high breadth of expression, high level of expression, AT-rich coding sequences, housekeeping functions, and in gene ontology categories with the largest deviation from expected stop codon usage. Selective forces associated with the breadth and the level of expression favoured AT-rich sequences in the mRNA region including the stop site and its proximal 3'-UTR, but acted with scarce effects on sense codons, generating two regions, upstream and downstream of the stop codon, with strongly different base composition. By favouring low levels of GC-content, selection promoted labile local secondary structures at the stop site and its proximal 3'-UTR. The compositional and structural context favoured by selection was surprisingly emphasized in the class of ribosomal proteins and was consistent with sequence elements that increase the efficiency of translational termination. Stop codons were also heterogeneously distributed among chromosomes by a mechanism that was strongly correlated with the GC-content of coding sequences.

CONCLUSIONS

In human genome, the nucleotide composition and the thermodynamic stability of stop codon site and its proximal 3'-UTR are correlated with the GC-content of coding sequences and with the breadth and the level of gene expression. In highly expressed genes stop codon usage is compositionally and structurally consistent with highly efficient translation termination signals.

Clinical potential of ataluren in the treatment of Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is an autosomal dominant, X-linked neuromuscular disorder caused by mutations in dystrophin, one of the largest genes known to date. Dystrophin gene mutations are generally transmitted from the mother to male offspring and can occur throughout the coding length of the gene. The majority of the methodologies aimed at treating the disorder have focused on restoring a shorter, although partially functional, dystrophin protein. The approach has the potential of converting a severe DMD phenotype into a milder form of the disease known as Becker muscular dystrophy. Others have focused on ameliorating the disease by targeting secondary pathologies such as inflammation or loss of regeneration. Of great potential is the development of strategies that are capable of restoring full-length dystrophin expression due to their ability to produce a normal, fully functional protein. Among these strategies, the use of read-through compounds (RTCs) that could be administered orally represents an ideal option. Gentamicin has been previously tested in clinical trials for DMD with limited or no success, and its use in the clinic has been dismissed due to issues of toxicity and lack of clear benefits to patients. More recently, new RTCs have been identified and tested in animal models for DMD. This review will focus on one of those RTCs known as ataluren that has now completed Phase III clinical studies for DMD and at providing an overview of the different stages that have led to its clinical development for the disease. The impact that this new drug may have on DMD and its future perspectives will also be described, with an emphasis on the importance of further assessing the clinical benefits of this molecule in patients as it becomes available on the market in different countries.

A flow cytometry-based reporter assay identifies macrolide antibiotics as nonsense mutation read-through agents

A large number of human diseases are caused by nonsense mutations. These mutations result in premature protein termination and the expression of truncated, usually nonfunctional products. A promising therapeutic strategy for patients suffering from premature termination codon (PTC)-mediated disorders is to suppress the nonsense mutation and restore the expression of the affected protein. Such a suppression approach using specific antibiotics and other read-through promoting agents has been shown to suppress PTCs and restore the production of several important proteins. Here, we report the establishment of a novel, rapid, and very efficient method for screening stop-codon read-through agents. We also show that, in both mammalian cells and in a transgenic mouse model, distinct members of the macrolide antibiotic family can induce read-through of disease-causing stop codons leading to re-expression of several key proteins and to reduced disease phenotypes. Taken together, our results may help in the identification and characterization of well-needed customized pharmaceutical PTC suppression agents.

KEY MESSAGES

Establishment of a flow cytometry-based reporter assay to identify nonsense mutation read-through agents. Macrolide antibiotics can induce read-through of disease-causing stop codons. Macrolide-induced protein restoration can alleviate disease-like phenotypes.

High-level expression of a novel liver-targeting fusion interferon with preferred Escherichia coli codon preference and its anti-hepatitis B virus activity in vivo

Background

In our previous study, a novel liver-targeting fusion interferon (IFN-CSP) combining IFN α2b with plasmodium region I peptide was successfully constructed. IFN-CSP has significant inhibition effects on HBV-DNA replication in HepG2.2.15 cells. The aim of the present investigation was focused on how to produce high levels of recombinant IFN-CSP and its in vivo anti-hepatitis B virus (HBV) activity.

Methods

A modified DNA fragment encoding IFN-CSP was synthesized according to Escherichia coli (E. coli) preferred codon usage and transformed into E. coli BL21 (DE3) for protein expression. The induction conditions were systematically examined by combining one-factor experiments with an orthogonal test (L(9)(3)(4)). The antigenicity of the purified protein was characterized by western blot analysis. The in vivo tissue distribution were assayed and compared with native IFN α2b. HBV-transgenic mice were used as in vivo model to evaluate the anti-HBV effect of the recombinant IFN-CSP.

Results

The results showed that the E. coli expression system was very efficient to produce target protein.

Conclusion

Our current research demonstrates for the first time that IFN-CSP gene can be expressed at high levels in E. coli through codon and expression conditions optimization. The purified recombinant IFN-CSP showed liver-targeting potentiality and anti-HBV activity in vivo. The present study further supported the application of IFN-CSP in liver-targeting anti-HBV medicines.

A Comprehensive Analysis of Codon Usage Patterns in Blunt Snout Bream (Megalobrama amblycephala) Based on RNA-Seq Data

Blunt snout bream (Megalobrama amblycephala) is an important fish species for its delicacy and high economic value in China. Codon usage analysis could be helpful to understand its codon biology, mRNA translation and vertebrate evolution. Based on RNA-Seq data for M. amblycephala, high-frequency codons (CUG, AGA, GUG, CAG and GAG), as well as low-frequency ones (NUA and NCG codons) were identified. A total of 724 high-frequency codon pairs were observed. Meanwhile, 14 preferred and 199 avoided neighboring codon pairs were also identified, but bias was almost not shown with one or more intervening codons inserted between the same pairs. Codon usage bias in the regions close to start and stop codons indicated apparent heterogeneity, which even occurs in the flanking nucleotide sequence. Codon usage bias (RSCU and SCUO) was related to GC₃ (GC content of 3rd nucleotide in codon) bias. Six GO (Gene ontology) categories and the number of methylation targets were influenced by GC₃. Codon usage patterns comparison among 23 vertebrates showed species specificities by using GC contents, codon usage and codon context analysis. This work provided new insights into fish biology and new information for breeding projects.

Comprehensive analysis of stop codon usage in bacteria and its correlation with release factor abundance

We present a comprehensive analysis of stop codon usage in bacteria by analyzing over eight million coding sequences of 4684 bacterial sequences. Using a newly developed program called "stop codon counter," the frequencies of the three classical stop codons TAA, TAG, and TGA were analyzed, and a publicly available stop codon database was built. Our analysis shows that with increasing genomic GC content the frequency of the TAA codon decreases and that of the TGA codon increases in a reciprocal manner. Interestingly, the release factor 1-specific codon TAG maintains a more or less uniform frequency (∼20%) irrespective of the GC content. The low abundance of TAG is also valid with respect to expression level of the genes ending with different stop codons. In contrast, the highly expressed genes predominantly end with TAA, ensuring termination with either of the two release factors. Using three model bacteria with different stop codon usage (Escherichia coli, Mycobacterium smegmatis, and Bacillus subtilis), we show that the frequency of TAG and TGA codons correlates well with the relative steady state amount of mRNA and protein for release factors RF1 and RF2 during exponential growth. Furthermore, using available microarray data for gene expression, we show that in both fast growing and contrasting biofilm formation conditions, the relative level of RF1 is nicely correlated with the expression level of the genes ending with TAG.

Nonsense-mediated decay in genetic disease: friend or foe?

Eukaryotic cells utilize various RNA quality control mechanisms to ensure high fidelity of gene expression, thus protecting against the accumulation of nonfunctional RNA and the subsequent production of abnormal peptides. Messenger RNAs (mRNAs) are largely responsible for protein production, and mRNA quality control is particularly important for protecting the cell against the downstream effects of genetic mutations. Nonsense-mediated decay (NMD) is an evolutionarily conserved mRNA quality control system in all eukaryotes that degrades transcripts containing premature termination codons (PTCs). By degrading these aberrant transcripts, NMD acts to prevent the production of truncated proteins that could otherwise harm the cell through various insults, such as dominant negative effects or the ER stress response. Although NMD functions to protect the cell against the deleterious effects of aberrant mRNA, there is a growing body of evidence that mutation-, codon-, gene-, cell-, and tissue-specific differences in NMD efficiency can alter the underlying pathology of genetic disease. In addition, the protective role that NMD plays in genetic disease can undermine current therapeutic strategies aimed at increasing the production of full-length functional protein from genes harboring nonsense mutations. Here, we review the normal function of this RNA surveillance pathway and how it is regulated, provide current evidence for the role that it plays in modulating genetic disease phenotypes, and how NMD can be used as a therapeutic target.

Codon usage patterns in Chinese bayberry (Myrica rubra) based on RNA-Seq data

BACKGROUND

Codon usage analysis has been a classical topic for decades and has significances for studies of evolution, mRNA translation, and new gene discovery, etc. While the codon usage varies among different members of the plant kingdom, indicating the necessity for species-specific study, this work has mostly been limited to model organisms. Recently, the development of deep sequencing, especial RNA-Seq, has made it possible to carry out studies in non-model species.

RESULT

RNA-Seq data of Chinese bayberry was analyzed to investigate the bias of codon usage and codon pairs. High frequency codons (AGG, GCU, AAG and GAU), as well as low frequency ones (NCG and NUA codons) were identified, and 397 high frequency codon pairs were observed. Meanwhile, 26 preferred and 141 avoided neighboring codon pairs were also identified, which showed more significant bias than the same pairs with one or more intervening codons. Codon patterns were also analyzed at the plant kingdom, organism and gene levels. Changes during plant evolution were evident using RSCU (relative synonymous codon usage), which was even more significant than GC3s (GC content of 3rd synonymous codons). Nine GO categories were differentially and independently influenced by CAI (codon adaptation index) or GC3s, especially in 'Molecular function' category. Within a gene, the average CAI increased from 0.720 to 0.785 in the first 50 codons, and then more slowly thereafter. Furthermore, the preferred as well as avoided codons at the position just following the start codon AUG were identified and discussed in relation to the key positions in Kozak sequences.

CONCLUSION

A comprehensive codon usage Table and number of high-frequency codon pairs were established. Bias in codon usage as well as in neighboring codon pairs was observed, and the significance of this in avoiding DNA mutation, increasing protein production and regulating protein synthesis rate was proposed. Codon usage patterns at three levels were revealed and the significance in plant evolution analysis, gene function classification, and protein translation start site predication were discussed. This work promotes the study of codon biology, and provides some reference for analysis and comprehensive application of RNA-Seq data from other non-model species.

Pharmaceutical therapies to recode nonsense mutations in inherited diseases

Nonsense codons, generated from nonsense mutations or frameshifts, contribute significantly to the spectrum of inherited human diseases such as cystic fibrosis, Duchenne muscular dystrophy, hemophilia, spinal muscular atrophy, and many forms of cancer. The presence of a mutant nonsense codon results in premature termination to preclude the synthesis of a full-length protein and leads to aberrations in gene expression. Suppression therapy to recode a premature termination codon with an amino acid allowing readthrough to rescue the production of a full-length protein presents a promising strategy for treatment of patients suffering from debilitating nonsense-mediated disorders. Suppression therapy using aminoglycosides to promote readthrough in vitro have been known since the sixties. Recent progress in the field of recoding via pharmaceuticals has led to the continuous discovery and development of several pharmacological agents with nonsense suppression activities. Here, we review the mechanisms that are involved in discriminating normal versus premature termination codons, the factors involved in readthrough efficiency, the epidemiology of several well-known nonsense-mediated diseases, and the various pharmacological agents (aminoglycoside and non-aminoglycoside compounds) that are currently being employed in nonsense suppression therapy studies. We also discuss how these therapeutic agents can be used to regulate gene expression for gene therapy applications.

Bacterial phylogenetic tree construction based on genomic translation stop signals

Background

The efficiencies of the stop codons TAA, TAG, and TGA in protein synthesis termination are not the same. These variations could allow many genes to be regulated. There are many similar nucleotide trimers found on the second and third reading-frames of a gene. They are called premature stop codons (PSC). Like stop codons, the PSC in bacterial genomes are also highly bias in terms of their quantities and qualities on the genes. Phylogenetically related species often share a similar PSC profile. We want to know whether the selective forces that influence the stop codons and the PSC usage biases in a genome are related. We also wish to know how strong these trimers in a genome are related to the natural history of the bacterium. Knowing these relations may provide better knowledge in the phylogeny of bacteria

Results

A 16SrRNA-alignment tree of 19 well-studied α-, β- and γ-Proteobacteria Type species is used as standard reference for bacterial phylogeny. The genomes of sixty-one bacteria, belonging to the α-, β- and γ-Proteobacteria subphyla, are used for this study. The stop codons and PSC are collectively termed “Translation Stop Signals” (TSS). A gene is represented by nine scalars corresponding to the numbers of counts of TAA, TAG, and TGA on each of the three reading-frames of that gene. “Translation Stop Signals Ratio” (TSSR) is the ratio between the TSS counts. Four types of TSSR are investigated. The TSSR-1, TSSR-2 and TSSR-3 are each a 3-scalar series corresponding respectively to the average ratio of TAA: TAG: TGA on the first, second, and third reading-frames of all genes in a genome. The Genomic-TSSR is a 9-scalar series representing the ratio of distribution of all TSS on the three reading-frames of all genes in a genome. Results show that bacteria grouped by their similarities based on TSSR-1, TSSR-2, or TSSR-3 values could only partially resolve the phylogeny of the species. However, grouping bacteria based on thier Genomic-TSSR values resulted in clusters of bacteria identical to those bacterial clusters of the reference tree. Unlike the 16SrRNA method, the Genomic-TSSR tree is also able to separate closely related species/strains at high resolution. Species and strains separated by the Genomic-TSSR grouping method are often in good agreement with those classified by other taxonomic methods. Correspondence analysis of individual genes shows that most genes in a bacterial genome share a similar TSSR value. However, within a chromosome, the Genic-TSSR values of genes near the replication origin region (Ori) are more similar to each other than those genes near the terminus region (Ter).

Conclusion

The translation stop signals on the three reading-frames of the genes on a bacterial genome are interrelated, possibly due to frequent off-frame recombination facilitated by translational-associated recombination (TSR). However, TSR may not occur randomly in a bacterial chromosome. Genes near the Ori region are often highly expressed and a bacterium always maintains multiple copies of Ori. Frequent collisions between DNA- polymerase and RNA-polymerase would create many DNA strand-breaks on the genes; whereas DNA strand-break induced homologues-recombination is more likely to take place between genes with similar sequence. Thus, localized recombination could explain why the TSSR of genes near the Ori region are more similar to each other. The quantity and quality of these TSS in a genome strongly reflect the natural history of a bacterium. We propose that the Genomic- TSSR can be used as a subjective biomarker to represent the phyletic status of a bacterium.

Comparative genomic and transcriptomic characterization of the toxigenic marine dinoflagellate Alexandrium ostenfeldii

Many dinoflagellate species are notorious for the toxins they produce and ecological and human health consequences associated with harmful algal blooms (HABs). Dinoflagellates are particularly refractory to genomic analysis due to the enormous genome size, lack of knowledge about their DNA composition and structure, and peculiarities of gene regulation, such as spliced leader (SL) trans-splicing and mRNA transposition mechanisms. Alexandrium ostenfeldii is known to produce macrocyclic imine toxins, described as spirolides. We characterized the genome of A. ostenfeldii using a combination of transcriptomic data and random genomic clones for comparison with other dinoflagellates, particularly Alexandrium species. Examination of SL sequences revealed similar features as in other dinoflagellates, including Alexandrium species. SL sequences in decay indicate frequent retro-transposition of mRNA species. This probably contributes to overall genome complexity by generating additional gene copies. Sequencing of several thousand fosmid and bacterial artificial chromosome (BAC) ends yielded a wealth of simple repeats and tandemly repeated longer sequence stretches which we estimated to comprise more than half of the whole genome. Surprisingly, the repeats comprise a very limited set of 79–97 bp sequences; in part the genome is thus a relatively uniform sequence space interrupted by coding sequences. Our genomic sequence survey (GSS) represents the largest genomic data set of a dinoflagellate to date. Alexandrium ostenfeldii is a typical dinoflagellate with respect to its transcriptome and mRNA transposition but demonstrates Alexandrium-like stop codon usage. The large portion of repetitive sequences and the organization within the genome is in agreement with several other studies on dinoflagellates using different approaches. It remains to be determined whether this unusual composition is directly correlated to the exceptionally genome organization of dinoflagellates with a low amount of histones and histone-like proteins.

Transposable elements in disease-associated cryptic exons

Transposable elements (TEs) make up a half of the human genome, but the extent of their contribution to cryptic exon activation that results in genetic disease is unknown. Here, a comprehensive survey of 78 mutation-induced cryptic exons previously identified in 51 disease genes revealed the presence of TEs in 40 cases (51%). Most TE-containing exons were derived from short interspersed nuclear elements (SINEs), with Alus and mammalian interspersed repeats (MIRs) covering >18 and >16% of the exonized sequences, respectively. The majority of SINE-derived cryptic exons had splice sites at the same positions of the Alu/MIR consensus as existing SINE exons and their inclusion in the mRNA was facilitated by phylogenetically conserved changes that improved both traditional and auxiliary splicing signals, thus marking intronic TEs amenable for pathogenic exonization. The overrepresentation of MIRs among TE exons is likely to result from their high average exon inclusion levels, which reflect their strong splice sites, a lack of splicing silencers and a high density of enhancers, particularly (G)AA(G) motifs. These elements were markedly depleted in antisense Alu exons, had the most prominent position on the exon-intron gradient scale and are proposed to promote exon definition through enhanced tertiary RNA interactions involving unpaired (di)adenosines. The identification of common mechanisms by which the most dynamic parts of the genome contribute both to new exon creation and genetic disease will facilitate detection of intronic mutations and the development of computational tools that predict TE hot-spots of cryptic exon activation.

Evolution and differentiation of the prion protein gene (PRNP) among species

A total of 937 prion protein gene (PRNP) sequences belonging to 83 species in 56 genera of 26 families were analyzed in order to investigate its evolution and differentiation among species. The length of PRNP coding sequence for all species analyzed varied from 567 to 825 bp, which was mainly because of insertion or deletion in the repeat region within and among the species. TGA and TAG are the main stop codons in the PRNP gene. Bos taurus had the smallest variation in terms of average number of nucleotide differences (0.811), nucleotide diversity (0.0011), and nonsynonymous nucleotide diversity (0.0002) among all the ruminants. The reconstructed phylogenetic tree of PRNP of families and species was basically consistent with the taxonomy of National Center for Biotechnology Information except for Felidae (Felis catus), which was initially clustered with Moschidae rather than Mustelidae or Canidae.

A stochastic model of gene evolution with chaotic mutations

We develop here a new class of stochastic models of gene evolution in which the mutations are chaotic, i.e. a random subset of the 64 possible trinucleotides mutates at each evolutionary time t according to some substitution probabilities. Therefore, at each time t, the numbers and the types of mutable trinucleotides are unknown. Thus, the mutation matrix changes at each time t. The chaotic model developed generalizes the standard model in which all the trinucleotides mutate at each time t. It determines the occurrence probabilities at time t of trinucleotides which chaotically mutate according to three substitution parameters associated with the three trinucleotide sites. Two theorems prove that this chaotic model has a probability vector at each time t and that it converges to a uniform probability vector identical to that of the standard model. Furthermore, four applications of this chaotic model (with a uniform random strategy for the 64 trinucleotides and with a particular strategy for the three stop codons) allow an evolutionary study of the three circular codes identified in both eukaryotic and prokaryotic genes. A circular code is a particular set of trinucleotides whose main property is the retrieval of the frames in genes locally, i.e. anywhere in genes and particularly without start codons, and automatically with a window of a few nucleotides. After a certain evolutionary time and with particular values for the three substitution parameters, the chaotic models retrieve the main statistical properties of the three circular codes observed in genes. These applications also allow an evolutionary comparison between the standard and chaotic models.

Correlation of mRNA expression and protein abundance affected by multiple sequence features related to translational efficiency in Desulfovibrio vulgaris: a quantitative analysis

The modest correlation between mRNA expression and protein abundance in large-scale data sets is explained in part by experimental challenges, such as technological limitations, and in part by fundamental biological factors in the transcription and translation processes. Among various factors affecting the mRNA-protein correlation, the roles of biological factors related to translation are poorly understood. In this study, using experimental mRNA expression and protein abundance data collected from Desulfovibrio vulgaris by DNA microarray and liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) proteomic analysis, we quantitatively examined the effects of several translational-efficiency-related sequence features on mRNA-protein correlation. Three classes of sequence features were investigated according to different translational stages: (i) initiation, Shine-Dalgarno sequences, start codon identity, and start codon context; (ii) elongation, codon usage and amino acid usage; and (iii) termination, stop codon identity and stop codon context. Surprisingly, although it is widely accepted that translation initiation is the rate-limiting step for translation, our results showed that the mRNA-protein correlation was affected the most by the features at elongation stages, i.e., codon usage and amino acid composition (5.3-15.7% and 5.8-11.9% of the total variation of mRNA-protein correlation, respectively), followed by stop codon context and the Shine-Dalgarno sequence (3.7-5.1% and 1.9-3.8%, respectively). Taken together, all sequence features contributed to 15.2-26.2% of the total variation of mRNA-protein correlation. This study provides the first comprehensive quantitative analysis of the mRNA-protein correlation in bacterial D. vulgaris and adds new insights into the relative importance of various sequence features in prokaryotic protein translation.

Comparison of codon usage and tRNAs in mitochondrial genomes of Candida species

To gain insight into the nature of the mitochondrial genomes (mtDNA) of different Candida species, the synonymous codon usage bias of mitochondrial protein coding genes and the tRNAs in C. albicans, C. parapsilosis, C. stellata, C. glabrata and the closely related yeast Saccharomyces cerevisiae were analyzed. Common features of the mtDNA in Candida species are a strong A+T pressure on protein coding genes, and insufficient mitochondrial tRNA species are encoded to perform protein synthesis. The wobble site of the anticodon is always U for the NNR (NNA and NNG) codon families, which are dominated by A-ending codons, and always G for the NNY (NNC and NNU) codon families, which is dominated by U-ending codons, and always U for the NNN (NNA, NNU, NNC and NNG) codon families, which are dominated by A-ending codons and U-ending codons. Patterns of synonymous codon usage of Candida species can be classified into three groups: (1) optimal codon-anticodon usage, Glu, Lys, Leu (translated by anti-codon UAA), Gln, Arg (translated by anti-codon UCU) and Trp are containing NNR codons. NNA, whose corresponding tRNA is encoded in the mtDNA, is used preferentially. (2) Non-optimal codon-anticodon usage, Cys, Asp, Phe, His, Asn, Ser (translated by anti-codon GCU) and Tyr are containing NNY codons. The NNU codon, whose corresponding tRNA is not encoded in the mtDNA, is used preferentially. (3) Combined codon-anticodon usage, Ala, Gly, Leu (translated by anti-codon UAG), Pro, Ser (translated by anti-codon UGA), Thr and Val are containing NNN codons. NNA (tRNA encoded in the mtDNA) and NNU (tRNA not encoded in the mtDNA) are used preferentially. In conclusion, we propose that in Candida species, codons containing A or U at third position are used preferentially, regardless of whether corresponding tRNAs are encoded in the mtDNA. These results might be useful in understanding the common features of the mtDNA in Candida species and patterns of synonymous codon usage.

Frequent appearance of novel protein-coding sequences by frameshift translation

Genomic duplication, followed by divergence, contributes to organismal evolution. Several mechanisms, such as exon shuffling and alternative splicing, are responsible for novel gene functions, but they generate homologous domains and do not usually lead to drastic innovation. Major novelties can potentially be introduced by frameshift mutations and this idea can explain the creation of novel proteins. Here, we employ a strategy using simulated protein sequences and identify 470 human and 108 mouse frameshift events that originate new gene segments. No obvious interspecies overlap was observed, suggesting high rates of acquisition of evolutionary events. This inference is supported by a deficiency of TpA dinucleotides in the protein-coding sequences, which decreases the occurrence of translational termination, even on the complementary strand. Increased usage of the TGA codon as the termination signal in newer genes also supports our inference. This suggests that tolerated frameshift changes are a prevalent mechanism for the rapid emergence of new genes and that protein-coding sequences can be derived from existing or ancestral exons rather than from events that result in noncoding sequences becoming exons.