Comparative analysis of nuclear tRNA genes of Nasonia vitripennis and other arthropods, and relationships to codon usage bias

Using bioinformatics methods, we identified a total of 221 and 199 tRNA genes in the nuclear genomes of Nasonia vitripennis and honey bee (Apis mellifera), respectively. We performed comparative analyses of Nasonia tRNA genes with honey bee and other selected insects to understand genomic distribution, sequence evolution and relationship of tRNA copy number with codon usage patterns. Many tRNA genes are located physically close to each other in the form of small clusters in the Nasonia genome. However, the number of clusters and the tRNA genes that form such clusters vary from species to species. In particular, the Ala-, Pro-, Tyr- and His-tRNA genes tend to accumulate in clusters in Nasonia but not in honey bee, whereas the bee contains a long cluster of 15 tRNA genes (of which 13 are Gln-tRNAs) that is absent in Nasonia. Though tRNA genes are highly conserved, contrasting patterns of nucleotide diversity are observed among the arm and loop regions of tRNAs between Nasonia and honey bee. Also, the sequence convergence between the reconstructed ancestral tRNAs and the present day tRNAs suggests a common ancestral origin of Nasonia and honey bee tRNAs. Furthermore, we also present evidence that the copy number of isoacceptor tRNAs (those having a different anticodon but charge the same amino acid) is correlated with codon usage patterns of highly expressed genes in Nasonia.