Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440

Pseudomonas putida is a metabolically versatile saprophytic soil bacterium that has been certified as a biosafety host for the cloning of foreign genes. The bacterium also has considerable potential for biotechnological applications. Sequence analysis of the 6.18 Mb genome of strain KT2440 reveals diverse transport and metabolic systems. Although there is a high level of genome conservation with the pathogenic Pseudomonad Pseudomonas aeruginosa (85% of the predicted coding regions are shared), key virulence factors including exotoxin A and type III secretion systems are absent. Analysis of the genome gives insight into the non-pathogenic nature of P. putida and points to potential new applications in agriculture, biocatalysis, bioremediation and bioplastic production.

The genome sequence of Schizosaccharomyces pombe

We have sequenced and annotated the genome of fission yeast (Schizosaccharomyces pombe), which contains the smallest number of protein-coding genes yet recorded for a eukaryote: 4,824. The centromeres are between 35 and 110 kilobases (kb) and contain related repeats including a highly conserved 1.8-kb element. Regions upstream of genes are longer than in budding yeast (Saccharomyces cerevisiae), possibly reflecting more-extended control regions. Some 43% of the genes contain introns, of which there are 4,730. Fifty genes have significant similarity with human disease genes; half of these are cancer related. We identify highly conserved genes important for eukaryotic cell organization including those required for the cytoskeleton, compartmentation, cell-cycle control, proteolysis, protein phosphorylation and RNA splicing. These genes may have originated with the appearance of eukaryotic life. Few similarly conserved genes that are important for multicellular organization were identified, suggesting that the transition from prokaryotes to eukaryotes required more new genes than did the transition from unicellular to multicellular organization.

General method of rapid Smith/Birnstiel mapping adds for gap closure in shotgun microbial genome sequencing projects: application to Pseudomonas putida KT2440

A physical mapping strategy has been developed to verify and accelerate the assembly and gap closure phase of a microbial genome shotgun-sequencing project. The protocol was worked out during the ongoing Pseudomonas putida KT2440 genome project. A macro-restriction map was constructed by linking probe hybridisation of SwaI- or I-CeuI-restricted chromosomes to serve as a backbone for the quick quality control of sequence and contig assemblies. The library of PCR-generated SwaI linking probes was derived from the sequence assembly after 3- and 6-fold genome coverage. In order to support gap closure in regions with ambiguous assemblies such as the repetitive sequence of the seven ribosomal operons, high-resolution Smith/Birnstiel maps were generated by Southern hybridisation of pulsed-field gel electrophoresis-separated rare-cutter complete/frequent-cutter partial digestions with rare-cutter fragment end probes. Overall 1.5 Mb of the 6.1 Mb P.putida KT2440 genome has been subjected to high-resolution physical mapping in order to align assemblies generated from shotgun sequencing.

Automated sample-preparation technologies in genome sequencing projects

A robotic workstation system (BioRobot 96OO, QIAGEN) and a 96-well UV spectrophotometer (Spectramax 250, Molecular Devices) were integrated in to the process of high-throughput automated sequencing of double-stranded plasmid DNA templates. An automated 96-well miniprep kit protocol (QIAprep Turbo, QIAGEN) provided high-quality plasmid DNA from shotgun clones. The DNA prepared by this procedure was used to generate more than two mega bases of final sequence data for two genomic projects (Arabidopsis thaliana and Schizosaccharomyces pombe), three thousand expressed sequence tags (ESTs) plus half a mega base of human full-length cDNA clones, and approximately 53,000 single reads for a whole genome shotgun project (Pseudomonas putida).

Sequence and analysis of chromosome 4 of the plant Arabidopsis thaliana

The higher plant Arabidopsis thaliana (Arabidopsis) is an important model for identifying plant genes and determining their function. To assist biological investigations and to define chromosome structure, a coordinated effort to sequence the Arabidopsis genome was initiated in late 1996. Here we report one of the first milestones of this project, the sequence of chromosome 4. Analysis of 17.38 megabases of unique sequence, representing about 17% of the genome, reveals 3,744 protein coding genes, 81 transfer RNAs and numerous repeat elements. Heterochromatic regions surrounding the putative centromere, which has not yet been completely sequenced, are characterized by an increased frequency of a variety of repeats, new repeats, reduced recombination, lowered gene density and lowered gene expression. Roughly 60% of the predicted protein-coding genes have been functionally characterized on the basis of their homology to known genes. Many genes encode predicted proteins that are homologous to human and Caenorhabditis elegans proteins.

High-throughput robotic system for sequencing of microbial genomes

