Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Paralogous Genes”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33Linked to original sources

Concerted and nonconcerted evolution of the Hsp70 gene superfamily in two sibling species of nematodes.

We have identified the Hsp70 gene superfamily of the nematode Caenorhabditis briggsae and investigated the evolution of these genes in comparison with Hsp70 genes from C. elegans, Drosophila, and yeast. The Hsp70 genes are classified into three monophyletic groups according to their subcellular localization, namely, cytoplasm (CYT), endoplasmic reticulum (ER), and mitochondria (MT). The Hsp110 genes can be classified into the polyphyletic CYT group and the monophyletic ER group. The different Hsp70 and Hsp110 groups appeared to evolve following the model of divergent evolution. This model can also explain the evolution of the ER and MT genes. On the other hand, the CYT genes are divided into heat-inducible and constitutively expressed genes. The constitutively expressed genes have evolved more or less following the birth-and-death process, and the rates of gene birth and gene death are different between the two nematode species. By contrast, some heat-inducible genes show an intraspecies phylogenetic clustering. This suggests that they are subject to sequence homogenization resulting from gene conversion-like events. In addition, the heat-inducible genes show high levels of sequence conservation in both intra-species and inter-species comparisons, and in most cases, amino acid sequence similarity is higher than nucleotide sequence similarity. This indicates that purifying selection also plays an important role in maintaining high sequence similarity among paralogous Hsp70 genes. Therefore, we suggest that the CYT heat-inducible genes have been subjected to a combination of purifying selection, birth-and-death process, and gene conversion-like events.

Animals↗

Deletions of 2q14 that include the homeobox engrailed 1 (EN1) transcription factor are compatible with a normal phenotype.

A novel transmitted 2-3 Mb deletion of 2q14.1-q14.2 was found in an affected boy from a consanguineous family with a possible diagnosis of PEHO syndrome (OMIM 260565). BAC FISH showed that the deletion included a minimum of 20 genes including the homeobox engrailed 1 gene (EN1). However, the same deletion was also found in his phenotypically normal father and brother (family 1). The phenotype of the proband may, therefore, have been coincidental to the deletion, a result of a recessive condition within or outside the deleted segment or possibly due to variable dosage compensation of EN1 by the paralogous EN2 gene at 7q36. BAC FISH also showed that this deletion overlapped with a previously reported transmitted deletion of 2q13-q14.1 that had no phenotypic consequences (family 2). The deleted regions contained a total of 32 genes and comprise the final 5.25 Mb of the ancestral chromosome 2B from which chromosome 2 was formed in man. These families provide further evidence that heterozygous deletions of regions of low gene density are compatible with a normal phenotype.

Abnormalities, Multiple↗

Region-specific expression of murine Hox genes implies the Hox code-mediated patterning of the digestive tract.

BACKGROUND: Hox genes encode transcription factors which are involved in the establishment of regional identities along the anteroposterior (AP) body axis. To elucidate the AP patterning of the digestive tract, we have systematically examined the expression patterns of Hox genes belonging to paralogue groups 6, 7, 8 and 9 by whole-mount in situ hybridization and by section in situ hybridization analyses. RESULTS: The expression patterns of these genes showed co-linearity along the wall of the digestive tract, thereby yielding the Hox code of the gut. The expression boundaries of the Hox genes at later stages (12.5 d.p.c.) corresponded to the morphological boundaries of individual gut subdomains. CONCLUSIONS: The visceral mesoderm-restricted expression suggested that the Hox code primarily functions in the mesenchymal specification which eventually leads to the regional differentiation of gut subdomains as the result of epithelial-mesenchymal interactions. Overlapping expression patterns were found among the paralogous Hox genes, indicating that the paralogues may have redundant functions in the specification of the gut.

Animals↗

A new osteochondrodysplasia with severe osteopenia, preaxial polydactyly, clefting and dysmorphic features resembling filamin-related disorders.

