Independent genome reduction and phylogenetic reclassification of the oceanic SAR11 clade

The SAR11 clade, here represented by Candidatus Pelagibacter ubique, is the most successful group of bacteria in the upper surface waters of the oceans. In contrast to previous studies that have associated the 1.3 Mb genome of Ca. Pelagibacter ubique with the less than 1.5 Mb genomes of the Rickettsiales, our phylogenetic analysis suggests that Ca. Pelagibacter ubique is most closely related to soil and aquatic Alphaproteobacteria with large genomes. This implies that the SAR11 clade and the Rickettsiales have undergone genome reduction independently. A gene flux analysis of 46 representative alphaproteobacterial genomes indicates the loss of more than 800 genes in each of Ca. Pelagibacter ubique and the Rickettsiales. Consistent with their different phylogenetic affiliations, the pattern of gene loss differs with a higher loss of genes for repair and recombination processes in Ca. Pelagibacter ubique as compared with a more extensive loss of genes for biosynthetic functions in the Rickettsiales. Some of the lost genes in Ca. Pelagibacter ubique, such as mutLS, recFN, and ruvABC, are conserved in all other alphaproteobacterial genomes including the small genomes of the Rickettsiales. The mismatch repair genes mutLS are absent from all currently sequenced SAR11 genomes and also underrepresented in the global ocean metagenome data set. We hypothesize that the unique loss of genes involved in repair and recombination processes in Ca. Pelagibacter ubique has been driven by selection and that this helps explain many of the characteristics of the SAR11 population, such as the streamlined genomes, the long branch lengths, the high recombination frequencies, and the extensive sequence divergence within the population.     

TFAM detects co-evolution of tRNA identity rules with lateral transfer of histidyl-tRNA synthetase

We present TFAM, an automated, statistical method to classify the identity of tRNAs. TFAM, currently optimized for bacteria, classifies initiator tRNAs and predicts the charging identity of both typical and atypical tRNAs such as suppressors with high confidence. We show statistical evidence for extensive variation in tRNA identity determinants among bacterial genomes due to variation in overall tDNA base content. With TFAM we have detected the first case of eukaryotic-like tRNA identity rules in bacteria. An alpha-proteobacterial clade encompassing Rhizobiales, Caulobacter crescentus and Silicibacter pomeroyi, unlike a sister clade containing the Rickettsiales, Zymomonas mobilis and Gluconobacter oxydans, uses the eukaryotic identity element A73 instead of the highly conserved prokaryotic element C73. We confirm divergence of bacterial histidylation rules by demonstrating perfect covariation of alpha-proteobacterial tRNA(His) acceptor stems and residues in the motif IIb tRNA-binding pocket of their histidyl-tRNA synthetases (HisRS). Phylogenomic analysis supports lateral transfer of a eukaryotic-like HisRS into the alpha-proteobacteria followed by in situ adaptation of the bacterial tDNA(His) and identity rule divergence. Our results demonstrate that TFAM is an effective tool for the bioinformatics, comparative genomics and evolutionary study of tRNA identity.     

Comparative and Phylogenomic Evidence That the Alphaproteobacterium HIMB59 Is Not a Member of the Oceanic SAR11 Clade

SAR11 is a globally abundant group of Alphaproteobacteria in the oceans that is taxonomically not well defined. It has been suggested SAR11 should be classified into the novel order Pelagibacterales. Features such as conservation of gene content and synteny have been taken as evidence that also the divergent member HIMB59 should be included in the order. However, this proposition is controversial since phylogenetic analyses have questioned the monophyly of this grouping. Here, we performed phylogenetic analyses and reinvestigated the genomic similarity of SAR11 and HIMB59. Our phylogenetic analysis confirmed that HIMB59 is not a sister group to the other SAR11 strains. By placing the comparison in the context of the evolution of the Alphaproteobacteria, we found that none of the measures of genomic similarity supports a clustering of HIMB59 and SAR11 to the exclusion of other Alphaproteobacteria. First, pairwise sequence similarity measures for the SAR11 and HIMB59 genomes were within the range observed for unrelated pairs of Alphaproteobacteria. Second, pairwise comparisons of gene contents revealed a higher similarity of SAR11 to several other alphaproteobacterial genomes than to HIMB59. Third, the SAR11 genomes are not more similar in gene order to the HIMB59 genome than what they are to several other alphaproteobacterial genomes. Finally, in contrast to earlier reports, we observed no sequence similarity between the hypervariable region HVR2 in the SAR11 genomes and the region located at the corresponding position in the HIMB59 genome. Based on these observations, we conclude that the alphaproteobacterium HIMB59 is not monophyletic with the SAR11 strains and that genome streamlining has evolved multiple times independently in Alphaproteobacteria adapted to the upper surface waters of the oceans.     