A high-throughput robotic workstation system was used for double-stranded plasmid DNA template preparation and sequencing reaction setup to streamline the sequencing process in genome projects. All 96-well miniprep kits that were tested provided high quality plasmid DNA suitable for fluorescent DNA sequencing. After quantitation in a 96-well UV spectrophotometer, the plasmid DNA was used as template to automatically set up sequencing reactions. The setup was controlled by spread sheets that were imported into the robotic system. We utilized this integrated system to prepare all necessary shotgun templates for our contributions to a number of large-scale genome projects as well as a full-length cDNA sequencing project.

Introns and intein coding sequence in the ribonucleotide reductase genes of Bacillus subtilis temperate bacteriophage SPbeta

The two putative ribonucleotide reductase subunits of the Bacillus subtilis bacteriophage SPbeta are encoded by the bnrdE and bnrdF genes that are highly similar to corresponding host paralogs, located on the opposite replication arm. In contrast to their bacterial counterparts, bnrdE and bnrdF each are interrupted by a group I intron, efficiently removed in vivo by mRNA processing. The bnrdF intron contains an ORF encoding a polypeptide similar to homing endonucleases responsible for intron mobility, whereas the bnrdE intron has no obvious trace of coding sequence. The downstream bnrdE exon harbors an intervening sequence not excised at the level of the primary transcript, which encodes an in-frame polypeptide displaying all the features of an intein. Presently, this is the only intein identified in bacteriophages. In addition, bnrdE provides an example of a group I intron and an intein coding sequence within the same gene.

Analysis of 1.9 Mb of contiguous sequence from chromosome 4 of Arabidopsis thaliana

The plant Arabidopsis thaliana (Arabidopsis) has become an important model species for the study of many aspects of plant biology. The relatively small size of the nuclear genome and the availability of extensive physical maps of the five chromosomes provide a feasible basis for initiating sequencing of the five chromosomes. The YAC (yeast artificial chromosome)-based physical map of chromosome 4 was used to construct a sequence-ready map of cosmid and BAC (bacterial artificial chromosome) clones covering a 1.9-megabase (Mb) contiguous region, and the sequence of this region is reported here. Analysis of the sequence revealed an average gene density of one gene every 4.8 kilobases (kb), and 54% of the predicted genes had significant similarity to known genes. Other interesting features were found, such as the sequence of a disease-resistance gene locus, the distribution of retroelements, the frequent occurrence of clustered gene families, and the sequence of several classes of genes not previously encountered in plants.

The complete genome sequence of the gram-positive bacterium Bacillus subtilis

Bacillus subtilis is the best-characterized member of the Gram-positive bacteria. Its genome of 4,214,810 base pairs comprises 4,100 protein-coding genes. Of these protein-coding genes, 53% are represented once, while a quarter of the genome corresponds to several gene families that have been greatly expanded by gene duplication, the largest family containing 77 putative ATP-binding transport proteins. In addition, a large proportion of the genetic capacity is devoted to the utilization of a variety of carbon sources, including many plant-derived molecules. The identification of five signal peptidase genes, as well as several genes for components of the secretion apparatus, is important given the capacity of Bacillus strains to secrete large amounts of industrially important enzymes. Many of the genes are involved in the synthesis of secondary metabolites, including antibiotics, that are more typically associated with Streptomyces species. The genome contains at least ten prophages or remnants of prophages, indicating that bacteriophage infection has played an important evolutionary role in horizontal gene transfer, in particular in the propagation of bacterial pathogenesis.

The nucleotide sequence of Saccharomyces cerevisiae chromosome XII

The yeast Saccharomyces cerevisiae is the pre-eminent organism for the study of basic functions of eukaryotic cells. All of the genes of this simple eukaryotic cell have recently been revealed by an international collaborative effort to determine the complete DNA sequence of its nuclear genome. Here we describe some of the features of chromosome XII.

The nucleotide sequence of Saccharomyces cerevisiae chromosome XIV and its evolutionary implications

In 1992 we started assembling an ordered library of cosmid clones from chromosome XIV of the yeast Saccharomyces cerevisiae. At that time, only 49 genes were known to be located on this chromosome and we estimated that 80% to 90% of its genes were yet to be discovered. In 1993, a team of 20 European laboratories began the systematic sequence analysis of chromosome XIV. The completed and intensively checked final sequence of 784,328 base pairs was released in April, 1996. Substantial parts had been published before or had previously been made available on request. The sequence contained 419 known or presumptive protein-coding genes, including two pseudogenes and three retrotransposons, 14 tRNA genes, and three small nuclear RNA genes. For 116 (30%) protein-coding sequences, one or more structural homologues were identified elsewhere in the yeast genome. Half of them belong to duplicated groups of 6-14 loosely linked genes, in most cases with conserved gene order and orientation (relaxed interchromosomal synteny). We have considered the possible evolutionary origins of this unexpected feature of yeast genome organization.

Transposon mutagenesis reinforces the correlation between Mycoplasma pneumoniae cytoskeletal protein HMW2 and cytadherence

