The Fitness Effects of Codon Composition of the Horizontally Transferred Antibiotic Resistance Genes Intensify at Sub-lethal Antibiotic Levels.

The rampant variability in codon bias existing between bacterial genomes is expected to interfere with horizontal gene transfer (HGT), a phenomenon that drives bacterial adaptation. However, delineating the constraints imposed by codon bias on functional integration of the transferred genes is complicated by multiple genomic and functional barriers controlling HGT, and by the dependence of the evolutionary outcomes of HGT on the host's environment. Here, we designed an experimental system in which codon composition of the transferred genes is the only variable triggering fitness change of the host. We replaced Escherichia coli's chromosomal folA gene encoding dihydrofolate reductase, an essential enzyme that constitutes a target for trimethoprim, with combinatorial libraries of synonymous codons of folA genes from trimethoprim-sensitive Listeria grayi and trimethoprim-resistant Neisseria sicca. The resulting populations underwent selection at a range of trimethoprim concentrations, and the ensuing changes in variant frequencies were used to infer the fitness effects of the individual combinations of codons. We found that when HGT causes overstabilization of the 5'-end mRNA, the fitness contribution of mRNA folding stability dominates over that of codon optimality. The 5'-end overstabilization can also lead to mRNA accumulation outside of the polysome, thus preventing the decay of the foreign transcripts despite the codon composition-driven reduction in translation efficiency. Importantly, the fitness effects of mRNA stability or codon optimality become apparent only at sub-lethal levels of trimethoprim individually tailored for each library, emphasizing the central role of the host's environment in shaping the codon bias compatibility of horizontally transferred genes.

Genomic imprints of unparalleled growth.

Chlorella ohadii was isolated from desert biological soil crusts, one of the harshest habitats on Earth, and is emerging as an exciting new green model for studying growth, photosynthesis and metabolism under a wide range of conditions. Here, we compared the genome of C. ohadii, the fastest growing alga on record, to that of other green algae, to reveal the genomic imprints empowering its unparalleled growth rate and resistance to various stressors, including extreme illumination. This included the genome of its close relative, but slower growing and photodamage sensitive, C. sorokiniana UTEX 1663. A larger number of ribosome-encoding genes, high intron abundance, increased codon bias and unique genes potentially involved in metabolic flexibility and resistance to photodamage are all consistent with the faster growth of C. ohadii. Some of these characteristics highlight general trends in Chlorophyta and Chlorella spp. evolution, and others open new broad avenues for mechanistic exploration of their relationship with growth. This work entails a unique case study for the genomic adaptations and costs of exceptionally fast growth and sheds light on the genomic signatures of fast growth in photosynthetic cells. It also provides an important resource for future studies leveraging the unique properties of C. ohadii for photosynthesis and stress response research alongside their utilization for synthetic biology and biotechnology aims.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages.

BACKGROUND: In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages. RESULTS: Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations. CONCLUSIONS: This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

The Codon Statistics Database: A Database of Codon Usage Bias.

We present the Codon Statistics Database, an online database that contains codon usage statistics for all the species with reference or representative genomes in RefSeq (over 15,000). The user can search for any species and access two sets of tables. One set lists, for each codon, the frequency, the Relative Synonymous Codon Usage, and whether the codon is preferred. Another set of tables lists, for each gene, its GC content, Effective Number of Codons, Codon Adaptation Index, and frequency of optimal codons. Equivalent tables can be accessed for (1) all nuclear genes, (2) nuclear genes encoding ribosomal proteins, (3) mitochondrial genes, and (4) chloroplast genes (if available in the relevant assembly). The user can also search for any taxonomic group (e.g., "primates") and obtain a table comparing all the species in the group. The database is free to access without registration at http://codonstatsdb.unr.edu.

Evolution of the Tóxicos en Levadura 63 (TL63) gene family in plants and functional characterization of Arabidopsis thaliana TL63 under oxidative stress.

