A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences

Summary Some simple formulae were obtained which enable us to estimate evolutionary distances in terms of the number of nucleotide substitutions (and, also, the evolutionary rates when the divergence times are known). In comparing a pair of nucleotide sequences, we distinguish two types of differences; if homologous sites are occupied by different nucleotide bases but both are purines or both pyrimidines, the difference is called type I (or transition type), while, if one of the two is a purine and the other is a pyrimidine, the difference is called type II (or transversion type). Letting P and Q be respectively the fractions of nucleotide sites showing type I and type II differences between two sequences compared, then the evolutionary distance per site is K = — (1/2) ln {(1 — 2P — Q) }. The evolutionary rate per year is then given by k = K/(2T), where T is the time since the divergence of the two sequences. If only the third codon positions are compared, the synonymous component of the evolutionary base substitutions per site is estimated by K'_S = — (1/2) ln (1 — 2P — Q). Also, formulae for standard errors were obtained. Some examples were worked out using reported globin sequences to show that synonymous substitutions occur at much higher rates than amino acid-altering substitutions in evolution.Contribution No. 1330 from the National Institute of Genetics, Mishima, 411 Japan     

Molecular evolution of bacteriophages: Discrete patterns of codon usage in T4 genes are related to the time of gene expression

Summary Patterns of codon usage in certain coliphages are adapted to expression inEscherichia coli. Bacteriophage T4 may be an exception to test the rule, as it produces eight tRNAs with specificities that are otherwise rare inE. coli. A database of all known T4 DNA sequences has been compiled, comprising 174 genes and a total of 115 kb (approximately 70% of the T4 genome). Codon usage has been examined in all T4 genes; some of these are known to be expressed before, and some after, the production of phage tRNAs. The results show two different patterns of codon usage: by comparison with the early genes, the late genes exhibit a shift in preference toward those codons recognized by the phage-encoded tRNAs. The T4 tRNAs translate A-ending codons, and it is possible that the phage acquired the tRNA genes because the mutation bias of the T4 DNA polymerase forces the T4 genome toward A+T-richness.Presented at the NATO Advanced Workshop on Genome Organization and Evolution, held in Spetses, Greece, September 1990     

Insight in to the phylogeny of polyhydroxyalkanoate biosynthesis: horizontal gene transfer

Polyhydroxyalkanoates (PHAs) are gaining more and more importance the world over due to their structural diversity and close analogy to plastics. Their biodegradability makes them extremely desirable substitutes for synthetic plastics. PHAs are produced in organisms under certain stress conditions. Here, we investigated 253 sequenced (completely and unfinished) genomes for the diversity and phylogenetics of the PHA biosynthesis. Discrepancies in the phylogenetic trees for phaA, phaB and phaC genes of the PHA biosynthesis have led to the suggestion that horizontal gene transfer (HGT) may be a major contributor for its evolution. Twenty four organisms belonging to diverse taxa were found to be involved in HGT. Among these, Bacillus cereus ATCC 14579 and Xanthomonas axonopodis pv. citri str. 306 seem to have acquired all the three genes through HGT events and have not been characterized so far as PHA producers. This study also revealed certain potential organisms such as Streptomyces coelicolor A3(2), Staphylococcus epidermidis ATCC 12228, Brucella suis 1330, Burkholderia sp., DSMZ 9242 and Leptospira interrogans serovar lai str. 56601, which can be transformed into novel PHA producers through recombinant DNA technology.     

A site-directed integration system for the nonuniversal CUGSer codon usage species Pichia farinosa by electroporation

Halotolerant yeast, Pichia farinosa, is a valuable yeast strain in fermentation industry because it produces high yield of glycerol and xylitol, and can tolerate both contamination and high-density growth during fermentation. However, the lack of genetic manipulation tools makes it less popular as a gene engineering strain. Expression systems commonly used in other yeast systems, such as Saccharomyces cerevisiae and Pichia pastoris cannot be used in P. farinosa because it translates universal Leu codon CUG as Ser. Here we reported a modified expression vector and a transformation system with enhanced efficiency in P. farinosa. The results showed that cells of OD₆₀₀ 0.8–1.0 with DTT treatment can obtain high transformation efficiency. The optimized electroporation condition was 900 V, 25 μF, and 200 Ω. The DNA concentration did not influence the transformation. Our system provides the potential not only for applying P. farinosa as an industrial strain of gene engineering, but also for studying gene function in its native host.     

