Mutation accumulation in nuclear, organelle, and prokaryotic transfer RNA genes

A comparative analysis of the transfer RNA genes in the genomes of the major kingdoms of eukaryotes and prokaryotes leads to the general conclusion that the rate of evolution of organelle tRNA genes is typically equal to of greater than that of their nuclear counterparts. Situations where this is not the case, most notably in vascular plants, are attributable to an elevated mutation rate in the nuclear genome. Through a comparison of rates of mutation with rates of nucleotide substitution, it is shown that there is a reduction in the efficiency of selection on new mutations in organelle genes. Numerous lines of evidence, including observed reductions in stem duplex stability and changes in loop sizes, suggest that the excess changes observed in the organelle genes are mildly deleterious. Uniparental inheritance of organelles causes a reduction in the efficiency of selection through the joint effects of an increase in linkage disequilibrium and a decrease in effective population size. These results provide molecular support for the idea that asexually propagating genomes are subject to long-term, gradual fitness loss and raise questions about the role of organelle mutations in the long-term survival of major phylogenetic lineages.      

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Fragile T-stem in disease-associated human mitochondrial tRNA sensitizes structure to local and distant mutations

Mutations in human mitochondrial isoleucine tRNA (hs mt tRNA(Ile)) are associated with cardiomyopathy and opthalmoplegia. A recent study showed that opthalmoplegia-related mutations gave rise to severe decreases in aminoacylation efficiencies and that the defective mutant tRNAs were effective inhibitors of aminoacylation of the wild-type substrate. The results suggested that the effectiveness of the mutations was due in large part to an inherently fragile mitochondrial tRNA structure. Here, we investigate mutant tRNAs associated with cardiomyopathy, and a series of rationally designed second-site substitutions introduced into both opthalmoplegia- and cardiomyopathy-related mutant tRNAs. A source of structural fragility was uncovered. An inherently unstable T-stem appears susceptible to misalignments. This susceptibility sensitizes both domains of the L-shaped tRNA structure to base substitutions that are deleterious. Thus, the fragile T-stem makes the structure of this human mitochondrial tRNA particularly vulnerable to local and distant mutations.     

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 (ΔG) was calculated using the free energy rules as implemented in mfold version 3.1. The ΔG’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 ΔG 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.Reviewing Editor: Dr. Rüdiger Cerff     

Estimating the distribution of selection coefficients from phylogenetic data using sitewise mutation-selection models

Estimation of the distribution of selection coefficients of mutations is a long-standing issue in molecular evolution. In addition to population-based methods, the distribution can be estimated from DNA sequence data by phylogenetic-based models. Previous models have generally found unimodal distributions where the probability mass is concentrated between mildly deleterious and nearly neutral mutations. Here we use a sitewise mutation-selection phylogenetic model to estimate the distribution of selection coefficients among novel and fixed mutations (substitutions) in a data set of 244 mammalian mitochondrial genomes and a set of 401 PB2 proteins from influenza. We find a bimodal distribution of selection coefficients for novel mutations in both the mitochondrial data set and for the influenza protein evolving in its natural reservoir, birds. Most of the mutations are strongly deleterious with the rest of the probability mass concentrated around mildly deleterious to neutral mutations. The distribution of the coefficients among substitutions is unimodal and symmetrical around nearly neutral substitutions for both data sets at adaptive equilibrium. About 0.5% of the nonsynonymous mutations and 14% of the nonsynonymous substitutions in the mitochondrial proteins are advantageous, with 0.5% and 24% observed for the influenza protein. Following a host shift of influenza from birds to humans, however, we find among novel mutations in PB2 a trimodal distribution with a small mode of advantageous mutations.     

Avoidance of antisense, antiterminator tRNA anticodons in vertebrate mitochondria

Protein synthesis (translation) stops at stop codons, codons not complemented by tRNA anticodons. tRNAs matching stops, antitermination (Ter) tRNAs, prevent translational termination, producing dysfunctional proteins. Genomes avoid tRNAs with anticodons whose complement (the anticodon of the ‘antisense’ tRNA) matches stops. This suggests that antisense tRNAs, which also form cloverleaves, are occasionally expressed. Mitochondrial antisense tRNA expression is plausible, because both DNA strands are transcribed as single RNAs, and tRNA structures signal RNA maturation. Results describe potential antisense Ter tRNAs in mammalian mitochondrial genomes detected by tRNAscan-SE, and evidence for adaptations preventing translational antitermination: genomes possessing Ter tRNAs use less corresponding stop codons; antisense Ter tRNAs form weaker cloverleaves than homologuous non-Ter antisense tRNAs; and genomic stop codon usages decrease with stabilities of codon-anticodon interactions and of Ter tRNA cloverleaves. This suggests that antisense tRNAs frequently function in translation. Results suggest that opposite strand coding is exceptional in modern genes, yet might be frequent for mitochondrial tRNAs. This adds antisense tRNA templating to other mitochondrial tRNA functions: sense tRNA templating, formation and regulation of secondary (light strand DNA) replication origins. Antitermination probably affects mitochondrial degenerative diseases and ageing: pathogenic mutations are twice as frequent in tRNAs with antisense Ter anticodons than in other tRNAs, and species lacking mitochondrial antisense Ter tRNAs have longer mean maximal lifespans than those possessing antisense Ter tRNAs.     