BACKGROUND: We report a 19-week gestation female foetus with a new syndrome characterised by increased nuchal translucency and severe micromelia with campomelia evident from the early second trimester. METHODS AND RESULTS: Cytogenetic studies performed on amniocytes revealed a normal female karyotype. Autopsy after termination of pregnancy showed facial dysmorphism, cleft palate, bowed, shortened limbs, hypoplasia of the preaxial elements in all four limbs with accompanying accessory ossification centres in the feet, and severe calvarial underossification. A diagnosis of otopalatodigital syndrome type 2, associated with mutations in FLNA, a gene encoding the cytoskeletal protein filamin A, was considered but discarded due to the severity of micromelia, early lethality, and the presence of generalised osteopenia instead of hyperostosis. The degree of undermodelling and campomelia was reminiscent of another group of conditions that include atelosteogenesis types 1 and 3, caused by mutations in FLNB. Sequencing analysis did not reveal any pathogenic mutation in the three paralogous filamin genes: FLNA, FLNB and FLNC. CONCLUSION: Clinical, radiological and cytogenetic findings suggest that this phenotype is a new entity whose aetiopathogenesis may be functionally related to the filaminopathies.

Abnormalities, Multiple↗

Coordinated expression of 3' hox genes during murine embryonal gut development: an enteric Hox code.

BACKGROUND & AIMS: Hox genes are highly conserved developmental control genes that may be organized and expressed in the form of a code required for correct morphogenesis. Little is known about their control of the embryonal gut. However, Hox paralogues 4 and 5, which are expressed at the sites of origin of vagal neural crest cells and splanchnic mesoderm, are likely to be important. We have studied the expression domains of these genes in the gut both spatially and temporally. METHODS: CD1 mice embryos of embryonic days E8.5-E17.5 were studied. The spatial and temporal expression patterns of messenger RNA of Hoxa4, b4, c4, d4, a5, c5, and b5 homeoprotein were determined by in situ hybridization and immunohistochemistry in whole embryos, whole gastrointestinal tracts, and vibratome sections. RESULTS: There were different spatial, temporal, and combinatorial expression patterns in different morphological regions: foregut, prececal gut, cecum, and postcecal gut. Two dynamic gradients, rostral and caudal, were coordinated with nested expression domains along the gut primordium. Region-specific domains were present in the stomach and cecum. CONCLUSIONS: The expression patterns of genes in paralogous groups 4 and 5 suggest that they are organized to form a specific enteric Hox code required for correct enteric development.

Animals↗

Paralog-aware assembly and filtering strategies reveal minimal nucleotide variation on the macro germline-restricted chromosome of the zebra finch.

The germline-restricted chromosome (GRC) of passerines is a remarkable tissue-specific chromosome that accumulated paralogs of genes from the regular "A chromosomes" over millions of years, often amplified into dozens of gene copies. In addition to its repetitive content, typically uniparental inheritance, and lack of recombination, the GRC resembles non-recombining sex chromosomes and some B chromosomes, for all of which assembly and single-nucleotide polymorphisms (SNPs) calling are difficult. Here, we first show that much of the Australian zebra finch macro-GRC can be assembled using accurate long reads. We then describe a paralog-aware Snakemake pipeline, ParaVar, to map short reads from the GRC to retrieve GRC regions suitable for haplotype-based analysis. ParaVar reliably calls hundreds of SNPs across the GRC, thereby providing an estimate of nucleotide diversity on the highly repetitive zebra finch macro-GRC. Our results show significantly lower nucleotide diversity (20- to 50-fold lower) on the GRC compared to the mitogenome and autosomes, and a strong phylogenetic discordance between the GRC and the mitochondrial genome. Beyond the contribution of background selection, our results suggest that a single GRC haplotype recently spread through the populations while jumping across matrilines via occasional paternal inheritance. We anticipate that our paralog-aware pipeline will be useful for SNP calling and population genetics analyses of repetitive GRCs, sex chromosomes, and B chromosomes.

Animals↗

Double-stranded RNA-mediated silencing of genomic tandem repeats and transposable elements in the D. melanogaster germline.