Association between translation efficiency and horizontal gene transfer within microbial communities

Horizontal gene transfer (HGT) is a major force in microbial evolution. Previous studies have suggested that a variety of factors, including restricted recombination and toxicity of foreign gene products, may act as barriers to the successful integration of horizontally transferred genes. This study identifies an additional central barrier to HGT-the lack of co-adaptation between the codon usage of the transferred gene and the tRNA pool of the recipient organism. Analyzing the genomic sequences of more than 190 microorganisms and the HGT events that have occurred between them, we show that the number of genes that were horizontally transferred between organisms is positively correlated with the similarity between their tRNA pools. Those genes that are better adapted to the tRNA pools of the target genomes tend to undergo more frequent HGT. At the community (or environment) level, organisms that share a common ecological niche tend to have similar tRNA pools. These results remain significant after controlling for diverse ecological and evolutionary parameters. Our analysis demonstrates that there are bi-directional associations between the similarity in the tRNA pools of organisms and the number of HGT events occurring between them. Similar tRNA pools between a donor and a host tend to increase the probability that a horizontally acquired gene will become fixed in its new genome. Our results also suggest that frequent HGT may be a homogenizing force that increases the similarity in the tRNA pools of organisms within the same community.     

The SAR11 group of alpha-proteobacteria is not related to the origin of mitochondria

Although free living, members of the successful SAR11 group of marine alpha-proteobacteria contain a very small and A+T rich genome, two features that are typical of mitochondria and related obligate intracellular parasites such as the Rickettsiales. Previous phylogenetic analyses have suggested that Candidatus Pelagibacter ubique, the first cultured member of this group, is related to the Rickettsiales+mitochondria clade whereas others disagree with this conclusion. In order to determine the evolutionary position of the SAR11 group and its relationship to the origin of mitochondria, we have performed phylogenetic analyses on the concatenation of 24 proteins from 5 mitochondria and 71 proteobacteria. Our results support that SAR11 group is not the sistergroup of the Rickettsiales+mitochondria clade and confirm that the position of this group in the alpha-proteobacterial tree is strongly affected by tree reconstruction artefacts due to compositional bias. As a consequence, genome reduction and bias toward a high A+T content may have evolved independently in the SAR11 species, which points to a different direction in the quest for the closest relatives to mitochondria and Rickettsiales. In addition, our analyses raise doubts about the monophyly of the newly proposed Pelagibacteraceae family.     

The leaky UGA termination codon of tobacco rattle virus RNA is suppressed by tobacco chloroplast and cytoplasmic tRNAs(Trp) with CmCA anticodon.

RNA-1 molecules from tobacco rattle virus (TRV) and pea early-browning virus (PEBV), two members of the tobravirus group, have recently been shown to contain internal, in-frame UGA termination codons which are suppressed in vitro. Our results suggest that a UGA stop codon also exists in RNA-1 of pepper ringspot virus (PRV), another tobravirus. UGA suppression may therefore be a universal feature of the expression of tobravirus genomes. We have isolated two natural suppressor tRNAs from uninfected tobacco plants on the basis of their ability to promote readthrough over the leaky UGA codon of TRV RNA-1 in a wheat germ extract depleted of endogenous mRNAs and tRNAs. Their amino acid acceptance and nucleotide sequences identify the two UGA-suppressor tRNAs as chloroplast (chl) and cytoplasmic (cyt) tryptophan-specific tRNAs with the anticodon CmCA. These are the first UGA suppressor tRNAs to be identified in plants. They have several interesting features. (i) Chl tRNA(Trp) suppresses the UGA stop codon more efficiently than cyt tRNA(Trp). (ii) Chl tRNA(Trp) contains an A24:U11 pair in the D-stem as does the mutated Escherichia coli UGA-suppressor tRNA(Trp) which is a more active suppressor than wild-type tRNA(Trp). (iii) The suppressor activity of chl tRNA(Trp) is dependent on the nucleotides surrounding the stop codon because it recognizes UGA in the TRV context but not the UGA in the beta-globin context.     

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.     

Comparison of base composition and codon usage in insect mitochondrial genomes

