- Research article
- Open Access
Codon usage and codon context bias in Xanthophyllomyces dendrorhous
BMC Genomics volume 16, Article number: 293 (2015)
Synonymous codons are used differentially in organisms from the three domains of life, a phenomenon referred to as codon usage bias. In addition, codon pair bias, particularly in the 3’ codon context, has also been described in several organisms and is associated with the accuracy and rate of translation. An improved understanding of both types of bias is important for the optimization of heterologous protein expression, particularly in biotechnologically important organisms, such as the yeast Xanthophyllomyces dendrorhous, a promising bioresource for the carotenoid astaxanthin. Using genomic and transcriptomic data, the codon usage and codon context biases of X. dendrorhous open reading frames (ORFs) were analyzed to determine their expression levels, GC% and sequence lengths. X. dendrorhous totiviral ORFs were also included in these analyses.
A total of 1,695 X. dendrorhous ORFs were identified through comparison with sequences in multiple databases, and the intron-exon structures of these sequences were determined. Although there were important expression variations among the ORFs under the studied conditions (different phases of growth and available carbon sources), most of these sequences were highly expressed under at least one of the analyzed conditions. Independent of the culture conditions, the highly expressed genes showed a strong bias in both codon usage and the 3’ context, with a minor association with the GC% and no relationship to the sequence length. The codon usage and codon-pair bias of the totiviral ORFs were highly variable with no similarities to the host ORFs.
There is a direct relation between the level of gene expression and codon usage and 3′ context bias in X. dendrorhous, which is more evident for ORFs that are expressed at the highest levels under the studied conditions. However, there is no direct relation between the totiviral ORF biases and the host ORFs.
With the exception of methionine and tryptophan, amino acids are encoded by two to six synonymous codons according to the standard genetic code, and degenerate codons are used at different frequencies, a phenomenon known as codon usage bias (CUB) . Several biological factors, such as the gene GC composition and length, mutation frequency and pattern, gene expression level, tRNA abundance, gene translation initiation signals and protein structure, influence the CUB [2-9].
The existence of CUB has been described in metazoans , D. melanogaster , bacteria [12,13], insects , archaea  and viruses [16-18]. It has been proposed that viral genomes adapt to the host codon usage to efficiently use the host’s translational resources [19,20]. Previous studies have reported interspecies or even intraspecies differences between highly and poorly expressed genes, likely associated with translational efficiency [21,22]. Highly expressed genes typically exhibit higher bias in synonymous codon usage, and it has been proposed that mutation pressure and natural selection are the major forces influencing this phenomenon, favoring translationally superior codons [23-27]. Thus, the most optimal codons are significantly more represented in highly expressed genes than in poorly expressed genes [28,29]. In addition to CUB, codon context bias reflects preferences related to the sequentiality of a pair of codons (codon pair). Codon context bias is likely associated with the accuracy of decoding, indicating the ability of the translational machinery to detect codon pairs present at ribosomal decoding sites [30-33]. One hypothesis is that translation rates are influenced by the compatibilities of adjacent tRNAs at the A- and P-sites on the surface of translating ribosomes. The results of a recent in vivo study suggested that the codon context primarily influences the speed at which proteins are synthesized in E. coli . Preferred and avoided codon pairs have been observed in the three domains of life, and it has been reported that 3′ codons primarily show selective effects on the codon context .
Both CUB and codon context bias analyses have been recommended for the optimization of heterologous gene expression, as parameters that significantly favor gene expression . Thus, knowledge of the CUB and codon context bias is of critical interest for genetic improvement when heterologous expression is used to favor the productivity of biotechnologically important microorganisms. The basidiomycetous yeast Xanthophyllomyces dendrorhous is relevant to biotechnology, as this microorganism synthesizes the carotenoid astaxanthin. This pigment has strong antioxidant properties beneficial for human health, including potential benefits for the treatment of degenerative diseases . In addition, astaxanthin is commonly used in aquaculture for the pigmentation of the flesh of salmonid fishes, which is a considerably important factor in this industry. Although X. dendrorhous is a promising source of natural astaxanthin, natural production in wild-type strains is not sufficient to be economically competitive against the chemical synthesis of this pigment. Therefore, considerable effort has been made to improve the production of carotenoids in X. dendrorhous, including culture optimization, classical random mutagenesis and metabolic engineering approaches (reviewed in ). Unfortunately, the molecular tools to genetically modify this yeast remain scarce , limiting the number of potential modifications that may be of interest. Thus, knowledge of the CUB and codon context bias for this yeast would be a pivotal contribution to the design of new metabolic engineering strategies to improve astaxanthin biosynthesis in this organism. In addition, totiviruses have recently been identified in X. dendrorhous strain UCD 67–385 ; unlike mammalian viruses, these viruses lack an extracellular infection route and are cytoplasmically transmitted.
Although the codon usage of X. dendrorhous has been previously described, the analysis was performed using only ten ribosomal genes . However, the current application of next-generation technologies has provided additional information to conduct more representative studies. In the present study, we evaluated the codon usage and codon context bias of multiple X. dendrorhous genes using genomic and transcriptomic data obtained from the yeast cultured with two different carbon sources (glucose and succinate) during two different phases of growth (exponential and stationary). The level of gene expression was included as a parameter in these analyses for the comparison of codon usage and codon context biases among highly and lowly expressed genes, and the gene expression was also compared against totivirus genes resident in this yeast.
Open reading frame (ORF) identification and expression analysis
The X. dendrorhous strain UCD 67–385 was grown in minimal media supplemented with glucose or succinate as the sole carbon source, and the cells were collected at the early exponential (~18 h) and initial stationary (~72 h) phases of growth, generating a total of four different conditions (G18, S18, G72 and S72: Glucose or Succinate and 18 or 72 h of culture). Total RNA was purified from the yeast pellets, and the quality of samples was assessed and sequenced using the Illumina GAII and HiSeq platforms. Open reading frames (ORFs) of at least 300 bp in length were predicted using transcriptome contigs, and subsequently these sequences were mapped to five genomic scaffolds of 1.1 to 2.4 Mbp in length (approximately 8.1 Mbp in total). Only the mapped ORFs identified under the four conditions were analyzed, and ORFs showing 100% identity with genome sequences, including a well-defined exon-intron structure, were selected and compared with the database using the Blast2GO server . Among the 2,434 sequences analyzed, 1,695 sequences showed positive Blastx hits to at least one conserved protein domain in the InterPro database  (maximum e-value 10−85) (Additional file 1). The remaining 739 sequences with no Blastx hits were not included in the following analyses. In each of the four conditions, the transcriptional levels of each ORF were quantified as reads per kilobases per million mapped reads (RPKM) as previously described . In general, the analyzed ORFs were highly expressed (Figure 1A), and among the four conditions, the percentages of ORFs with RPKM values considered as low- to moderate- (1–30 RPKM), quite high- to high- (31–100 RPKM) and over- (>100 RPKM) expressed, ranged from 2.9 to 10.7, 14.7 to 31 and 58.3 to 82.4%, respectively. The major percentages of over-expressed ORFs were observed after culturing X. dendrorhous in both carbon sources for 72 h, with 82.4% for succinate and 75.8% for glucose. Considering the highest RPKM value for each ORF observed among the four conditions, the percentages of low- to moderate-, quite high- to high- and over-expressed ORFs were 1.6, 10 and 88.4%, respectively. Variations in the expression levels of each ORF were determined by normalizing the RPKM value of each ORF in the reference condition to the lowest RPKM value of the respective ORF among the four conditions. The majority of the ORFs showed considerable variations in expression among the four conditions, although most of the genes were over-expressed (Figure 1B). Smaller differences in the RPKM values were observed after 18 h of culture, and the lowest values were observed using succinate as the sole carbon source. Taking the ratio between the highest and the lowest RPKM value of each ORF among the four conditions as a fold-change in expression, the percentages of ORFs with 1–2, 2.1-5, 5.1-10, 10.1-50 and >50-fold-changes, were 21, 59, 14, 5 and 1%, respectively. The ten ORFs showing the highest expression levels and the highest fold-changes, without considering the ribosomal genes, are listed in Table 1.
Codon usage bias analysis
To analyze the X. dendrorhous CUB, the ORFs were classified according to their expression level under Each Condition (EC grouping), with RPKM values ranging from i) 1–30, ii) 31–70, iii) 71–100, iv) 101–999 and v) ≥1,000. The ribosomal ORFs commonly used as references for highly expressed genes were grouped separately (R grouping). However, as an ORF can be poorly expressed under one condition but highly expressed under another, the ORFs were also classified using the same RPKM value ranges but only considering the Highest RPKM Value observed among the four conditions (HV grouping). In addition, the ORFs were also classified according to the Average RPKM Value among the four conditions (AV grouping). The analysis of relative synonymous codon usage was performed for each group within a classification using the CodonW program/server/software (http://mobyle.pasteur.fr/cgi-bin/portal.py#forms::CodonW, ), and the results are illustrated in Figure 2. Although a direct relation between the expression level and the codon usage was observed in the EC grouping, some variations were observed, depending on the condition (Figure 2A). However, a clearer tendency was observed in the HV and AV grouping, where ORFs with higher expression levels showed a greater preference for some codons including the ribosomal ORFs. Using the codon bias of ribosomal ORFs as a reference, a pattern similar to that of the highly expressed genes in the three different groupings was observed, but this tendency was clearly detected when the HV and AV grouping was compared (Figure 2B). This finding might reflect the differential expression of the analyzed ORFs under different conditions, affecting the number of ORFs in each group (Table 2). The relative synonymous codon usage (RSCU) for the ribosomal ORFs, the highly expressed genes defined in the HV grouping and the ribosomal ORFs are shown in Table 3.
In addition, the ORFs were classified according to sequence length and GC% to analyze the CUB. When the ORFs were grouped according to sequence length, all groups showed similar codon usage, and only the shorter sequences, ranging from 300 to 499 bp, showed some differences with the larger ORFs (Additional file 2 A and B). A direct relation between the CUB and the GC% was observed, with a greater bias in ORFs with a higher GC content (Additional file 2 C). Greater differences in the RSCU ratios between the data for each GC% group and that for the group with 54% GC were observed, whereas the differences in GC% in the ORFs increased (Additional file 2 D).
We also specifically analyzed the codon usage of the X. dendrorhous viral ORFs from totiviruses XdV-L1A and XdV-L1B . As shown in Figure 3A, the codon usage for a majority of the amino acids was quite different among the totiviral ORFs. Compared with the host, only the XdV-L1B totiviral polymerase ORF was similar to the highly expressed X. dendrorhous ORFs, whereas the remaining totiviral ORFs did not show similarities with the lowly or highly expressed ORFs from X. dendrorhous (Figure 3B).
Codon context bias analysis
The 3′ codon context analysis was performed using the Anaconda software  and the HV grouping of ORFs. A 3′ codon context bias was observed in all groups, differing according to the expression level (Figure 4). When the 3′ codon context was compared among the groups, a direct relationship between the expression level and differences in the codon context was observed: ORFs with greater differences in expression level showed greater differences in the 3′ codon context, and ORFs with RPKM values of 101–999 and ≥1,000 showed more similarities (Figure 5). The top five preferred and non-preferred codon pairs in each HV group are listed in Table 4. The non-preferred codon pairs CTT-AAG and CCT-AAG appeared in five and four groups, respectively, whereas the preferred codon pairs, TCA-TCC, AAG-AAG and GAA-GAA, appeared in three groups.
A codon context bias in all groups classified according to the ORFs GC% was observed (Additional file 3 A). The comparison analysis revealed that groups having nearly 50% GC content showed a similar 3′ codon bias. For example, groups having 53 and 54% GC content were more different than groups with 49 and 50% or 50 and 51% GC content (Additional file 3 B). In the case of the ORF length, a 3′ codon context bias was observed in all groups (Additional file 4 A), and these findings were similar among groups of ORFs with 500 or more bases. Shorter ORFs of 300 to 499 bases showed a different 3′ codon context bias (Additional file 4 B).
The analysis of the 3′ codon context of the totiviral ORFs showed differences between the capsid protein ORFs of XdV-L1A and XdV-L1B and the viral polymerase ORFs from both totiviruses (Additional file 5). When the totiviral ORFs were compared with the cellular ORFs, no similarities were observed between poorly or highly expressed ORFs (Additional file 5).
In the present study, based on genomic and transcriptomic data, 1,695 ORFs were selected from the X. dendrorhous strain UCD 67–385, represented in two phases of growth for the yeast cultured using two different carbon sources, glucose or succinate. Furthermore, these ORFs encode a polypeptide with conserved domains listed in the InterPro database and showed a well-defined exon-intron structure. The expression analysis indicated that a majority of the analyzed genes were highly expressed under the four conditions, particularly in older cultures (glucose and succinate, 72 h). Among the four conditions, lower expression levels were observed for the 18-h cultures containing succinate as the sole carbon source, and the expression of a majority of the genes increased from 18 to 72 h of culture in both carbon sources. As several factors influence the CUB, in the present study, we examined the X. dendrorhous CUB according to the gene expression level, sequence length and GC% of the ORF. Clearly defined ORFs were classified according to these parameters, and the 3′ codon context was also analyzed.
In the first analyses, the ORFs were classified into groups according to the expression levels observed under each of the four conditions, independently exhibiting CUB differences among the groups without a clear relationship between the CUB and the expression level. The highest or average expression levels among the four conditions were used to group the ORFs according to their expression level. In these cases, a direct relationship between the CUB and the expression level was observed, and the highly expressed genes showed a major bias, with the exception of the Asp, Cys and His codons. A comparison of the CUB among all groups, based on the expression levels, revealed that the codon usage was similar among genes with similar levels of expression. Although this finding seems rather obvious, the gene expression varied under different conditions; therefore, to classify a gene as lowly or highly expressed based on only one culture condition and state of growth could lead to errors in gene classification and analysis.
Previously, the CUB for the X. dendrorhous strain CBS 6938 was described using ten ribosomal genes . However, when considering a higher number of highly expressed genes from another strain (UCD 67–385), important differences were detected. In the previous study, the usage of the codons GCG (Ala), CGC (Arg), GGG (Gly), AUA (Ile), CUA and CUG (Leu), UUU (Phe), UGA (TER) and UAU (Tyr) in X. dendrorhous was not observed; however, in the present study, we observed that although these codons are not the preferred codons for each amino acid or for a stop codon, these codons are indeed used in X. dendrorhous. In addition, there are more than ten-fold differences in the RSCU values determined for codons GCA (Ala), AGA (Arg), CAA (Gln), GAA (Glu), CAU (His), AGC (Ser) and GUA (Val) between the previous and present studies.
The results obtained from transcriptomic analyses are consistent with the results of previous expression studies. For example, the expression of the genes encoding astaxanthin synthase and phytoene-beta carotene synthase was quantified using RT-qPCR in X. dendrorhous cultured in glucose and succinate as the sole carbon sources , and similar results were observed. The direct relation between the gene expression levels and codon usage biases observed in X. dendrorhous was also consistent with that of other organisms in which highly expressed genes generally show a higher synonymous codon usage bias attributed to selection for efficient translation [24,25,27,47]. Other factors, including gene length  and GC% [48,49], might also influence codon bias. Therefore, we analyzed the codon bias according to these parameters, but no relation was observed for gene length, and although a direct relation regarding the GC% was observed, this association was not as evident as for the gene expression level.
Information regarding the CUB is important in the field of heterologous gene expression to achieve the efficient production of recombinant proteins, for example, enzymes relevant in the biotechnology industry . In recent years, it has been suggested that the codon context or codon pairs might influence translational accuracy and speed, as preferences for specific codon pairs are observed in the three domains of life, referred to as codon context bias [34,35]. Actually, the codon context bias, particularly the 3′ codon context, has been proposed to have as much or even more influence on heterologous gene expression than the CUB [36,51,52]. We observed variations in the 3′ codon context among the groups of genes with different expression levels, and we detected major differences between genes with different expression levels. When we analyzed the genes based on the GC% or sequence length, a codon context bias was observed in which genes with nearly 50% CG content had similar biases. In the case of the gene length, genes of 500 or more bp showed a similar codon context bias. In the three domains of life, the codon pairs with nnUAnn, nnGGnn, nnGnnC, nnCGCn, GUCCnn, CUCCnn, nnCnnA or UUCGnn patterns are most frequently avoided, and codon pairs with nnGCnn, nnCAnn or nnUnCn patterns are most frequently preferred . In the present study, the most avoided codon pairs in X. dendrorhous were consistent with the described patterns, i.e., CCUAAG, GAGUCC and AACCGA, and the most preferred codon pairs were GCGCTC, GUUCCC, ACUCCU, UCUUCU and UCUUCC (the most conserved nucleotides in each pattern are in italics).
In general, viruses do not encode tRNAs, and the synthesis of viral proteins is dependent on the host translational machinery. Thus, several virus sequences have adapted to the host codon usage, including viruses that infect humans and other mammals, particularly for highly expressed genes [19,53]. Two totivirus genomes are present in X. dendrorhous strain UCD 67–385 ; thus, we analyzed the codon usage and the 3′ codon context bias of four totiviral genes with observed variations in both types of bias in all the analyzed genes. Compared with the cellular genes, no similarities with any group classified according to expression level were observed.
In general, the identified X. dendrorhous ORFs are highly expressed, particularly during the stationary phase of growth using succinate or glucose as the sole carbon source, and the majority of the ORFs showed considerable variations in expression under the conditions studied. The codon usage bias and the 3′ codon context bias showed a clear direct relation with the expression levels and GC% of the ORFs, but not the sequence length. However, no similarities among the totiviral and host ORFs were observed for either codon usage or 3′ codon context biases.
X. dendrorhous cultivation conditions and nucleic acid purification
The wild-type X. dendrorhous strain UCD 67–385 (ATCC 24230) was used for next-generation whole genome and transcriptome sequencing and analysis. The strain was cultured at 22°C with constant agitation in YM medium (1% glucose, 0.3% yeast extract, 0.3% malt extract and 0.5% peptone) for DNA extraction or in Vogel minimal medium (MMv) supplemented with 2% glucose or 2% succinate for RNA extraction.
The yeast RNA was purified from the early exponential (18 h) and initial stationary (72 h) phases of growth from cultures grown in MMv medium supplemented with 2% glucose or 2% succinate. After 18 h of culture in MMv medium supplemented with 2% glucose, 1% glucose remained in the medium (confirmed using the DNS method ).
Purification of genomic DNA
X. dendrorhous DNA was isolated from protoplasts as previously described , resulting in a high yield of chromosomal DNA fragments larger than 50 kb. The DNA was purified using phenolic extraction (pH 8.0), including three washes with saturated phenol, three washes with phenol: chloroform: isoamyl alcohol (25: 24: 1) and one wash with chloroform: isoamyl alcohol (24: 1). Subsequently, the DNA was precipitated with 98% ethanol and washed with 70% ethanol. The dried DNA was suspended in Tris: EDTA (10: 1; pH 8.0) with 40 μg/ml of RNase A and incubated for 30 min at 37°C. The DNA was diluted five times with sterile water, and the described phenol extraction protocol was repeated. DNA samples at a 260/280 ratio of 1.7 to 1.9 and a 260/230 ratio >2, measured using a V-630 UV–vis Spectrophotometer (JASCO), were used for next-generation sequencing.
Purification of total RNA
Total RNA was extracted from the cell pellets via mechanical rupture with 0.5 mm glass beads (BioSpec) by vortexing for 10 min, followed by the addition of Tri-Reagent (Ambion). The lysate was incubated for 10 min at room temperature, and subsequently 200 μl of chloroform per ml of Tri-Reagent was added, mixed, and centrifuged for 5 min at 4,000 x g. The aqueous phase was recovered, and two consecutive extractions with acidic phenol: chloroform (1: 1) were performed. The RNA was precipitated with two volumes of isopropanol for 10 min at room temperature, and the RNA was washed with 75% ethanol and suspended in RNase-free water. RNA samples at a 260/280 ratio >1.9, measured using a V-630 UV–vis Spectrophotometer, were used for next-generation sequencing.
Next-generation Sequencing (NGS)
The genome of X. dendrorhous strain UCD 67– 385 was sequenced using the Illumina GAII Sequencing System at Amplicon Express Inc. (http://ampliconexpress.com/, Pullman, Washington, USA) and the Illumina HiSeq2000 System at Macrogen Inc. (http://dna.macrogen.com/eng/index.jsp, Seoul, Republic of Korea). Read assembly and genome and transcriptome analyses were performed using the CLC Genomics Workbench 5. We estimated that the current collection of genomic scaffolds and contigs should cover approximately 95% of the haploid genome of the yeast. For Illumina GAII genome sequencing, a 250-350-bp paired-end library and a 2,500-3,500-bp mate pair library were constructed and sequenced. In addition, 48 primer pairs across 48 gaps were designed for bulk gap closure by sequencing the PCR products with 96 primers. For Illumina HiSeq2000 genome sequencing, a 100-bp paired-end library was constructed and sequenced. The RNA samples from the 72-h culture were sequenced using a Illumina GAII, including a 250-350-bp paired-end library, and the RNA samples from the 18-h culture were sequenced using Illumina HiSeq2000, including a 100-bp paired-end library.
ORFs and gene prediction, annotation and expression level analysis
Using the transcriptome data obtained under each condition, the open reading frames (ORFs) of at least 300 bp in length were predicted using the standard genetic code and the software Geneious® 8.0.2. ORFs that were present in yeast cultured under the four conditions were selected and mapped to five genomic scaffolds of 1,116,253; 1,334,503; 1,461,881; 1,770,274; and 2,396,803 bp. The mapped ORF sequences showing 100% identity with genome sequences, including a correct exon-intron structure, were selected, compared with the database and annotated using the Blast2GO  server: i) the sequences were compared against the National Center for Bioinformatics (NCBI) using the Blastx tool, with an E-value cut off of 10−3; ii) the blast hits of each sequence were mapped using the Gene Ontology Consortium (functional information of known gene products); iii) the GO functional annotation was completed using a cutoff value of 10−6; and iv) functional annotation was performed using InterPro annotations. The expression level of each ORF was calculated under each condition as reads per kilobase per million mapped reads (RPKM), as previously described . The results in which the percentage of coverage of each sequence was at least 90% were used.
Availability of supporting data
The X. dendrorhous ORF names, sequences, RPKM values and the viral ORF sequences used in this work are included in the Additional file 1. The Genbank accession numbers of the XdV-L1A and XdV-L1B viral genomes are [NC_020903 and JN997473], respectively.
Ikemura T. Codon usage and tRNA content in unicellular and multicellular organisms. Mol Biol Evol. 1985;2:13–34.
Bernardi G, Bernardi G. Compositional constraints and genome evolution. J Mol Evol. 1986;24:1–11.
D’Onofrio G, Ghosh TC, Bernardi G. The base composition of the genes is correlated with the secondary structures of the encoded proteins. Gene. 2002;300:179–87.
Hooper SD, Berg OG. Gradients in nucleotide and codon usage along Escherichia coli genes. Nucleic Acids Res. 2000;28:3517–23.
Bulmer M. The selection-mutation-drift theory of synonymous codon usage. Genetics. 1991;129:897–907.
Akashi H. Codon bias evolution Drosophila. Population genetics of mutation-selection drift. Gene. 1997;205:269–78.
Chen SL, Lee W, Hottes AK, Shapiro L, McAdams HH. Codon usage between genomes is constrained by genome-wide mutational processes. Proc Natl Acad Sci U S A. 2004;101:3480–5.
Hershberg R, Petrov DA. Selection on codon bias. Annu Rev Genet. 2008;42:287–99.
Plotkin JB, Kudla G. Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet. 2010;12:32–42.
Duret L. Evolution of synonymous codon usage in metazoans. Curr Opin Genet Dev. 2002;12:640–9.
Akashi H. Synonymous codon usage in Drosophila melanogaster: natural selection and translational accuracy. Genetics. 1994;136:927–35.
Ermolaeva MD. Synonymous codon usage in bacteria. Curr Issues Mol Biol. 2001;3:91–7.
Kanaya S, Kinouchi M, Abe T, Kudo Y, Yamada Y, Nishi T, et al. Analysis of codon usage diversity of bacterial genes with a self-organizing map (SOM): characterization of horizontally transferred genes with emphasis on the E. coli O157 genome. Gene. 2001;276:89–99.
Behura SK, Severson DW. Codon usage bias: causative factors, quantification methods and genome-wide patterns: with emphasis on insect genomes. Biol Rev Camb Philos Soc. 2013;88:49–61.
Khrustalev VV, Barkovsky EV. Study of completed archaeal genomes and proteomes: hypothesis of strong mutational AT pressure existed in their common predecessor. Genom Prot Bioinformatics. 2010;8:22–32.
Cao HW, Zhang H, Cui YD. Synonymous codon usage bias of E2 genes of classical swine fever virus. Israel J Vet Med. 2012;67:253–8.
Hu C, Chen J, Ye L, Chen R, Zhang L, Xue X. Codon usage bias in human cytomegalovirus and its biological implication. Gene. 2014;545:5–14.
Li M, Zhao Z, Chen J, Wang B, Li Z, Li J, et al. Characterization of synonymous codon usage bias in the pseudorabies virus US1 gene. Virol Sin. 2012;27:303–15.
Bahir I, Fromer M, Prat Y, Linial M. Viral adaptation to host: a proteome-based analysis of codon usage and amino acid preferences. Mol Syst Biol. 2009;5:311.
Zhou JH, Zhang J, Chen HT, Ma LN, Ding YZ, Pejsak Z, et al. The codon usage model of the context flanking each cleavage site in the polyprotein of foot-and-mouth disease virus. Infect Genet Evol. 2011;11:1815–9.
dos Reis M, Wernisch L, Savva R. Unexpected correlations between gene expression and codon usage bias from microarray data for the whole Escherichia coli K 12 genome. Nucleic Acids Res. 2003;31:6976–85.
Goetz RM, Fuglsang A. Correlation of codon bias measures with mRNA levels: analysis of transcriptome data from Escherichia coli. Biochem Biophys Res Commun. 2005;327:4–7.
Karlin S, Mrazek J. What drives codon choices in human genes? J Mol Biol. 1996;262:459–72.
Bennetzen JL, Hall BD. Codon selection in yeast. J Biol Chem. 1982;257:3026–31.
Duret L, Mouchiroud D. Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis. Proc Natl Acad Sci U S A. 1999;96:4482–7.
Akashi H. Translational selection and yeast proteome evolution. Genetics. 2003;164:1291–303.
Behura SK, Stanke M, Desjardins CA, Werren JH, Severson DW. Comparative analysis of nuclear tRNA genes of Nasonia vitripennis and other arthropods, and relationships to codon usage bias. Insect Mol Biol. 2010;19 Suppl 1:49–58.
Novoa EM, Ribas de Pouplana L. Speeding with control: codon usage, tRNAs, and ribosomes. Trends Genet. 2012;28:574–81.
Klumpp S, Dong J, Hwa T. On ribosome load, codon bias and protein abundance. PLoS One. 2012;7, e48542.
Irwin B, Heck JD, Hatfield GW. Codon pair utilization biases influence translational elongation step times. J Biol Chem. 1995;270:22801–6.
Parker J. Errors and alternatives in reading the universal genetic code. Microbiol Rev. 1989;53:273–98.
Precup J, Parker J. Missense misreading of asparagine codons as a function of codon identity and context. J Biol Chem. 1987;262:11351–5.
Gutman GA, Hatfield GW. Nonrandom utilization of codon pairs in Escherichia coli. Proc Natl Acad Sci U S A. 1989;86:3699–703.
Chevance FF, Le Guyon S, Hughes KT. The effects of codon context on in vivo translation speed. PLoS Genet. 2014;10, e1004392.
Tats A, Tenson T, Remm M. Preferred and avoided codon pairs in three domains of life. BMC Genomics. 2008;9:463.
Lanza AM, Curran KA, Rey LG, Alper HS. A condition-specific codon optimization approach for improved heterologous gene expression in Saccharomyces cerevisiae. BMC Syst Biol. 2014;8:33.
Higuera-Ciapara I, Felix-Valenzuela L, Goycoolea FM. Astaxanthin: a review of its chemistry and applications. Crit Rev Food Sci Nutr. 2006;46:185–96.
Schmidt I, Schewe H, Gassel S, Jin C, Buckingham J, Humbelin M, et al. Biotechnological production of astaxanthin with Phaffia rhodozyma/Xanthophyllomyces dendrorhous. Appl Microbiol Biotechnol. 2011;89:555–71.
Contreras G, Barahona S, Rojas MC, Baeza M, Cifuentes V, Alcaino J. Increase in the astaxanthin synthase gene (crtS) dose by in vivo DNA fragment assembly in Xanthophyllomyces dendrorhous. BMC Biotechnol. 2013;13:84.
Baeza M, Bravo N, Sanhueza M, Flores O, Villarreal P, Cifuentes V. Molecular characterization of totiviruses in Xanthophyllomyces dendrorhous. Virol J. 2012;9:140.
Verdoes JC, Ooyen AJJV. Codon usage in Xanthophyllomyces dendrorhous (formerly Phaffia rhodozyma). Biotechnol Lett. 2000;22:9–13.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.
Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, et al. InterPro in 2011: new developments in the family and domain prediction database. Nucleic Acids Res. 2012;40:D306–12.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Moura G, Pinheiro M, Silva R, Miranda I, Afreixo V, Dias G, et al. Comparative context analysis of codon pairs on an ORFeome scale. Genome Biol. 2005;6:R28.
Wozniak A, Lozano C, Barahona S, Niklitschek M, Marcoleta A, Alcaino J, et al. Differential carotenoid production and gene expression in Xanthophyllomyces dendrorhous grown in a nonfermentable carbon source. FEMS Yeast Res. 2011;11:252–62.
Sharp PM, Tuohy TMF, Mosurski KR. Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986;14:5125–43.
Comeron JM, Kreitman M, Aguade M. Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila. Genetics. 1999;151:239–49.
Marais G, Mouchiroud D, Duret L. Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes. Proc Natl Acad Sci U S A. 2001;98:5688–92.
Elena C, Ravasi P, Castelli ME, Peiru S, Menzella HG. Expression of codon optimized genes in microbial systems: current industrial applications and perspectives. Front Microbiol. 2014;5:21.
Chung BK, Yusufi FN, Mariati, Yang Y, Lee DY. Enhanced expression of codon optimized interferon gamma in CHO cells. J Biotechnol. 2013;167:326–33.
Chung BK, Lee DY. Computational codon optimization of synthetic gene for protein expression. BMC Syst Biol. 2012;6:134.
Lobo FP, Mota BE, Pena SD, Azevedo V, Macedo AM, Tauch A, et al. Virus-host coevolution: common patterns of nucleotide motif usage in Flaviviridae and their hosts. PLoS One. 2009;4, e6282.
Miller GL. Use of dinitrosalicylic acid reagent for determination of reducing sugar. Anal Chem. 1959;31:426–8.
Cifuentes V, Hermosilla G, Martinez C, Leon R, Pincheira G, Jimenez A. Genetics and electrophoretic karyotyping of wild-type and astaxanthin mutant strains of Phaffia rhodozyma. A Van Leeuw. 1997;72:111–7.
This study was financially supported through FONDECYT grants (nos. 1140504, 1130333 and 11121200).
The authors declare that they have no competing interests.
MB designed and coordinated the study, JA and VC participated in the experimental design and analyses. SB and DS performed the X. dendrorhous cultures and nucleic acid extraction. JA, MB, and VC drafted the manuscript. All authors read and approved the final manuscript.
Sequences and RPKM values of ORFs used in this work.
Codon usage in X. dendrorhous ORFs grouped according to their length and GC%. ORFs were grouped according to sequence length (A and B) and GC% (C and D); parameter ranges are indicated at the top of the figure. A and C: Codon usage for each amino acid in which the most common, intermediary and less preferred codons are represented in red, green and black, respectively, in the continuous color scale. B and D: Graphical representation of the RSCU ratios between data from each parameter ranging from 2,000 – 9,000 bp ORFs (B) or 54% GC (D). Ratios between 0.9 and 1.1 are represented in blue and indicate similarity, whereas ratios out of this range are represented in black.
3′ codon context in X. dendrorhous ORFs classified according to their GC content. ORFs were classified into nine groups (44–46, 47, 48, 49, 50, 51, 52, 53 and 54) according to their GC%. The upper bar indicates the amino acid codons in the 3′ position, and the left bar indicates the reference amino acid codons for each analysis. A: 3′ codon context in ORFs classified according to their GC%. The color coding scale is indicated at the bottom of the figure, where red represents the avoided codons (negative values: bad context), and green indicates the preferred codons (positive values: good context). B: Relationship between the 3′ codon context among X. dendrorhous ORFs according to the GC%. The color coding scale is indicated in the figure, where black denotes similarities, and blue indicates differences.
3′ codon context in X. dendrorhous ORFs classified according to their length. ORFs were classified into five groups (300–499, 500–999, 1,000-1,499, 1,500-1,999 and 2,000-9,000) according to their length in bases. A: 3′ codon context in ORFs classified according to length. The upper bar indicates the amino acid codons in the 3′ position, and the left bar indicates the reference amino acid codons. The color coding scale is indicated at the bottom of the figure, where red represent the avoided codons (negative values: bad context), and green indicates the preferred codons (positive values: good context). B: Relationship between the 3′ codon context among X. dendrorhous ORFs according to length. The color coding scale is indicated in the figure, where black denotes similarities, and blue indicates differences.
Relationships of the 3′ codon context of X. dendrorhous totiviruses. The ORFs L1A-CP and L1B-CP encode the capsid protein, and L1A-Pol and L1B-Pol encode the polymerases from the XdV-L1A and XdV-L1B totiviruses, respectively. The upper bar indicates the amino acid codons in the 3′ position, and the left bar indicates the reference amino acid codons for each analysis. Upper panels: relationship between the 3′ codon context among four the totiviral ORFs. Middle and lower panels: relationship between the 3′ codon context with the totiviral ORFs and the ORFs with the lowest (1–30) and highest (≥1,000) expression levels in X. dendrorhous. The color coding scale is indicated in the figure, where black denotes similarity, and blue indicates differences.
About this article
Cite this article
Baeza, M., Alcaíno, J., Barahona, S. et al. Codon usage and codon context bias in Xanthophyllomyces dendrorhous . BMC Genomics 16, 293 (2015) doi:10.1186/s12864-015-1493-5
- Codon usage bias
- X. dendrorhous
- Codon context bias