A new genetic locus associated with Mycoplasma pneumoniae cytadherence was previously identified by transposon mutagenesis with Tn4001. This locus maps approximately 160 kbp from the genes encoding cytadherence-associated proteins HMW1 and HMW3, and yet insertions therein result in loss of these proteins and a hemadsorption-negative (HA-) phenotype, prompting the designation cytadherence-regulatory locus (crl). In the current study, passage of transformants in the absence of antibiotic selection resulted in loss of the transposon, a wild-type protein profile, and a HA+ phenotype, underscoring the correlation between crl and M. pneumoniae cytadherence. Nucleotide sequence analysis of crl revealed open reading frames (ORFs) orfp65, orfp216, orfp41, and orfp24, arranged in tandem and flanked by a promoter-like and a terminator-like sequence, suggesting a single transcriptional unit, the P65 operon. The 5' end of orfp65 mRNA was mapped by primer extension, and a likely promoter was identified just upstream. The product of each ORF was identified by using antisera prepared against fusion proteins. The previously characterized surface protein P65 is encoded by orfp65, while the 190,000 Mr cytadherence-associated protein HMW2 is a product of orfp216. Proteins with sizes of 47,000 and 41,000 Mr and unknown function were identified for orfp41 and orfp24, respectively. Structural analyses of HMW2 predict a periodicity highly characteristic of a coiled-coil conformation and five leucine zipper motifs, indicating that HMW2 probably forms dimers in vivo, which is consistent with a structural role in cytadherence. Each transposon insertion mapped to orfp216 but affected the levels of all products of the P65 operon. HMW2 is thought to form a disulfide-linked dimer, formerly designated HMW5, and examination of an hmw2 deletion mutant confirms that HMW5 is a product of the hmw2 gene.

Comparative analysis of the genomes of the bacteria Mycoplasma pneumoniae and Mycoplasma genitalium

The sequenced genomes of the two closely related bacteria Mycoplasma genitalium and Mycoplasma pneumoniae were compared with emphasis on genome organization and coding capacity. All the 470 proposed open reading frames (ORFs) of the smaller M.genitalium genome (580 kb) were contained in the larger genome (816 kb) of M.pneumoniae. There were some discrepancies in annotation, but inspection of the DNA sequences showed that the corresponding DNA was always present in M. pneumoniae. The two genomes could be subdivided into six segments. The order of orthologous genes was well conserved within individual segments but the order of these segments in both bacteria was different. We explain the different organization of the segments by translocation via homologous recombination. The translocations did not disturb the continuous bidirectional course of transcription in both genomes, starting at the proposed origin of replication. The additional 236 kb in M.pneumoniae,compared with theM.genitalium genome, were coding for 209 proposed ORFs not identified in M.genitalium. Of these ORFs, 110 were specific to M.pneumoniae exhibiting no significant similarity to M.genitalium ORFs, while 76 ORFs were amplifications of ORFs existing mainly as single copies in M. genitalium. In addition, 23 ORFs containing a copy of either one of the three repetitive DNA sequences RepMP2/3, RepMP4 and RepMP5 were annotated in M.pneumoniae but not in M.genitalium,although similar DNA sequences were present. TheM.pneumoniae-specific genes included a restriction-modification system, two transport systems for carbohydrates, the complete set of three genes coding for the arginine dihydrolase pathway and 14 copies of the repetitive DNA sequence RepMP1 which were part of several different translated genes with unknown function.

Complete sequence analysis of the genome of the bacterium Mycoplasma pneumoniae

The entire genome of the bacterium Mycoplasma pneumoniae M129 has been sequenced. It has a size of 816,394 base pairs with an average G+C content of 40.0 mol%. We predict 677 open reading frames (ORFs) and 39 genes coding for various RNA species. Of the predicted ORFs, 75.9% showed significant similarity to genes/proteins of other organisms while only 9.9% did not reveal any significant similarity to gene sequences in databases. This permitted us tentatively to assign a functional classification to a large number of ORFs and to deduce the biochemical and physiological properties of this bacterium. The reduction of the genome size of M. pneumoniae during its reductive evolution from ancestral bacteria can be explained by the loss of complete anabolic (e.g. no amino acid synthesis) and metabolic pathways. Therefore, M. pneumoniae depends in nature on an obligate parasitic lifestyle which requires the provision of exogenous essential metabolites. All the major classes of cellular processes and metabolic pathways are briefly described. For a number of activities/functions present in M. pneumoniae according to experimental evidence, the corresponding genes could not be identified by similarity search. For instance we failed to identify genes/proteins involved in motility, chemotaxis and management of oxidative stress.

Sequence analysis and characterization of the hmw gene cluster of Mycoplasma pneumoniae

