Specific correlations between relative synonymous codon usage and protein secondary structure

The first complete chloroplast genome of Thalictrum fargesii: insights into phylogeny and species identification

Introduction

Thalictrum fargesii is a medicinal plant belonging to the genus Thalictrum of the Ranunculaceae family and has been used in herbal medicine in the Himalayan regions of China and India. This species is taxonomically challenging because of its morphological similarities to other species within the genus. Thus, herbal drugs from this species are frequently adulterated, substituted, or mixed with other species, thereby endangering consumer safety.

Methods

The present study aimed to sequence and assemble the entire chloroplast (cp) genome of T. fargesii using the Illumina HiSeq 2500 platform to better understand the genomic architecture, gene composition, and phylogenetic relationships within the Thalictrum.

Results and discussion

The cp genome was 155,929 bp long and contained large single-copy (85,395 bp) and small single-copy (17,576 bp) regions that were segregated by a pair of inverted repeat regions (26,479 bp) to form a quadripartite structure. The cp genome contains 133 genes, including 88 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. Additionally, this genome contains 64 codons that encode 20 amino acids, the most preferred of which are alanine and leucine. We identified 68 SSRs, 27 long repeats, and 242 high-confidence C-to-U RNA-editing sites in the cp genome. Moreover, we discovered seven divergent hotspot regions in the cp genome of T. fargesii, among which ndhD-psaC and rpl16-rps3 may be useful for developing molecular markers for identifying ethnodrug species and their contaminants. A comparative study with eight other species in the genus revealed that pafI and rps19 had highly variable sites in the cp genome of T. fargesii. Additionally, two special features, (i) the shortest length of the ycf1 gene at the IRA-SSC boundary and (ii) the distance between the rps19 fragment and trnH at the IRA-LSC junction, distinguish the cp genome of T. fargesii from those of other species within the genus. Furthermore, phylogenetic analysis revealed that T. fargesii was closely related to T. tenue and T. petaloidium.

Conclusion

Considering all these lines of evidence, our findings offer crucial molecular and evolutionary information that could play a significant role in further species identification, evolution, and phylogenetic studies on T. fargesii.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

Expression of blue pigment synthetase a from Streptomyces lavenduale reveals insights on the effects of refactoring biosynthetic megasynthases for heterologous expression in Escherichia coli

High GC bacteria from the genus Streptomyces harbor expansive secondary metabolism. The expression of biosynthetic proteins and the characterization and identification of biological "parts" for synthetic biology purposes from such pathways are of interest. However, the high GC content of proteins from actinomycetes in addition to the large size and multi-domain architecture of many biosynthetic proteins (such as non-ribosomal peptide synthetases; NRPSs, and polyketide synthases; PKSs often called "megasynthases") often presents issues with full-length translation and folding. Here we evaluate a non-ribosomal peptide synthetase (NRPS) from Streptomyces lavenduale, a multidomain "megasynthase" gene that comes from a high GC (72.5%) genome. While a preliminary step in revealing differences, to our knowledge this presents the first head-to-head comparison of codon-optimized sequences versus a native sequence of proteins of streptomycete origin heterologously expressed in E. coli. We found that any disruption in co-translational folding from codon mismatch that reduces the titer of indigoidine is explainable via the formation of more inclusion bodies as opposed to compromising folding or posttranslational modification in the soluble fraction. This result supports that one could apply any refactoring strategies that improve soluble expression in E. coli without concern that the protein that reaches the soluble fraction is differentially folded.

The Influence of the Selection at the Amino Acid Level on Synonymous Codon Usage from the Viewpoint of Alternative Genetic Codes

Synonymous codon usage can be influenced by mutations and/or selection, e.g., for speed of protein translation and correct folding. However, this codon bias can also be affected by a general selection at the amino acid level due to differences in the acceptance of the loss and generation of these codons. To assess the importance of this effect, we constructed a mutation-selection model model, in which we generated almost 90,000 stationary nucleotide distributions produced by mutational processes and applied a selection based on differences in physicochemical properties of amino acids. Under these conditions, we calculated the usage of fourfold degenerated (4FD) codons and compared it with the usage characteristic of the pure mutations. We considered both the standard genetic code (SGC) and alternative genetic codes (AGCs). The analyses showed that a majority of AGCs produced a greater 4FD codon bias than the SGC. The mutations producing more thymine or adenine than guanine and cytosine increased the differences in usage. On the other hand, the mutational pressures generating a lot of cytosine or guanine with a low content of adenine and thymine decreased this bias because the nucleotide content of most 4FD codons stayed in the compositional equilibrium with these pressures. The comparison of the theoretical results with those for real protein coding sequences showed that the influence of selection at the amino acid level on the synonymous codon usage cannot be neglected. The analyses indicate that the effect of amino acid selection cannot be disregarded and that it can interfere with other selection factors influencing codon usage, especially in AT-rich genomes, in which AGCs are usually used.

Comparative phylogenetic analysis of oolong tea (Phoenix Dancong tea) using complete chloroplast genome sequences

Phoenix Dancong tea, a variety of oolong tea, is produced in Chaozhou, Guangdong Province, China, and is characterized by numerous hybridizations and polyploidization. To assess the genetic diversity and phylogenetic relationships among Phoenix Dancong tea and other oolong teas, an integrated circular chloroplast genome was constructed for thirty species of Phoenix Dancong tea from Chaozhou. The genome of Phoenix dancong tea is a circular molecule of 157,041-157,137 bp, with a pair of inverted repeats (26,072-26,610 bp each) separated by a large single copy (86,615-86,658 bp) and small single copy (18,264-18,284 bp). A total of 135 unique genes were encoded, including 90 protein coding genes, 37 tRNAs and 8 rRNAs. A comparative analysis with the other seven species in the oolong tea family that have been sequenced to date revealed similarities in structural organization, gene content and arrangement. Repeated sequence analysis identified 17-23 tandem repeats, 20-24 forward repeats and 25-27 palindromic repeats. Additionally, a total of 65-70 simple sequence repeats were detected, with mononucleotide repeats being the most common. Phylogenetic analyses showed that Phoenix Dancong tea and Fujian oolong tea were clustered with other cultivated Camellia sinensis in the genus Camellia of the family Theaceae, while the two oolong tea species were relatively independently cross-embedded in the genus, Camellia. Close genetic relationships were observed between Phoenix Dancong tea and other oolong teas, and the overall chloroplast genomes of oolong tea showed patterns with low variations and conserved evolution. The availability of Phoenix Dancong tea chloroplast genomes not only elucidated the relationship among oolong teas from different origins in Guangdong and Fujian but also provided valuable genetic resources to assist further molecular studies on the taxonomic and phylogenomic resolution of the genus Camellia.