BACKGROUND: The injection of double-stranded RNA (dsRNA) has been shown to induce a potent sequence-specific inhibition of gene function in diverse invertebrate and vertebrate species. The homology-dependent posttranscriptional gene silencing (PTGS) caused by the introduction of transgenes in plants may be mediated by dsRNA. The analysis of Caenorhabditis elegans mutants impaired with dsRNA-mediated silencing and studies in plants implicate a biological role of dsRNA-mediated silencing as a transposon-repression and antiviral mechanism. RESULTS: We investigated the silencing of testis-expressed Stellate genes by paralogous Su(Ste) tandem repeats, which are known to be involved in the maintenance of male fertility in Drosophila melanogaster. We found that both strands of repressor Su(Ste) repeats are transcribed, producing sense and antisense RNA. The Stellate silencing is associated with the presence of short Su(Ste) RNAs. Cotransfection experiments revealed that Su(Ste) dsRNA can target and eliminate Stellate transcripts in Drosophila cell culture. The short fragment of Stellate gene that is homologous to Su(Ste) was shown to be sufficient to confer Su(Ste)-dependent silencing of a reporter construct in testes. We demonstrated that Su(Ste) dsRNA-mediated silencing affects not only Stellate expression but also the level of sense Su(Ste) RNA providing a negative autogenous regulation of Su(Ste) expression. Mutation in the spindle-E gene relieving Stellate silencing also leads to a derepression of the other genomic tandem repeats and retrotransposons in the germline. CONCLUSIONS: Homology-dependent gene silencing was shown to be used to inhibit Stellate gene expression in the D. melanogaster germline, ensuring male fertility. dsRNA-mediated silencing may provide a basis for negative autogenous control of gene expression. The related surveillance system is implicated to control expression of retrotransposons in the germline.

Adenosine Triphosphatases↗

Dcs2, a novel stress-induced modulator of m7GpppX pyrophosphatase activity that locates to P bodies.

The eukaryotic "scavenger" type decapping enzyme, an m(7)GpppX pyrophosphatase, is active in cellular mRNA metabolism and thereby influences posttranscriptional gene expression. The yeast version of this enzyme, Dcs1, catalyses cleavage of 5'end m(7)G-oligoribonucleotide fragments generated by 3'-->5' exonucleolytic decay, and cleavage of m(7)GDP generated by Dcp1/Dcp2-mediated decapping in the 5'-->3' decay pathway. We show that Dcs1 is active as a homodimer with low KM values for cleavage of m(7)GpppG (0.14 microM) and m(7)GDP (0.26 microM). Previous work showed that the paralogous DCS2 gene is transcriptionally induced via the amp-PKA pathway as yeast enters diauxie. The resulting Dcs2 protein forms a heterodimer together with Dcs1, both modulating Dcs1 substrate specificity and suppressing its k(cat). Since Dcs2 is recruited into cytoplasmic P bodies, its inhibitory function may be focused in these centres of mRNA storage/turnover. Dcs2 is therefore a novel type of stress-induced regulatory protein that modulates m(7)GpppX pyrophosphatase activity. Moreover, inhibition of Dcs1 activity by Dcs2, like depletion of Dcs1, reduces chronological life span, possibly by modulating m(7)G misincorporation into nucleic acids. This could potentially link control of mRNA metabolism with senescence.

Adaptor Proteins, Signal Transducing↗

Caenorhabditis elegans embryonic axial patterning requires two recently discovered posterior-group Hox genes.

Hox genes encode highly conserved transcription factors that control regional identities of cells and tissues along the developing anterior-posterior axis, probably in all bilaterian metazoans. However, in invertebrate embryos other than Drosophila, Hox gene functions remain largely unknown except by inference from sequence similarities and expression patterns. Recent genomic sequencing has shown that Caenorhabditis elegans has three Hox genes of the posterior paralog group [Ruvkun, G. & Hobert, O. (1998) Science 282, 2033-2041]. However, only one has been previously identified genetically, and it is not required for embryonic development [Chisholm, A. (1991) Development (Cambridge, U.K.) 111, 921-932]. Herein, we report identification of the remaining two posterior paralogs as the nob-1 gene and the neighboring php-3 gene. Elimination of nob-1 and php-3 functions causes gross embryonic defects in both posterior patterning and morphogenetic movements of the posterior hypodermis, as well as posterior-to-anterior cell fate transformations and lethality. The only other Hox gene essential for embryogenesis is the labial/Hox1 homolog ceh-13, required for more anterior patterning [Brunschwig, K., Wittmann, C., Schnabel, R., Burglin, T. R., Tobler, H. & Muller, F. (1999) Development (Cambridge, U.K.) 126, 1537-1546]. Therefore, essential embryonic patterning in C. elegans requires only Hox genes of the anterior and posterior paralog groups, raising interesting questions about evolution of the medial-group genes.

