The non-redundant Bacillus subtilis (NRSub) database: update 1998.

The non-redundant Bacillus subtilis database (NRSub) has been developed in the context of the sequencing project devoted to this bacterium. As this project has reached completion, the whole genome is now available as a single contig. Thanks to the ACNUC database management system and its associated retrieval system Query_win, each functional region of the genome can be accessed individually. Extra annotations have been added such as accession numbers for the genes, locations on the genetic map, codon adaptation index values, as well as cross-references with other collections. NRSub is distributed through anonymous FTP as a text file in EMBL format and as an ACNUC database. It is also possible to access NRSub through two dedicated World Wide Web servers located in France (http://acnuc. univ-lyon1.fr/nrsub/nrsub.html ) and in Japan (http://ddbjs4h.genes. nig.ac.jp/ ).     

Detecting lateral gene transfers by statistical reconciliation of phylogenetic forests

Background  

To understand the evolutionary role of Lateral Gene Transfer (LGT), accurate methods are needed to identify transferred genes and infer their timing of acquisition. Phylogenetic methods are particularly promising for this purpose, but the reconciliation of a gene tree with a reference (species) tree is computationally hard. In addition, the application of these methods to real data raises the problem of sorting out real and artifactual phylogenetic conflict.     

Efficient likelihood computations with nonreversible models of evolution

Recent advances in heuristics have made maximum likelihood phylogenetic tree estimation tractable for hundreds of sequences. Noticeably, these algorithms are currently limited to reversible models of evolution, in which Felsenstein's pulley principle applies. In this paper we show that by reorganizing the way likelihood is computed, one can efficiently compute the likelihood of a tree from any of its nodes with a nonreversible model of DNA sequence evolution, and hence benefit from cutting-edge heuristics. This computational trick can be used with reversible models of evolution without any extra cost. We then introduce nhPhyML, the adaptation of the nonhomogeneous nonstationary model of Galtier and Gouy (1998; Mol. Biol. Evol. 15:871-879) to the structure of PhyML, as well as an approximation of the model in which the set of equilibrium frequencies is limited. This new version shows good results both in terms of exploration of the space of tree topologies and ancestral G+C content estimation. We eventually apply it to rRNA sequences slowly evolving sites and conclude that the model and a wider taxonomic sampling still do not plead for a hyperthermophilic last universal common ancestor.     

Tree pattern matching in phylogenetic trees: automatic search for orthologs or paralogs in homologous gene sequence databases

MOTIVATION: Comparative sequence analysis is widely used to study genome function and evolution. This approach first requires the identification of homologous genes and then the interpretation of their homology relationships (orthology or paralogy). To provide help in this complex task, we developed three databases of homologous genes containing sequences, multiple alignments and phylogenetic trees: HOBACGEN, HOVERGEN and HOGENOM. In this paper, we present two new tools for automating the search for orthologs or paralogs in these databases. RESULTS: First, we have developed and implemented an algorithm to infer speciation and duplication events by comparison of gene and species trees (tree reconciliation). Second, we have developed a general method to search in our databases the gene families for which the tree topology matches a peculiar tree pattern. This algorithm of unordered tree pattern matching has been implemented in the FamFetch graphical interface. With the help of a graphical editor, the user can specify the topology of the tree pattern, and set constraints on its nodes and leaves. Then, this pattern is compared with all the phylogenetic trees of the database, to retrieve the families in which one or several occurrences of this pattern are found. By specifying ad hoc patterns, it is therefore possible to identify orthologs in our databases.     

Synthetic approaches to a mononucleotide prodrug of cytarabine

Synthetic pathways to a mononucleotide prodrug of cytarabine (Ara-C) bearing S-pivaloyl-2-thioethyl (tBuSATE) groups, as biolabile phosphate protections, are reported. Using a common phosphoramidite approach, two different kinds of nucleoside protecting groups have been investigated. During this study, we observed an intermolecular migration of the Boc protecting group in the course of the tert-butyldimethylsilyl ether cleavage using tetrabutyl ammonium fluoride.     

Lipid phosphate phosphatases in Arabidopsis. Regulation of the AtLPP1 gene in response to stress

An Arabidopsis thaliana gene (AtLPP1) was isolated on the basis that it was transiently induced by ionizing radiation. The putative AtLPP1 gene product showed homology to the yeast and mammalian lipid phosphate phosphatase enzymes and possessed a phosphatase signature sequence motif. Heterologous expression and biochemical characterization of the AtLPP1 gene in yeast showed that it encoded an enzyme (AtLpp1p) that exhibited both diacylglycerol pyrophosphate phosphatase and phosphatidate phosphatase activities. Kinetic analysis indicated that diacylglycerol pyrophosphate was the preferred substrate for AtLpp1p in vitro. A second Arabidopsis gene (AtLPP2) was identified based on sequence homology to AtLPP1 that was also heterologously expressed in yeast. The AtLpp2p enzyme also utilized diacylglycerol pyrophosphate and phosphatidate but with no preference for either substrate. The AtLpp1p and AtLpp2p enzymes showed differences in their apparent affinities for diacylglycerol pyrophosphate and phosphatidate as well as other enzymological properties. Northern blot analyses showed that the AtLPP1 gene was preferentially expressed in leaves and roots, whereas the AtLPP2 gene was expressed in all tissues examined. AtLPP1, but not AtLPP2, was regulated in response to various stress conditions. The AtLPP1 gene was transiently induced by genotoxic stress (gamma ray or UV-B) and elicitor treatments with mastoparan and harpin. The regulation of the AtLPP1 gene in response to stress was consistent with the hypothesis that its encoded lipid phosphate phosphatase enzyme may attenuate the signaling functions of phosphatidate and/or diacylglycerol pyrophosphate that form in response to stress in plants.     

HoSeqI: automated homologous sequence identification in gene family databases

We present a web service allowing to automatically assign sequences to homologous gene families from a set of databases. After identification of the most similar gene family to the query sequence, this sequence is added to the whole alignment and the phylogenetic tree of the family is rebuilt. Thus, the phylogenetic position of the query sequence in its gene family can be easily identified. AVAILABILITY: http://pbil.univ-lyon1.fr/software/HoSeqI/.     

Early origin of foraminifera suggested by SSU rRNA gene sequences

Foraminifera are one of the largest groups of unicellular eukaryotes with probably the best known fossil record. However, the origin of foraminifera and their phylogenetic relationships with other eukaryotes are not well established. In particular, two recent reports, based on ribosomal RNA gene sequences, have reached strikingly different conclusions about foraminifera's evolutionary position within eukaryotes. Here, we present the complete small subunit (SSU) rRNA gene sequences of three species of foraminifera. Phylogenetic analysis of these sequences indicates that they branch very deeply in the eukaryotic evolutionary tree: later than those of the amitochondrial Archezoa, but earlier than those of the Euglenozoa and other mitochondria-bearing phyla. Foraminifera are clearly among the earliest eukaryotes with mitochondria, but because of the peculiar nature of their SSU genes we cannot be certain that they diverged first, as our data suggest.      

Codon contexts in enterobacterial and coliphage genes

This investigation of the codon context of enterobacteria, plasmid, and phage protein genes was based on a search for correlations between the presence of one base type at codon position III and the presence of another base type at some other position in adjacent codons. Enterobacterial genes were compared with eukaryotic sequences for codon context effects. In enterobacterial genes, base usage at codon position III is correlated with the third position of the upstream adjacent codon and with all three positions of the downstream codon. Plasmid genes are free of context biases. Phage genes are heterogeneous: MS2 codons have no biased context, whereas lambda genes partly follow the trends of the host bacterium, and T7 genes have biased codon contexts that differ from those of the host. It has been reported that two successive third-codon positions tend to be occupied by two purines or two pyrimidines in Escherichia coli genes of low expression level. Here, the extent to which highly expressed protein genes can modulate base usage at two successive codon positions III, given the constraints on codon usage and protein sequence that act on them, was quantified. This demonstrates that the above-mentioned favored patterns are not a characteristic of weakly expressed genes but occur in all genes in which codon context can vary appreciably. The correlation between successive third-codon positions is a distinct feature of enterobacteria and of some phages, one that may result from adaptation of gene structure to translational efficiency. Conversely, codon context in yeast and human genes is biased--but for reasons unrelated to translation.