Allelic variation at tRNA genes in three nematode species indicates mutation load despite strong purifying selection

Cytosolic transfer RNAs, which are encoded as hundreds of genes in nuclear genomes, experience exceptionally high mutation rates and have been hypothesized to confer substantial mutational load in natural populations. Although this phenomenon appears universal across multicellular eukaryotes, a comprehensive characterization of standing variation in tRNA repertoires is still lacking in any system. Here, we resolve within-species allelic variation in nuclear-encoded tRNAs in three nematode species: Caenorhabditis elegans, C. briggsae, and C. tropicalis. We show that these genes carry signatures of high rates of historical transcription-associated mutagenesis and of purifying selection, resulting in allelic variation that includes pervasive instances of within-gene mismatches between the amino acid recognized by the tRNA backbone and that indicated by the anticodon. Furthermore, patterns of tRNA genomic organization and variation differ markedly from those of protein-coding regions. Individual genomes harbor distinct complements of tRNA genes with predicted functional differences, an observation that coincides with recent evidence that variation in tRNA expression and regulation contributes to human disease. Our findings offer an entry point for identifying the micro-evolutionary processes that act on tRNA repertoires, and in turn connecting those processes to the macro-evolutionary patterns that have more frequently been the focus of study. Significance statementTransfer RNAs (tRNAs) are ancient and essential molecules that deliver amino acids to the growing polypeptide chain during protein synthesis. Despite their importance, little is known about how the genes encoding tRNAs vary among individuals, even though they accumulate large numbers of mutations due to their high rates of expression. Here, we report that three nematode species, including the model organism C. elegans, show within-species variation in their tRNA gene repertoires, which may cause differences in protein synthesis across individuals. Our findings may represent patterns common across multicellular species and they set the stage for future investigation into how this diversity influences cellular function and organismal fitness.