Insects, the most biodiverse taxonomic group, have high AT content in their mitochondrial genomes. Although codon usage tends to be AT-rich, base composition and codon usage of mitochondrial genomes may vary among taxa. Thus, we compare base composition and codon usage patterns of 49 insect mitochondrial genomes. For protein coding genes, AT content is as high as 80% in the Hymenoptera and Lepidoptera and as low as 72% in the Orthopotera. The AT content is high at positions 1 and 3, but A content is low at position 2. A close correlation occurs between codon usage and tRNA abundance in nuclear genomes. Optimal codons can pair well with the antr codons of the most abundant tRNAs. One tRNA gene translates a synonymous codon family in vertebrate mitochondrial genomes and these tRNA anticodons can pair with optimal codons. However, optimal codons cannot pair with anticodons in mtDNA ofCochiiomyia hominivorax (Dipteral: CaLliphoridae). Ten optimal codons cannot pair with tRNA anticodons in all 49 insect mitochondrial genomes; non-optimal codon-anticodon usage is common and codon usage is not influenced by tRNA abundance.     

Transfer RNA structure and coding specificity. II. A D-arm tertiary interaction that restricts coding range

We investigated the structural basis of the kinetic effect on coding specificity by the D-arm mutant (G24 to A) of Escherichia coli tRNATrp. A set of tRNA genes with structural alterations in the D-arm was constructed by site-directed mutagenesis in vitro, and we determined the in vivo translational activities of these tRNAs. Our results suggest that a hydrogen-bond donor in the major groove of the D-helix at position 24 is required for the expansion of tRNA wobble coding specificity. From inspection of tRNA crystal structure, we identified a potential new tertiary pairing of base 24 with the base at position 9 (this base links the acceptor and D-stems). We constructed tRNAs with mutations at position 9 and showed that the phenotypes of position 11-24 D-arm mutants are indeed dependent on the identity of base 9. Our analysis of the effects of these mutations on the interactions of tRNA with the ribosome and with aminoacyl-tRNA synthetase suggests that the conformation or conformational dynamics of the middle of the tRNA molecule alters the kinetics of the interaction with the ribosomal coding site. The 9-23 and putative 9-24 tertiaries, and perhaps other normal tertiary interactions in this region, modulate these kinetics to increase or decrease coding specificity.     

Identification and characterization of a tRNA decoding the rare AUA codon in Haloarcula marismortui

Annotation of the complete genome of the extreme halophilic archaeon Haloarcula marismortui does not include a tRNA for translation of AUA, the rare codon for isoleucine. This is a situation typical for most archaeal genomes sequenced to date. Based on computational analysis, it has been proposed recently that a single intron-containing tRNA gene produces two very similar but functionally different tRNAs by means of alternative splicing; a UGG-decoding tRNA(TrpCCA) and an AUA-decoding tRNA(IleUAU). Through analysis of tRNAs from H. marismortui, we have confirmed the presence of tRNA(TrpCCA), but found no evidence for the presence of tRNA(IleUAU). Instead, we have shown that a tRNA, currently annotated as elongator methionine tRNA and containing CAU as the anticodon, is aminoacylated with isoleucine in vivo and that this tRNA represents the missing isoleucine tRNA. Interestingly, this tRNA carries a base modification of C34 in the anticodon different from the well-known lysidine found in eubacteria, which switches the amino acid identity of the tRNA from methionine to isoleucine and its decoding specificity from AUG to AUA. The methods described in this work for the identification of individual tRNAs present in H. marismortui provide the tools necessary for experimentally confirming the presence of any tRNA in a cell and, thereby, to test computational predictions of tRNA genes.     

A model of horizontal gene transfer and the bacterial phylogeny problem

How much horizontal gene transfer (HGT) between species influences bacterial phylogenomics is a controversial issue. This debate, however, lacks any quantitative assessment of the impact of HGT on phylogenies and of the ability of tree-building methods to cope with such events. I introduce a Markov model of genome evolution with HGT, accounting for the constraints on time -- an HGT event can only occur between concomitantly living species. This model is used to simulate multigene sequence data sets with or without HGT. The consequences of HGT on phylogenomic inference are analyzed and compared to other well-known phylogenetic artefacts. It is found that supertree methods are quite robust to HGT, keeping high levels of performance even when gene trees are largely incongruent with each other. Gene tree incongruence per se is not indicative of HGT. HGT, however, removes the (otherwise observed) positive relationship between sequence length and gene tree congruence to the estimated species tree. Surprisingly, when applied to a bacterial and a eukaryotic multigene data set, this criterion rejects the HGT hypothesis for the former, but not the latter data set.     

