Organisation of the S10, spc and alpha ribosomal protein gene clusters in prokaryotic genomes

Although it is well known that there is no long range colinearity in gene order in bacterial genomes, it is thought that there are several regions that are under strong structural constraints during evolution, in which gene order is extremely conserved. One such region is the str locus, containing the S10-spc-alpha operons. These operons contain genes coding for ribosomal proteins and for a number of housekeeping genes. We compared the organisation of these gene clusters in 111 sequenced prokaryotic genomes (99 bacterial and 12 archaeal genomes). We also compared the organisation to the phylogeny based on 16S ribosomal RNA gene sequences and the sequences of the ribosomal proteins L22, L16 and S14. Our data indicate that there is much variation in gene order and content in these gene clusters, both in bacterial as well as in archaeal genomes. Our data indicate that differential gene loss has occurred on multiple occasions during evolution. We also noted several discrepancies between phylogenetic trees based on 16S rRNA gene sequences and sequences of ribosomal proteins L16, L22 and S14, suggesting that horizontal gene transfer did play a significant role in the evolution of the S10-spc-alpha gene clusters.     

Evolution of paralogous genes: Reconstruction of genome rearrangements through comparison of multiple genomes within Staphylococcus aureus

Analysis of evolution of paralogous genes in a genome is central to our understanding of genome evolution. Comparison of closely related bacterial genomes, which has provided clues as to how genome sequences evolve under natural conditions, would help in such an analysis. With species Staphylococcus aureus, whole-genome sequences have been decoded for seven strains. We compared their DNA sequences to detect large genome polymorphisms and to deduce mechanisms of genome rearrangements that have formed each of them. We first compared strains N315 and Mu50, which make one of the most closely related strain pairs, at the single-nucleotide resolution to catalogue all the middle-sized (more than 10 bp) to large genome polymorphisms such as indels and substitutions. These polymorphisms include two paralogous gene sets, one in a tandem paralogue gene cluster for toxins in a genomic island and the other in a ribosomal RNA operon. We also focused on two other tandem paralogue gene clusters and type I restriction-modification (RM) genes on the genomic islands. Then we reconstructed rearrangement events responsible for these polymorphisms, in the paralogous genes and the others, with reference to the other five genomes. For the tandem paralogue gene clusters, we were able to infer sequences for homologous recombination generating the change in the repeat number. These sequences were conserved among the repeated paralogous units likely because of their functional importance. The sequence specificity (S) subunit of type I RM systems showed recombination, likely at the homology of a conserved region, between the two variable regions for sequence specificity. We also noticed novel alleles in the ribosomal RNA operons and suggested a role for illegitimate recombination in their formation. These results revealed importance of recombination involving long conserved sequence in the evolution of paralogous genes in the genome.     

Comparative study of human mitochondrial proteome reveals extensive protein subcellular relocalization after gene duplications

Background  

Gene and genome duplication is the principle creative force in evolution. Recently, protein subcellular relocalization, or neolocalization was proposed as one of the mechanisms responsible for the retention of duplicated genes. This hypothesis received support from the analysis of yeast genomes, but has not been tested thoroughly on animal genomes. In order to evaluate the importance of subcellular relocalizations for retention of duplicated genes in animal genomes, we systematically analyzed nuclear encoded mitochondrial proteins in the human genome by reconstructing phylogenies of mitochondrial multigene families.     

Intrastrain and interstrain genetic variation within a paralogous gene family in Chlamydia pneumoniae

Background

Chlamydia pneumoniae causes human respiratory diseases and has recently been associated with atherosclerosis. Analysis of the three recently published C. pneumoniae genomes has led to the identification of a new gene family (the Cpn 1054 family) that consists of 11 predicted genes and gene fragments. Each member encodes a polypeptide with a hydrophobic domain characteristic of proteins localized to the inclusion membrane.

Results