MAIN CONCLUSION: TL63 orthologs were angiosperm specific and had undergone motifs loss and gain, and increased purifying selection. AtTL63 was involved in the response of yeast and Arabidopsis plants to oxidative stress. The Tóxicos en Levadura (TL) family, a class of E3 ubiquitin ligases with typical RING-H2 type zinc finger structure, plays a pivotal role in mediating physiological processes and responding to stress in plants. However, the evolution and function of TL63 remain unclear. In this study, TL63 homologs were dated roughly back to the origin of land plants and confirmed to have subjected to the gain and loss of motifs and increased purifying selection. Phylogenetic analysis displayed that 279 TL63s could be divided into four main clades (Clade A-D). Notably, the ancestral tandem TL40/41 cluster contributed to the expansion of modern Brassicaceae TL40/41. The substitution rate tests revealed that the TL63 lineage was evidently different from other lineages. The codon usage index exhibited that monocotyledons preferred to use not A3s and T3s, but C3s, G3s, CAI, CBI and Fop. Sequence analysis showed that the TL63 homologs had conserved TM and GLD motifs and RING-H2 domain whose key amino acid residues accounted for the high average abundance. Particularly, Arabidopsis thaliana TL63 (AtTL63) was located in the nuclei, cell membranes and peroxisomes and expressed universally and significantly throughout A. thaliana development. Under H2O2 treatment, low or moderate expression of the AtTL63 held beneficial effects on the growth and viability of yeast cells and the mutation or overexpression of the AtTL63 positively affected the growth of A. thaliana plants. In brief, this study could supply useful insight into the evolution of the plant TL63s and the AtTL63 functions under oxidative stress.

[Complete chloroplast genome sequencing and phylogeny of wild Atractylodes lancea from Yuexi, Anhui province].

This study investigated the choroplast genome sequence of wild Atractylodes lancea from Yuexi in Anhui province by high-throughput sequencing, followed by characterization of the genome structure, which laid a foundation for the species identification, analysis of genetic diversity, and resource conservation of A. lancea. To be specific, the total genomic DNA was extracted from the leaves of A. lancea with the improved CTAB method. The chloroplast genome of A. lancea was sequenced by the high-throughput sequencing technology, followed by assembling by metaSPAdes and annotation by CPGAVAS2. Bioiformatics methods were employed for the analysis of simple sequence repeats(SSRs), inverted repeat(IR) border, codon bias, and phylogeny. The results showed that the whole chloroplast genome of A. lancea was 153 178 bp, with an 84 226 bp large single copy(LSC) and a 18 658 bp small single copy(SSC) separated by a pair of IRs(25 147 bp). The genome had the GC content of 37.7% and 124 genes: 87 protein-coding genes, 8 rRNA genes, and 29 tRNA genes. It had 26 287 codons and encoded 20 amino acids. Phylogenetic analysis showed that Atractylodes species clustered into one clade and that A. lancea had close genetic relationship with A. koreana. This study established a method for sequencing the chloroplast genome of A. lancea and enriched the genetic resources of Compositae. The findings are expected to lay a foundation for species identification, analysis of genetic diversity, and resource conservation of A. lancea.

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.

On the Base Composition of Transposable Elements.

Transposable elements exhibit a base composition that is often different from the genomic average and from hosts' genes. The most common compositional bias is towards Adenosine and Thymine, although this bias is not universal, and elements with drastically different base composition can coexist within the same genome. The AT-richness of transposable elements is apparently maladaptive because it results in poor transcription and sub-optimal translation of proteins encoded by the elements. The cause(s) of this unusual base composition remain unclear and have yet to be investigated. Here, I review what is known about the nucleotide content of transposable elements and how this content can affect the genome of their host as well as their own replication. The compositional bias of transposable elements could result from several non-exclusive processes including horizontal transfer, mutational bias, and selection. It appears that mutation alone cannot explain the high AT-content of transposons and that selection plays a major role in the evolution of the compositional bias. The reason why selection would favor a maladaptive nucleotide content remains however unexplained and is an area of investigation that clearly deserves attention.

Comprehensive Identification and Analyses of the GRF Gene Family in the Whole-Genome of Four Juglandaceae Species.

