Evolution of gene orders in genomes of cyanobacteria

Genomes of 23 strains of cyanobacteria were comparatively analyzed using quantitative methods of estimation of gene order similarity. It has been found that reconstructions of phylogenesis of cyanobacteria based on the comparison of the orders of genes in chromosomes and nucleotide sequences appear to be similar. This confirms the applicability of quantitative measures of similarity of gene orders for phylogenetic reconstructions. In the evolution of marine unicellular planktonic cyanobacteria, genome rearrangements are fixed with a low rate (about 3% of gene order changes per 1% of 16S rRNA changes), whereas in other groups of cyanobacteria the gene order can change several times more rapidly. The gene orders in genomes of cyanobacteria and chloroplasts preserve a considerable degree of similarity. The closest relatives of chloroplasts among the analyzed cyanobacteria are likely to be strains from hot springs belonging to the genus Synechococcus. Comparative analysis of gene orders and nucleotide sequences strongly suggests that Synechococcus strains from different environments (sea, fresh waters, hot springs) are not related and belong to evolutionally distant lines.     

Gene orders in genomes of alpha-proteobacteria: similarity and evolution

The gene orders in the genomes of nine alpha-proteobacteria were compared using quantitative indices S (the relative number of common pairs of adjacent genes) and L (the mean difference between intergenic distances). A sample of 200 homologous genes, occurring in all 11 strains, was studied. In all of the genomes examined, 20 conserved, "uninterrupted" regions, including in total 63 out of 200 genes, were found. The rate of evolutionary change in the gene order widely varied in different evolutionary lineages. The highest rate (40 to 60 genome rearrangements per 100 Myr) was characteristic of the intercellular parasite Wolbachia (Rickettsiales). Computer simulation has showed that the S to L ratio observed in the sample testified that the probability of large genome rearrangements was somewhat lower than that of small ones.     

Gene Orders in Genomes of Alpha-Proteobacteria: Similarity and Evolution

The gene orders in the genomes of nine alpha-proteobacteria were compared using quantitative indices S (the relative number of common pairs of adjacent genes) and L (the mean difference between intergenic distances). A sample of 200 homologous genes, occurring in all 11 strains, was studied. In all of the genomes examined, 20 conserved, “uninterrupted” regions, including in total 63 out of 200 genes, were found. The rate of evolutionary change in the gene order widely varied in different evolutionary lineages. The highest rate (40 to 60 genome rearrangements per 100 Myr) was characteristic of the intercellular parasite Wolbachia (Rickettsiales). Computer simulation has showed that the S to L ratio observed in the sample testified that the probability of large genome rearrangements was somewhat lower than that of small ones.     

The evolutionary fate of nonfunctional DNA in the bacterial endosymbiont Buchnera aphidicola

Reduction of the genome size in endosymbiotic bacteria is the main feature linked to the adaptation to a host-associated lifestyle. We have analyzed the fate of the nonfunctional DNA in Buchnera aphidicola, the primary endosymbiont of aphids. At least 164 gene losses took place during the recent evolution of three B. aphidicola strains, symbionts of the aphids Acyrthosiphon pisum (BAp), Schizaphis graminum (BSg), and Baizongia pistacia (BBp). A typical pattern starts with the inactivation of a gene, which produces a pseudogene, and is followed by the progressive loss of its DNA. Our results show that during the period from the separation of the Aphidinae and Pemphiginae lineages (86-164 MYA) to the divergence of BAp and BSg (50-70 MYA) the half-life of a pseudogene was 23.9 Myr. For the remaining periods of evolution, the ranges of values obtained for this parameter are of the same order of magnitude. These results have revealed that a gene inactivated during B. aphidicola evolution requires 40-60 Myr to become almost completely disintegrated. Moreover, we have shown a positive correlation between the decrease in the GC content and the DNA loss for these nonfunctional DNA regions. When gene losses are classified, based on the detection of a pseudogene or otherwise of an absent gene in the modern B. aphidicola genomes, we have observed a drastic reduction of DNA length in the latter versus the former relative to the functional gene. Finally, we have also detected a slight reduction in size of the intergenic regions in the three B. aphidicola strains, when they are compared with the size of the close relative Escherichia coli.     

The mitochondrial genome sequence of Enterobius vermicularis (Nematoda: Oxyurida)--an idiosyncratic gene order and phylogenetic information for chromadorean nematodes

