Non-random arrangement of synonymous codons in archaea coding sequences

• The order of synonymous codons is not randomly arranged in most archaea.
• Identical codon pairing (ICP) is the major pattern for synonymous codon ordering.
• Selection has driven synonymous codons to be arranged in a beneficial order.
• Translational selection influences both frequencies and orders of synonymous codons.

Non-random arrangement of synonymous codons in coding sequences has been recently reported in eukaryotic and bacterial genomes, but the case in archaeal genomes is largely undetermined. Here, we systematically investigated 122 archaeal genomes for their synonymous codon co-occurrence patterns. We found that in most archaeal coding sequences, the order of synonymous codons is not arranged randomly, but rather some successive codon pairs appear significantly more often than expected. Importantly, such codon pairing bias (CPB) pattern in archaea does not seem to completely follow the co-tRNA codon pairing (CCP) rule previously reported for eukaryotes, but largely obeys an identical codon pairing (ICP) rule. Further, synonymous codon permutation test demonstrated that in many archaeal genomes, random mutation alone is unable to cause the observed high level of ICP bias, which strongly indicates that selection force has been involved to shape synonymous codon orders, potentially meeting a global requirement to optimize translation rate.