An investigation of codon usage pattern analysis in pancreatitis associated genes

BACKGROUND

Pancreatitis is an inflammatory disorder resulting from the autoactivation of trypsinogen in the pancreas. The genetic basis of the disease is an old phenomenon, and evidence is accumulating for the involvement of synonymous/non-synonymous codon variants in disease initiation and progression.

RESULTS

The present study envisaged a panel of 26 genes involved in pancreatitis for their codon choices, compositional analysis, relative dinucleotide frequency, nucleotide disproportion, protein physical properties, gene expression, codon bias, and interrelated of all these factors. In this set of genes, gene length was positively correlated with nucleotide skews and codon usage bias. Codon usage of any gene is dependent upon its AT and GC component; however, AGG, CGT, and CGA encoding for Arg, TCG for Ser, GTC for Val, and CCA for Pro were independent of nucleotide compositions. In addition, Codon GTC showed a correlation with protein properties, isoelectric point, instability index, and frequency of basic amino acids. We also investigated the effect of various evolutionary forces in shaping the codon usage choices of genes.

CONCLUSIONS

This study will enable us to gain insight into the molecular signatures associated with the disease that might help identify more potential genes contributing to enhanced risk for pancreatitis. All the genes associated with pancreatitis are generally associated with physiological function, and mutations causing loss of function, over or under expression leads to an ailment. Therefore, the present study attempts to envisage the molecular signature in a group of genes that lead to pancreatitis in case of malfunction.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

Biosynthetic Protein Folding and Molecular Chaperons

The problem of linear polypeptide chain folding into a unique tertiary structure is one of the fundamental scientific challenges. The process of folding cannot be fully understood without its biological context, especially for big multidomain and multisubunit proteins. The principal features of biosynthetic folding are co-translational folding of growing nascent polypeptide chains and involvement of molecular chaperones in the process. The review summarizes available data on the early events of nascent chain folding, as well as on later advanced steps, including formation of elements of native structure. The relationship between the non-uniformity of translation rate and folding of the growing polypeptide is discussed. The results of studies on the effect of biosynthetic folding features on the parameters of folding as a physical process, its kinetics and mechanisms, are presented. Current understanding and hypotheses on the relationship of biosynthetic folding with the fundamental physical parameters and current views on polypeptide folding in the context of energy landscapes are discussed.

Complete Chloroplast Genome Sequence of Fagus longipetiolata Seemen (Fagaceae): Genome Structure, Adaptive Evolution, and Phylogenetic Relationships

Fagus longipetiolata Seemen is a deciduous tree of the Fagus genus in Fagaceae, which is endemic to China. In this study, we successfully sequenced the cp genome of F. longipetiolata, compared the cp genomes of the Fagus genus, and reconstructed the phylogeny of Fagaceae. The results showed that the cp genome of F. longipetiolata was 158,350 bp, including a pair of inverted repeat (IRA and IRB) regions with a length of 25,894 bp each, a large single-copy (LSC) region of 87,671 bp, and a small single-copy (SSC) region of 18,891 bp. The genome encoded 131 unique genes, including 81 protein-coding genes, 37 transfer RNA genes (tRNAs), 8 ribosomal RNA genes (rRNAs), and 5 pseudogenes. In addition, 33 codons and 258 simple sequence repeats (SSRs) were identified. The cp genomes of Fagus were relatively conserved, especially the IR regions, which showed the best conservation, and no inversions or rearrangements were found. The five regions with the largest variations were the rps12, rpl32, ccsA, trnW-CCA, and rps3 genes, which spread over in LSC and SSC. The comparison of gene selection pressure indicated that purifying selection was the main selective pattern maintaining important biological functions in Fagus cp genomes. However, the ndhD, rpoA, and ndhF genes of F. longipetiolata were affected by positive selection. Phylogenetic analysis revealed that F. longipetiolata and F. engleriana formed a close relationship, which partially overlapped in their distribution in China. Our analysis of the cp genome of F. longipetiolata would provide important genetic information for further research into the classification, phylogeny and evolution of Fagus.

Glue Genes Are Subjected to Diverse Selective Forces during Drosophila Development

Molecular evolutionary studies usually focus on genes with clear roles in adult fitness or on developmental genes expressed at multiple time points during the life of the organism. Here, we examine the evolutionary dynamics of Drosophila glue genes, a set of eight genes tasked with a singular primary function during a specific developmental stage: the production of glue that allows animal pupa to attach to a substrate for several days during metamorphosis. Using phenotypic assays and available data from transcriptomics, PacBio genomes, and sequence variation from global populations, we explore the selective forces acting on glue genes within the cosmopolitan Drosophila melanogaster species and its five closely related species, D. simulans, D. sechellia, D. mauritiana, D. yakuba, and D. teissieri. We observe a three-fold difference in glue adhesion between the least and the most adhesive D. melanogaster strain, indicating a strong genetic component to phenotypic variation. These eight glue genes are among the most highly expressed genes in salivary glands yet they display no notable codon bias. New copies of Sgs3 and Sgs7 are found in D. yakuba and D. teissieri with the Sgs3 coding sequence evolving rapidly after duplication in the D. yakuba branch. Multiple sites along the various glue genes appear to be constrained. Our population genetics analysis in D. melanogaster suggests signals of local adaptive evolution for Sgs3, Sgs5, and Sgs5bis and traces of selective sweeps for Sgs1, Sgs3, Sgs7, and Sgs8. Our work shows that stage-specific genes can be subjected to various dynamic evolutionary forces.