The complete mitochondrial genome sequence was determined for the human pinworm Enterobius vermicularis (Oxyurida: Nematoda) and used to infer its phylogenetic relationship to other major groups of chromadorean nematodes. The E. vermicularis genome is a 14,010-bp circular DNA molecule that encodes 36 genes (12 proteins, 22 tRNAs, and 2 rRNAs). This mtDNA genome lacks atp8, as reported for almost all other nematode species investigated. Phylogenetic analyses (maximum parsimony, maximum likelihood, neighbor joining, and Bayesian inference) of nucleotide sequences for the 12 protein-coding genes of 25 nematode species placed E. vermicularis, a representative of the order Oxyurida, as sister to the main Ascaridida+Rhabditida group. Tree topology comparisons using statistical tests rejected an alternative hypothesis favoring a closer relationship among Ascaridida, Spirurida, and Oxyurida, which has been supported from most studies based on nuclear ribosomal DNA sequences. Unlike the relatively conserved gene arrangement found for most chromadorean taxa, E. vermicularis mtDNA gene order is very unique, not sharing similarity to any other nematode species reported to date. This lack of gene order similarity may represent idiosyncratic gene rearrangements unique to this specific lineage of the oxyurids. To more fully understand the extent of gene rearrangement and its evolutionary significance within the nematode phylogenetic framework, additional mitochondrial genomes representing a greater evolutionary diversity of species must be characterized.     

Complete mitochondrial genome sequences of the three pelagic chaetognaths Sagitta nagae,Sagitta decipiens and Sagitta enflata

The complete nucleotide sequences of the mitochondrial genomes were determined for the three pelagic chaetognaths, Sagitta nagae, Sagitta decipiens, and Sagitta enflata. The mitochondrial genomes of these species which were 11,459, 11,121, and 12,631 bp in length, respectively, contained 14 genes (11 protein-coding genes, one transfer RNA gene, and two ribosomal RNA genes), and were found to have lost 23 genes that are present in the typical metazoan mitochondrial genome. The same mitochondrial genome contents have been reported from the benthic chaetognaths belonging to the family Spadellidae, Paraspadella gotoi and Spadella cephaloptera. Within the phylum Chaetognatha, Sagitta and Spadellidae are distantly related, suggesting that the gene loss occurred in the ancestral species of the phylum. The gene orders of the three Sagitta species are markedly different from those of the other non-Chaetognatha metazoans. In contrast to the region with frequent gene rearrangements, no gene rearrangements were observed in the gene cluster encoding COII–III, ND1–3, srRNA, and tRNA^met. Within this conserved gene cluster, gene rearrangements were not observed in the three Sagitta species or between the Sagitta and Spadellidae species. The gene order of this cluster was also assumed to be the ancestral state of the phylum.     

Inference and analysis of the relative stability of bacterial chromosomes

The stability of genomes is highly variable, both in terms of gene content and gene order. Here I calibrate the loss of gene order conservation (GOC) through time by fitting a simple probabilistic model on pairwise comparisons involving 126 bacterial genomes. The model computes the probability of separation of pairs of contiguous genes per unit of time and fits the data better than previous ones while allowing a mechanistic interpretation for the loss of GOC with time. Although the information on operons is not used in the model, I observe, as expected, that most highly conserved pairs of genes are indeed within operons. However, even the other pairs are much more conserved than expected given the observed experimental rearrangement rates. After 500 Myr, about 50% of the originally contiguous orthologues remain so in the average genome. Hence, the large majority of rearrangements must be deleterious and random genome rearrangements are unlikely to provide for positively selected structural changes. I then use the deviations from the model to define an intrinsic measure of genome stability that allowed the comparison of distantly related genomes and the inference of ancestral states. This shows that clades differ in genome stability, with cyanobacteria being the least stable and gamma-proteobacteria the most stable. Without correction for phylogeny, free-living bacteria are the least stable group of genomes, followed by pathogens, and then endomutualists. However, after correction for phylogenetic inertia (or the removal of cyanobacteria from the analysis), there is no significant association between genome stability and lifestyle or genome size. Hence, although this method has allowed uncovering some of mechanisms leading to rearrangements, we still ignore the forces that differentially shape selection upon genome stability in different species.     

A new genomic evolutionary model for rearrangements, duplications, and losses that applies across eukaryotes and prokaryotes