Synonymous and Nonsynonymous Substitution Rates in Diatoms: A Comparison Between Chloroplast and Nuclear Genes

Rates of synonymous and nonsynonymous nucleotide substitutions and codon usage bias (ENC) were estimated for a number of nuclear and chloroplast genes in a sample of centric and pennate diatoms. The results suggest that DNA evolution has taken place, on an average, at a slower rate in the chloroplast genes than in the nuclear genes: a rate variation pattern similar to that observed in land plants. Synonymous substitution rates in the chloroplast genes show a negative association with the degree of codon usage bias, suggesting that genes with a higher degree of codon usage bias have evolved at a slower rate. While this relationship has been shown in both prokaryotes and multicellular eukaryotes, it has not been demonstrated before in diatoms. Received: 3 June 1998 / Accepted: 11 August 1998     

A vector for the site-directed,genomic integration of foreign DNA into soybean root-nodule bacteria

A non-essential DNA region carrying two different repeated sequences (RS3 and RS9) adjacent to a nitrogen fixation (nif) gene cluster has been identified previously in Bradyrhizobium japonicum strain 110. In closely related B. japonicum strains a similar genomic arrangement was found. We constructed a mobilizable plasmid vector carrying RS3 and RS9, and a kanamycin resistance cassette (nptII gene) plus suitable cloning sites inserted between the two repeated sequences. Using this vector (pRJ1035), stable integration of a lacZ gene fusion into the B. japonicum genomic RS region was achieved. The resulting strain yielded more than 10-fold higher -galactosidase activity in soybean root nodules as compared to a B. japonicum strain carrying the same lacZ fusion on a pRK290-based plasmid.     

The evolutionary landscape of the Mycobacterium tuberculosis genome

Mycobacterium tuberculosis is one of the most deadly human pathogens. The major mechanism for the adaptations of M. tuberculosis is nucleotide substitution. Previous studies have relied on the nonsynonymous-to-synonymous substitution rate (d_N/d_S) ratio as a measurement of selective constraint based on the assumed selective neutrality of synonymous substitutions. However, this assumption has been shown to be untrue in many cases. In this study, we used the substitution rate in intergenic regions (d_i) of the M. tuberculosis genome as the neutral reference, and conducted a genome-wide profiling for d_i, d_S, and the rate of insertions/deletions (indel rate) as compared with the genome of M. canettii using a 50 kb sliding window. We demonstrate significant variations in all of the three evolutionary measurements across the M. tuberculosis genome, even for regions in close vicinity. Furthermore, we identified a total of 233 genes with their d_S deviating significantly from d_i within the same window. Interestingly, d_S also varies significantly in some of the windows, indicating drastic changes in mutation rate and/or selection pressure within relatively short distances in the M. tuberculosis genome. Importantly, our results indicate that selection on synonymous substitutions is common in the M. tuberculosis genome. Therefore, the d_N/d_S ratio test must be applied carefully for measuring selection pressure on M. tuberculosis genes.     

P-related sequences inDrosophila bifasciata: A molecular clue to the understanding of P-element evolution in the genusDrosophila

Summary Two P-elements (bif1 and bif2) were isolated from a genomic library ofDrosophila bifasciata. Both elements are internally deleted and have lost the coding capacity for a functional transposase. One of the elements (bif2) contains an insert consisting of a repetitive sequence. The terminal inverted repeats and the segments necessary for passive mobility are well conserved. Element bif2 has retained rudiments of the coding sequence of exon 0 and exon 3, but the reading frame is destroyed by insertions and deletions. The comparison of theD. bifasciata P-elements with P-elements ofDrosophila melanogaster andDrosophila nebulosa reveals that the two latter sequences are more similar to each other than either of them is to theD. bifasciata elements. This finding contradicts the phylogenetic relationship of the species and can be taken as an indirect but unequivocal evidence for recent horizontal gene transfer from a relative ofD. nebulosa to the gene pool ofD. melanogaster. The P-elements ofD. bifasciata are phylogenetically ancient and have evolved independently for about 50 million years. A higher substitution rate at the third codon position as well as a predominance of conservative replacements at the amino acid level indicates that the P-elements ofD. bifasciata have been under selective constraint over a long period and that immobilization has occurred only recently.     