Evidence for inefficient selection against deleterious mutations in cytochrome oxidase I of asexual bdelloid rotifers

Evolutionary theory predicts that natural selection should be less efficient in asexually than in sexually reproducing organisms. Obligate asexuals are expected to adapt slowly to changing environments and to accumulate mildly deleterious mutations to their genomes, potentially explaining their typically short evolutionary lifespans. One group of animals that appear to challenge these ideas is the bdelloid rotifers, a large and ancient clade of obligate asexuals. Previous work has found no evidence for inefficient selection against deleterious mutations in protein-coding genes of bdelloids. However, these studies relied mostly on between-species comparisons and were therefore unable to detect mildly deleterious mutations that persist within populations but are removed by selection over longer time periods. Here, we test for inefficient purifying selection acting on the cytochrome oxidase I (cox1) mitochondrial gene in 3 clades of bdelloids. Patterns of variation are compared to those of two facultatively sexual clades: a monogonont rotifer (Brachionus) and a branchiopod crustacean (Daphnia). As predicted due to the strict linkage between mitochondrial and nuclear genomes, bdelloids exhibit higher frequencies of putatively deleterious amino acid polymorphism within populations than the two facultatively sexual clades. While the monophyly and age of bdelloids makes it hard to rule out other explanations for the observed differences, several possible confounding factors, such as differences in effective population size or patterns of codon usage, are shown not to explain the observed differences. We therefore conclude that bdelloid mitochondrial DNA variation does display the signature of inefficient selection expected of obligate asexuals.     

Increased variation in mtDNA in patients with familial sensorineural hearing impairment

Analyses of mitochondrial DNA (mtDNA) sequences have revealed non-neutral patterns, suggesting that many amino acid mutations in animal mtDNA may be mildly deleterious, but this has not been verified in human clinical series. Since sensorineural hearing impairment (SNHI) is a common manifestation in many of the syndromes caused by mutations in mtDNA, this may be regarded as the phenotype of choice in attempts to detect mutations that may have a mildly deleterious effect on mitochondrial function. We selected 32 subjects from among 117 unrelated SNHI patients with SNHI in maternal relatives by means of family history, determined the entire coding region sequence of mtDNA and compared the sequence variation with that in 32 haplogroup-matched controls taken at random from 192 Finnish sequences. The 32 control sequences differed from the remaining 160 sequences by 36±9 substitutions (mean ± SD), while the difference for the 32 patients was 58±4 substitutions (P=0.005 for difference; Wilcoxon signed rank test). Differences were also found in the number of new haplotypes and new non-synonymous mutations or mutations in tRNA or rRNA genes. A total of 12 rare mtDNA variants were detected in the patients, and only 3 of these were considered to be neutral in effect. It is proposed that increased sequence variation in mtDNA may be a genetic risk factor for SNHI, and the increased frequency of rare haplotypes in these patients points to the presence of mildly deleterious mutations in mtDNA.     

High proportions of deleterious polymorphisms in constrained human genes

Previous studies on human mitochondrial genomes showed that the ratio of intra-specific diversities at nonsynonymous-to-synonymous positions was two to ten times higher than the ratio of interspecific divergences at these positions, suggesting an excess of slightly deleterious nonsynonymous polymorphisms. However, such an overabundance of nonsynonymous single nucleotide polymorphisms (SNPs) was not found in human nuclear genomes. Here, genome-wide estimates using >14,000 human-chimp nuclear genes and 1 million SNPs from four human genomes showed a significant proportion of deleterious nonsynonymous SNPs (～ 15%). Importantly, this study reveals a negative correlation between the magnitude of selection pressure and the proportion of deleterious SNPs on human genes. The proportion of deleterious amino acid replacement polymorphisms is 3.5 times higher in genes under high purifying selection compared with that in less constrained genes (28% vs. 8%). These results are explained by differences in the extent of contribution of mildly deleterious mutations to diversity and substitution.     

Mitochondrial translation in absence of local tRNA aminoacylation and methionyl tRNAMet formylation in Apicomplexa