Genomic rearrangements have been studied since the beginnings of modern genetics and models for such rearrangements have been the subject of many papers over the last 10 years. However, none of the extant models can predict the evolution of genomic organization into circular unichromosomal genomes (as in most prokaryotes) and linear multichromosomal genomes (as in most eukaryotes). Very few of these models support gene duplications and losses--yet these events may be more common in evolutionary history than rearrangements and themselves cause apparent rearrangements. We propose a new evolutionary model that integrates gene duplications and losses with genome rearrangements and that leads to genomes with either one (or a very few) circular chromosome or a collection of linear chromosomes. Our model is based on existing rearrangement models and inherits their linear-time algorithms for pairwise distance computation (for rearrangement only). Moreover, our model predictions fit observations about the evolution of gene family sizes and agree with the existing predictions about the growth in the number of chromosomes in eukaryotic genomes.     

Use of quantitative measures of gene order similarity to phylogenetic reconstructions (exemplified by bacteria of the genus Rickettsia)

Data reflecting evolutionary changes in chromosomal gene order can be used for phylogenetic reconstructions along with the results of nucleotide sequence comparison. By the example of bacteria of the genus Rickettsia, we have shown that phylogenetic reconstructions based on quantitative estimates of the similarity and cladistic analysis of gene order data, may, in some cases, amend and fill up classical phylogenetic trees. When applied, these approaches enabled us to substantiate the hypothesis that Rickettsia felis species had split before the typhus (R. typhi, R. prowazekii) and spotted fever (R. connorii) group divergence and thus R. felis does not belong to the latter group. In general, rickettsias evolved towards increasing intracellular parasitic specialization. Five Rickettsia species whose genomes have been sequenced and annotated completely actually form an evolutionary series R. hellii-R. felis-R. connorii-R. prowazekii-R. typhi. Within this series, a reduction in genome size and rapid decrease of genome rearrangement rates (genome plasticity loss) gradually occur.     

Specific evolution of F1-like ATPases in mycoplasmas

F(1)F(0) ATPases have been identified in most bacteria, including mycoplasmas which have very small genomes associated with a host-dependent lifestyle. In addition to the typical operon of eight genes encoding genuine F(1)F(0) ATPase (Type 1), we identified related clusters of seven genes in many mycoplasma species. Four of the encoded proteins have predicted structures similar to the α, β, γ and ε subunits of F(1) ATPases and could form an F(1)-like ATPase. The other three proteins display no similarity to any other known proteins. Two of these proteins are probably located in the membrane, as they have three and twelve predicted transmembrane helices. Phylogenomic studies identified two types of F(1)-like ATPase clusters, Type 2 and Type 3, characterized by a rapid evolution of sequences with the conservation of structural features. Clusters encoding Type 2 and Type 3 ATPases were assumed to originate from the Hominis group of mycoplasmas. We suggest that Type 3 ATPase clusters may spread to other phylogenetic groups by horizontal gene transfer between mycoplasmas in the same host, based on phylogeny and genomic context. Functional analyses in the ruminant pathogen Mycoplasma mycoides subsp. mycoides showed that the Type 3 cluster genes were organized into an operon. Proteomic analyses demonstrated that the seven encoded proteins were produced during growth in axenic media. Mutagenesis and complementation studies demonstrated an association of the Type 3 cluster with a major ATPase activity of membrane fractions. Thus, despite their tendency toward genome reduction, mycoplasmas have evolved and exchanged specific F(1)-like ATPases with no known equivalent in other bacteria. We propose a model, in which the F(1)-like structure is associated with a hypothetical X(0) sector located in the membrane of mycoplasma cells.     

Genome rearrangement distances and gene order phylogeny in gamma-Proteobacteria

Genome rearrangements have been studied in 30 gamma-proteobacterial complete genomes by comparing the order of a reduced set of genes on the chromosome. This set included those genes fulfilling several characteristics, the main ones being that an ortholog was present in every genome and that none of them had been acquired by horizontal gene transfer. Genome rearrangement distances were estimated based on either the number of breakpoints or the minimal number of inversions separating two genomes. Breakpoint and inversion distances were highly correlated, indicating that inversions were the main type of rearrangement event in gamma-Proteobacteria. In general, the progressive increase in sequence-based distances between genome pairs was associated with the increase in their rearrangement-based distances but with several groups of distances not following this pattern. Compared with free-living enteric bacteria, the lineages of Pasteurellaceae were evolving, on average, to relatively higher rates of between 2.02 and 1.64, while the endosymbiotic bacterial lineages of Buchnera aphidicola and Wigglesworthia glossinidia were evolving at moderately higher rates of 1.38 and 1.35, respectively. Because we know that the rearrangement rate in the Bu. aphidicola lineage was close to zero during the last 100-150 Myr of evolution, we deduced that a much higher rate took place in the first period of lineage evolution after the divergence of the Escherichia coli lineage. On the other hand, the lineage of the endosymbiont Blochmannia floridanus did present an almost identical rate to free-living enteric bacteria, indicating that the increase in the genome rearrangement rate is not a general change associated with bacterial endosymbiosis. Phylogenetic reconstruction based on rearrangement distances showed a different topology from the one inferred by sequence information. This topology broke the proposed monophyly of the three endosymbiotic lineages and placed Bl. floridanus as a closer relative to E. coli than Yersinia pestis. These results indicate that the phylogeny of these insect endosymbionts is still an open question that will require the development of specific phylogenetic methods to confirm whether the sisterhood of the three endosymbiotic lineages is real or a consequence of a long-branch attraction phenomenon.     

