Estimate of human gene number provided by genome-wide analysis using Tetraodon nigroviridis DNA sequence.

The number of genes in the human genome is unknown, with estimates ranging from 50,000 to 90,000 (refs 1, 2), and to more than 140,000 according to unpublished sources. We have developed 'Exofish', a procedure based on homology searches, to identify human genes quickly and reliably. This method relies on the sequence of another vertebrate, the pufferfish Tetraodon nigroviridis, to detect conserved sequences with a very low background. Similar to Fugu rubripes, a marine pufferfish proposed by Brenner et al. as a model for genomic studies, T. nigroviridis is a more practical alternative with a genome also eight times more compact than that of human. Many comparisons have been made between F. rubripes and human DNA that demonstrate the potential of comparative genomics using the pufferfish genome. Application of Exofish to the December version of the working draft sequence of the human genome and to Unigene showed that the human genome contains 28,000-34,000 genes, and that Unigene contains less than 40% of the protein-coding fraction of the human genome.

Bacterial mode of replication with eukaryotic-like machinery in a hyperthermophilic archaeon.

Despite a rapid increase in the amount of available archaeal sequence information, little is known about the duplication of genetic material in the third domain of life. We identified a single origin of bidirectional replication in Pyrococcus abyssi by means of in silico analyses of cumulative oligomer skew and the identification of an early replicating chromosomal segment. The replication origin in three Pyrococcus species was found to be highly conserved, and several eukaryotic-like DNA replication genes were clustered around it. As in Bacteria, the chromosomal region containing the replication terminus was a hot spot of genome shuffling. Thus, although bacterial and archaeal replication proteins differ profoundly, they are used to replicate chromosomes in a similar manner in both prokaryotic domains.

A survey of polypeptide deformylase function throughout the eubacterial lineage.

N-terminal formylation of ribosome-synthesized polypeptides is assumed to be among the most conserved features that distinguish the eubacterial line of descent from other living phyla. In order to assess the ancientness of this trait, def genes encoding polypeptide deformylase were characterized from four eubacterial species, Lactococcus lactis, Bacillus subtilis, Calothrix PCC7601 and Thermotoga maritima, taking advantage of the conditional viability of the def mutants of Escherichia coli. Altogether, eight sequences of polypeptide deformylase have been obtained from all the eubacterial sources which were investigated, either through systematic genome sequence analysis or through genetic screening, yielding a highly homologous family. A gene putatively encoding Met-tRNAi formyltransferase, fmt, was found downstream of the deformylase gene except in L. lactis, Mycoplasma genitalium, Calothrix PCC7601 and T. maritima. These results argue strongly for the ancestral character of N-terminal formylation in eubacteria. Most of the wide deviations of amino acid usage observed in def- and fmt-encoded proteins among species is best accounted for by the nucleotide composition of genomes. Furthermore, the species of origin of each protein appears to be more recognizable than its function, considering only its amino acid composition.

Sequence relationships between integral inner membrane proteins of binding protein-dependent transport systems: evolution by recurrent gene duplications.

Periplasmic binding protein-dependent transport systems are composed of a periplasmic substrate-binding protein, a set of 2 (sometimes 1) very hydrophobic integral membrane proteins, and 1 (sometimes 2) hydrophilic peripheral membrane protein that binds and hydrolyzes ATP. These systems are members of the superfamily of ABC transporters. We performed a molecular phylogenetic analysis of the sequences of 70 hydrophobic membrane proteins of these transport systems in order to investigate their evolutionary history. Proteins were grouped into 8 clusters. Within each cluster, protein sequences displayed significant similarities, suggesting that they derive from a common ancestor. Most clusters contained proteins from systems transporting analogous substrates such as monosaccharides, oligopeptides, or hydrophobic amino acids, but this was not a general rule. Proteins from diverse bacteria are found within each cluster, suggesting that the ancestors of current clusters were present before the divergence of bacterial groups. The phylogenetic trees computed for hydrophobic membrane proteins of these permeases are similar to those described for the periplasmic substrate-binding proteins. This result suggests that the genetic regions encoding binding protein-dependent permeases evolved as whole units. Based on the results of the classification of the proteins and on the reconstructed phylogenetic trees, we propose an evolutionary scheme for periplasmic permeases. According to this model, it is probable that these transport systems derive from an ancestral system having only 1 hydrophobic membrane protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Versatility of a vector for expressing foreign polypeptides at the surface of gram-negative bacteria.