Isolation of novel tRNA(Ala) mutants by library selection in a tRNA(Ala) knockout strain

The relationship between tRNA structure and function has been widely investigated by site-directed mutagenesis. This method has been a very useful tool to reveal the critical bases in tRNAs that are important for recognition and aminoacylation, but has been limited by the large number of possible base combinations in tRNA molecules. We have devised a new method that uses tRNA knockout cells for selection of functional tRNAs from a mutant tRNA gene library to overcome this limitation. To explore the mechanism of tRNA(Ala) recognition, the bases of the acceptor-stem region were randomized and active mutants were selected in a tRNA(Ala) knockout strain. Mutants of tRNA(Ala) having diverse sequence combinations in the acceptor-stem region and a broad range of functional activity to support knockout cell growth were isolated. The mutant tRNAs selected by the method included molecules containing novel base substitutions as well as extensively altered base combinations that would not be readily generated by rationally designed site-directed mutagenesis. Our results emphasize the importance of the acceptor stem as a structural unit in which some nucleotides may carry more weight than others, but in summation every nucleotide contributes to the interaction with the enzyme.     

The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Determining the influence of horizontal gene transfer (HGT) on phylogenomic analyses and the retrieval of a tree of life is relevant for our understanding of microbial genome evolution. It is particularly difficult to differentiate between phylogenetic incongruence due to noise and that resulting from HGT. We have performed a large-scale, detailed evolutionary analysis of the different phylogenetic signals present in the genomes of Xanthomonadales, a group of Proteobacteria. We show that the presence of phylogenetic noise is not an obstacle to infer past and present HGTs during their evolution. The scenario derived from this analysis and other recently published reports reflect the confounding effects on bacterial phylogenomics of past and present HGT. Although transfers between closely related species are difficult to detect in genome-scale phylogenetic analyses, past transfers to the ancestor of extant groups appear as conflicting signals that occasionally might make impossible to determine the evolutionary origin of the whole genome.     

Another cut for lysine tRNA: application of the hyperprocessing reaction reveals another stabilization strategy in metazoan lysine tRNAs

Recently, we revealed that the cloverleaf structure of some eukaryotic tRNAs is not always stable in vitro, and the denatured structures of these tRNAs are sometimes detected in bacterial RNase P reactions. We have designated the unusual internal cleavage reaction of these tRNAs as hyperprocessing. We have developed this hyperprocessing strategy as a useful tool for examining the stability of the tRNA cloverleaf structure. There are some common features in such unstable, hyperprocessible tRNAs, and the criteria for the hyperprocessing reaction of tRNA are extracted. Metazoan initiator methionine tRNAs and lysine tRNAs commonly fit the criteria, and are predicted to be hyperprocessible. The RNase P reactions of two metazoan lysine tRNAs from Homo sapiens and Caenorhabditis elegans, which fit the criteria, resulted in resistance to the internal cleavage reaction, while one bacterial lysine tRNA from Acholeplasma laidlawii, which also fits the criteria, was internally cleaved by the RNase P. The results showed that the metazoan lysine tRNAs examined are very stable without base modifications even under in vitro conditions. We also examined the 3'-half short construct of the human lysine tRNA, and the results showed that this RNA was internally cleaved by the enzyme. The results indicated that the human lysine tRNA has the ability to be hyperprocessed but is structurally stabilized in spite of lacking base modifications. A comparative study suggested, moreover, that the acceptor-stem bases should take part in the stabilization of metazoan lysine tRNAs. Our data strongly suggest that the cloverleaf shape of other metazoan lysine tRNAs should also be stabilized by means of similar strategies to in the case of human tRNA(Lys3).     

Tissue-specific differences in human transfer RNA expression

Over 450 transfer RNA (tRNA) genes have been annotated in the human genome. Reliable quantitation of tRNA levels in human samples using microarray methods presents a technical challenge. We have developed a microarray method to quantify tRNAs based on a fluorescent dye-labeling technique. The first-generation tRNA microarray consists of 42 probes for nuclear encoded tRNAs and 21 probes for mitochondrial encoded tRNAs. These probes cover tRNAs for all 20 amino acids and 11 isoacceptor families. Using this array, we report that the amounts of tRNA within the total cellular RNA vary widely among eight different human tissues. The brain expresses higher overall levels of nuclear encoded tRNAs than every tissue examined but one and higher levels of mitochondrial encoded tRNAs than every tissue examined. We found tissue-specific differences in the expression of individual tRNA species, and tRNAs decoding amino acids with similar chemical properties exhibited coordinated expression in distinct tissue types. Relative tRNA abundance exhibits a statistically significant correlation to the codon usage of a collection of highly expressed, tissue-specific genes in a subset of tissues or tRNA isoacceptors. Our findings demonstrate the existence of tissue-specific expression of tRNA species that strongly implicates a role for tRNA heterogeneity in regulating translation and possibly additional processes in vertebrate organisms.     