A Code Within a Code: How Codons Fine-Tune Protein Folding in the Cell

The genetic code sets the correspondence between the sequence of a given nucleotide triplet in an mRNA molecule, called a codon, and the amino acid that is added to the growing polypeptide chain during protein synthesis. With four bases (A, G, U, and C), there are 64 possible triplet codons: 61 sense codons (encoding amino acids) and 3 nonsense codons (so-called, stop codons that define termination of translation). In most organisms, there are 20 common/standard amino acids used in protein synthesis; thus, the genetic code is redundant with most amino acids (with the exception of Met and Trp) are being encoded by more than one (synonymous) codon. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in mRNA suggested that the specific codon choice might have functional implications beyond coding for amino acid. Observation of nonequivalent use of codons in mRNAs implied a possibility of the existence of auxiliary information in the genetic code. Indeed, it has been found that genetic code contains several layers of such additional information and that synonymous codons are strategically placed within mRNAs to ensure a particular translation kinetics facilitating and fine-tuning co-translational protein folding in the cell via step-wise/sequential structuring of distinct regions of the polypeptide chain emerging from the ribosome at different points in time. This review summarizes key findings in the field that have identified the role of synonymous codons and their usage in protein folding in the cell.

Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB

Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite-induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duodenalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G. duodenalis CUB formation. This study systematically analyzes G. duodenalis codon usage pattern and clarifies the mechanisms of G. duodenalis CUB. These results will be very useful to identify new genes, molecular genetic manipulation, and study of G. duodenalis evolution.

Are Nonsynonymous Transversions Generally More Deleterious than Nonsynonymous Transitions?

It has been suggested that, due to the structure of the genetic code, nonsynonymous transitions are less likely than transversions to cause radical changes in amino acid physicochemical properties so are on average less deleterious. This view was supported by some but not all mutagenesis experiments. Because laboratory measures of fitness effects have limited sensitivities and relative frequencies of different mutations in mutagenesis studies may not match those in nature, we here revisit this issue using comparative genomics. We extend the standard codon model of sequence evolution by adding the parameter η that quantifies the ratio of the fixation probability of transitional nonsynonymous mutations to that of transversional nonsynonymous mutations. We then estimate η from the concatenated alignment of all protein-coding DNA sequences of two closely related genomes. Surprisingly, η ranges from 0.13 to 2.0 across 90 species pairs sampled from the tree of life, with 51 incidences of η < 1 and 30 incidences of η >1 that are statistically significant. Hence, whether nonsynonymous transversions are overall more deleterious than nonsynonymous transitions is species-dependent. Because the corresponding groups of amino acid replacements differ between nonsynonymous transitions and transversions, η is influenced by the relative exchangeabilities of amino acid pairs. Indeed, an extensive search reveals that the large variation in η is primarily explainable by the recently reported among-species disparity in amino acid exchangeabilities. These findings demonstrate that genome-wide nucleotide substitution patterns in coding sequences have species-specific features and are more variable among evolutionary lineages than are currently thought.

Some theoretical aspects of reprogramming the standard genetic code

Reprogramming of the standard genetic code to include non-canonical amino acids (ncAAs) opens new prospects for medicine, industry, and biotechnology. There are several methods of code engineering, which allow us for storing new genetic information in DNA sequences and producing proteins with new properties. Here, we provided a theoretical background for the optimal genetic code expansion, which may find application in the experimental design of the genetic code. We assumed that the expanded genetic code includes both canonical and non-canonical information stored in 64 classical codons. What is more, the new coding system is robust to point mutations and minimizes the possibility of reversion from the new to old information. In order to find such codes, we applied graph theory to analyze the properties of optimal codon sets. We presented the formal procedure in finding the optimal codes with various number of vacant codons that could be assigned to new amino acids. Finally, we discussed the optimal number of the newly incorporated ncAAs and also the optimal size of codon groups that can be assigned to ncAAs.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.

Network analysis of synonymous codon usage

MOTIVATION

Most amino acids are encoded by multiple synonymous codons, some of which are used more rarely than others. Analyses of positions of such rare codons in protein sequences revealed that rare codons can impact co-translational protein folding and that positions of some rare codons are evolutionarily conserved. Analyses of their positions in protein 3-dimensional structures, which are richer in biochemical information than sequences alone, might further explain the role of rare codons in protein folding.

RESULTS

We model protein structures as networks and use network centrality to measure the structural position of an amino acid. We first validate that amino acids buried within the structural core are network-central, and those on the surface are not. Then, we study potential differences between network centralities and thus structural positions of amino acids encoded by conserved rare, non-conserved rare and commonly used codons. We find that in 84% of proteins, the three codon categories occupy significantly different structural positions. We examine protein groups showing different codon centrality trends, i.e. different relationships between structural positions of the three codon categories. We see several cases of all proteins from our data with some structural or functional property being in the same group. Also, we see a case of all proteins in some group having the same property. Our work shows that codon usage is linked to the final protein structure and thus possibly to co-translational protein folding.

AVAILABILITY AND IMPLEMENTATION

https://nd.edu/∼cone/CodonUsage/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Plastid Phylogenomics of Dendroseris (Cichorieae; Asteraceae): Insights Into Structural Organization and Molecular Evolution of an Endemic Lineage From the Juan Fernández Islands