A wide variety of peptides in terms of length and sequence can be expressed at the surface of the bacterium Escherichia coli by genetic insertion into a 'permissive' site of the outer membrane protein LamB, used as a carrier. The resulting hybrid proteins essentially keep their biological activities with inserts of up to about 60 amino acid residues, and of a large range of predicted structures or hydrophobicities. This reflects a remarkable flexibility in the organization of the protein, but also in the export machinery. The method used to select such a permissive site is quite general and its potential to generate applications, including a versatile type of live bacterial vaccine, are discussed.

Location of tolerated insertions/deletions in the structure of the maltose binding protein.

In a previous study [(1987) J. Mol. Biol. 194, 663-673], we isolated ten insertion/deletion mutants (indels) of the maltose binding protein for which the maltose binding constant was only a little or not at all affected. In this paper, we have localized these mutations in the recently solved three-dimensional structure. Contrary to the general expectation, most of the insertion/deletion modifications occurred within elements of secondary structure. An analysis of the inserted residues for three indels found within alpha helices allowed an interpretation regarding protein structure accommodation to such modifications.

Genomic exploration of the hemiascomycetous yeasts: 2. Data generation and processing.

The generation of sequencing data for the hemiascomycetous yeast random sequence tag project was performed using the procedures established at GENOSCOPE. These procedures include a series of protocols for the sequencing reactions, using infra-red labelled primers, performed on both ends of the plasmid inserts in the same reaction tube, and their analysis on automated DNA sequencers. They also include a package of computer programs aimed at detecting potential assignation errors, selecting good quality sequences and estimating their useful length.

Genomic exploration of the hemiascomycetous yeasts: 1. A set of yeast species for molecular evolution studies.

The identification of molecular evolutionary mechanisms in eukaryotes is approached by a comparative genomics study of a homogeneous group of species classified as Hemiascomycetes. This group includes Saccharomyces cerevisiae, the first eukaryotic genome entirely sequenced, back in 1996. A random sequencing analysis has been performed on 13 different species sharing a small genome size and a low frequency of introns. Detailed information is provided in the 20 following papers. Additional tables available on websites describe the ca. 20000 newly identified genes. This wealth of data, so far unique among eukaryotes, allowed us to examine the conservation of chromosome maps, to identify the 'yeast-specific' genes, and to review the distribution of gene families into functional classes. This project conducted by a network of seven French laboratories has been designated 'Génolevures'.

Genomic exploration of the hemiascomycetous yeasts: 18. Comparative analysis of chromosome maps and synteny with Saccharomyces cerevisiae.

We have analyzed the evolution of chromosome maps of Hemiascomycetes by comparing gene order and orientation of the 13 yeast species partially sequenced in this program with the genome map of Saccharomyces cerevisiae. From the analysis of nearly 8000 situations in which two distinct genes having homologs in S. cerevisiae could be identified on the sequenced inserts of another yeast species, we have quantified the loss of synteny, the frequency of single gene deletion and the occurrence of gene inversion. Traces of ancestral duplications in the genome of S. cerevisiae could be identified from the comparison with the other species that do not entirely coincide with those identified from the comparison of S. cerevisiae with itself. From such duplications and from the correlation observed between gene inversion and loss of synteny, a model is proposed for the molecular evolution of Hemiascomycetes. This model, which can possibly be extended to other eukaryotes, is based on the reiteration of events of duplication of chromosome segments, creating transient merodiploids that are subsequently resolved by single gene deletion events.

Genomic exploration of the hemiascomycetous yeasts: 19. Ascomycetes-specific genes.

Comparisons of the 6213 predicted Saccharomyces cerevisiae open reading frame (ORF) products with sequences from organisms of other biological phyla differentiate genes commonly conserved in evolution from 'maverick' genes which have no homologue in phyla other than the Ascomycetes. We show that a majority of the 'maverick' genes have homologues among other yeast species and thus define a set of 1892 genes that, from sequence comparisons, appear 'Ascomycetes-specific'. We estimate, retrospectively, that the S. cerevisiae genome contains 5651 actual protein-coding genes, 50 of which were identified for the first time in this work, and that the present public databases contain 612 predicted ORFs that are not real genes. Interestingly, the sequences of the 'Ascomycetes-specific' genes tend to diverge more rapidly in evolution than that of other genes. Half of the 'Ascomycetes-specific' genes are functionally characterized in S. cerevisiae, and a few functional categories are over-represented in them.

