Two C or not two C: recurrent disruption of Zn-ribbons, gene duplication, lineage-specific gene loss, and horizontal gene transfer in evolution of bacterial ribosomal proteins

Background

Ribosomal proteins are encoded in all genomes of cellular life forms and are, generally, well conserved during evolution. In prokaryotes, the genes for most ribosomal proteins are clustered in several highly conserved operons, which ensures efficient co-regulation of their expression. Duplications of ribosomal-protein genes are infrequent, and given their coordinated expression and functioning, it is generally assumed that ribosomal-protein genes are unlikely to undergo horizontal transfer. However, with the accumulation of numerous complete genome sequences of prokaryotes, several paralogous pairs of ribosomal protein genes have been identified. Here we analyze all such cases and attempt to reconstruct the evolutionary history of these ribosomal proteins.

Results

Complete bacterial genomes were searched for duplications of ribosomal proteins. Ribosomal proteins L36, L33, L31, S14 are each duplicated in several bacterial genomes and ribosomal proteins L11, L28, L7/L12, S1, S15, S18 are so far duplicated in only one genome each. Sequence analysis of the four ribosomal proteins, for which paralogs were detected in several genomes, two of the ribosomal proteins duplicated in one genome (L28 and S18), and the ribosomal protein L32 showed that each of them comes in two distinct versions. One form contains a predicted metal-binding Zn-ribbon that consists of four conserved cysteines (in some cases replaced by histidines), whereas, in the second form, these metal-chelating residues are completely or partially replaced. Typically, genomes containing paralogous genes for these ribosomal proteins encode both versions, designated C+ and C-, respectively. Analysis of phylogenetic trees for these seven ribosomal proteins, combined with comparison of genomic contexts for the respective genes, indicates that in most, if not all cases, their evolution involved a duplication of the ancestral C+ form early in bacterial evolution, with subsequent alternative loss of the C+ and C- forms in different lineages. Additionally, evidence was obtained for a role of horizontal gene transfer in the evolution of these ribosomal proteins, with multiple cases of gene displacement 'in situ', that is, without a change of the gene order in the recipient genome.

Conclusions

A more complex picture of evolution of bacterial ribosomal proteins than previously suspected is emerging from these results, with major contributions of lineage-specific gene loss and horizontal gene transfer. The recurrent theme of emergence and disruption of Zn-ribbons in bacterial ribosomal proteins awaits a functional interpretation.     

Molecular evolution of rDNA in early diverging Metazoa: First comparative analysis and phylogenetic application of complete SSU rRNA secondary structures in Porifera

Background  

The cytoplasmic ribosomal small subunit (SSU, 18S) ribosomal RNA (rRNA) is the most frequently-used gene for molecular phylogenetic studies. However, information regarding its secondary structure is neglected in most phylogenetic analyses. Incorporation of this information is essential in order to apply specific rRNA evolutionary models to overcome the problem of co-evolution of paired sites, which violates the basic assumption of the independent evolution of sites made by most phylogenetic methods. Information about secondary structure also supports the process of aligning rRNA sequences across taxa. Both aspects have been shown to increase the accuracy of phylogenetic reconstructions within various taxa.     

RibAlign: a software tool and database for eubacterial phylogeny based on concatenated ribosomal protein subunits

Background  

Until today, analysis of 16S ribosomal RNA (rRNA) sequences has been the de-facto gold standard for the assessment of phylogenetic relationships among prokaryotes. However, the branching order of the individual phlya is not well-resolved in 16S rRNA-based trees. In search of an improvement, new phylogenetic methods have been developed alongside with the growing availability of complete genome sequences. Unfortunately, only a few genes in prokaryotic genomes qualify as universal phylogenetic markers and almost all of them have a lower information content than the 16S rRNA gene. Therefore, emphasis has been placed on methods that are based on multiple genes or even entire genomes. The concatenation of ribosomal protein sequences is one method which has been ascribed an improved resolution. Since there is neither a comprehensive database for ribosomal protein sequences nor a tool that assists in sequence retrieval and generation of respective input files for phylogenetic reconstruction programs, RibAlign has been developed to fill this gap.     