Dendroseris D. Don comprises 11 species endemic to the Juan Fernández islands in Chile. They demonstrate spectacular and unusual growth forms of rosette trees with extremely variable morphology and occupy wide ecological ranges on the islands. These unique plants are now highly threatened with extinction with very small population sizes, typically consisting of 10 or fewer individuals in wild. Despite morphological and ecological divergence among species of Dendroseris, their monophyly has been supported in previous studies, but with little resolution among subgeneric groups. We assembled seven complete plastome sequences from seven species of Dendroseris, including representatives from three subgenera, and carried out comparative phylogenomic analyses. The plastomes are highly conserved in gene content and order, with size ranging from 152,199 to 152,619 bp and containing 130 genes (87 coding genes, 6 rRNA genes, and 37 tRNA genes). Plastid phylogenomic analyses based on both the complete plastome sequences and 81 concatenated coding genes only show Dendroseris nested within Sonchus sensu lato, and also that inter-subgeneric relationships are fully resolved. Subg. Phoenicoseris is resolved as sister to the remaining species of the genus and a sister relationship between the two subgenera Dendroseris and Rea. Ten mutation hotspots from LSC and SSC regions and variable SSRs are identified as potential chloroplast markers for future phylogenetic and phylogeographic studies of Sonchus and related groups.

A code within the genetic code: codon usage regulates co-translational protein folding

The genetic code is degenerate, and most amino acids are encoded by two to six synonymous codons. Codon usage bias, the preference for certain synonymous codons, is a universal feature of all genomes examined. Synonymous codon mutations were previously thought to be silent; however, a growing body evidence now shows that codon usage regulates protein structure and gene expression through effects on co-translational protein folding, translation efficiency and accuracy, mRNA stability, and transcription. Codon usage regulates the speed of translation elongation, resulting in non-uniform ribosome decoding rates on mRNAs during translation that is adapted to co-translational protein folding process. Biochemical and genetic evidence demonstrate that codon usage plays an important role in regulating protein folding and function in both prokaryotic and eukaryotic organisms. Certain protein structural types are more sensitive than others to the effects of codon usage on protein folding, and predicted intrinsically disordered domains are more prone to misfolding caused by codon usage changes than other domain types. Bioinformatic analyses revealed that gene codon usage correlates with different protein structures in diverse organisms, indicating the existence of a codon usage code for co-translational protein folding. This review focuses on recent literature on the role and mechanism of codon usage in regulating translation kinetics and co-translational protein folding.

Synthetic Biology: A Solution for Tackling Nanomaterial Challenges

Bioengineers have mastered practical techniques for tuning a biomaterial's properties with only limited information on the relationship between the material's structure and function. These techniques have been quintessential to engineering proteins, which are most often riddled with ill-defined structure-function relationships. In this Perspective, we review bioengineering approaches aimed at overcoming the elusive protein structure-function relation. We extend these principles to engineering synthetic nanomaterials, specifically applying the underlying theory to optical sensors based on single-stranded DNA-wrapped single-walled carbon nanotubes (ssDNA-SWCNTs). Bioengineering techniques such as directed evolution, computational design, and noncanonical synthesis are reviewed in the broader context of nanomaterials engineering. We further provide an order-of-magnitude analysis of empirical approaches that rely on random or guided searches for designing new nanomaterials. The underlying concepts presented in these approaches can be further extended to a broad range of engineering fields confronted with empirical design strategies, including catalysis, metal-organic frameworks (MOFs), pharmaceutical dosing, and optimization algorithms.

The Complete Chloroplast Genome Sequencing and Comparative Analysis of Reed Canary Grass (Phalaris arundinacea) and Hardinggrass (P. aquatica)

There are 22 species in the Phalaris genera that distribute almost all over the temperate regions of the world. Among them, reed canary grass (Phalaris arundinacea, tetraploid and hexaploid) and hardinggrass (P. aquatica, tetraploid) have been long cultivated as forage grass and have received attention as bio-energy materials in recent years. We aimed to facilitate inter-species/ploidies comparisons, and to illuminate the degree of sequence variation within existing gene pools, chloroplast (cp) genomes of three Phalaris cytotypes (P. aquatica/4x, P. arundinacea/4x and P. arundinacea/6x) were sequenced and assembled. The result indicated that certain sequence variations existed between the cp genomes of P. arundinacea and P. aquatica. Several hotspot regions (atpI~atpH, trnT-UGU~ndhJ, rbcL~psaI, and ndhF~rpl32) were found, and variable genes (infA, psaI, psbK, etc.) were detected. SNPs (single nucleotide polymorphisms) and/or indels (insertions and deletions) were confirmed by the high Ka/Ks and Pi value. Furthermore, distribution and presence of cp simple sequence repeats (cpSSRs) were identified in the three Phalaris cp genomes, although little difference was found between hexaploid and tetraploid P. arundinacea, and no rearrangement was detected among the three Phalaris cp genomes. The evolutionary relationship and divergent time among these species were discussed. The RNA-seq revealed several differentially expressed genes (DEGs), among which psaA, psaB, and psbB related to leaf color were further verified by leaf color differences.

Differences in Codon Usage Bias between Photosynthesis-Related Genes and Genetic System-Related Genes of Chloroplast Genomes in Cultivated and Wild Solanum Species