Genomic exploration of the hemiascomycetous yeasts: 20. Evolution of gene redundancy compared to Saccharomyces cerevisiae.

We have evaluated the degree of gene redundancy in the nuclear genomes of 13 hemiascomycetous yeast species. Saccharomyces cerevisiae singletons and gene families appear generally conserved in these species as singletons and families of similar size, respectively. Variations of the number of homologues with respect to that expected affect from 7 to less than 24% of each genome. Since S. cerevisiae homologues represent the majority of the genes identified in the genomes studied, the overall degree of gene redundancy seems conserved across all species. This is best explained by a dynamic equilibrium resulting from numerous events of gene duplication and deletion rather than by a massive duplication event occurring in some lineages and not in others.

Genomic exploration of the hemiascomycetous yeasts: 21. Comparative functional classification of genes.

We explored the biological diversity of hemiascomycetous yeasts using a set of 22000 newly identified genes in 13 species through BLASTX searches. Genes without clear homologue in Saccharomyces cerevisiae appeared to be conserved in several species, suggesting that they were recently lost by S. cerevisiae. They often identified well-known species-specific traits. Cases of gene acquisition through horizontal transfer appeared to occur very rarely if at all. All identified genes were ascribed to functional classes. Functional classes were differently represented among species. Species classification by functional clustering roughly paralleled rDNA phylogeny. Unequal distribution of rapidly evolving, ascomycete-specific, genes among species and functions was shown to contribute strongly to this clustering. A few cases of gene family amplification were documented, but no general correlation could be observed between functional differentiation of yeast species and variations of gene family sizes. Yeast biological diversity seems thus to result from limited species-specific gene losses or duplications, and for a large part from rapid evolution of genes and regulatory factors dedicated to specific functions.

HIV-1 isolates are rapidly evolving quasispecies: evidence for viral mixtures and preferred nucleotide substitutions.

RNA viruses are renowed for their genetic variability. The human immunodeficiency viruses (HIV) are no exception. A rapid method has been established for the genetic identification and differentiation of viral strains based on the sequencing of many M13 clones of gene-amplified products. Some isolates are internally relatively homogeneous while others are heterogeneous. There was no correlation between virus complexity and disease stage. One isolate was in fact a mixture of two distinct strains. A strong preference for G----A base substitutions was observed. These data indicate that HIV isolates cannot be described in simple molecular terms and should rather be considered as quasispecies.

[Extraction of symbolic determinants common to a family of biological sequences]. / Extraction de déterminants symboliques communs à une famille de bioséquences.

A set of sequences can be defined by their common subsequences, and the length of these is a measure of the overall resemblance of the set. Each subsequence corresponds to a succession of symbols embedded in every sequence, following the same order but not necessarily contiguous. Determining the longest common subsequence (LCS) requires the exhaustive testing of all possible common subsequences, which sum up to about 2L, if L is the length of the shortest sequence. We present a polynomial algorithm (O(n X L4), where n is the number of sequences) for generating strings related to the LCS and constructed with the sequence alphabet and an indetermination symbol. Such strings are iteratively improved by deleting indetermination symbols and concomitantly introducing the greatest number of alphabet symbols. Processed accordingly, nucleic acid and protein sequences lead to key-words encompassing the salient positions of homologous chains, which can be used for aligning or classifying them, as well as for finding related sequences in data banks.

Progress in the identification of interaction sites on the periplasmic maltose binding protein from E coli.