Comparative analysis of this gene family within the published genome sequences provided evidence that multiple levels of genetic variation are evident within this single collection of paralogous genes. Frameshift mutations are found that result in both truncated gene products and pseudogenes that vary among isolates. Several genes in this family contain polycytosine (polyC) tracts either upstream or within the terminal 5' end of the predicted coding sequence. The length of the polyC stretch varies between paralogous genes and within single genes in the three genomes. Sequence analysis of genomic DNA from a collection of 12 C. pneumoniae clinical isolates was used to determine the extent of the variation in the Cpn 1054 gene family.

Conclusions

These studies demonstrate that sequence variability is present both among strains and within strains at several of the loci. In particular, changes in the length of the polyC tract associated with the different Cpn 1054 gene family members are common within each tested C. pneumoniae isolate. The variability identified within this newly described gene family may modulate either phase or antigenic variation and subsequent physiologic diversity within a C. pneumoniae population.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and archaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacteria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variations in some operons across different organisms within each domain, and these variations are informative on the evolutionary relations among the organisms. This method provides a new potential for studying the origin and evolution of old species.     

After the duplication: gene loss and adaptation in Saccharomyces genomes

The ancient duplication of the Saccharomyces cerevisiae genome and subsequent massive loss of duplicated genes is apparent when it is compared to the genomes of related species that diverged before the duplication event. To learn more about the evolutionary effects of the duplication event, we compared the S. cerevisiae genome to other Saccharomyces genomes. We demonstrate that the whole genome duplication occurred before S. castellii diverged from S. cerevisiae. In addition to more accurately dating the duplication event, this finding allowed us to study the effects of the duplication on two separate lineages. Analyses of the duplication regions of the genomes indicate that most of the duplicated genes (approximately 85%) were lost before the speciation. Only a small amount of paralogous gene loss (4-6%) occurred after speciation. On the other hand, S. castellii appears to have lost several hundred genes that were not retained as duplicated paralogs. These losses could be related to genomic rearrangements that reduced the number of chromosomes from 16 to 9. In addition to S. castellii, other Saccharomyces sensu lato species likely diverged from S. cerevisiae after the duplication. A thorough analysis of these species will likely reveal other important outcomes of the whole genome duplication.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and ar-chaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacte-ria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variat     

Simultaneous horizontal gene transfer of a gene coding for ribosomal protein l27 and operational genes in Arthrobacter sp

Phylogenetic analysis of bacterial L27 ribosomal proteins showed that, against taxonomy, the L27 protein from the Actinobacteria Arthrobacter sp. clusters with protein sequences from the Bacillus group. The L27 gene clusters in the Arthrobacter sp. genome with six genes responsible for creatinine and sarcosine degradation. Phylogenetic analyses of orthologue proteins encoded by three of these genes also showed a phylogenetic relationship with Bacillus species. Comparisons between the synonymous codon usage of the Arthrobacter sp. genes and those from complete genomes showed that Arthrobacter genes encoding the L27 ribosomal protein and the proteins responsible for the degradation of creatinine and sarcosine have a codon usage that is more similar to that of Bacillus species than that of Arthrobacter. We suggest that the Arthrobacter sp. genes encoding the L27 ribosomal protein and the proteins responsible for the degradation of creatinine and sarcosine were acquired simultaneously through horizontal gene transfer from an unknown Bacillus species.     

Genomic, regulatory and epigenetic mechanisms underlying duplicated gene evolution in the natural allotetraploid Oryza minuta

Background

Polyploid species contribute to Oryza diversity. However, the mechanisms underlying gene and genome evolution in Oryza polyploids remain largely unknown. The allotetraploid Oryza minuta, which is estimated to have formed less than one million years ago, along with its putative diploid progenitors (O. punctata and O. officinalis), are quite suitable for the study of polyploid genome evolution using a comparative genomics approach.

Results