Solanum is one of the largest genera, including two important crops-potato (Solanum tuberosum) and tomato (Solanum lycopersicum). In this study we compared the chloroplast codon usage bias (CUB) among 12 Solanum species, between photosynthesis-related genes (Photo-genes) and genetic system-related genes (Genet-genes), and between cultivated species and wild relatives. The Photo-genes encode proteins for photosystems, the photosynthetic electron transport chain, and RuBisCO, while the Genet-genes encode proteins for ribosomal subunits, RNA polymerases, and maturases. The following findings about the Solanum chloroplast genome CUB were obtained: (1) the nucleotide composition, gene expression, and selective pressure are identified as the main factors affecting chloroplast CUB; (2) all these 12 chloroplast genomes prefer A/U over G/C and pyrimidines over purines at the third-base of codons; (3) Photo-genes have higher codon adaptation indexes than Genet-genes, indicative of a higher gene expression level and a stronger adaptation of Photo-genes; (4) gene function is the primary factor affecting CUB of Photo-genes but not Genet-genes; (5) Photo-genes prefer pyrimidine over purine, whereas Genet-genes favor purine over pyrimidine, at the third position of codons; (6) Photo-genes are mainly affected by the selective pressure, whereas Genet-genes are under the underlying mutational bias; (7) S. tuberosum is more similar with Solanum commersonii than with Solanum bulbocastanum; (8) S. lycopersicum is greatly different from the analyzed seven wild relatives; (9) the CUB in codons for valine, aspartic acid, and threonine are the same between the two crop species, S. tuberosum and S. lycopersicum. These findings suggest that the chloroplast CUB contributed to the differential requirement of gene expression activity and function between Photo-genes and Genet-genes and to the performance of cultivated potato and tomato.

Genome-Wide Quantification of the Effect of Gene Overexpression on Escherichia coli Growth

Recombinant protein production plays an essential role in both biological studies and pharmaceutical production. Escherichia coli is one of the most favorable hosts for this purpose. Although a number of strategies for optimizing protein production have been developed, the effect of gene overexpression on host cell growth has been much less studied. Here, we performed high-throughput tests on the E. coli a complete set of E. coli K-12 ORF archive (ASKA) collection to quantify the effects of overexpressing individual E. coli genes on its growth. The results indicated that overexpressing membrane-associated proteins or proteins with high abundances of branched-chain amino acids tended to impair cell growth, the latter of which could be remedied by amino acid supplementation. Through this study, we expect to provide an index for a fast pre-study estimate of host cell growth in order to choose proper rescuing approaches when working with different proteins.

Evolution of Synonymous Codon Usage Bias in West African and Central African Strains of Monkeypox Virus

The evolution of bias in synonymous codon usage in chosen monkeypox viral genomes and the factors influencing its diversification have not been reported so far. In this study, various trends associated with synonymous codon usage in chosen monkeypox viral genomes were investigated, and the results are reported. Identification of factors that influence codon usage in chosen monkeypox viral genomes was done using various codon usage indices, such as the relative synonymous codon usage, the effective number of codons, and the codon adaptation index. The Spearman rank correlation analysis and a correspondence analysis were used for correlating various factors with codon usage. The results revealed that mutational pressure due to compositional constraints, gene expression level, and selection at the codon level for utilization of putative optimal codons are major factors influencing synonymous codon usage bias in monkeypox viral genomes. A cluster analysis of relative synonymous codon usage values revealed a grouping of more virulent strains as one major cluster (Central African strains) and a grouping of less virulent strains (West African strains) as another major cluster, indicating a relationship between virulence and synonymous codon usage bias. This study concluded that a balance between the mutational pressure acting at the base composition level and the selection pressure acting at the amino acid level frames synonymous codon usage bias in the chosen monkeypox viruses. The natural selection from the host does not seem to have influenced the synonymous codon usage bias in the analyzed monkeypox viral genomes.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

Genome-wide codon usage profiling of ocular infective Chlamydia trachomatis serovars and drug target identification

Chlamydia trachomatis (C.t) is a Gram-negative obligate intracellular bacteria and is a major causative of infectious blindness and sexually transmitted diseases. Among the varied serovars of this organism, A, B and C are reported as prominent ocular pathogens. Genomic studies of these strains shall aid in deciphering potential drug targets and genomic influence on pathogenesis. Hence, in this study we performed deep statistical profiling of codon usage in these serovars. The overall base composition analysis reveals that these serovars are over biased to AU than GC. Similarly, relative synonymous codon usage also showed preference towards A/U ending codons. Parity Rule 2 analysis inferred unequal distribution of AT and GC, indicative of other unknown factors acting along with mutational pressure to influence codon usage bias (CUB). Moreover, absolute quantification of CUB also revealed lower bias across these serovars. The effect of natural selection on CUB was also confirmed by neutrality plot, reinforcing natural selection under mutational pressure turned to be a pivotal role in shaping the CUB in the strains studied. Correspondence analysis (COA) clarified that, C.t C/TW-3 to show a unique trend in codon usage variation. Host influence analysis on shaping the codon usage pattern also inferred some speculative relativity. In a nutshell, our finding suggests that mutational pressure is the dominating factor in shaping CUB in the strains studied, followed by natural selection. We also propose potential drug targets based on cumulative analysis of strand bias, CUB and human non-homologue screening.

Equalizing the information amounts of protein and mRNA by information theory

Based on the Shannon's information communication theory, information amount of the entire length of a polymeric macromolecule can be calculated in bits through adding the entropies of each building block. Proteins, DNA and RNA are such macromolecules. When only the building blocks' variation is considered as the source of entropy, there is seemingly lower information in case of the protein if this approach is applied directly on a protein of specific size and the coding sequence size of the mRNA corresponding to the particular length of the protein. This decrease in the information amount seems contradictory but this apparent conflict is resolved by considering the conformational variations in proteins as a new variable in the calculation and balancing the approximated entropy of the coding part of the mRNA and the protein. Probabilities can change therefore we also assigned hypothetical probabilities to the conformational states, which represent the uneven distribution as the time spent in one conformation, providing the probability of the presence in either or one of the possible conformations. Results that are obtained by using hypothetical probabilities are in line with the experimental values of variations in the conformational-state of protein populations. This equalization approach has further biological relevance that it compensates for the degeneracy in the codon usage during protein translation and it leads to the conclusion that the alphabet size for the protein is rather optimal for the proper protein functioning within the thermodynamic milieu of the cell. The findings were also discussed in relation to the codon bias and have implications in relation to the codon evolution concept. Eventually, this work brings the fields of protein structural studies and molecular protein translation processes together with a novel approach.