Molecular adaptation and expression evolution following duplication of genes for organellar ribosomal protein S13 in rosids

Background  

Gene duplication has been a fundamental process in the evolution of eukaryotic genomes. After duplication one copy (or both) can undergo divergence in sequence, expression pattern, and function. Two divergent copies of the ribosomal protein S13 gene (rps13) of chloroplast origin are found in the nucleus of the rosids Arabidopsis, Gossypium, and Glycine. One encodes chloroplast-imported RPS13 (nucp rps13), while the other encodes mitochondria-imported RPS13 (numit rps13). The rps13 gene has been lost from mitochondrial DNA (mt rps13) of many rosids.     

A genome-wide survey of changes in protein evolutionary rates across four closely related species of <Emphasis Type="Italic">Saccharomyces</Emphasis> sensu stricto group

Background  

Changes in protein evolutionary rates among lineages have been frequently observed during periods of notable phenotypic evolution. It is also known that, following gene duplication and loss, the protein evolutionary rates of genes involved in such events changed because of changes in functional constraints acting on the genes. However, in the evolution of closely related species, excluding the aforementioned situations, the frequency of changes in protein evolutionary rates is still not clear at the genome-wide level. Here we examine the constancy of protein evolutionary rates in the evolution of four closely related species of the Saccharomyces sensu stricto group (S. cerevisiae, S. paradoxus, S. mikatae and S. bayanus).     

Introducing W.A.T.E.R.S.: a Workflow for the Alignment,Taxonomy, and Ecology of Ribosomal Sequences

Background  

For more than two decades microbiologists have used a highly conserved microbial gene as a phylogenetic marker for bacteria and archaea. The small-subunit ribosomal RNA gene, also known as 16 S rRNA, is encoded by ribosomal DNA, 16 S rDNA, and has provided a powerful comparative tool to microbial ecologists. Over time, the microbial ecology field has matured from small-scale studies in a select number of environments to massive collections of sequence data that are paired with dozens of corresponding collection variables. As the complexity of data and tool sets have grown, the need for flexible automation and maintenance of the core processes of 16 S rDNA sequence analysis has increased correspondingly.     

Supernumerary proteins of mitochondrial ribosomes

Background

Messenger RNAs encoded by mitochondrial genomes are translated on mitochondrial ribosomes that have unique structure and protein composition compared to prokaryotic and cytoplasmic ribosomes. Mitochondrial ribosomes are a patchwork of core proteins that share homology with prokaryotic ribosomal proteins and new, supernumerary proteins that can be unique to different organisms. In mammals, there are specific supernumerary ribosomal proteins that are not present in other eukaryotes.

Scope of review

Here we discuss the roles of supernumerary proteins in the regulation of mitochondrial gene expression and compare them among different eukaryotic systems. Furthermore, we consider if differences in the structure and organization of mitochondrial genomes may have contributed to the acquisition of mitochondrial ribosomal proteins with new functions.

Major conclusions

The distinct and diverse compositions of mitochondrial ribosomes illustrate the high evolutionary divergence found between mitochondrial genetic systems.

General significance

Elucidating the role of the organism-specific supernumerary proteins may provide a window into the regulation of mitochondrial gene expression through evolution in response to distinct evolutionary paths taken by mitochondria in different organisms. This article is part of a Special Issue entitled Frontiers of Mitochondrial Research.     

A short-oligonucleotide microarray that allows improved detection of gastrointestinal tract microbial communities

Background  

The human gastrointestinal (GI) tract contains a diverse collection of bacteria, most of which are unculturable by conventional microbiological methods. Increasingly molecular profiling techniques are being employed to examine this complex microbial community. The purpose of this study was to develop a microarray technique based on 16S ribosomal gene sequences for rapidly monitoring the microbial population of the GI tract.     

