Analysis of codon usage pattern in Taenia saginata based on a transcriptome dataset
© Yang et al.; licensee BioMed Central. 2014
Received: 11 April 2014
Accepted: 6 November 2014
Published: 2 December 2014
Codon usage bias is an important evolutionary feature in a genome and has been widely documented in many genomes. Analysis of codon usage bias has significance for mRNA translation, design of transgenes, new gene discovery, and studies of molecular biology and evolution, etc. However, the information about synonymous codon usage pattern of T. saginata genome remains unclear. T. saginata is a food-borne zoonotic cestode which infects approximataely 50 million humans worldwide, and causes significant health problems to the host and considerable socio-economic losses as a consequence. In this study, synonymous codon usage in T. saginata were examined.
Total RNA was isolated from T. saginata cysticerci and 91,487 unigenes were generated using Illumina sequencing technology. After filtering, the final sequence collection containing 11,399 CDSs was used for our analysis.
Neutrality analysis showed that the T. saginata had a wide GC3 distribution and a significant correlation was observed between GC12 and GC3. NC-plot showed most of genes on or close to the expected curve, but only a few points with low-ENC values were below it, suggesting that mutational bias plays a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that GC and AT were not used proportionally. We also identified twenty-three optimal codons in the T. saginata genome, all of which were ended with a G or C residue. These results suggest that mutational and selection forces are probably driving factors of codon usage bias in T. saginata genome. Meanwhile, other factors such as protein length, gene expression, GC content of genes, the hydropathicity of each protein also influence codon usage.
Here, we systematically analyzed the codon usage pattern and identified factors shaping in codon usage bias in T. saginata. Currently, no complete nuclear genome is available for codon usage analysis at the genome level in T. saginata. This is the first report to investigate codon biology in T. sagninata. Such information does not only bring about a new perspective for understanding the mechanisms of biased usage of synonymous codons but also provide useful clues for molecular genetic engineering and evolutionary studies.
KeywordsTaenia saginata Codon usage bias Trancriptome Optimal codon
Codon usage bias (CUB) refers to the phenomenon where synonymous codons are not used with equal frequencies during translation of genes. CUB is a common phenomenon in a wide variety of organisms, including prokaryotes and eukaryotes –. Many factors have been reported to influence codon usage in various organisms. Weak natural selection and mutational pressure are thought to be the main factors that account for the codon usage variation among the genes in these organisms . Genome-wide investigations of codon usage patterns has an immense importance in understanding the basic features of molecular organization of a genome. In addition, analysis of CUB has many other important applied aspects, such as heterologous gene expression , the determining of the origins of species , the design of degenerate primers , the prediction of expression level of genes ,, as well as the prediction of gene functions . However, most of numerous reports on CUB have focused on model organisms and many microorganisms, such as Caenorhabditis, Drosophila, Arabidopsis, yeast, Giardia lamblia, Entamoeba histolytica, Streptomyces, Borrelia burgdorferi, and Saccharomyces cerevisiae. For example, in C. elegans it is observed that most favored codons are ended with G and/or C (majority are C ending) . In contrast, there are few studies on tapeworms. T. saginata is an important parasitic tapeworm which is widely distributed in the world . The adult worms mainly parasitize in the small intestines of humans ,. T. saginata can cause great economic losses and endangers public health ,. However, the information about synonymous codon usage pattern of T. saginata remains unclear. In this study, we investigated the codon usage profile of T. saginata through transcriptome data using a multivariate statistical analysis. Analysis of codon usage pattern in T. saginata would provide a basis for understanding the related mechanism for biased usage of synonymous codons and for choosing appropriate host expression systems for an optimized expression of target genes.
This study was approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Approval No. LVRIAEC2009-2012). The cattle from which Taenia saginata cysticerci were collected for transcriptome sequencing, were handled in accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People's Republic of China.
RNA extraction, cDNA library preparation and Illumina sequencing
Total RNA was extracted from cysticerci using the Trizol reagent (Invitrogen, Carlsbad,CA), following the manufacturer’s instructions. The quantity and quality of total RNA was analyzed using Agilent 2100 RNA Nanochip (Agilent, Santa Clara, CA, USA) and gel electrophoresis. A total of 16.1 μg of RNA was pooled for the preparation of the cDNA library.
The OligoTex mRNA mini kit (Qiagen) was used to poly-T+ RNA after total RNA was collected according to the manufacturer’s protocol. The mRNA was mixed with fragmentation buffer and fragmented into short fragments. cDNA was synthesized using the mRNA fragments as templates. Short fragment (200 ± 25 bp) were gel extracted from an agarose gel and PCR amplified for 15 cycles. Finally, the library was sequenced using the Illumina HiSeq 2000 sequencer (Beijing Genomics Institute, BGI, Shenzhen, Guangdong, China).
De novo assembly
Using Solexa/Illumina RNA-seq deep sequencing technology, we obtained a total of 55.49 million raw reads (4.99 Gb). Further, raw reads were filtered to remove the low-quality reads. The filtration steps were as follows: 1) remove adaptor sequence; 2) remove reads containing the unknown nucleotide “N” over 10%; 3) remove low quality reads containing more than 10 bases with Q-value ≤ 20. Then, the remaining high-quality reads were used for further analysis. Transcriptome raw reads dataset has been submitted to the NCBI (http://www.ncbi.nlm.nih.gov/bioproject/PRJNA260140).
In this study, a total of 91,487 T. saginata unigenes were obtained. Based on a sequence similarity with known proteins, a total of 59,262 unigenes were annotated. Up to 57,607 of which were annotated against the NCBI non-redundant (Nr) protein database, 24,860 were assigned to the protein database Clusters of Orthologous Groups (COG), 26,476 were assigned to the term annotation database of Gene Ontology (GO), and 43,575 were assigned to 200 pathways in the database of Kyoto Encyclopedia of Genes and Genomes (KEGG). Among the annotated unigenes, 61,941 coding sequences (CDS) were obtained by the BLASTx algorithm . All CDSs were analyzed using the FrameDP software , which has the ability to self-train directly on EST clusters instead of requiring curated cDNA sets to train the underlying ESTScan and DECODER software .
To minimise the sampling error, only CDS sequences longer than 300 bp were used for this study. The final sequence collection containing 11,399 CDSs was used for our analyses.
Indices of codon usage
where s represents the given (G + C)3 % value .
where gij is the observed number of the i th codon for j th amino acid which has n i type of synonymous codons. The codon with RSCU value more than 1.0 has positive codon usage bias, while the value <1.0 has relative negative codon usage bias. When RSCU value is equal to 1.0, it means that this codon is chosen equally and randomly.
The GC content of first, second and third codon position (GC1, GC2 and GC3 respectively) were then calculated. GC12 is the average of GC 1 and GC2, and was used for analysis of neutrality plots (GC12vsGC3) . The codon adaptation index (CAI) was used to estimate the extent of bias toward codons that were known to be preferred in highly expressed genes. A CAI value is between 0 and 1.0, and a higher value means a likely stronger codon usage bias and a potential higher expression level .
Correspondence analysis (CA) has been widely used to explore codon usage variation among genes. CA is a sophisticated multivariate statistical technique in which the codon usage data (59 codons) are plotted in a multidimensional space of 59 axes (excluding Met, Trp and stop codons) and then it identifies the axes which represent the most prominent factors contributing to variation among genes ,.
Determination of optimal codons
We selected 5% of the total genes with extremely high and low CAI values which were regarded as the high and low expression genes datasets, respectively. Codon usage was compared using Chi squared contingency test of the two groups, and codons whose frequency of usage were significantly higher (P < 0.01) in highly expressed genes than in genes with low level of expression would be defined as the optimal codons .
CodonW 1.4.4 software was used to analyze the indices of codon usage. Correlation analysis was carried out using the Spearman’s rank correlation analysis method wrapped in the multianalysis software SPSS version 19.0.
Codon usage in T. saginata
Codon usage in T. saginata and T. pisiformis
In general, the pattern of codon usage is similar among closely related organisms, but differs significantly among distantly related species, such as Escherichia coli, Saccharomyces cerevisiae and Drosophila melanogaster. In this study, patterns of codon usage are compared in T. saginata and T. pisiformis (Table 1) , and we found that there are high similarities between them. With the exception of UCA and GGA, the two species have the same preferred codon for all amino acids.
Nucleotide content of genes
Relation between ENC and GC3
In order to analyze the codon usage of different kinds of gene, we selected the hydrophobic genes with gene scores >5, the aromatic genes with gene scores ≥0.15,ribosomal genes and other genes from 11399 genes. The distribution of the four types of genes were shown in Figure 5C. We employed a multivariate analysis of variance (MANOVA) and found that there was a statistically significant difference among four types of genes in codon usage (p < 0.01).
Gene expression level and synonymous codon usage bias
Protein length and synonymous codon usage bias
Effect of the hydrophobicity and aromaticity of encoded protein on codon bias
Numerous studies have shown that hydrophobicity and aromaticity of encoded protein play important roles in shaping codon usage of some species. In order to investigate if the same thing is happening to T. saginata, we performed a correlation analysis to evaluate whether Gravy and Aromo values were related to ENC values. The correlation analyses between the hydrophobicity of each protein and ENC value showed that the correlation coefficients (r = −0.0883, P <0.001) were significantly correlated. The aromaticity of each protein was not significantly correlated with ENC (r =0.0097, P > 0.05). The analysis results indicated that variation in codon usage were associated with the degree of hydrophobicity, but not with the aromatic amino acids .
Translational optimal codons of T. saginata
Codon usage bias is an important and complex evolutionary phenomenon, and it exists in a wide variety of organisms, from prokaryotes, to unicellular and multicellular eukaryotes. Some hypotheses are proposed to explain the origin of codon usage bias, among which neutral theory  and the selection-mutation-drift balance model , are the most representative ones. According to neutral theory, mutations at degenerate coding positions should be selectively neutral, thus resulting in random synonymous codon choice. In the selection-mutation-drift model, codon bias is thought to be determined by a balance between mutation pressure, genetic drift, and weak selection. In other words, if a gene experiences a highly selective pressure, such as high expression, it may be inclined to stronger codon usage bias. However, in recent years, with the completion of genome projects of many organisms, the two hypotheses are not sufficient to explain codon usage anymore. Many other factors have been reported to influence CUB, including gene length , GC-content ,, recombination rate ,,, gene expression level ,,, RNA structure –, protein structure , intron length , population size , evolutionary age of the genes , environmental stress , the hydrophobicity and the aromaticity of the encoded proteins ,, and so on. In this study, the factors involved in shaping codon usage of the Taenia saginata genome at least includes gene expression level, gene compositional constraint, protein length, as well as the hydrophobicity of each protein (slightly).
Nucleotide composition could be one of the most important factors in shaping codon usage among genes and genomes. GC-rich organisms, such as Bacteria, Archea, Fungi. Triticum Aestivum, Hordium vulgare and Oryza sativa,, tend to use G or C in the third position. And meanwhile, AT-rich organisms show a preference for A or T in third position, such as Onchocerca volvulus, Mycoplasma capricolum and Plasmodium falciparum–. The genomic G + C content for T. saginata is 43.61%. Although the genome would thus appear to be slightly A + T rich, overall codon usage is biased toward C- and G-ending codons (Table 2), this is similar to that in Giardia lamblia.
Previous studies have found significant negative correlations between protein length and CUB in variety of organsims, such as yeast, Caenorhabditis elegans, Drosophila melanogaster], Arabidopsis thaliana and Silene latifolia. Similar results have also been found in T. saginata. There is an explanation proposed by Moriyama and Powell for this phenomenon: namely, if shorter proteins could perform similar functions to those of longer ones, longer proteins become energy-expensive and disadvantageous, thus the selection constraint acts to reduce the size of highly expressed genes, dominantly determines the relationship between codon bias and gene length .
As we know, it is difficult to quantify the expression level of genes in a differentiated multicellular eukaryote, where genes are expressed at different levels in different tissues and at different developmental stages. In the T. saginata genome, the expression level of an individual gene is lacking. It is known that EST counting is efficient for assessing gene expression level. Nevertheless, due to the limitation of EST numbers and rough prediction of gene expression level by counting ESTs, so we use the “Codon Adaptation Index” to evaluate the expression level of examined genes. CAI has been widely used to examine the expressivities of genes by many researchers and has now been considered as a well-accepted measure of gene expression ,.
In this study, we identified 23 codons as the optimal codons. Most of all optimal codons in the T. saginata genome end with G or C. This is very similar to the pattern observed in other eukaryotic genomes, such as Dictyostelium discoideum, D. melanogaster, C. elegans, Giardia lambliaand Schizosaccharomyces pombe. The identification of optimal codons may provide useful clues for molecular genetic engineering and evolutionary studying.
For the first time, we have reported the pattern of codon usage bias in the T. saginata genome and its causative factors. Evidence suggests that the codon usage pattern in T. saginata appears to be the result of a complex equilibrium between different forces, namely mutation bias, natural selection, the GC content of genes, protein length, gene expression level and hydropathicity. Meanwhile, 23 optimal codons were identified, all of which ended with either a G or C residue, this will be useful for cloning and expression of foreign genes in the organism. Such information from this study will provide a better understanding of the characteristics of synonymous codon usage in T. saginata and its molecular evolution, and provide a new resource to underpin the development of urgently needed treatments and control.
We thank several reviewers for helpful comments on the work presented here. Project support was provided by the Science Fund for Creative Research Groups of Gansu Province (Grant No. 1210RJIA006) and opening projects of National Key Laboratory of Veterinary Etiological Biology at Lanzhou Veterinary Research Institute of Chinese Academy of Agricultural Sciences (Grant Number: 201001).
- Akashi H, Eyre-Walker A: Translational selection and molecular evolution. Curr Opin Genet Dev. 1998, 8 (6): 688-693. 10.1016/S0959-437X(98)80038-5.View ArticlePubMedGoogle Scholar
- Akashi H: Gene expression and molecular evolution. Curr Opin Genet Dev. 2001, 11 (6): 660-666. 10.1016/S0959-437X(00)00250-1.View ArticlePubMedGoogle Scholar
- Duret L: Evolution of synonymous codon usage in metazoans. Curr Opin Genet Dev. 2002, 12 (6): 640-649. 10.1016/S0959-437X(02)00353-2.View ArticlePubMedGoogle Scholar
- Hershberg R, Petrov DA: Selection on codon bias. Annu Rev Genet. 2008, 42: 287-299. 10.1146/annurev.genet.42.110807.091442.View ArticlePubMedGoogle Scholar
- Kane JF: Effects of rare codon clusters on high-level expression of heterologous proteins in Escherichia coli. Curr Opin Biotechnol. 1995, 6 (5): 494-500. 10.1016/0958-1669(95)80082-4.View ArticlePubMedGoogle Scholar
- Ahn I, Jeong B-J, Bae S-E, Jung J, Son HS: Genomic analysis of influenza A viruses, including avian flu (H5N1) strains. Eur J Epidemiol. 2006, 21 (7): 511-519. 10.1007/s10654-006-9031-z.View ArticlePubMedGoogle Scholar
- Zheng Y, Zhao WM, Wang H, Zhou YB, Luan Y, Qi M, Cheng YZ, Tang W, Liu J, Yu H, Yu XP, Fan YZ, Yang X: Codon usage bias in Chlamydia trachomatis and the effect of codon modification in the MOMP gene on immune responses to vaccination. Biochem Cell Biol. 2007, 85 (2): 218-226. 10.1139/o06-211.View ArticlePubMedGoogle Scholar
- Naya H, Romero H, Carels N, Zavala A, Musto H: Translational selection shapes codon usage in the GC-rich genome of Chlamydomonas reinhardtii. FEBS Lett. 2001, 501 (2): 127-130. 10.1016/S0014-5793(01)02644-8.View ArticlePubMedGoogle Scholar
- Gupta S, Bhattacharyya T, Ghosh TC: Synonymous codon usage in Lactococcus lactis: mutational bias versus translational selection. J Biomol Struct Dyn. 2004, 21 (4): 527-535. 10.1080/07391102.2004.10506946.View ArticlePubMedGoogle Scholar
- Lin K, Kuang Y, Joseph JS, Kolatkar PR: Conserved codon composition of ribosomal protein coding genes in Escherichia coli, Mycobacterium tuberculosis and Saccharomyces cerevisiae: lessons from supervised machine learning in functional genomics. Nucleic Acids Res. 2002, 30 (11): 2599-2607. 10.1093/nar/30.11.2599.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (8): 4482-4487. 10.1073/pnas.96.8.4482.PubMed CentralView ArticlePubMedGoogle Scholar
- Kliman RM, Irving N, Santiago M: Selection conflicts, gene expression, and codon usage trends in yeast. J Mol Evol. 2003, 57 (1): 98-109. 10.1007/s00239-003-2459-9.View ArticlePubMedGoogle Scholar
- Lafay B, Sharp PM: Synonymous codon usage variation among Giardia lamblia genes and isolates. Mol Biol Evol. 1999, 16 (11): 1484-1495. 10.1093/oxfordjournals.molbev.a026060.View ArticlePubMedGoogle Scholar
- Ghosh TC, Gupta SK, Majumdar S: Studies on codon usage in Entamoeba histolytica. Int J Parasitol. 2000, 30 (6): 715-722. 10.1016/S0020-7519(00)00042-4.View ArticlePubMedGoogle Scholar
- Wright F, Bibb MJ: Codon usage in the G + C-rich Streptomyces genome. Gene. 1992, 113 (1): 55-65. 10.1016/0378-1119(92)90669-G.View ArticlePubMedGoogle Scholar
- McInerney JO: Replicational and transcriptional selection on codon usage in Borrelia burgdorferi. Proc Natl Acad Sci U S A. 1998, 95 (18): 10698-10703. 10.1073/pnas.95.18.10698.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharp PM, Cowe E: Synonymous codon usage in Saccharomyces cerevisiae. Yeast. 1991, 7 (7): 657-678. 10.1002/yea.320070702.View ArticlePubMedGoogle Scholar
- Stenico M, Lloyd AT, Sharp PM: Codon usage in Caenorhabditis elegans: delineation of translational selection and mutational biases. Nucleic Acids Res. 1994, 22 (13): 2437-2446. 10.1093/nar/22.13.2437.PubMed CentralView ArticlePubMedGoogle Scholar
- Wanzala W, Onyango-Abuje JA, Kang'ethe EK, Zessin KH, Kyule NM, Baumann MP, Ochanda H, Harrison LJ: Control of Taenia saginata by post-mortem examination of carcasses. Afr Health Sci. 2003, 3 (2): 68-76.PubMed CentralPubMedGoogle Scholar
- Dorny P, Vercammen F, Brandt J, Vansteenkiste W, Berkvens D, Geerts S: Sero-epidemiological study of Taenia saginata cysticercosis in Belgian cattle. Vet Parasitol. 2000, 88 (1): 43-49. 10.1016/S0304-4017(99)00196-X.View ArticlePubMedGoogle Scholar
- Lightowlers MW, Rolfe R, Gauci CG:Taenia saginata: Vaccination against Cysticercosis in Cattle with Recombinant Oncosphere Antigens. Exp Parasitol. 1996, 84 (3): 330-338. 10.1006/expr.1996.0121.View ArticlePubMedGoogle Scholar
- Matuchansky C, Lenormand Y: Images in clinical medicine. Taenia saginata N Engl J Med. 1999, 341 (23): 1737-10.1056/NEJM199912023412305.View ArticlePubMedGoogle Scholar
- Lees W, Nightingale J, Brown D, Scandrett B, Gajadhar A: Outbreak of Cysticercus bovis (Taenia saginata) in feedlot cattle in Alberta. Can Vet J. 2002, 43 (3): 227-228.PubMed CentralPubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Journet EP, van Tuinen D, Gouzy J, Crespeau H, Carreau V, Farmer MJ, Niebel A, Schiex T, Jaillon O, Chatagnier O, Godiard L, Micheli F, Kahn D, Gianinazzi-Pearson V, Gamas P: Exploring root symbiotic programs in the model legume Medicago truncatula using EST analysis. Nucleic Acids Res. 2002, 30 (24): 5579-5592. 10.1093/nar/gkf685.PubMed CentralView ArticlePubMedGoogle Scholar
- Fukunishi Y, Hayashizaki Y: Amino acid translation program for full-length cDNA sequences with frameshift errors. Physiol Genomics. 2001, 5 (2): 81-87.PubMedGoogle Scholar
- Sharp PM, Li W-H: An evolutionary perspective on synonymous codon usage in unicellular organisms. J Mol Evol. 1986, 24 (1–2): 28-38. 10.1007/BF02099948.View ArticlePubMedGoogle Scholar
- Wright F: The 'effective number of codons' used in a gene. Gene. 1990, 87 (1): 23-29. 10.1016/0378-1119(90)90491-9.View ArticlePubMedGoogle Scholar
- Sueoka N: Directional mutation pressure and neutral molecular evolution. Proc Natl Acad Sci U S A. 1988, 85 (8): 2653-2657. 10.1073/pnas.85.8.2653.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharp PM, Li W-H: The codon adaptation index-a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987, 15 (3): 1281-1295. 10.1093/nar/15.3.1281.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang H-C, Hickey DA: Rapid divergence of codon usage patterns within the rice genome. BMC Evol Biol. 2007, 7 (Suppl 1): S6-10.1186/1471-2148-7-S1-S6.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Q, Feng Y, Zhao X, Dong H, Xue Q: Synonymous codon usage bias in Oryza sativa. Plant Sci. 2004, 167 (1): 101-105. 10.1016/j.plantsci.2004.03.003.View ArticleGoogle Scholar
- Liu Q: Analysis of codon usage pattern in the radioresistant bacterium Deinococcus radiodurans. Biosystems. 2006, 85 (2): 99-106. 10.1016/j.biosystems.2005.12.003.View ArticlePubMedGoogle Scholar
- Sharp PM, Cowe E, Higgins DG, Shields DC, Wolfe KH, Wright F: Codon usage patterns in Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens; a review of the considerable within-species diversity. Nucleic Acids Res. 1988, 16 (17): 8207-8211. 10.1093/nar/16.17.8207.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen L, Liu T, Yang D, Nong X, Xie Y, Fu Y, Wu X, Huang X, Gu X, Wang S, Peng X, Yang G: Analysis of codon usage patterns in Taenia pisiformis through annotated transcriptome data. Biochem Biophys Res Commun. 2013, 430 (4): 1344-1348. 10.1016/j.bbrc.2012.12.078.View ArticlePubMedGoogle Scholar
- Kawabe A, Miyashita NT: Patterns of codon usage bias in three dicot and four monocot plant species. Genes Genet Syst. 2003, 78 (5): 343-352. 10.1266/ggs.78.343.View ArticlePubMedGoogle Scholar
- Sueoka N, Kawanishi Y: DNA G+ C content of the third codon position and codon usage biases of human genes. Gene. 2000, 261 (1): 53-62. 10.1016/S0378-1119(00)00480-7.View ArticlePubMedGoogle Scholar
- Nakamura Y, Gojobori T, Ikemura T: Codon usage tabulated from the international DNA sequence databases. Nucleic Acids Res. 1997, 25 (1): 244-245. 10.1093/nar/25.1.244.PubMed CentralView ArticlePubMedGoogle Scholar
- Bulmer M: Are codon usage patterns in unicellular organisms determined by selection‐mutation balance?. J Evol Biol. 1988, 1 (1): 15-26. 10.1046/j.1420-9101.1988.1010015.x.View ArticleGoogle Scholar
- Comeron JM, Kreitman M, Aguade M: Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila. Genetics. 1999, 151 (1): 239-249.PubMed CentralPubMedGoogle Scholar
- 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 (10): 5688-5692. 10.1073/pnas.091427698.PubMed CentralView ArticlePubMedGoogle Scholar
- Hey J, Kliman RM: Interactions between natural selection, recombination and gene density in the genes of Drosophila. Genetics. 2002, 160 (2): 595-608.PubMed CentralPubMedGoogle Scholar
- Kliman RM, Hey J: Hill-Robertson interference in Drosophila melanogaster: reply to Marais, Mouchiroud and Duret. Genet Res. 2003, 81 (2): 89-90. 10.1017/S0016672302006067.View ArticlePubMedGoogle Scholar
- Hartl DL, Moriyama EN, Sawyer SA: Selection intensity for codon bias. Genetics. 1994, 138 (1): 227-234.PubMed CentralPubMedGoogle Scholar
- Chen Y, Carlini DB, Baines JF, Parsch J, Braverman JM, Tanda S, Stephan W: RNA secondary structure and compensatory evolution. Genes Genet Syst. 1999, 74 (6): 271-286. 10.1266/ggs.74.271.View ArticlePubMedGoogle Scholar
- Carlini DB, Chen Y, Stephan W: The relationship between third-codon position nucleotide content, codon bias, mRNA secondary structure and gene expression in the drosophilid alcohol dehydrogenase genes Adh and Adhr. Genetics. 2001, 159 (2): 623-633.PubMed CentralPubMedGoogle Scholar
- Oresic M, Dehn M, Korenblum D, Shalloway D: Tracing specific synonymous codon-secondary structure correlations through evolution. J Mol Evol. 2003, 56 (4): 473-484. 10.1007/s00239-002-2418-x.View ArticlePubMedGoogle Scholar
- Vinogradov AE: Intron length and codon usage. J Mol Evol. 2001, 52 (1): 2-5. 10.1007/s002390010128.View ArticlePubMedGoogle Scholar
- Berg OG: Selection intensity for codon bias and the effective population size of Escherichia coli. Genetics. 1996, 142 (4): 1379-1382.PubMed CentralPubMedGoogle Scholar
- Prat Y, Fromer M, Linial N, Linial M: Codon usage is associated with the evolutionary age of genes in metazoan genomes. BMC Evol Biol. 2009, 9: 285-10.1186/1471-2148-9-285.PubMed CentralView ArticlePubMedGoogle Scholar
- Goodarzi H, Torabi N, Najafabadi HS, Archetti M: Amino acid and codon usage profiles: adaptive changes in the frequency of amino acids and codons. Gene. 2008, 407 (1–2): 30-41. 10.1016/j.gene.2007.09.020.View ArticlePubMedGoogle Scholar
- Romero H, Zavala A, Musto H: Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces. Nucleic Acids Res. 2000, 28 (10): 2084-2090. 10.1093/nar/28.10.2084.PubMed CentralView ArticlePubMedGoogle Scholar
- Rispe C, Delmotte F, van Ham RC, Moya A: Mutational and selective pressures on codon and amino acid usage in Buchnera, endosymbiotic bacteria of aphids. Genome Res. 2004, 14 (1): 44-53. 10.1101/gr.1358104.PubMed CentralView ArticlePubMedGoogle Scholar
- Hershberg R, Petrov DA: General rules for optimal codon choice. PLoS Genet. 2009, 5 (7): e1000556-10.1371/journal.pgen.1000556.PubMed CentralView ArticlePubMedGoogle Scholar
- Saul A, Battistutta D: Codon usage in Plasmodium falciparum. Mol Biochem Parasitol. 1988, 27 (1): 35-42. 10.1016/0166-6851(88)90022-9.View ArticlePubMedGoogle Scholar
- Milhon JL, Tracy JW: Updated codon usage in Schistosoma. Exp Parasitol. 1995, 80 (2): 353-356. 10.1006/expr.1995.1046.View ArticlePubMedGoogle Scholar
- Muto A, Yamao F, Osawa S: The genome of Mycoplasma capricolum. Prog Nucleic Acid Res Mol Biol. 1987, 34: 29-58. 10.1016/S0079-6603(08)60492-4.View ArticlePubMedGoogle Scholar
- Ingvarsson PK: Gene expression and protein length influence codon usage and rates of sequence evolution in Populus tremula. Mol Biol Evol. 2007, 24 (3): 836-844. 10.1093/molbev/msl212.View ArticlePubMedGoogle Scholar
- Qiu S, Bergero R, Zeng K, Charlesworth D: Patterns of codon usage bias in Silene latifolia. Mol Biol Evol. 2011, 28 (1): 771-780. 10.1093/molbev/msq251.View ArticlePubMedGoogle Scholar
- Moriyama EN, Powell JR: Codon usage bias and tRNA abundance in Drosophila. J Mol Evol. 1997, 45 (5): 514-523. 10.1007/PL00006256.View ArticlePubMedGoogle Scholar
- Sharp PM, Li W-H: On the rate of DNA sequence evolution in Drosophila. J Mol Evol. 1989, 28 (5): 398-402. 10.1007/BF02603075.View ArticlePubMedGoogle Scholar
- Shields DC, Sharp PM, Higgins DG, Wright F: " Silent" sites in Drosophila genes are not neutral: evidence of selection among synonymous codons. Mol Biol Evol. 1988, 5 (6): 704-716.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.