Here, we performed a comparative study of a large genomic region surrounding the Shattering4 locus in O. minuta, as well as in O. punctata and O. officinalis. Duplicated genomes in O. minuta have maintained the diploid genome organization, except for several structural variations mediated by transposon movement. Tandem duplicated gene clusters are prevalent in the Sh4 region, and segmental duplication followed by random deletion is illustrated to explain the gene gain-and-loss process. Both copies of most duplicated genes still persist in O. minuta. Molecular evolution analysis suggested that these duplicated genes are equally evolved and mostly manipulated by purifying selection. However, cDNA-SSCP analysis revealed that the expression patterns were dramatically altered between duplicated genes: nine of 29 duplicated genes exhibited expression divergence in O. minuta. We further detected one gene silencing event that was attributed to gene structural variation, but most gene silencing could not be related to sequence changes. We identified one case in which DNA methylation differences within promoter regions that were associated with the insertion of one hAT element were probably responsible for gene silencing, suggesting a potential epigenetic gene silencing pathway triggered by TE movement.

Conclusions

Our study revealed both genetic and epigenetic mechanisms involved in duplicated gene silencing in the allotetraploid O. minuta.     

Evolution by leaps: gene duplication in bacteria

Background  

Sequence related families of genes and proteins are common in bacterial genomes. In Escherichia coli they constitute over half of the genome. The presence of families and superfamilies of proteins suggest a history of gene duplication and divergence during evolution. Genome encoded protein families, their size and functional composition, reflect metabolic potentials of the organisms they are found in. Comparing protein families of different organisms give insight into functional differences and similarities.     

Diversity of 5S rRNA genes within individual prokaryotic genomes

We examined intragenomic variation of paralogous 5S rRNA genes to evaluate the concept of ribosomal constraints. In a dataset containing 1161 genomes from 779 unique species, 96 species exhibited >?3% diversity. Twenty-seven species with >?10% diversity contained a total of 421 mismatches between all pairs of the most dissimilar copies of 5S rRNA genes. The large majority (401 of 421) of the diversified positions were conserved at the secondary structure level. The high diversity was associated with partial rRNA operon, split operon, or spacer length-related divergence. In total, these findings indicated that there are tight ribosomal constraints on paralogous 5S rRNA genes in a genome despite of the high degree of diversity at the primary structure level.     

The Mammalian 2′-5′ Oligoadenylate Synthetase Gene Family: Evidence for Concerted Evolution of Paralogous Oas1 Genes in Rodentia and Artiodactyla

Multiple 2′-5′ oligoadenylate (2-5A) synthetases are important components of innate immunity in mammals. Gene families encoding these proteins have previously been studied mainly in humans and mice. To reconstruct the evolution of this gene family in mammals, a search for additional 2-5A synthetase genes was performed in rat, cattle, pig, and dog. Twelve 2′-5′ oligoadenylate synthetase (Oas) genes were identified in the rat genome, including eight Oas1 genes, two Oas1 pseudogenes, single copies of Oas2 and Oas3, and two Oas-like genes, Oasl1 and Oasl2. Four OAS genes were detected in the pig genome and five OAS genes were found in both the cattle and dog genomes. An OAS3 gene was not found in either the cattle or the pig genome. While two tandemly duplicated OAS-like (OASL) genes were identified in the dog genome, only a single OASL orthologue was found in both the cattle and the pig genomes. The bovine and porcine OASL genes contain premature stop codons and encode truncated proteins, which lack the typical C-terminal double ubiquitin domains. The cDNA sequences of the rat, cattle, pig, and dog OAS genes were amplified, sequenced and compared with each other and with those in the human, mouse, horse, and chicken genomes. Evidence of concerted evolution of paralogous 2′-5′ oligoadenylate synthetase 1 genes was obtained in rodents (Rodentia) and even-toed ungulates (Artiodactyla). Calculations using the nonparametric Kolmogorov-Smirnov test suggested that the homogenization of paralogous OAS1 sequences was due to gene conversion rather than stabilizing selection. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. Martin Kreitman     

Ancient gene duplications and the root(s) of the tree of life