Analysis of synonymous codon usage bias and phylogeny of coat protein gene in banana bract mosaic virus isolates

The banana is one of the world's most important livelihood crops. Banana plants are principally infected by four virus species, Banana bunchy top virus (genus Babuvirus), Cucumber mosaic virus (genus Cucumovirus), Banana streak virus (genus Badnavirus) and Banana bract mosaic virus (genus Potyvirus). The objective of this study is to understand the codon usage pattern and phylogeny of coat protein gene in different banana bract mosaic virus (BBrMV) isolates. The BBrMV Coat Protein (CP) gene was amplified from BBrMV infected banana plant samples collected from different districts of Tamil Nadu and Karnataka, India. Six new BBrMV isolates were submitted to National Center for Biotechnology Information. Phylogenetic analysis and codon usage indices were studied along with other isolates of BBrMV. Phylogenetic analysis of CP genes shows that most of BBrMV isolates are closely related to each other except KF385484.1 and KF385478.1. Relative codon usage patterns among different BBrMV isolates were calculated by software CodonW version 1.4.2. In BBrMV, codons with A-ended or U ended are the most preferential except the Leu and Gln whose optimized codons are CAG and UUG ending by G. The codon usage patterns of BBrMV isolates are principally influenced by mutational bias; however, compositional constraints along with mutational bias also play a major role.

The Impact of Selection at the Amino Acid Level on the Usage of Synonymous Codons

There are two main forces that affect usage of synonymous codons: directional mutational pressure and selection. The effectiveness of protein translation is usually considered as the main selectional factor. However, biased codon usage can also be a byproduct of a general selection at the amino acid level interacting with nucleotide replacements. To evaluate the validity and strength of such an effect, we superimposed >3.5 billion unrestricted mutational processes on the selection of nonsynonymous substitutions based on the differences in physicochemical properties of the coded amino acids. Using a modified evolutionary optimization algorithm, we determined the conditions in which the effect on the relative codon usage is maximized. We found that the effect is enhanced by mutational processes generating more adenine and thymine than guanine and cytosine, as well as more purines than pyrimidines. Interestingly, this effect is observed only under an unrestricted model of nucleotide substitution, and disappears when the mutational process is time-reversible. Comparison of the simulation results with data for real protein coding sequences indicates that the impact of selection at the amino acid level on synonymous codon usage cannot be neglected. Furthermore, it can considerably interfere, especially in AT-rich genomes, with other selections on codon usage, e.g., translational efficiency. It may also lead to difficulties in the recognition of other effects influencing codon bias, and an overestimation of protein coding sequences whose codon usage is subjected to adaptational selection.

Role of mRNA structure in the control of protein folding

Specific structures in mRNA modulate translation rate and thus can affect protein folding. Using the protein structures from two eukaryotes and three prokaryotes, we explore the connections between the protein compactness, inferred from solvent accessibility, and mRNA structure, inferred from mRNA folding energy (ΔG). In both prokaryotes and eukaryotes, the ΔG value of the most stable 30 nucleotide segment of the mRNA (ΔGmin) strongly, positively correlates with protein solvent accessibility. Thus, mRNAs containing exceptionally stable secondary structure elements typically encode compact proteins. The correlations between ΔG and protein compactness are much more pronounced in predicted ordered parts of proteins compared to the predicted disordered parts, indicative of an important role of mRNA secondary structure elements in the control of protein folding. Additionally, ΔG correlates with the mRNA length and the evolutionary rate of synonymous positions. The correlations are partially independent and were used to construct multiple regression models which explain about half of the variance of protein solvent accessibility. These findings suggest a model in which the mRNA structure, particularly exceptionally stable RNA structural elements, act as gauges of protein co-translational folding by reducing ribosome speed when the nascent peptide needs time to form and optimize the core structure.

Comparative Analysis of Codon Usage Bias Patterns in Microsporidian Genomes

The sub-3 Mbp genomes from microsporidian species of the Encephalitozoon genus are the smallest known among eukaryotes and paragons of genomic reduction and compaction in parasites. However, their diminutive stature is not characteristic of all Microsporidia, whose genome sizes vary by an order of magnitude. This large variability suggests that different evolutionary forces are applied on the group as a whole. In this study, we have compared the codon usage bias (CUB) between eight taxonomically distinct microsporidian genomes: Encephalitozoon intestinalis, Encephalitozoon cuniculi, Spraguea lophii, Trachipleistophora hominis, Enterocytozoon bieneusi, Nematocida parisii, Nosema bombycis and Nosema ceranae. While the CUB was found to be weak in all eight Microsporidia, nearly all (98%) of the optimal codons in S. lophii, T. hominis, E. bieneusi, N. parisii, N. bombycis and N. ceranae are fond of A/U in third position whereas most (64.6%) optimal codons in the Encephalitozoon species E. intestinalis and E. cuniculi are biased towards G/C. Although nucleotide composition biases are likely the main factor driving the CUB in Microsporidia according to correlation analyses, directed mutational pressure also likely affects the CUB as suggested by ENc-plots, correspondence and neutrality analyses. Overall, the Encephalitozoon genomes were found to be markedly different from the other microsporidians and, despite being the first sequenced representatives of this lineage, are uncharacteristic of the group as a whole. The disparities observed cannot be attributed solely to differences in host specificity and we hypothesize that other forces are at play in the lineage leading to Encephalitozoon species.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

Mass spectrometry approach and ELISA reveal the effect of codon optimization on N-linked glycosylation of HIV-1 gp120

