Takashi Kunisawa

Gene Order Breakpoint Evidence in Animal Mitochondrial Phylogeny

Abstract. Multiple genome rearrangement methodology facilitates the inference of animal phylogeny from gene orders on the mitochondrial genome. The breakpoint distance is preferable to other, highly correlated but computationally more difficult, genomic distances when applied to these data. A number of theories of metazoan evolution are compared to phylogenies reconstructed by ancestral genome optimization, using a minimal total breakpoints criterion. The notion of unambiguously reconstructed segments is introduced as a way of extracting the invariant aspects of multiple solutions for a given ancestral genome; this enables a detailed reconstruction of the evolution of non-tRNA mitochondrial gene order.

Parametric genome rearrangement

Algorithms inspired by comparative genomics calculate an edit distance between two linear orders based on elementary edit operations such as inversion, transposition and reciprocal translocation. All operations are generally assigned the same weight, simply by default, because no systematic empirical studies exist verifying whether algorithmic outputs involve realistic proportion of each. Nor do we have data on how weights should vary with the length of the inverted or transposed segment of the chromosome. In this paper, we present a rapid algorithm that allows each operation to take on a range of weights, producing an relatively tight upper bound on the distance between single-chromosome genomes, by means of a greedy search with look-ahead. The efficiency of this algorithm allows us to test random genomes for each parameter setting, to detect gene order similarity and to infer the parameter values most appropriate to the phylogenetic domain under study. We apply this method to genome segments in which the same gene order is conserved in Escherichia coli and Bacillus subtilis, as well as to the gene order in human versus Drosophila mitochondrial genomes. In both cases, we conclude that it is most appropriate to assign somewhat more than twice the weight to transpositions and inverted transpositions than to inversions. We also explore segment-length weighting for fungal mitochondrial gene orders.

Gene arrangements and branching orders of gram-positive bacteria

The availability of complete genomic sequence data allows one to develop new methods of reconstructing phylogenetic trees. A simple method of reconstructing branching orders based on gene transposition (or lateral transfer) is presented. It is argued that specific gene arrangements on four different genomes could determine a branching order. A computer search for such gene arrangements was carried out against gene order data of completely sequenced Gram-positive bacteria. Gene arrangements around ribosomal protein S4 gene, murC (UDP-N-acetylmuramate:alanine ligase) gene and dnaE (DNA polymerase III alpha chain) gene each suggest a branching order in which actinobacteria with a high genomic G+C content first branched off from other Gram-positives with a low G+C content and then a split occurred between Mycoplasma species and a group closely related to Bacillus subtilis. A recently sequenced thermophilic bacterium Thermoanaerobacter tengcongensis is suggested to have branched off from the lineage leading to the low G+C Gram-positives prior to the split between the Mycoplasma and Bacillus groups. By contrast to the indel analysis in which a single evolutionary event of insertion or deletion of a signature sequence is assumed, the present method does not necessarily require such a parsimonious assumption of gene transposition.

Identification and chromosomal distribution of DNA sequence segments conserved since divergence of Escherichia coli and Bacillus subtilis

DNA sequence segments conserved since divergence of Escherichia coli and Bacillus subtilis were identified, using the GenBank sequence database. Chromosomal locations of the conserved segments were compared between the two bacteria, and the following three features were observed. (1) Although the two genomes are nearly identical in size, chromosomal arrangements of the conserved segments are considerably different from each other. (2) In many cases, chromosomal locations of a conserved segment in the two species have deviated from each other by a multiple of 60°. (3) There are many instances in which a contiguous segment in one genome is split into two or more segments located at distinct positions in the other genome, and these split segments were found to tend to lie on the E. coli or B. subtilis genome separated by distances of multiples of 60°. On the basis of these observations, genome organizations of the two bacteria were discussed in terms of genome doublings as well as random chromosomal rearrangements.

Gene arrangements characteristic of the phylum Actinobacteria

The availability of genomic sequence data allows new challenges to various biological problems. One of such attempts is the extraction of phylogenetic information from gene order data of genomes. Phylogenetic inferences are most commonly carried out on the basis of 16S rRNA trees, which sometimes produce unresolved or unreliable branching orders. One example for such a low resolution is recognized in the branching pattern among the phylum Actinobacteria. Here, gene arrangements characteristic of the Actinobacteria were identified, based on which Symbiobacterium thermophilum is phylogenetically placed outside that phylum, this being in contrast to 16S rRNA trees and to the current taxonomy in GenBank. Three transposition suggestive arrangements were found which support a notion that Rubrobacter xylanophilus is the earliest diverging species among the completely sequenced Actinobacteria. The gene arrangements identified here serve as a complement to previously reported indels and proteins characteristic of this phylum.

