Mitochondrial tRNA sequences as unusual replication origins: pathogenic implications for Homo sapiens

The heavy strand of vertebrate mitochondrial genomes accumulates deaminations proportionally to the time it spends single-stranded during replication. A previous study showed that the strength of genome-wide deamination gradients originating from tRNA gene's locations increases with their capacities to form secondary structures resembling mitochondrial origins of light strand replication (OL), suggesting an alternative function for tRNA sequences. We hypothesize that this function is frequently pathogenic for those tRNA genes that normally do not form OL-like structures, because this could cause excess mutations in genome regions unadapted to tolerate them. In human mitochondrial genomes, pathogenic tRNA variants usually form less OL-like structures than non-pathogenic ones in cases where the normal non-pathogenic tRNA variant can function as OL, as evolutionary analyses reveal. For tRNAs lacking the putative OL-like functioning capacity, pathogenic variants form more OL-like secondary structures, particularly structures that might invoke bi-directional replication (true for 14 among 21 tRNA species, p<0.05, sign test; significantly at p<0.05 (1 tailed test) for 7 tRNA species), but not more unidirectional replication invoking structures. Accounting for the functional cloverleaf-like structure-forming capacities of tRNAs yields similar results. Rare, non-pathogenic tRNA mutants tend to form more OL-like structures than the common, non-pathogenic ones, suggesting weak directional selection also among non-pathogenic variants. The duration spent single stranded by a region of the heavy strand (D(ssH)) during replication, estimated by integrating over all regions that can function as OL in Homo sapiens mitochondrial genomes, increases with distance of that region from the Dloop. This suggests convergence of single-strandedness during replication and transcription, and explains conserved locations of tRNA species in mitochondrial genomes and bacterial operons. These locations minimize deamination costs only in anticodons and not in other tRNA regions, during replication and transcription. Therefore, putative functioning as OLs by tRNA sequences is normal at some locations and pathogenic at others.