The genes encoding many viral proteins such as HIV-1 envelope glycoprotein gp120 have a tendency for codons that are poorly used by the human genome. Why these codons are frequently present in the HIV genome is not known. The presence of these codons limits expression of HIV-1 gp120 for biochemical studies. The poor codons are replaced by synonymous codons that are frequently present in the highly expressed human genes to overexpress this protein. Whether this codon optimization affects functional properties of gp120 such as its N-linked glycosylation is unknown. We applied a bottom-up mass-spectrometry-based workflow for the direct measurement of deglycosylated and unglycosylated peptides with putative N-linked glycosylation sites, that is, NxS/T motifs. Using this mass-spectrometry approach in combination with ELISA, it is found that codon optimization significantly reduces the frequency with which the dolichol pyrophosphate-linked oligosaccharide is added by the catalytic subunits of oligosaccharide transferase complex to the glycosylation sites. This reduction affects binding of glycan-dependent broadly neutralizing antibodies. These data are essential for biochemical studies of gp120 and successful development of a vaccine against HIV-1. Furthermore, they demonstrate a mass-spectrometry approach for studying the site-specific N-linked glycosylation efficiency of glycoproteins.

Different types of secondary information in the genetic code

Whole-genome and functional analyses suggest a wealth of secondary or auxiliary genetic information (AGI) within the redundancy component of the genetic code. Although there are multiple aspects of biased codon use, we focus on two types of auxiliary information: codon-specific translational pauses that can be used by particular proteins toward their unique folding and biased codon patterns shared by groups of functionally related mRNAs with coordinate regulation. AGI is important to genetics in general and to human disease; here, we consider influences of its three major components, biased codon use itself, variations in the tRNAome, and anticodon modifications that distinguish synonymous decoding. AGI is plastic and can be used by different species to different extents, with tissue-specificity and in stress responses. Because AGI is species-specific, it is important to consider codon-sensitive experiments when using heterologous systems; for this we focus on the tRNA anticodon loop modification enzyme, CDKAL1, and its link to type 2 diabetes. Newly uncovered tRNAome variability among humans suggests roles in penetrance and as a genetic modifier and disease modifier. Development of experimental and bioinformatics methods are needed to uncover additional means of auxiliary genetic information.

Unraveling patterns of site-to-site synonymous rates variation and associated gene properties of protein domains and families

In protein-coding genes, synonymous mutations are often thought not to affect fitness and therefore are not subject to natural selection. Yet increasingly, cases of non-neutral evolution at certain synonymous sites were reported over the last decade. To evaluate the extent and the nature of site-specific selection on synonymous codons, we computed the site-to-site synonymous rate variation (SRV) and identified gene properties that make SRV more likely in a large database of protein-coding gene families and protein domains. To our knowledge, this is the first study that explores the determinants and patterns of the SRV in real data. We show that the SRV is widespread in the evolution of protein-coding sequences, putting in doubt the validity of the synonymous rate as a standard neutral proxy. While protein domains rarely undergo adaptive evolution, the SRV appears to play important role in optimizing the domain function at the level of DNA. In contrast, protein families are more likely to evolve by positive selection, but are less likely to exhibit SRV. Stronger SRV was detected in genes with stronger codon bias and tRNA reusage, those coding for proteins with larger number of interactions or forming larger number of structures, located in intracellular components and those involved in typically conserved complex processes and functions. Genes with extreme SRV show higher expression levels in nearly all tissues. This indicates that codon bias in a gene, which often correlates with gene expression, may often be a site-specific phenomenon regulating the speed of translation along the sequence, consistent with the co-translational folding hypothesis. Strikingly, genes with SRV were strongly overrepresented for metabolic pathways and those associated with several genetic diseases, particularly cancers and diabetes.

Recombinant protein expression in Escherichia coli: advances and challenges

Escherichia coli is one of the organisms of choice for the production of recombinant proteins. Its use as a cell factory is well-established and it has become the most popular expression platform. For this reason, there are many molecular tools and protocols at hand for the high-level production of heterologous proteins, such as a vast catalog of expression plasmids, a great number of engineered strains and many cultivation strategies. We review the different approaches for the synthesis of recombinant proteins in E. coli and discuss recent progress in this ever-growing field.

Characterisation of a type 1 Avian Paramyxovirus belonging to a divergent group

Newcastle disease, induced by a type 1 Avian Paramyxovirus (APMV-1), is one of the most serious poultry diseases. APMV-1 are divided into two classes based on genetic analysis: class II strains have been recovered from wild or domestic birds and include virulent and avirulent isolates whereas class I strains have been mainly isolated from wild birds and are avirulent. Within class I, a new proposed genotype has recently been reported. The only full genome strain of this group is presently characterised from the point of view of codon usage with reference to class I and class II specificities. Class-specific residues were identified on HN and F proteins that are the two major proteins involved in cell attachment and pathogenicity. Comparison of protein patterns and codon usage for this newly identified APMV-1 strain indicates it is similar to class I viruses but contains a few characteristics close to the viruses of class II. Transmission of viruses from this recently identified divergent group from wild birds to domestic birds could have a major impact on the domestic poultry industry.

Analysis of codon use features of stearoyl-acyl carrier protein desaturase gene in Camellia sinensis

The stearoyl-acyl carrier protein desaturase (SAD) gene widely exists in all kinds of plants. In this paper, the Camellia sinensis SAD gene (CsSAD) sequence was firstly analyzed by Codon W, CHIPS, and CUSP programs online, and then compared with genomes of the tea plant, other species and SAD genes from 11 plant species. The results show that the CsSAD gene and the selected 73 of C. sinensis genes have similar codon usage bias. The CsSAD gene has a bias toward the synonymous codons with A and T at the third codon position, the same as the 73 of C. sinensis genes. Compared with monocotyledons such as Triticum aestivum and Zea mays, the differences in codon usage frequency between the CsSAD gene and dicotyledons such as Arabidopsis thaliana and Nicotiana tobacum are less. Therefore, A. thaliana and N. tobacum expression systems may be more suitable for the expression of the CsSAD gene. The analysis result of SAD genes from 12 plant species also shows that most of the SAD genes are biased toward the synonymous codons with G and C at the third codon position. We believe that the codon usage bias analysis presented in this study will be essential for providing a theoretical basis for discussing the structure and function of the CsSAD gene.