Apicomplexans possess three translationally active compartments: the cytosol, a single tubular mitochondrion, and a vestigial plastid organelle called apicoplast. Mitochondrion and apicoplast are of bacterial evolutionary origin and therefore depend on a bacterial‐like translation machinery. The minimal mitochondrial genome contains only three ORFs, and in Toxoplasma gondii the absence of mitochondrial tRNA genes is compensated for by the import of cytosolic eukaryotic tRNAs. Although all compartments require a complete set of charged tRNAs, the apicomplexan nuclear genomes do not hold sufficient aminoacyl‐tRNA synthetase (aaRSs) genes to be targeted individually to each compartment. This study reveals that aaRSs are either cytosolic, apicoplastic or shared between the two compartments by dual targeting but are absent from the mitochondrion. Consequently, tRNAs are very likely imported in their aminoacylated form. Furthermore, the unexpected absence of tRNA^Met formyltransferase and peptide deformylase implies that the requirement for a specialized formylmethionyl‐tRNA^Met for translation initiation is bypassed in the mitochondrion of Apicomplexa.     

Human mitochondrial tRNA quality control in health and disease: a channelling mechanism?

Mutations in human mitochondrial tRNA genes are associated with a number of multisystemic disorders. These single nucleotide substitutions in various domains of tRNA molecules may affect different steps of tRNA biogenesis. Often, the prominent decrease of aminoacylation and/or steady-state levels of affected mitochondrial tRNA have been demonstrated in patients' tissues and in cultured cells. Similar effect has been observed for pathogenic mutations in nuclear genes encoding mitochondrial aminoacyl-tRNA-synthetases, while over-expression of mitochondrial aminoacyl-tRNA synthetases or elongation factor EF-Tu rescued mutated tRNAs from degradation. In this review we summarize experimental data concerning the possible regulatory mechanisms governing mitochondrial tRNA steady-state levels, and propose a hypothesis based on the tRNA channelling principle. According to this hypothesis, interaction of mitochondrial tRNA with proteins ensures not only tRNA synthesis, maturation and function, but also protection from degradation. Mutations perturbing this interaction lead to decreased tRNA stability.     

Mamit-tRNA, a database of mammalian mitochondrial tRNA primary and secondary structures