The RIN: an RNA integrity number for assigning integrity values to RNA measurements

Background  

The integrity of RNA molecules is of paramount importance for experiments that try to reflect the snapshot of gene expression at the moment of RNA extraction. Until recently, there has been no reliable standard for estimating the integrity of RNA samples and the ratio of 28S:18S ribosomal RNA, the common measure for this purpose, has been shown to be inconsistent. The advent of microcapillary electrophoretic RNA separation provides the basis for an automated high-throughput approach, in order to estimate the integrity of RNA samples in an unambiguous way.     

Duplicability of self-interacting human genes

Background  

There is increasing interest in the evolution of protein-protein interactions because this should ultimately be informative of the patterns of evolution of new protein functions within the cell. One model proposes that the evolution of new protein-protein interactions and protein complexes proceeds through the duplication of self-interacting genes. This model is supported by data from yeast. We examined the relationship between gene duplication and self-interaction in the human genome.     

Phylogenetic relationships and biogeography of the genus Algansea Girard (Cypriniformes: Cyprinidae) of central Mexico inferred from molecular data

Background  

The genus Algansea is one of the most representative freshwater fish groups in central Mexico due to its wide geographic distribution and unusual level of endemicity. Despite the small number of species, this genus has had an unsettled taxonomic history due to high levels of intraspecific morphological variation. Moreover, several phylogenetic hypotheses among congeners have been proposed but have had the following shortcomings: the use of homoplasious morphological characters, the use of character codification and polarisation methods that lacked objectivity, and incomplete taxonomic sampling. In this study, a phylogenetic analysis among species of Algansea is presented. This analysis is based upon two molecular markers, the mitochondrial gene cytochrome b and the first intron of the ribosomal protein S7 gene.     

Genome factor and gene pleiotropy hypotheses in protein evolution

   

The debate of genomic correlations between sequence conservation, protein connectivity, gene essentiality and gene expression, has generated a number of new hypotheses that are challenging the classical framework of molecular evolution. For instance, the translational selection hypothesis claims that the determination of the rate of protein evolution is the protein stability to avoid the misfolding toxicity. In this short article, we propose that gene pleiotropy, the capacity for affecting multiple phenotypes, may play a vital role in molecular evolution. We discuss several approaches to testing this hypothesis.     

Generalized linear mixed model for segregation distortion analysis

Background

Concerted evolution refers to the pattern in which copies of multigene families show high intraspecific sequence homogeneity but high interspecific sequence diversity. Sequence homogeneity of these copies depends on relative rates of mutation and recombination, including gene conversion and unequal crossing over, between misaligned copies. The internally repetitive intergenic spacer (IGS) is located between the genes for the 28S and 18S ribosomal RNAs. To identify patterns of recombination and/or homogenization within IGS repeat arrays, and to identify regions of the IGS that are under functional constraint, we analyzed 13 complete IGS sequences from 10 individuals representing four species in the Daphnia pulex complex.

Results

Gene conversion and unequal crossing over between misaligned IGS repeats generates variation in copy number between arrays, as has been observed in previous studies. Moreover, terminal repeats are rarely involved in these events. Despite the occurrence of recombination, orthologous repeats in different species are more similar to one another than are paralogous repeats within species that diverged less than 4 million years ago. Patterns consistent with concerted evolution of these repeats were observed between species that diverged 8-10 million years ago. Sequence homogeneity varies along the IGS; the most homogeneous regions are downstream of the 28S rRNA gene and in the region containing the core promoter. The inadvertent inclusion of interspecific hybrids in our analysis uncovered evidence of both inter- and intrachromosomal recombination in the nonrepetitive regions of the IGS.

Conclusions