The GRF gene family plays an important role in plant growth and development as regulators involved in plant hormone signaling and metabolism. However, the Juglandaceae GRF gene family remains to be studied. Here, we identified 15, 15, 19, and 20 GRF genes in J. regia, C. illinoinensis, J. sigillata, and J. mandshurica, respectively. The phylogeny shows that the Juglandaceae family GRF is divided into two subfamilies, the ε-group and the non-ε-group, and that selection pressure analysis did not detect amino acid loci subject to positive selection pressure. In addition, we found that the duplications of the Juglandaceae family GRF genes were all segmental duplication events, and a total of 79 orthologous gene pairs and one paralogous homologous gene pair were identified in four Juglandaceae families. The Ka/KS ratios between these homologous gene pairs were further analyzed, and the Ka/KS values were all less than 1, indicating that purifying selection plays an important role in the evolution of the Juglandaceae family GRF genes. The codon bias of genes in the GRF family of Juglandaceae species is weak, and is affected by both natural selection pressure and base mutation, and translation selection plays a dominant role in the mutation pressure in codon usage. Finally, expression analysis showed that GRF genes play important roles in pecan embryo development and walnut male and female flower bud development, but with different expression patterns. In conclusion, this study will serve as a rich genetic resource for exploring the molecular mechanisms of flower bud differentiation and embryo development in Juglandaceae. In addition, this is the first study to report the GRF gene family in the Juglandaceae family; therefore, our study will provide guidance for future comparative and functional genomic studies of the GRF gene family in the Juglandaceae specie.

Effect of genome composition and codon bias on infectious bronchitis virus evolution and adaptation to target tissues.

BACKGROUND: Infectious bronchitis virus (IBV) is one of the most relevant viruses affecting the poultry industry, and several studies have investigated the factors involved in its biological cycle and evolution. However, very few of those studies focused on the effect of genome composition and the codon bias of different IBV proteins, despite the remarkable increase in available complete genomes. In the present study, all IBV complete genomes were downloaded (n = 383), and several statistics representative of genome composition and codon bias were calculated for each protein-coding sequence, including but not limited to, the nucleotide odds ratio, relative synonymous codon usage and effective number of codons. Additionally, viral codon usage was compared to host codon usage based on a collection of highly expressed genes in IBV target and nontarget tissues. RESULTS: The results obtained demonstrated a significant difference among structural, non-structural and accessory proteins, especially regarding dinucleotide composition, which appears under strong selective forces. In particular, some dinucleotide pairs, such as CpG, a probable target of the host innate immune response, are underrepresented in genes coding for pp1a, pp1ab, S and N. Although genome composition and dinucleotide bias appear to affect codon usage, additional selective forces may act directly on codon bias. Variability in relative synonymous codon usage and effective number of codons was found for different proteins, with structural proteins and polyproteins being more adapted to the codon bias of host target tissues. In contrast, accessory proteins had a more biased codon usage (i.e., lower number of preferred codons), which might contribute to the regulation of their expression level and timing throughout the cell cycle. CONCLUSIONS: The present study confirms the existence of selective forces acting directly on the genome and not only indirectly through phenotype selection. This evidence might help understanding IBV biology and in developing attenuated strains without affecting the protein phenotype and therefore immunogenicity.

Codon usage bias and environmental adaptation in microbial organisms.

In each genome, synonymous codons are used with different frequencies; this general phenomenon is known as codon usage bias. It has been previously recognised that codon usage bias could affect the cellular fitness and might be associated with the ecology of microbial organisms. In this exploratory study, we investigated the relationship between codon usage bias, lifestyles (thermophiles vs. mesophiles; pathogenic vs. non-pathogenic; halophilic vs. non-halophilic; aerobic vs. anaerobic and facultative) and habitats (aquatic, terrestrial, host-associated, specialised, multiple) of 615 microbial organisms (544 bacteria and 71 archaea). Principal component analysis revealed that species with given phenotypic traits and living in similar environmental conditions have similar codon preferences, as represented by the relative synonymous codon usage (RSCU) index, and similar spectra of tRNA availability, as gauged by the tRNA gene copy number (tGCN). Moreover, by measuring the average tRNA adaptation index (tAI) for each genome, an index that can be associated with translational efficiency, we observed that organisms able to live in multiple habitats, including facultative organisms, mesophiles and pathogenic bacteria, are characterised by a reduced translational efficiency, consistently with their need to adapt to different environments. Our results show that synonymous codon choices might be under strong translational selection, which modulates the choice of the codons to differently match tRNA availability, depending on the organism's lifestyle needs. To our knowledge, this is the first large-scale study that examines the role of codon bias and translational efficiency in the adaptation of microbial organisms to the environment in which they live.