Analysis of the genome and proteome composition of Bdellovibrio bacteriovorus: indication for recent prey-derived horizontal gene transfer

The genome/proteome composition of Bdellovibrio bacteriovorus, the predatory microorganism that preys on other Gram-negative bacteria, has been analyzed. The study elucidates that translational selection plays a major role in genome compositional variation with higher intensity compared to other deltaproteobacteria. Other sources of variations having relatively minor contributions are local GC-bias, horizontal gene transfer and strand-specific mutational bias. The study identifies a group of AT-rich genes with distinct codon composition that is presumably acquired by Bdellovibrio recently from Gram-negative prey-bacteria other than deltaproteobacteria. The proteome composition of this species is influenced by various physico-chemical factors, viz, alcoholicity, residue-charge, aromaticity and hydropathy. Cell-wall-surface-anchor-family (CSAPs) and transporter proteins with distinct amino acid composition and specific secondary-structure also contribute notably to proteome compositional variation. CSAPs, which are low molecular-weight, outer-membrane proteins with highly disordered secondary-structure, have preference toward polar-uncharged residues and cysteine that presumably help in prey-predator interaction by providing particular bonds of attachment.     

Phylogenetic relationships of the Cobitoidea (Teleostei: Cypriniformes) inferred from mitochondrial and nuclear genes with analyses of gene evolution

The superfamily Cobitoidea of the order Cypriniformes is a diverse group of fishes, inhabiting freshwater ecosystems across Eurasia and North Africa. The phylogenetic relationships of this well-corroborated natural group and diverse clade are critical to not only informing scientific communities of the phylogeny of the order Cypriniformes, the world's largest freshwater fish order, but are key to every area of comparative biology examining the evolution of traits, functional structures, and breeding behaviors to their biogeographic histories, speciation, anagenetic divergence, and divergence time estimates. In the present study, two mitochondrial gene sequences (COI, ND4+5) and four single-copy nuclear gene segments (RH1, RAG1, EGR2B, IRBP) were used to infer the phylogenetic relationships of the Cobitoidea as reconstructed from maximum likelihood (ML) and partitioned Bayesian Analysis (BA). Analyses of the combined mitochondrial/nuclear gene datasets revealed five strongly supported monophyletic Cobitoidea families and their sister-group relationships: Botiidae+(Vaillantellidae+(Cobitidae+(Nemacheilidae+Balitoridae))). These recovered relationships are in agreement with previous systematic studies on the order Cypriniformes and/or those focusing on the superfamily Cobitoidea. Using these relationships, our analyses revealed pattern lineage- or ecological-group-specific evolution of these genes for the Cobitoidea. These observations and results corroborate the hypothesis that these group-specific-ancestral ecological characters have contributed in the diversification and/or adaptations within these groups. Positive selections were detected in RH1 of nemacheilids and in RAG1 of nemacheilids and genus Vaillantella, which indicated that evolution of RH1 (related to eye's optic sense) and RAG1 (related to immunity) genes appeared to be important for the diversification of these groups. The balitorid lineage (those species inhabiting fast-flowing riverine habitats) had, as compared with other cobitoid lineages, significantly different dN/dS, dN and dS values for ND4 and IRBP genes. These significant differences are usually indicative of weaker selection pressure, and lineage-specific evolution on genes along the balitorid lineage. Furthermore, within Cobitoidea, excluding balitorids, species living in subtropics had significantly higher dN/dS values in RAG1 and IRBP genes than those living in temperate and tropical zones. Among tropical cobitoids, genes COI, ND5, EGR2B, IRBP and RH1, had a significantly higher mean dS value than those species in subtropical and temperate groups. These findings suggest that the evolution of these genes could also be ecological-group-specific and may have played an important role in the adaptive evolution and diversification of these groups. Thus, we hypothesize that the genes included in the present study were actively involved in lineage- and/or ecological-group-specific evolutionary processes of the highly diverse Cobitoidea. These two evolutionary patterns, both subject to further testing, are hypothesized as integral in the diversification with this major clade of the world's most diverse group of freshwater fishes.     