Nucleotide Composition of CO1 Sequences in Chelicerata (Arthropoda): Detecting New Mitogenomic Rearrangements

Here we study the evolution of nucleotide composition in third codon-positions of CO1 sequences of Chelicerata, using a phylogenetic framework, based on 180 taxa and three markers (CO1, 18S, and 28S rRNA; 5,218?nt). The analyses of nucleotide composition were also extended to all CO1 sequences of Chelicerata found in GenBank (1,701 taxa). The results show that most species of Chelicerata have a positive strand bias in CO1, i.e., in favor of C nucleotides, including all Amblypygi, Palpigradi, Ricinulei, Solifugae, Uropygi, and Xiphosura. However, several taxa show a negative strand bias, i.e., in favor of G nucleotides: all Scorpiones, Opisthothelae spiders and several taxa within Acari, Opiliones, Pseudoscorpiones, and Pycnogonida. Several reversals of strand-specific bias can be attributed to either a rearrangement of the control region or an inversion of a fragment containing the CO1 gene. Key taxa for which sequencing of complete mitochondrial genomes will be necessary to determine the origin and nature of mtDNA rearrangements involved in the reversals are identified. Acari, Opiliones, Pseudoscorpiones, and Pycnogonida were found to show a strong variability in nucleotide composition. In addition, both mitochondrial and nuclear genomes have been affected by higher substitution rates in Acari and Pseudoscorpiones. The results therefore indicate that these two orders are more liable to fix mutations of all types, including base substitutions, indels, and genomic rearrangements.     

Microbial genome evolution: sources of variability

Comparative genome analyses of close relatives have yielded exciting insight into the sources of microbial genome variability with respect to gene content, gene order and evolution of genes with unknown functions. The genomes of free-living bacteria often carry phages and repetitive sequences that mediate genomic rearrangements in contrast to the small genomes of obligate host-associated bacteria. This suggests that genomic stability correlates with the genomic content of repeated sequences and movable genetic elements, and thereby with bacterial lifestyle. Genes with unknown functions present in a single species tend to be shorter than conserved, functional genes, indicating that the fraction of unique genes in microbial genomes has been overestimated.     

Complete mitochondrial genome DNA sequence for two ophiuroids and a holothuroid: the utility of protein gene sequence and gene maps in the analyses of deep deuterostome phylogeny

The complete mitochondrial genome sequences have been determined for the holothuroid Cucumaria miniata and two ophiuroid species Ophiopholis aculeata and Ophiura lütkeni. In addition, the nucleotide sequence of the mitochondrial protein-coding genes for the asteroid Pisaster ochraceus has been completed. Maximum-likelihood and LogDet distance analyses of concatenated protein-coding sequences produced a series of trees that did not conclusively support generally accepted models of echinoderm phylogeny. The ophiuroid data consistently demonstrated accelerated nucleotide divergence rates and lack of stationarity. This confounds the phylogenetic analyses. Molecular investigations using individual protein-coding gene alignments demonstrated that the cytochrome b gene exhibits the least deviation in rate and stationarity and generated some trees consistent with proposed echinoderm phylogenies. Phylogenies based on echinoderm mitochondrial gene rearrangements also proved problematic because of extensive variation in gene order between and within classes. A comparison of the two distinctive ophiuroid mitochondrial gene orders supports the hypothesis that O. lütkeni has a more derived mitochondrial gene order versus O. aculeata. The variation in the echinoderm mitochondrial gene maps reinforces the limitations of the application of mitochondrial gene rearrangements as a global phylogenetic tool.     

Adaptive evolution of chloroplast genome structure inferred using a parametric bootstrap approach