Amino Acid Sequence↗

N-linked glycosylation and sequence changes in a critical negative control region of the ASCT1 and ASCT2 neutral amino acid transporters determine their retroviral receptor functions.

A widely dispersed interference group of retroviruses that includes the feline endogenous virus (RD114), baboon endogenous virus (BaEV), human endogenous virus type W (HERV-W), and type D primate retroviruses uses the human Na(+)-dependent neutral amino acid transporter type 2 (hASCT2; gene name, SLC1A5) as a common cell surface receptor. Although hamster cells are fully resistant to these viruses and murine cells are susceptible only to BaEV and HERV-W pseudotype viruses, these rodent cells both become highly susceptible to all of the viruses after treatment with tunicamycin, an inhibitor of protein N-linked glycosylation. A partial explanation for these results was recently provided by findings that the orthologous murine transporter mASCT2 is inactive as a viral receptor, that a related (ca. 55% identity) murine paralog (mASCT1; gene name, SLC1A4) mediates infections specifically of BaEV and HERV-W, and that N-deglycosylation of mASCT1 activates it as a receptor for all viruses of this interference group. Because the only two N-linked oligosaccharides in mASCT1 occur in the carboxyl-terminal region of extracellular loop 2 (ECL2), it was inferred that this region contributes in an inhibitory manner to infections by RD114 and type D primate viruses. To directly and more thoroughly investigate the receptor active sites, we constructed and analyzed a series of hASCT2/mASCT2 chimeras and site-directed mutants. Our results suggest that a hypervariable sequence of 21 amino acids in the carboxyl-terminal portion of ECL2 plays a critical role in determining the receptor properties of ASCT2 proteins for all viruses in this interference group. In addition, we analyzed the tunicamycin-dependent viral susceptibility of hamster cells. In contrast to mASCT1, which contains two N-linked oligosaccharides that partially restrict viral infections, hamster ASCT1 contains an additional N-linked oligosaccharide clustered close to the others in the carboxyl-terminal region of ECL2. Removal of this N-linked oligosaccharide by mutagenesis enabled hamster ASCT1 to function as a receptor for all viruses of this interference group. These results strongly suggest that combinations of amino acid sequence changes and N-linked oligosaccharides in a critical carboxyl-terminal region of ECL2 control retroviral utilization of both the ASCT1 and ASCT2 receptors.

Amino Acid Sequence↗

Ancestral paralogs and pseudoparalogs and their role in the emergence of the eukaryotic cell.

Gene duplication is a crucial mechanism of evolutionary innovation. A substantial fraction of eukaryotic genomes consists of paralogous gene families. We assess the extent of ancestral paralogy, which dates back to the last common ancestor of all eukaryotes, and examine the origins of the ancestral paralogs and their potential roles in the emergence of the eukaryotic cell complexity. A parsimonious reconstruction of ancestral gene repertoires shows that 4137 orthologous gene sets in the last eukaryotic common ancestor (LECA) map back to 2150 orthologous sets in the hypothetical first eukaryotic common ancestor (FECA) [paralogy quotient (PQ) of 1.92]. Analogous reconstructions show significantly lower levels of paralogy in prokaryotes, 1.19 for archaea and 1.25 for bacteria. The only functional class of eukaryotic proteins with a significant excess of paralogous clusters over the mean includes molecular chaperones and proteins with related functions. Almost all genes in this category underwent multiple duplications during early eukaryotic evolution. In structural terms, the most prominent sets of paralogs are superstructure-forming proteins with repetitive domains, such as WD-40 and TPR. In addition to the true ancestral paralogs which evolved via duplication at the onset of eukaryotic evolution, numerous pseudoparalogs were detected, i.e. homologous genes that apparently were acquired by early eukaryotes via different routes, including horizontal gene transfer (HGT) from diverse bacteria. The results of this study demonstrate a major increase in the level of gene paralogy as a hallmark of the early evolution of eukaryotes.