Evaluation of the phylogenetic position of the sulfate-reducing bacterium Thermodesulfovibrio yellowstonii (phylum Nitrospirae) by means of gene order data from completely sequenced genomes.

The phylogenetic placement of Thermodesulfovibrio yellowstonii was investigated on the basis of gene order data from completely sequenced bacterial genomes. T. yellowstonii was found to share four gene arrangements characteristic of the Proteobacteria, Aquificae, Planctomycetes, Spirochaetes, Bacteroidetes, Chlorobi, Acidobacteria, Verrucomicrobia and termite group 1, the presence of which defines superphylum 2. The remaining phyla show sets of alternative gene arrangements and form superphylum 1. An analysis of conserved gene pairs showed that the overall genome organization of T. yellowstonii is most similar to that of deltaproteobacteria. Three arrangements that suggest gene translocations were identified that are likely to have occurred in a common ancestor of T. yellowstonii and the Proteobacteria exclusive of virtually all other major bacterial phyla. The translocation events suggest the closest evolutionary relationship between T. yellowstonii and the Proteobacteria.

Inference of the phylogenetic position of the phylum Deferribacteres from gene order comparison

The phylogenetic placement of the phylum Deferribacteres was investigated on the basis of gene order comparisons of completely sequenced bacterial genomes. Two completely sequenced Deferribacteres species share five sets of gene arrangements with a group of phyla, Proteobacteria, Aquificae, Planctomycetes, Spirochaetes, Bacteroidetes, Chlorobi, Acidobacteria, Verrucomicrobia, Elusimicrobia and Nitrospirae, while the other group of phyla, Synergistetes, Firmicutes, Actinobacteria, Thermotogae, Chloroflexi and Deinococcus-Thermus, Fusobacteria, shares alternative sets of gene arrangements, suggesting that the Deferribacteres is classified in the former group of phyla. Gene transfers that are thought to have occurred in a common ancestor of the Deferribacteres, Deltaproteobacteria and Nitrospirae exclusive of virtually all other phyla were identified, which suggests that the Deferribacteres is phylogenetically proximal to the Proteobacteria and Nitrospirae.

A gene order database of plastid genomes

A gene order database of 32 completely sequenced plastid genomes was developed. The data structure is formally identical to that of the feature tables in the major GenBank/EMBL/DDBJ databases. The quality of annotations was largely improved. A normalizing gene-labeling system across the complete plastid genomes was developed so that comparative studies are made available without having to go back to sequence analysis. Many incorrect coordinates of tRNA-encoding regions found in the major databases were corrected. We attempted to distinctively label tRNA genes with the anticodon sequence CAT, which encodes either the initiator tRNA, elongator tRNA, or Ile-tRNA. The database is available at http://www.rs.noda.tus.ac.jp/~kunisawa.

Accumulation pattern of amino acid substitutions in protein evolution

SummaryA simple method for the evolutionary analysis of amino acid sequence data is presented and used to examine whether the number of variable sites (NVS) of a protein is constant during its evolution. The NVSs for hemoglobin and for mitochondrial cytochrome c are each found to be almost constant, and the ratio between the NVSs is close to the ratio between the unit evolutionary periods. This indicates that the substitution rate per variable site is almost uniform for these proteins, as the neutral theory claims. An advantage of the present analysis is that it can be done without knowledge of paleontological divergence times and can be extended to bacterial proteins such as bacterial c-type cytochromes. It is suggested that the NVS of cytochrome c has been almost constant even over the long period (ca. 3.0 billion years) of bacterial evolution but that at least two different substitution rates are necessary to describe the accumulated changes in the sequence. This “two clock” interpretation is consistent with fossil evidence for the appearance times of photosynthetic bacteria and eukaryotes.

Theoretical evaluation of transcriptional pausing effect on the attenuation in trp leader sequence

The effect of transcriptional pausing on attenuation is investigated theoretically on the basis of the attenuation control mechanism presented by Oxender et al. (Oxender, D. L., G. Zurawski, and C. Yanofsky, 1979, Proc. Natl. Acad. Sci. USA. 76:5524-5528). An extended stochastic model including the RNA polymerase pausing in the leader region is developed to calculate the probability of relative position between the RNA polymerase transcribing the trp leader sequence and the ribosome translating the transcript. The present study results in a new rationale that the transcriptional pausing site in the leader sequence makes the attenuation control both more sensitive as an on/off switch and less sensitive to variations in the concentration of cellular metabolites not connected with the need for expressing, or not expressing, the particular operon. It is also proposed that the transcriptional pausing diminishes the dependence of attenuation control characteristics on the number of nucleotides in the leader sequence. This result may be useful for understanding the attenuation control efficiencies of other amino acid leader sequences with different lengths of nucleotides.