Evolutionary pressures and codon bias in low complexity regions of plasmodia.

The biological meaning of low complexity regions in the proteins of Plasmodium species is a topic of discussion in evolutionary biology. There is a debate between selectionists and neutralists, who either attribute or do not attribute an effect of low-complexity regions on the fitness of these parasites, respectively. In this work, we comparatively study 22 Plasmodium species to understand whether their low complexity regions undergo a neutral or, rather, a selective and species-dependent evolution. The focus is on the connection between the codon repertoire of the genetic coding sequences and the occurrence of low complexity regions in the corresponding proteins. The first part of the work concerns the correlation between the length of plasmodial proteins and their propensity at embedding low complexity regions. Relative synonymous codon usage, entropy, and other indicators reveal that the incidence of low complexity regions and their codon bias is species-specific and subject to selective evolutionary pressure. We also observed that protein length, a relaxed selective pressure, and a broad repertoire of codons in proteins, are strongly correlated with the occurrence of low complexity regions. Overall, it seems plausible that the codon bias of low-complexity regions contributes to functional innovation and codon bias enhancement of proteins on which Plasmodium species rest as successful evolutionary parasites.

Codon bias confers stability to human mRNAs.

Codon bias has been implicated as one of the major factors contributing to mRNA stability in several model organisms. However, the molecular mechanisms of codon bias on mRNA stability remain unclear in humans. Here, we show that human cells possess a mechanism to modulate RNA stability through a unique codon bias. Bioinformatics analysis showed that codons could be clustered into two distinct groups-codons with G or C at the third base position (GC3) and codons with either A or T at the third base position (AT3): the former stabilizing while the latter destabilizing mRNA. Quantification of codon bias showed that increased GC3-content entails proportionately higher GC-content. Through bioinformatics, ribosome profiling, and in vitro analysis, we show that decoupling the effects of codon bias reveals two modes of mRNA regulation, one GC3- and one GC-content dependent. Employing an immunoprecipitation-based strategy, we identify ILF2 and ILF3 as RNA-binding proteins that differentially regulate global mRNA abundances based on codon bias. Our results demonstrate that codon bias is a two-pronged system that governs mRNA abundance.

Codon usage bias regulates gene expression and protein conformation in yeast expression system P. pastoris.

BACKGROUND: Protein synthesis is one of the extremely important anabolic pathways in the yeast expression system Pichia pastoris. Codon optimization is a commonly adopted strategy for improved protein expression, although unexpected failures did appear sometimes waiting for further exploration. Recently codon bias has been studied to regulate protein folding and activity in many other organisms. RESULTS: Here the codon bias profile of P. pastoris genome was examined first and a direct correlation between codon translation efficiency and usage frequency was identified. By manipulating the codon choices of both endogenous and heterologous signal peptides, secretion abilities of N-terminal signal peptides were shown to be tolerant towards codon changes. Then two gene candidates with different levels of structural disorder were studied, and full-length codon optimization was found to affect their expression profiles differentially. Finally, more evidences were provided to support possible protein conformation change brought by codon optimization in structurally disordered proteins. CONCLUSION: Our results suggest that codon bias regulates gene expression by modulating several factors including transcription and translation efficiency, protein folding and activity. Because of sequences difference, the extent of affection may be gene specific. For some genes, special codon optimization strategy should be adopted to ensure appropriate expression and conformation.