Amino Acid Sequence↗

Intermedin, a novel calcitonin family peptide that exists in teleosts as well as in mammals: a comparison with other calcitonin/intermedin family peptides in vertebrates.

Endocrine regulation in vertebrates is critical for the adaptation and regulation of homeostasis. The G protein-coupled receptor (GPCR) signaling transduction system represents one of the most ancient forms of cell surface signaling. Recently, comparative sequence analysis has aided in the identification and pairing of a variety of ligand/GPCR signaling systems. Among the ligands of type II GPCRs, the calcitonin family peptides including calcitonin, alpha-calcitonin gene-related peptide (alphaCGRP), betaCGRP, adrenomedullin, and amylin are among the best studied hormones, and the founding member, calcitonin, was originally identified and isolated from teleosts. This unique group of peptides shares a conserved tertiary structure with an N-terminal disulfide-bridged ring. In mammals, these peptides signal through two closely related type II GPCRs and three unique receptor activity-modifying proteins. Recently, based on the analysis of multiple vertebrate genomes, we identified a novel calcitonin/CGRP family peptide named intermedin. Here we show that in humans the five paralogous family genes, calcitonin, CGRP, amylin, adrenomedullin, and intermedin, evolved before the emergence of modern vertebrates, and that teleost genomes carry multiple copies of these co-evolved hormone genes. Sequence comparison showed that each of these genes is highly conserved in different vertebrates and that multiple copies of these peptides in teleosts could be derived from ancient genome duplication and/or lineage-specific intragenic duplications. The present article provides an overview of the calcitonin/intermedin family peptides found in teleost and mammalian genomes, and describes their putative functions. In addition, we demonstrate that one of the intermedin orthologs deduced from the pufferfish (Fugu rubripes) genome shares a conserved signaling activity with mammalian intermedin. The combined results indicate that the physiology associated with each of these family peptides likely evolved during early vertebrate evolution and diverged to serve select physiological functions in different vertebrates.

Adrenomedullin↗

Ectopic expression of Hox-2.3 induces craniofacial and skeletal malformations in transgenic mice.

To better understand the role of the Hox-2.3 murine homeobox gene during development, a dominant gain-of-function mutation was generated. The developmental malformations that resulted when the chicken beta-actin promoter was used to direct widespread expression of the Hox-2.3 gene in transgenic mice included early postnatal death as well as craniofacial abnormalities, including open eyes and cleft palate. Ventricular septal defects were also observed in the hearts of three transgenic mice. Skeletal malformations were seen in the bones of the craniocervical transition, with the occipital, basisphenoid, and atlas bones deficient or misshapen. Interestingly, one mutant exhibited an extra pair of ribs as well as alterations in cervical vertebrae identities. Some of the malformations observed in Hox-2.3 gain-of-function mutants overlap with those seen in Hox-1.1 and Hox-2.2 misexpression mutants which suggests functional similarities between paralogous homeobox genes. The results of these experiments are consistent with a role for Hox-2.3 in specifying positional information during development.

Abnormalities, Multiple↗

Tdrd15 is dispensable for male fertility and spermatogenesis in the golden hamster.

Tudor domain-containing proteins (TDRDs) constitute an evolutionarily conserved protein family and are critical for germline development and piRNA pathway regulation, with established roles in male fertility. While multiple TDRD family members have been functionally linked to spermatogenic impairment, the precise biological role of TDRD15 remains to be elucidated. We used CRISPR/Cas9-mediated gene editing to generate Tdrd15 knockout (KO) golden hamsters (Mesocricetus auratus), a model necessitated by the absence of a functional Tdrd15 ortholog in the mouse genome, to investigate its function in male reproduction. Phylogenetic analysis demonstrated that TDRD15 is strongly conserved among eutherian mammals, with testis-restricted expression patterns in hamsters. Despite the successful induction of frameshift mutations and significant transcriptional knockdown, Tdrd15 KO males maintained normal fertility parameters, including unaltered testicular architecture, spermatogenic progression (confirmed by periodic acidic-Schiff (PAS) staining and immunohistochemistry), and sperm quality metrics determined using a computer-assisted analysis. Quantitative polymerase chain reaction (qPCR) revealed compensatory overexpression of paralogous Tdrd genes in KO testes, implying functional redundancy within this protein family. This study provides the first experimental evidence that TDRD15 is dispensable for male fertility in golden hamsters under physiological conditions, thereby challenging the prevailing assumptions of its obligatory function in spermiogenesis. Altogether, these findings support a more targeted allocation of research efforts within the field of male reproductive biology.

