Comparative analysis of codon usage bias in Crenarchaea and Euryarchaea genome reveals differential preference of synonymous codons to encode highly expressed ribosomal and RNA polymerase proteins

The present study was undertaken to investigate the pattern of optimal codon usage in Archaea. Comparative analysis was executed to understand the pattern of codon usage bias between the high expression genes (HEG) and the whole genomes in two Archaeal phyla, Crenarchaea and Euryarchaea. The G + C% of the HEG was found to be less in comparison to the genome G + C% in Crenarchaea, whereas reverse was the case in Euryarchaea. The preponderance of U/A ending codons that code for HEG in Crenarchaea was in sharp contrast to the C/G ended ones in Euryarchaea. The analysis revealed prevalence of U-ending codons even within the WWY (nucleotide ambiguity code) families in Crenarchaea vis-à-vis Euryarchaea, bacteria and Eukarya. No plausible interpretation of the observed disparity could be made either in the context of tRNA gene composition or genome G + C%. The results in this study attested that the preferential biasness for codons in HEG of Crenarchaea might be different from Euryarchaea. The main highlights are (i) varied CUB in the HEG and in the whole genomes in Euryarchaea and Crenarchaea. (ii) Crenarchaea was found to have some unusual optimal codons (OCs) compared to other organisms. (iii) G + C% (and GC₃) of the HEG were different from the genome G + C% in the two phyla. (iv) Genome G + C% and tRNA gene number failed to explain CUB in Crenarchaea. (v) Translational selection is possibly responsible for A + T rich OCs in Crenarchaea.