Circadian Oscillators: Around the Transcription-Translation Feedback Loop and on to Output.

From cyanobacteria to mammals, organisms have evolved timing mechanisms to adapt to environmental changes in order to optimize survival and improve fitness. To anticipate these regular daily cycles, many organisms manifest ∼24h cell-autonomous oscillations that are sustained by transcription-translation-based or post-transcriptional negative-feedback loops that control a wide range of biological processes. With an eye to identifying emerging common themes among cyanobacterial, fungal, and animal clocks, some major recent developments in the understanding of the mechanisms that regulate these oscillators and their output are discussed. These include roles for antisense transcription, intrinsically disordered proteins, codon bias in clock genes, and a more focused discussion of post-transcriptional and translational regulation as a part of both the oscillator and output.

Selection of functional 2A sequences within foot-and-mouth disease virus; requirements for the NPGP motif with a distinct codon bias.

Foot-and-mouth disease virus (FMDV) has a positive-sense ssRNA genome including a single, large, open reading frame. Splitting of the encoded polyprotein at the 2A/2B junction is mediated by the 2A peptide (18 residues long), which induces a nonproteolytic, cotranslational "cleavage" at its own C terminus. A conserved feature among variants of 2A is the C-terminal motif N16P17G18/P19, where P19 is the first residue of 2B. It has been shown previously that certain amino acid substitutions can be tolerated at residues E14, S15, and N16 within the 2A sequence of infectious FMDVs, but no variants at residues P17, G18, or P19 have been identified. In this study, using highly degenerate primers, we analyzed if any other residues can be present at each position of the NPG/P motif within infectious FMDV. No alternative forms of this motif were found to be encoded by rescued FMDVs after two, three, or four passages. However, surprisingly, a clear codon preference for the wt nucleotide sequence encoding the NPGP motif within these viruses was observed. Indeed, the codons selected to code for P17 and P19 within this motif were distinct; thus the synonymous codons are not equivalent.

From SARS and MERS CoVs to SARS-CoV-2: Moving toward more biased codon usage in viral structural and nonstructural genes.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an emerging disease with fatal outcomes. In this study, a fundamental knowledge gap question is to be resolved by evaluating the differences in biological and pathogenic aspects of SARS-CoV-2 and the changes in SARS-CoV-2 in comparison with the two prior major COV epidemics, SARS and Middle East respiratory syndrome (MERS) coronaviruses. METHODS: The genome composition, nucleotide analysis, codon usage indices, relative synonymous codons usage, and effective number of codons (ENc) were analyzed in the four structural genes; Spike (S), Envelope (E), membrane (M), and Nucleocapsid (N) genes, and two of the most important nonstructural genes comprising RNA-dependent RNA polymerase and main protease (Mpro) of SARS-CoV-2, Beta-CoV from pangolins, bat SARS, MERS, and SARS CoVs. RESULTS: SARS-CoV-2 prefers pyrimidine rich codons to purines. Most high-frequency codons were ending with A or T, while the low frequency and rare codons were ending with G or C. SARS-CoV-2 structural proteins showed 5 to 20 lower ENc values, compared with SARS, bat SARS, and MERS CoVs. This implies higher codon bias and higher gene expression efficiency of SARS-CoV-2 structural proteins. SARS-CoV-2 encoded the highest number of over-biased and negatively biased codons. Pangolin Beta-CoV showed little differences with SARS-CoV-2 ENc values, compared with SARS, bat SARS, and MERS CoV. CONCLUSION: Extreme bias and lower ENc values of SARS-CoV-2, especially in Spike, Envelope, and Mpro genes, are suggestive for higher gene expression efficiency, compared with SARS, bat SARS, and MERS CoVs.

Conservation of the voltage-sensitive sodium channel protein within the Insecta.