The periplasmic maltose binding protein (MBP) is required for the high affinity transport of maltose and maltodextrins and for chemotaxis towards these sugars. In these functions, MBP interacts with proteins of the cytoplasmic membrane: MalF and MalG for transport, Tar for chemotaxis. A large number of MBP mutations have been isolated by us and other laboratories. We grouped these mutations into classes depending on the interactions affected and we represented the corresponding residues on the 3-D model for MBP so as to further identify the sites of MBP interacting with the MalF-MalG complex and with the Tar protein. MBP (like the other binding proteins) is composed of 2 lobes enclosing a cleft where the substrate binds. The face of the protein opposite the cleft seems to interact neither with MalF-MalG nor with Tar. The other face, corresponding to the cleft, contains sites for interactions with MalF-MalG and Tar. These sites appear to cover both sides of the cleft and may overlap in part. The present definition of the interaction sites suggests further that MBP has different in vivo orientations when it interacts with MalF-MalG or with Tar. This work constitutes an additional step in combining the use of genetic and structural analysis to define the interaction sites on MBP. Because of the structural similarities between periplasmic binding proteins, the regions of interaction defined could be relevant for other members of this family.

The maltoporin of Salmonella typhimurium: sequence and folding model.

The sequence of the lamB gene from Salmonella typhimurium was determined. It encodes the precursor to the LamB protein from S. typhimurium (pre-LamBS.t.; 452 residues) which presents extensive homologies with the pre-LamB protein from Escherichia coli (pre-LamBE.c.; 446 residues). The first third of pre-LamBS.t. is the most conserved, with 4% changes and strict identity between the signal peptides. The last two-third contains five "variable" segments where more than 50% of the residues are changed with respect to LamBE.c.. The three first variable segments are 8 to 14 residues long and contain only substitutions, while the two more distal ones are 24 and 29 residues long and also include insertions and deletions. It is remarkable that the variable segments correspond essentially to regions predicted to be extramembranous loops on our 2D folding model for LamBE.c.; they alternate with conserved predicted transmembranous segments. Four of the variable regions were predicted to be cell-surface-exposed loops on the basis of genetic and immunological data, while one of them (region II) was predicted to be periplasmic on the sole basis of folding rules. The LamB protein from S. typhimurium can substitute for the LamB protein from E. coli for maltodextrins binding and transport, but not for infection by any of the known E. coli phages using LamBE.c. for adsorption. A tetrapeptide, RGDS, assumed to be responsible for mammalian cell aggregation by LamBE.c. is conserved in LamBS.t., suggesting that it could have a functional role. The conservation of the binding and transport activity can be accounted for by the conservation of the regions known to be directly involved, namely the first third of the protein and a region corresponding to 352 to 374 of LamBS.t.. The phage resistance can be attributed to the variability of the four cell-surface-exposed loops previously identified as essential for phage adsorption. These results, together with those obtained with polyclonal and monoclonal antibodies directed against known LamB regions, strongly support the folding model presented for LamBE.c. and the idea that it can essentially be extended to LamBS.t., except perhaps for a region between residues 155 and 245. We propose that the existence of variable regions is due essentially, and perhaps only, to the local lack of structural constraints in the protein. The intergenic region between lamB and the following gene, malM, comprises conserved segments, including one palindromic unit.

The BIME family of bacterial highly repetitive sequences.

Palindromic units (PU or REP) were initially defined as a DNA sequence of 40 nucleotides which is highly repeated in the genome of several enterobacteria and found in clusters of up to six copies. It appears now that PU belong to a larger repeated DNA element, of up to 300 nucleotides, called BIME for bacterial interspersed mosaic element. BIME is a mosaic combination of ten small DNA motifs, including the PU sequence. A central question concerning BIME is to determine whether they play a critical role within the cell. BIME exhibit only limited effects on local gene expression; it seems unlikely that these weak effects alone can account for the high BIME sequence homogeneity. It has recently been shown that DNA gyrase and DNA polymerase I are able to specifically recognize BIME DNA in vitro. These findings suggest that BIME could play a role in the functional organization of the bacterial nucleoid. Hypotheses on their origin and evolution are discussed.

Phylogenetic analyses of the ATP-binding constituents of bacterial extracytoplasmic receptor-dependent ABC-type nutrient uptake permeases.

Thirty-eight ATP-binding cassette (ABC) protein constituents of bacterial extracytoplasmic receptor-dependent nutrient uptake systems, including one homologous chloroplast protein were analysed for sequence conservation and phylogenetic relatedness. The proteins were generally found to cluster in accordance with the clustering patterns previously observed for the extracytoplasmic receptors and the transmembrane channel-forming constituents of these permeases. The results suggest that these transport systems evolved from a single primordial system with minimal shuffling of the three dissimilar protein constituents of the systems.