Our analysis of variation in ribosomal IGS from Daphnia shows that levels of homogeneity within and between species result from the interaction between rates of recombination and selective constraint. Consequently, different regions of the IGS are on substantially different evolutionary trajectories.     

QPath: a method for querying pathways in a protein-protein interaction network

Background  

Sequence comparison is one of the most prominent tools in biological research, and is instrumental in studying gene function and evolution. The rapid development of high-throughput technologies for measuring protein interactions calls for extending this fundamental operation to the level of pathways in protein networks.     

Enzyme dimension of the ribosomal protein S4 across plant and animal kingdoms

Background

The protein S4 of the smaller ribosomal subunit is centrally important for its anchorage role in ribosome assembly and rRNA binding. Eubacterial S4 also facilitates synthesis of rRNA, and restrains translation of ribosomal proteins of the same polycistronic mRNA. Eukaryotic S4 has no homolog in eubacterial kingdom, nor are such extraribosomal functions of S4 known in plants and animals even as genetic evidence suggests that deficiency of S4X isoform in 46,XX human females may produce Turner syndrome (45,XO).

Methods

Recombinant human S4X and rice S4 were used to determine their enzymatic action in the cleavage of synthetic peptide substrates and natural proteins. We also studied autoproteolysis of the recombinant S4 proteins, and examined the growth and proliferation of S4-transfected human embryonic kidney cells.

Results

Extraribosomal enzyme nature of eukaryotic S4 is described. Both human S4X and rice S4 are cysteine proteases capable of hydrolyzing a wide spectrum of peptides and natural proteins of diverse origin. Whereas rice S4 also cleaves the -XXXD↓- consensus sequence assumed to be specific for caspase-9 and granzyme B, human S4 does not. Curiously, both human and rice S4 show multiple-site autoproteolysis leading to self-annihilation. Overexpression of human S4 blocks the growth and proliferation of transfected embryonic kidney cells, presumably due to the extraribosomal enzyme trait reported.

Conclusions

The S4 proteins of humans and rice, prototypes of eukaryota, are non-specific cysteine proteases in the extraribosomal milieu.

General significance

The enzyme nature of S4 is relevant toward understanding not only the origin of ribosomal proteins, but also processes in cell biology and diseases.     

Identification of a Novel G2073A Mutation in 23S rRNA in Amphenicol-Selected Mutants of Campylobacter jejuni

Objectives

This study was conducted to examine the development and molecular mechanisms of amphenicol resistance in Campylobacter jejuni by using in vitro selection with chloramphenicol and florfenicol. The impact of the resistance development on growth rates was also determined using in vitro culture.

Methods

Chloramphenicol and florfenicol were used as selection agents to perform in vitro stepwise selection. Mutants resistant to the selective agents were obtained from the selection process. The mutant strains were compared with the parent strain for changes in MICs and growth rates. The 23S rRNA gene and the L4 and L22 ribosomal protein genes in the mutant strains and the parent strain were amplified and sequenced to identify potential resistance-associated mutations.

Results

C. jejuni strains that were highly resistant to chloramphenicol and florfenicol were obtained from in vitro selection. A novel G2073A mutation in all three copies of the 23S rRNA gene was identified in all the resistant mutants examined, which showed resistance to both chloramphenicol and florfenicol. In addition, all the mutants selected by chloramphenicol also exhibited the G74D modification in ribosomal protein L4, which was previously shown to confer a low-level erythromycin resistance in Campylobacter species. The mutants selected by florfenicol did not have the G74D mutation in L4. Notably, the amphenicol-resistant mutants also exhibited reduced susceptibility to erythromycin, suggesting that the selection resulted in cross resistance to macrolides.

Conclusions

This study identifies a novel point mutation (G2073A) in 23S rRNA in amphenicol-selected mutants of C. jejuni. Development of amphenicol resistance in Campylobacter likely incurs a fitness cost as the mutant strains showed slower growth rates in antibiotic-free media.