Summary. Tracing organismal histories on the timescale of the tree of life remains one of the challenging tasks in evolutionary biology. The hotly debated questions include the evolutionary relationship between the three domains of life (e.g., which of the three domains are sister domains, are the domains para-, poly-, or monophyletic) and the location of the root within the universal tree of life. For the latter, many different points of view have been considered but so far no consensus has been reached. The only widely accepted rationale to root the universal tree of life is to use anciently duplicated paralogous genes that are present in all three domains of life. To date only few anciently duplicated gene families useful for phylogenetic reconstruction have been identified. Here we present results from a systematic search for ancient gene duplications using twelve representative, completely sequenced, archaeal and bacterial genomes. Phylogenetic analyses of identified cases show that the majority of datasets support a root between Archaea and Bacteria; however, some datasets support alternative hypotheses, and all of them suffer from a lack of strong phylogenetic signal. The results are discussed with respect to the impact of horizontal gene transfer on the ability to reconstruct organismal evolution. The exchange of genetic information between divergent organisms gives rise to mosaic genomes, where different genes in a genome have different histories. Simulations show that even low rates of horizontal gene transfer dramatically complicate the reconstruction of organismal evolution, and that the different most recent common molecular ancestors likely existed at different times and in different lineages. Correspondence and reprints: Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT 06269-3125, U.S.A. Present address: Genome Atlantic, Department of Biochemistry and Molecular Biology, Dalhousie University, Halifax, Nova Scotia, Canada.     

MetaHox gene clusters

Homeobox genes encode important developmental control proteins. The Drosophila fruit fly HOM complex genes are clustered in region 84-89 of chromosome 3. Probably due to large-scale genome duplication events, their human HOX orthologs belong to four paralogous regions. A series of 13 other homeobox genes are also clustered in region 88-94, on the same chromosome of Drosophila. We suggest that they also duplicated during vertebrate evolution and belong to paralogous regions in humans. These regions are on chromosome arms 4p, 5q, 10q, and 2p or 8p. We coined the term "paralogon" to designate paralogous regions in general. We propose to call these genes "meta Hox" genes. Like Hox genes, metaHox genes are present in one cluster in Drosophila and four clusters (metaHox A-D) in humans on the 4p/5q/10q paralogon.     

Chromosomal mapping,differential origin and evolution of the S100 gene family

S100 proteins are calcium-binding proteins, which exist only in vertebrates and which constitute a large protein family. The origin and evolution of the S100 family in vertebrate lineages remain a challenge. Here, we examined the synteny conservation of mammalian S100A genes by analysing the sequence of available vertebrate S100 genes in databases. Five S100A gene members, unknown previously, were identified by chromosome mapping analysis. Mammalian S100A genes are duplicated and clustered on a single chromosome while two S100A gene clusters are found on separate chromosomes in teleost fish, suggesting that S100A genes existed in fish before the fish-specific genome duplication took place. During speciation, tandem gene duplication events within the cluster of S100A genes of a given chromosome have probably led to the multiple members of the S100A gene family. These duplicated genes have been retained in the genome either by neofunctionalisation and/or subfunctionalisation or have evolved into non-coding sequences. However in vertebrate genomes, other S100 genes are also present i.e. S100P, S100B, S100G and S100Z, which exist as single copy genes distributed on different chromosomes, suggesting that they could have evolved from an ancestor different to that of the S100A genes.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.     

Gene duplication and the origins of morphological complexity in pancrustacean eyes,a genomic approach

Background  

Duplication and divergence of genes and genetic networks is hypothesized to be a major driver of the evolution of complexity and novel features. Here, we examine the history of genes and genetic networks in the context of eye evolution by using new approaches to understand patterns of gene duplication during the evolution of metazoan genomes. We hypothesize that 1) genes involved in eye development and phototransduction have duplicated and are retained at higher rates in animal clades that possess more distinct types of optical design; and 2) genes with functional relationships were duplicated and lost together, thereby preserving genetic networks. To test these hypotheses, we examine the rates and patterns of gene duplication and loss evident in 19 metazoan genomes, including that of Daphnia pulex - the first completely sequenced crustacean genome. This is of particular interest because the pancrustaceans (hexapods+crustaceans) have more optical designs than any other major clade of animals, allowing us to test specifically whether the high amount of disparity in pancrustacean eyes is correlated with a higher rate of duplication and retention of vision genes.