Mycoplasma pneumoniae (Mp) cytadherence requires the proper anchoring of cytadhesin proteins in the mycoplasmal membrane at an attachment organelle through their interaction with a cytoskeleton-like network of accessory proteins that includes HMW1 and HMW3. Approximately 8.25 kb of Mp DNA was sequenced, beginning at the 3' end of the hmw3 gene and continuing through hmw1. Comparison of the resulting deduced amino acid (aa) sequence with N terminus and internal peptide aa sequences from purified HMW1 permitted definitive identification of hmw1. HMW1 was characterized with respect to structure, hydrophobicity, possible phosphoacceptor sites and expression of the Mp recombinant protein in Escherichia coli. In addition, HMW1 membrane topography was examined for antibody accessibility on the mycoplasmal surface. hmw3 and hmw1 flank four open reading frames (ORFs) spanning approximately 4.3 kb and in the same orientation as the hmw genes. The sequences of their deduced products were evaluated for likely structural features and comparison with protein data banks. Finally, the Mp rpsD analog was identified immediately downstream from hmw1.

The P200 protein of Mycoplasma pneumoniae shows common features with the cytadherence-associated proteins HMW1 and HMW3

The gene coding for the P200 protein of the bacterium, Mycoplasma pneumoniae (Mp), was cloned and sequenced. The sequence-derived data and biochemical data indicated that P200 has several features in common with the well characterized cytadherence-associated proteins, HMW1 and HMW3. These features consist of abnormal migration in SDS-PAGE, a central acidic domain with a high Pro content, repeated peptide blocks within the Pro-rich domain and P200 partitioning similar to HMW1 and HMW3 in the insoluble fraction after extraction of Mp with the detergent Triton X-100.

Sequence analysis of 56 kb from the genome of the bacterium Mycoplasma pneumoniae comprising the dnaA region, the atp operon and a cluster of ribosomal protein genes

To sequence the entire 800 kilobase pair genome of the bacterium Mycoplasma pneumoniae, a plasmid library was established with contained the majority of the EcoR1 fragments from M.pneumoniae. The EcoR1 fragments were subcloned from an ordered cosmid library comprising the complete M.pneumoniae genome. Individual plasmid clones were sequenced in an ordered fashion mainly by primer walking. We report here the initial results from the sequence analysis of -56 kb comprising the dnaA region as a potential origin of replication, the ATPase operon and a region coding for a cluster of ribosomal protein genes. The data were compared with the corresponding genes/operons from Bacillus subtilis, Escherichia coli, Mycoplasma capricolum and Mycoplasma gallisepticum.

The proline-rich P65 protein of Mycoplasma pneumoniae is a component of the Triton X-100-insoluble fraction and exhibits size polymorphism in the strains M129 and FH

Previously, we described the identification of a novel Mycoplasma pneumoniae M129 protein, named P65 because of its apparent molecular mass of 65 kDa estimated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (T. Proft and R. Herrmann, Mol. Microbiol. 13:337-348, 1994). DNA sequence analysis of the P65 open reading frame (orfp65), however, revealed an ORF encoding a protein with a molecular weight of 47,034. This discrepancy can be explained by the unusual amino acid composition of this protein. According to the deduced amino acid sequence, the N-terminal half of P65 contains several penta- and hexapeptides (DPNAY and DPNQAY) forming a proline-rich acidic domain. Secondary-structure predictions indicated beta-sheets and turns within that region, suggesting an extended and rigid conformation. Near the C terminus of P65 the tripeptide Arg-Gly-Asp (RGD) was found. This motif is known to play an important role in binding of extracellular matrix proteins to integrins. P65 could be located exclusively to the Triton X-100-insoluble cell fraction. The results of immunofluorescence microscopy and of immunoadsorption experiments indicated that P65 carries surface-exposed regions. Mild treatment of whole cells with proteases resulted in cleavage of a limited amount of P65 molecules, suggesting either that only a small percentage of P65 molecules are exposed on the surface or that protease cleavage is hampered by a compact protein conformation or by binding of an unknown component to P65. P65 exhibits size polymorphism in M. pneumoniae M129 and FH. This is caused by an intragenetic duplication of a 54-bp sequence within the FH orfp65. As a consequence, the number of DPNAY pentapeptides increased from 9 to 12 repeats in the FH strain.

A spontaneous hemadsorption-negative mutant of Mycoplasma pneumoniae exhibits a truncated adhesin-related 30-kilodalton protein and lacks the cytadherence-accessory protein HMW1

A spontaneous, hemadsorption-negative mutant of Mycoplasma pneumoniae lacks the cytoskeleton-forming HMW1 protein and exhibits a truncated adhesin-related 30-kDa protein. Genetic analyses revealed deletion of one nucleotide in the hmw1 gene and loss of eight repeated sequences comprising 144 nucleotides in the gene for the adhesin-related 30-kDa protein.