Mamit-tRNA (http://mamit-tRNA.u-strasbg.fr), a database for mammalian mitochondrial genomes, has been developed for deciphering structural features of mammalian mitochondrial tRNAs and as a helpful tool in the frame of human diseases linked to point mutations in mitochondrial tRNA genes. To accommodate the rapid growing availability of fully sequenced mammalian mitochondrial genomes, Mamit-tRNA has implemented a relational database, and all annotated tRNA genes have been curated and aligned manually. System administrative tools have been integrated to improve efficiency and to allow real-time update (from GenBank Database at NCBI) of available mammalian mitochondrial genomes. More than 3000 tRNA gene sequences from 150 organisms are classified into 22 families according to the amino acid specificity as defined by the anticodon triplets and organized according to phylogeny. Each sequence is displayed linearly with color codes indicating secondary structural domains and can be converted into a printable two-dimensional (2D) cloverleaf structure. Consensus and typical 2D structures can be extracted for any combination of primary sequences within a given tRNA specificity on the basis of phylogenetic relationships or on the basis of structural peculiarities. Mamit-tRNA further displays static individual 2D structures of human mitochondrial tRNA genes with location of polymorphisms and pathology-related point mutations. The site offers also a table allowing for an easy conversion of human mitochondrial genome nucleotide numbering into conventional tRNA numbering. The database is expected to facilitate exploration of structure/function relationships of mitochondrial tRNAs and to assist clinicians in the frame of pathology-related mutation assignments.     

A nuclear encoded tRNA of Trypanosoma brucei is imported into mitochondria.

The mitochondrial genome of trypanosomes, unlike that of most other eukaryotes, does not appear to encode any tRNAs. Therefore, mitochondrial tRNAs must be either imported into the organelle or created through a novel mitochondrial process, such as RNA editing. Trypanosomal tRNA(Tyr), whose gene contains an 11-nucleotide intron, is present in both the cytosol and the mitochondrion and is encoded by a single-copy nuclear gene. By site-directed mutagenesis, point mutations were introduced into this tRNA gene, and the mutated gene was reintroduced into the trypanosomal nuclear genome by DNA transfection. Expression of the mutant tRNA led to the accumulation of unspliced tRNA(Tyr) (A. Schneider, K. P. McNally, and N. Agabian, J. Biol. Chem. 268:21868-21874, 1993). Cell fractionation revealed that a significant portion of the unspliced mutant tRNA(Tyr) was recovered in the mitochondrial fraction and was resistant to micrococcal nuclease treatment in the intact organelle. Expression of the nuclear integrated, mutated tRNA gene and recovery of its gene product in the mitochondrial fraction directly demonstrated import. In vitro experiments showed that the unspliced mutant tRNA(Tyr), in contrast to the spliced wild-type form, was no longer a substrate for the cognate aminoacyl synthetase. The presence of uncharged tRNA in the mitochondria demonstrated that aminoacylation was not coupled to import.     

The mitochondrial tRNAs of Trypanosoma brucei are nuclear encoded

The mitochondrial DNA of Trypanosoma brucei is organized as a catenated network of maxicircles and minicircles. The maxicircles are equivalent to the typical mitochondrial genome except that the genes for the mitochondrial tRNAs have not been identified by sequence analysis of the maxicircle DNA. The apparent absence of tRNA genes in the maxicircle DNA suggests that the mitochondrial tRNAs are encoded by either the minicircle or the nuclear DNA. In order to determine their genomic origin, we isolated and identified the mitochondrial tRNAs of T. brucei. We show that these mitochondrial tRNAs are truly mitochondrially located in vivo and that they are free from detectable contamination by cytosolic RNAs. By hybridization analysis, using mitochondrial tRNAs as the probe, we determined that the mitochondrial tRNAs are encoded by nuclear DNA. This implies that RNAs, like proteins, are imported into the mitochondria. We investigated the relationship between the cytosolic and the mitochondrial tRNA genes and show that there are unique cytosolic tRNA genes, unique mitochondrial tRNA genes, and tRNA genes which appear to be shared and whose products are therefore targeted to both the cytosol and the mitochondrion.     

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.     

The Rates of Molecular Evolution in Rodent and Primate Mitochondrial DNA

A higher rate of molecular evolution in rodents than in primates at synonymous sites and, to a lesser extent, at amino acid replacement sites has been reported previously for most nuclear genes examined. Thus in these genes the average ratio of amino acid replacement to synonymous substitution rates in rodents is lower than in primates, an observation at odds with the neutral model of molecular evolution. Under Ohta's mildly deleterious model of molecular evolution, these observations are seen as the consequence of the combined effects of a shorter generation time (driving a higher mutation rate) and a larger effective population size (resulting in more effective selection against mildly deleterious mutations) in rodents. The present study reports the results of a maximum-likelihood analysis of the ratio of amino acid replacements to synonymous substitutions for genes encoded in mitochondrial DNA (mtDNA) in these two lineages. A similar pattern is observed: in rodents this ratio is significantly lower than in primates, again consistent only with the mildly deleterious model. Interestingly the lineage-specific difference is much more pronounced in mtDNA-encoded than in nuclear-encoded proteins, an observation which is shown to run counter to expectation under Ohta's model. Finally, accepting certain fossil divergence dates, the lineage-specific difference in amino acid replacement-to-synonymous substitution ratio in mtDNA can be partitioned and is found to be entirely the consequence of a higher mutation rate in rodents. This conclusion is consistent with a replication-dependent model of mutation in mtDNA. Received: 24 September 1999 / Accepted: 18 September 2000     

Pathogenic mutations in antisense mitochondrial tRNAs

Pathogenic mutations in mitochondrial tRNAs are 6.5 times more frequent than in other mitochondrial genes. This suggests that tRNA mutations perturb more than one function. A potential additional tRNA gene function is that of templating for antisense tRNAs. Pathogenic mutations weaken cloverleaf secondary structures of sense tRNAs. Analyses here show similar effects for most antisense tRNAs, especially after adjusting for associations between sense and antisense cloverleaf stabilities. These results imply translational activity by antisense tRNAs. For sense tRNAs Ala and Ser UCN, pathogenicity associates as much with sense as with antisense cloverleaf formation. For tRNA Pro, pathogenicity seems associated only with antisense, not sense tRNA cloverleaf formation. Translational activity by antisense tRNAs is expected for the 11 antisense tRNAs processed by regular sense RNA maturation, those recognized by their cognate amino acid’s tRNA synthetase, and those forming relatively stable cloverleaves as compared to their sense counterpart. Most antisense tRNAs probably function routinely in translation and extend the tRNA pool (extension hypothesis); others do not (avoidance hypothesis). The greater the expected translational activity of an antisense tRNA, the more pathogenic mutations weaken its cloverleaf secondary structure. Some evidence for RNA interference, a more classical role for antisense tRNAs, exists only for tRNA Ser UCN. Mutation pathogenicity probably frequently results from a mixture of effects due to sense and antisense tRNA translational activity for many mitochondrial tRNAs. Genomic studies should routinely explore for translational activity by antisense tRNAs.