The effects of the synonymous codon usage and tRNA abundance on protein folding of the 3C protease of foot-and-mouth disease virus

The 3C protease of foot-and-mouth disease virus (FMDV) has a conserved amino acid sequence and is responsible for most cleavage in the viral polyprotein. The effects of the synonymous codon usage of FMDV 3C gene and tRNA abundance of the hosts on shaping different folding units (α-helix, β-strand and the coil) in the 3C protease were analyzed based on the structural information of the FMDV 3C protease from Protein Data Bank (PDB: 2BHG) and 210 genes of 3C for all serotypes of FMDV. The strong correlation between some codons usage and the specific folding unit in the FMDV 3C protease is found. As for the effect of translation speed caused by tRNA abundance on the formation of folding units, the codon positions with lowly abundant tRNA scatter in the 3C gene and there is the obvious fluctuation of tRNA abundance locating in the transition boundaries from the β-strand to the α-helix and the coil, respectively. However, codon positions with lowly abundant tRNA clustering into these boundaries are not found, suggesting that the adjustment of translation speed is likely also achieved by the single codon position with low tRNA abundance rather than a cluster. The observations can provide the information for insight into the role of the synonymous codon usage in the formation of 3C protease of FMDV and effect of the tRNA abundance of the hosts on this formation of protease.

Characterization of Synonymous Codon Usage Bias in the Riemerella Anatipestifer OmpA/MotB Gene

The analysis on codon usage bias of OmpA/MotB gene of Riemerella anatipestifer (RA) may provide a basis for understanding the evolution and pathogenesis of RA and for selecting appropriate host expression systems to improve the expression of target genes in vivo and in vitro. In our study, a comparative analysis of the codon usage bias in the newly discovered RA OmpA/MotB gene and the OmpA/MotB gene of 20 reference flavobacteriaceae was performed. The results of the codon adaptation indes (CAI), effective number of codon (ENC), and GC3s values indicated that synonymous codon usage bias in the OmpA/MotB gene of flavobacteriaceae. The results showed that codon usage bias in the RA OmpA/MotB gene was strong bias towards the synonymous codons with A and T at the third codon position. A high level of diversity in codon usage bias existed, and the effective number of codons used in a gene plot revealed that the G+C compositional constraint is the main factor that determines the codon usage bias in OmpA/MotB gene of flavobacteriaceae. Comparison of the codon usage in the OmpA/MotB gene of different organisms revealed that there were 31 codons showing distinct codon usage differences between the RA and E. coli, 41 between the RA and humans, but 29 between the RA and yeast. Therefore the yeast expression system may be more suitable for the expression of RA OmpA/MotB gene. These results may improve our understanding of the evolution, pathogenesis and functional studies of RA, as well as contribute significantly to the area of flavobacteriaceae research.

Are synonymous codons indeed synonymous?

It has long been known that the distribution and frequency of occurence of synonymous codons can vary greatly among different species, and that the abundance of isoaccepting tRNA species could also be very different. The interaction of these two factors may influence the rate and efficiency of protein synthesis and therefore synonymous mutations might influence the fitness of the organism and cannot be treated generally as ‘neutral’ in an evolutionary sense. These general effects of synonymous mutations, and their possible role in evolution, have been discussed in several recent papers. This review, however, will only deal with the influence of synonymous codon replacements on the expression of individual genes. It will describe the possible mechanisms of such effects and will present examples demonstrating the existence and effects of each of these mechanisms.

The complete mitochondrial genome analysis of the tiger (Panthera tigris)

The complete mitochondrial genomes of five tiger samples from three subspecies (P. t. sumatrae, P. t. altica, and P. t. tigris) were successfully obtained by using 26 specifically designed Panthera-specific primer sets. The genome organization and gene arrangement of the five tiger samples were similar to each other; however polymorphic tandem repeat sequences were observed in the control region (CR). This led to a difference in the genome lengths obtained from these five samples with an average size of 16,994 bp for the five tiger mitochondrial genomes. The nucleotide base composition was on average as follows: A, 31.8%; T, 27.0%; C, 26.6%; G, 14.6% and exhibited compositional asymmetry. Most of tiger mitochondrial genome characteristics are similar to those of other common vertebrate species; however, some distinctive features were observed in the CR. First, the repetitive sequence 2 (RS 2) contained two repeat units of 80 bp and the first 15 bp of what would be the third repeat motif. The repetitive sequence 3 (RS 3) contained 47-50 repeat motifs of a shorter 8 bp (ACGTAYAC)(n). Second, length heteroplasmy polycystosine (poly-C) stretches was observed at the end of the HV I locus in all tiger samples.

Mapping codon usage of the translation initiation region in porcine reproductive and respiratory syndrome virus genome

BACKGROUND

Porcine reproductive and respitatory syndrome virus (PRRSV) is a recently emerged pathogen and severely affects swine populations worldwide. The replication of PRRSV is tightly controlled by viral gene expression and the codon usage of translation initiation region within each gene could potentially regulate the translation rate. Therefore, a better understanding of the codon usage pattern of the initiation translation region would shed light on the regulation of PRRSV gene expression.

RESULTS

In this study, the codon usage in the translation initiation region and in the whole coding sequence was compared in PRRSV ORF1a and ORFs2-7. To investigate the potential role of codon usage in affecting the translation initiation rate, we established a codon usage model for PRRSV translation initiation region. We observed that some non-preferential codons are preferentially used in the translation initiation region in particular ORFs. Although some positions vary with codons, they intend to use codons with negative CUB. Furthermore, our model of codon usage showed that the conserved pattern of CUB is not directly consensus with the conserved sequence, but shaped under the translation selection.

CONCLUSIONS

The non-variation pattern with negative CUB in the PRRSV translation initiation region scanned by ribosomes is considered the rate-limiting step in the translation process.