Animals↗

Region-specific gastrointestinal Hox code during murine embryonal gut development.

Hox genes encode transcription factors, and they are involved in the specification of each body part along the anteroposterior (AP) body axis during embryogenesis. To clarify AP pattern formation of the digestive tract, the expression patterns of Hox genes belonging to paralogous groups 4 and 5, and parts of groups 6 and 7, were systematically examined by whole-mount and section in situ hybridization. The Hox gene expression pattern of paralogous groups 4-9 in the developing gut at 12.5 days post-coitum was fully examined. All HoxA and HoxB genes in paralogous groups 4-8 were expressed in the stomach, in contrast to the HoxC and HoxD genes. In the midgut region, all Hox cluster genes showed colinear expression within each cluster, yielding the Hox code; the more 3' located genes were expressed more rostrally and the 5' group genes more caudally. The colinear expression of HoxA and HoxB cluster genes started from the duodenum, that of HoxC cluster genes started from the jejunum, and HoxD cluster genes were expressed in the caudal part of the midgut, ileum and cecum. In the hindgut region, HoxD cluster genes and Abd-B family genes were expressed. Thus, a different Hox code seems to exist in each subdomain of developing gut (foregut, midgut and hindgut). The visceral mesoderm restricted expression also suggested that the Hox code primarily functions in mesenchymal specification, and then leads to the regional differentiation of gut subdomains as the result of epithelial-mesenchymal interactions.

Animals↗

Classification of sequences expressed during the primordial and basidiome stages of the cultivated mushroom Agaricus bisporus.

Two complementary DNA (cDNA) libraries were constructed from tissues isolated from primordia and basidiomes of Agaricus bisporus to characterize genes involved in mushroom development. Using single-pass sequencing of 869 cDNA clones, we found 477 expressed sequence tags (ESTs), including 466 not previously described in the databases for A. bisporus. A BLASTX search revealed that 374 ESTs had similarities with protein sequences available from databases; 193 of these ESTs were categorized according to their putative function. Most ESTs were assigned to one of four roles: metabolism (23%), cell structure (15%), cell growth and division (12%), and protein destination and storage (10%). The remaining ESTs with putative homologues were classified in 10 additional categories. Many ESTs could not be functionally assigned. Based on redundancy levels, at least 4 ESTs were preferentially expressed in each tissue type. Sequence analysis also suggested the presence of paralog tyrosinase genes in the A. bisporus genome.

Agaricus↗

Targeted mutations in hoxa-9 and hoxb-9 reveal synergistic interactions.

Mice were generated with a targeted disruption of the homeobox-containing gene hoxb-9. Mice homozygous for this mutation show defects in the development of the first and second ribs. In most cases the first and second ribs are fused near the point at which the first and second pairs of ribs normally attach to the sternum. Abnormalities of the sternum accompany the rib fusions. These include abnormal attachment of the ribs to the sternum, a reduction in the number of intercostal segments of the sternum, and abnormal growth of the intercostal segments. Over half of the homozygous mutants, as well as some heterozygotes, also have an eighth rib attached to the sternum. These results show that hoxb-9 plays a significant role in the specification of thoracic skeletal elements. To reveal potential interactions between the paralogous Hox genes hoxa-9 and hoxb-9, mice heterozygous for both mutations were intercrossed. Mice homozygous for both mutations show more severe phenotypes than predicted by the addition of the individual mutant phenotypes. Both the penetrance and the expressivity of the rib and sternal defects are increased, suggesting synergistic interactions between these genes. In particular, the sternum defects are greatly exacerbated. Interestingly, the defects in hoxb-9 and hoxa-9/ hoxb-9 mutant mice are concentrated along the axial column at points of transition between vertebral types.

Animals↗

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.

Bacterial Proteins↗