Sequence variation in the tRNA genes of human mitochondrial DNA

Recent analyses have shown that nonsynonymous variation in human mitochondrial DNA (mtDNA) contains nonneutral variants, suggesting the presence of mildly deleterious mutations. Many of the disease-causing mutations in mtDNA occur in the genes encoding the tRNAs. Nucleotide sequence variation in these genes has not been studied in human populations, nor have the structural consequences of nucleotide substitutions in tRNA molecules been examined. We therefore determined the nucleotide sequences of the 22 tRNA genes in the mtDNA of 477 Finns and, also, obtained 435 European sequences from the MitoKor database. No differences in population polymorphism indices were found between the two data sets. We assessed selective constraints against various tRNA domains by comparing allele frequencies between these domains and the synonymous and nonsynonymous sites, respectively. All tRNA domains except the variable loop were more conserved than synonymous sites, and T stem and D stem were more conserved than the respective loops. We also analyzed the energetic consequences of the 96 polymorphisms recovered in the two data sets or in the Mitomap database. The minimum free energy (DeltaG) was calculated using the free energy rules as implemented in mfold version 3.1. The DeltaG's were normally distributed among the 22 wild-type tRNA genes, whereas the 96 polymorphic tRNAs departed significantly from a normal distribution. The largest differences in DeltaG between the wild-type and the polymorphic tRNAs in the Finnish population tended to be in the polymorphisms that were present at low frequencies. Allele frequency distributions and minimum free energy calculations both suggested that some polymorphisms in tRNA genes are nonneutral.