A novel approach for gene therapy: engraftment of fibroblasts containing the artificial chromosome expression system at the site of inflammation

BACKGROUND: Rheumatoid arthritis is characterized by inflammation of the synovial tissue. High systemic doses are necessary to achieve therapeutic levels of anti-rheumatic drugs in the joints. Gene transfer might provide a more efficient delivery system for genes encoding therapeutic proteins. METHODS: The artificial chromosome expression system (ACE System) is a new non-integrating, non-viral gene expression system which functions like a natural chromosome. This technology offers advantages over current expression systems because it allows stable and predictable expression of proteins encoded by single or multiple genes over long periods of time. We are developing ex vivo gene therapy using murine artificial chromosomes containing a reporter gene (LacZ and red fluorescent protein (RFP)) for local delivery of genes in rats with adjuvant arthritis (AA). RESULTS: The delivery of the intact ACE System into rat fibroblast-like synoviocytes (FLS) and rat skin fibroblasts (RSF) was detected within 24 to 48 h post-transfection. After growing cells under selection, clones expressing LacZ and RFP were identified. Furthermore, we investigated the feasibility of local delivery of a reporter gene to the joints of rats with AA by ex vivo gene therapy. This resulted in engraftment of the injected cells in the synovial tissue microarchitecture and expression of the reporter gene. CONCLUSIONS: This work demonstrates the potential feasibility of treating arthritis and other inflammatory diseases using fibroblasts containing the ACE System as a non-viral vector for gene therapy.     

A novel approach for enhancing the catalytic efficiency of a protease at low temperature: Reduction in substrate inhibition by chemical modification

The alkaline protease, savinase was chemically modified to enhance the productivity of the enzyme at low temperatures on a complex polymeric protein (azocasein) substrate. At 5 and 15°C, savinase modified with ficol or dextran hydrolyzed fivefold more azocasein than the unmodified savinase. Kinetic studies showed that the catalytic improvements are associated with changes in uncompetitive substrate inhibition with K_i values of modified savinases sixfold higher than the unmodified savinase. Modeling of small‐angle scattering data indicates that two substrate molecules bind on opposing sides of the enzyme. The combined kinetic and structural data indicate that the polysaccharide modifier sterically blocks the allosteric site and reduces substrate inhibition. In contrast to the properties of cold‐active enzymes that generally manifest as low activation enthalpy and high flexibility, this study shows that increased activity and productivity at low temperature can be achieved by reducing uncompetitive substrate inhibition, and that this can be achieved using chemical modification with an enzyme in a commercial enzyme‐formulation. Biotechnol. Bioeng. 2009;103: 676–686. © 2009 Wiley Periodicals, Inc.     

A novel method for selective isolation of C-terminal peptides from tryptic digests of proteins by immobilized anhydrotrypsin: application to structural analyses of the tail sheath and tube proteins from bacteriophage T4

A novel method useful for selective isolation of the C-terminal peptide from a tryptic digestion mixture of a protein has been developed by taking advantage of a unique property of anhydrotrypsin, which has a strong specific affinity for the peptides containing arginine or lysine at their C-termini. Briefly, peptides produced by tryptic digestion of a protein are fractionated by affinity chromatography on a column of immobilized anhydrotrypsin. The C-terminal peptide is recovered in a breakthrough fraction, while the remainders are adsorbed on the column (unless the protein ends in arginine or lysine). The breakthrough fraction is then subjected to reversed-phase high-performance liquid chromatography in order to purify the C-terminal peptide. Using this method, we have successfully isolated the C-terminal peptides from tryptic digests of the sheath protein (gp 18) and the tube protein (gp 19) of bacteriophage T4. The analytical results on these peptides, together with the information on the N-terminal structures of the original proteins and on the nucleotide sequences of genes 18 and 19, allowed us to establish the complete primary structures of the two proteins.