The voltage-sensitive sodium channel (VSSC) is essential for the generation and propagation of action potentials. VSSC kinetics can be modified by producing different splice variants. The functionality of VSSC depends on features such as the voltage sensors, the selectivity filter and the inactivation loop. Mutations in Vssc conferring resistance to pyrethroid insecticides are known as knockdown resistance (kdr). We analysed the conservation of VSSC in both a broad scope and a narrow scope by three approaches: (1) we compared conservation of sequences and of differential exon use across orders of the Insecta; (2) we determined which kdr mutations were possible with a single nucleotide mutation in nine populations of Aedes aegypti; and (3) we examined the individual VSSC variation that exists within a population of Drosophila melanogaster. There is an increasing amount of transcript diversity possible from Diplura towards Diptera. The residues of the voltage sensors, selectivity filter and inactivation loop are highly conserved. The majority of exon sequences were >88.6% similar. Strain-specific differences in codon constraints exist for kdr mutations in nine strains of A. aegypti. Three Vssc mutations were found in one population of D. melanogaster. This study shows that, overall, Vssc is highly conserved across Insecta and within a population of an insect, but that important differences do exist.

Identification & characterization of leucine-rich repeat kinase 2 & parkin RBR E3 ubiquitin protein ligase variants in patients with Parkinson's disease.

BACKGROUND & OBJECTIVES: Parkinson's disease (PD) is a motor disorder that affects movement. More than 24 loci and 28 associated genes have been identified to be associated with this disease. The present study accounts for the contribution of two candidates, leucine-rich repeat kinase 2 ( LRRK2) and parkin RBR E3 ubiquitin protein ligase ( PRKN) in the PD patients, and their characterization in silico and in vitro. METHODS: A total of 145 sporadic PD cases and 120 ethnically matched healthy controls were enrolled with their informed consent. Mutation screening was performed by direct DNA sequencing of the targeted exons of LRRK2 and all exons flanking introns of PRKN. The effect of the pathogenic PRKN variants on a drug (MG-132) induced loss of mitochondrial membrane potential (△ΨM) was measured by a fluorescent dye tetramethylrhodamine methyl ester (TMRM). RESULTS: Twelve and 20 genetic variants were identified in LRRK2 and PRKN, respectively. Interestingly, five out of seven exonic LRRK2 variants were synonymous. Further assessment in controls confirmed the rarity of two such p.Y1527 and p.V1615. Among the pathogenic missense variations (as predicted in silico) in PRKN, two were selected (p.R42H and p.A82E) for their functional study in vitro, which revealed the reduced fluorescence intensity of TMRM as compared to wild type, in case of p.R42H but not the other. INTERPRETATION & CONCLUSIONS: About 6.2 per cent of the cases (9/145) in the studied patient cohort were found to carry pathogenic (as predicted in silico) missense variations in PRKN in heterozygous condition but not in case of LRRK2 which was rare. The presence of two rare synonymous variants of LRRK2 (p.Y1527 and p.V1615) may support the phenomenon of codon bias. Functional characterization of selected PRKN variations revealed p.R42H to cause disruption of mitochondrial membrane potential (△ΨM) rendering cells more susceptible to cellular stress.

An analysis of codon bias in six red yeast species.

Red yeasts, primarily species of Rhodotorula, Sporobolomyces, and other genera of Pucciniomycotina, are traditionally considered proficient systems for lipid and terpene production, and only recently have also gained consideration for the production of a wider range of molecules of biotechnological potential. Improvements of transgene delivery protocols and regulated gene expression systems have been proposed, but a dearth of information on compositional and/or structural features of genes has prevented transgene sequence optimization efforts for high expression levels. Here, the codon compositional features of genes in six red yeast species were characterized, and the impact that evolutionary forces may have played in shaping this compositional bias was dissected by using several computational approaches. Results obtained are compatible with the hypothesis that mutational bias, although playing a significant role, cannot alone explain synonymous codon usage bias of genes. Nevertheless, several lines of evidences indicated a role for translational selection in driving the synonymous codons that allow high expression efficiency. These optimal synonymous codons are identified for each of the six species analyzed. Moreover, the presence of intragenic patterns of codon usage, which are thought to facilitate polyribosome formation, was highlighted. The information presented should be taken into consideration for transgene design for optimal expression in red yeast species.