Background  

Genome rearrangements influence gene order and configuration of gene clusters in all genomes. Most land plant chloroplast DNAs (cpDNAs) share a highly conserved gene content and with notable exceptions, a largely co-linear gene order. Conserved gene orders may reflect a slow intrinsic rate of neutral chromosomal rearrangements, or selective constraint. It is unknown to what extent observed changes in gene order are random or adaptive. We investigate the influence of natural selection on gene order in association with increased rate of chromosomal rearrangement. We use a novel parametric bootstrap approach to test if directional selection is responsible for the clustering of functionally related genes observed in the highly rearranged chloroplast genome of the unicellular green alga Chlamydomonas reinhardtii, relative to ancestral chloroplast genomes.     

Block alignment: New representation and comparison method to study evolution of genomes

Genomes are not random sequences because natural selection has injected information in biological sequences for billions of years. Inspired by this idea, we developed a simple method to compare genomes considering nucleotide counts in subsequences (blocks) instead of their exact sequences.We introduce the Block Alignment method for comparing two genomes and based on this comparison method, define a similarity score and a distance. The presented model ignores nucleotide order in the sequence. On the other hand, in this block comparison method, due to exclusion of point mutations and small size variations, there is no need for high coverage sequencing which is responsible for the high costs of data production and storage; moreover, the sequence comparisons could be performed with higher speed.Phylogenetic trees of two sets of bacterial genomes were constructed and the results were in full agreement with their already constructed phylogenetic trees. Furthermore, a weighted and directed similarity network of each set of bacterial genomes was inferred ab initio by this model. Remarkably, the communities of these networks are in agreement with the clades of the corresponding phylogenetic trees which means these similarity networks also contain phylogenetic information about the genomes. Moreover, the block comparison method was used to distinguish rob(15;21)c-associated iAMP21 and sporadic iAMP21 rearrangements in subgroups of chromosome 21 in acute lymphoblastic leukemia. Our results show a meaningful difference between the number of contigs that mapped to chromosomes 15 and 21 in these cases. Furthermore, the presented block alignment model can select the candidate blocks to perform more accurate analysis and it is capable to find conserved blocks on a set of genomes.     

Cloning and gene map assignment of the Xiphophorus DNA ligase 1 gene

Fishes represent the stem vertebrate condition and have maintained several gene arrangements common to mammalian genomes throughout the 450 Myr of divergence from a common ancestor. One such syntenic arrangement includes the GPI-PEPD enzyme association on Xiphophorus linkage group IV and human chromosome 19. Previously we assigned the Xiphophorus homologue of the human ERCC2 gene to linkage group U5 in tight association with the CKM locus. CKM is also tightly linked to the ERCC2 locus on human chromosome 19, leading to speculation that human chromosome 19 may have arisen by fusion of two ancestral linkage groups which have been maintained in fishes. To investigate this hypothesis further, we isolated and sequenced Xiphophorus fish genomic regions exhibiting considerable sequence similarity to the human DNA ligase 1 amino acid sequence. Comparison of the fish DNA ligase sequence with those of other species suggests several modes of amino acid conservation in this gene. A 2.2-kb restriction fragment containing part of an X. maculatus DNA ligase 1 exon was used in backcross hybrid mapping with 12 enzyme or RFLP loci. Significant linkage was observed between the nucleoside phosphorylase (NP2) and the DNA ligase (LIG1) loci on Xiphophorus linkage group VI. This assignment suggests that the association of four DNA repair-related genes on human chromosome 19 may be the result of chance chromosomal rearrangements.      

Identification of functional open reading frames in chloroplast genomes

We have used a rapid computer dot-matrix comparison method to identify all DNA regions which have been evolutionarily conserved between the completely sequenced chloroplast genomes of tobacco and a liverwort. Analysis of these regions reveals 74 homologous open reading frames (ORFs) which have been conserved as to length and amino acid sequence; these ORFs also have an excess of nucleotide substitutions at silent sites of codons. Since the nonfunctional parts of these genomes have become saturated with mutations and show no sequence similarity whatsoever, the homologous ORFs are almost certainly functional. A further four pairs of ORFs show homology limited to only a short part of their putative gene products. Amino acid sequence identities range between 50 and 99%; some chloroplast proteins are seen to be among the most slowly evolving of all known proteins. A search of the nucleotide and amino acid sequence databanks has revealed several previously unidentified genes in chloroplast sequences from other species, but no new homologies to prokaryotic genes.