Identification of novel genomic islands associated with small RNAs

Genome evolution in prokaryotes is assisted by integration of gene pools from phages and plasmids. Regions downstream of tRNAs and tmRNAs are considered as hot spots for the integration of these gene pools or genomic islands. Till date, genomic islands have been identified only at tRNA/tmRNA genes in the enterobacterial genomes. Present work reports 10 distinct small RNAs as potent integration sites for genomic islands. A known tool tRNAcc 1.0 has been used to identify genomic islands associated with small RNAs c0362, oxyS, ryaA, rybB, rybD, ryeB, ryeE, rtT, sraE and tmRNA. The coordinates of 25 such small RNA associated genomic islands in three E. coli (strains: CFT073, EDL933 and K12) and Shigella flexneri (strain: 301) genomes are presented. Moreover cross-verification of the genomic sequences encoded within the identified genomic islands in horizontal gene transfer database, GenBank annotation features and atypical sequence compositions support our results. Again, all of the identified 25 genomic integration sites do exhibit genomic block rearrangements with respect to the associated small RNA. Similar to tRNAs/tmRNAs, the downstream regions of the small RNAs are found to be hotspots of integration.     

Characterization of a yeast mitochondrial locus necessary for tRNA biosynthesis. Deletion mapping and restriction mapping studies

Yeast mitochondrial DNA codes for a complete set of tRNAs. Although most components necessary for the biosynthesis of mitochondrial tRNA are coded by nuclear genes, there is one genetic locus on mitochondrial DNA necessary for the synthesis of mitochondrial tRNAs other than the mitochondrial tRNA genes themselves. Characterization of mutants by deletion mapping and restriction enzyme mapping studies has provided a precise location of this yeast mitochondrial tRNA synthesis locus. Deletion mutants retaining various segments of mitochondrial DNA were examined for their ability to synthesize tRNAs from the genes they retain. A subset of these strains was also tested for the ability to provide the tRNA synthesis function in complementation tests with deletion mutants unable to synthesize mature mitochondrial tRNAs. By correlating the tRNA synthetic ability with the presence or absence of certain wild-type restriction fragments, we have confined the locus to within 780 base pairs of DNA located between the tRNAMetf gene and tRNAPro gene, at 29 units on the wild-type map. Heretofore, no genetic function or gene product had been localized in this area of the yeast mitochondrial genome.     

A phylometagenomic exploration of oceanic alphaproteobacteria reveals mitochondrial relatives unrelated to the SAR11 clade

Background

According to the endosymbiont hypothesis, the mitochondrial system for aerobic respiration was derived from an ancestral Alphaproteobacterium. Phylogenetic studies indicate that the mitochondrial ancestor is most closely related to the Rickettsiales. Recently, it was suggested that Candidatus Pelagibacter ubique, a member of the SAR11 clade that is highly abundant in the oceans, is a sister taxon to the mitochondrial-Rickettsiales clade. The availability of ocean metagenome data substantially increases the sampling of Alphaproteobacteria inhabiting the oxygen-containing waters of the oceans that likely resemble the originating environment of mitochondria.

Methodology/Principal Findings

We present a phylogenetic study of the origin of mitochondria that incorporates metagenome data from the Global Ocean Sampling (GOS) expedition. We identify mitochondrially related sequences in the GOS dataset that represent a rare group of Alphaproteobacteria, designated OMAC (Oceanic Mitochondria Affiliated Clade) as the closest free-living relatives to mitochondria in the oceans. In addition, our analyses reject the hypothesis that the mitochondrial system for aerobic respiration is affiliated with that of the SAR11 clade.

Conclusions/Significance

Our results allude to the existence of an alphaproteobacterial clade in the oxygen-rich surface waters of the oceans that represents the closest free-living relative to mitochondria identified thus far. In addition, our findings underscore the importance of expanding the taxonomic diversity in phylogenetic analyses beyond that represented by cultivated bacteria to study the origin of mitochondria.