Systematic Sequencing of the 180 Kilobase Region of the Bacillus Subtilis Chromosome Containing the Replication Origin

We have determined a 180 kb contiguous sequence in the replicationorigin region of the Bacillus subtilis chromosome. Open readingframes (ORF) in this region were unambiguously identified fromthe determined sequence, using criteria characteristic for theB. subtilis gene structure, i.e., starting with an ATG, GTGor TTG codon preceded by sequences complementary to the 3' endof the 16S rRNA. Four rRNA gene sets, 7 individual tRNA genesand 1 scRNA gene were identified, occupying 20 kb in total.In the remaining 160 kb region, 158 ORFs were identified, suggestingthat 1 ORF is coded on average by 1 kb of DNA of the B. subtilisgenome. Among the 158 ORFs, the functions of 48 ORFs were assignedand those of 11 ORFs are suggested through significant similaritiesto known proteins present in data banks. However, the functionsof more than half of the ORFs (63%) remain to be determined.     

Sequence Analysis of a 25-kb Segment in the 17{degrees}-19{degrees} Region of the Bacillus subtilis Chromosome Containing ada Locus

As a part of the Bacillus subtilis genome sequencing project,we have determined a 25-kb sequence covering the 17°–19°region. This region contains 26 complete open reading frames(ORFs) including the alkA and adaA/B operon, which encode genesfor adaptive response to DNA alkylation. A homology search forthe newly identified 21 ORFs revealed that 4 of them exhibita significant similarity to known proteins, e.g., methicillin-resistantStaphylococcus aureus (MRSA) protein homolog, proteins involvedin chloramphenicol resistance, glucosamine synthase and an ABCtransporter protein. The remaining 17 ORFs did not show anysignificant sequence similarities to known gene products inthe database.     

Sequence Analysis of the groESL-cotA Region of the Bacillus subtilis Genome, Containing the Restriction/Modification System Genes

We have determined a 35-kb sequence of the groESL-gutR-cotA(45°–52°) region of the Bacillus subtilis genome.In addition to the groESL, gutRB and cotA genes reported previously,we have newly identified 24 ORFs including gutA and fruC genes,encoding glucitol permease and fructokinase, respectively. Theinherent restriction/modification system genes, hsdMR and hsdMM,were mapped between groESL and gutRB, and we have identifiedtwo open reading frames (ORFs) encoding 5-methylcytosine formingDNA methyl transferase and an operon probably encoding a restrictionenzyme complex. The unusual genome structure of few ORFs andlower GC content around the restriction/modification genes stronglysuggests that the region originated from a bacteriophage integratedduring evolution.     

An 8 kb Nucleotide Sequence at the 3' Flanking Region of the sspC Gene (184{degrees}) on the Bacillus subtilis 168 Chromosome Containing an Intein and an Intron

As part of the Bacillus subtilis genome sequencing project,we determined the complete nucleotide sequence of an 8000-bpfragment downstream of the sspC gene (184°) of the B. subtilis168 chromosome. The sequence analysis shows that the sspC geneis located inside of the SPß region, which differsfrom the current genetic map of B. subtilis 168. This regioncontains 12 putative ORFs (yojQ through yojZ and sspC). A homologysearch for the deduced products of the ORFs shows signi.cantsimilarities to enzymes involved in deoxyribonucleotide metabolism:ribonucleotide reductase (Nrd) E, NrdF, thioredoxinand dUTPase.Interestingly, this DNA fragment includes two split genes, yojPcontaining conserved motifs of an intein and yojQ and yojS withan 808-bp intervening sequence for a putative intron structure.In addition, the yojR gene includes a putative new DNA replicationterminator.     

Sequence Analysis of the Bacillus subtilis 168 Chromosome Region between the sspC and odhA Loci (184{degrees}-180{degrees})

The nucleotide sequence of 45,389 bp in the 184°-;180°region of the Bacillus subtilis chromosome, containing the cgecluster, which is controlled by the sporulation regulatory proteinGerE, was determined. Fifty-four putative ORFs with putativeribosome-binding sites were recognized. Seven of them correspondto previously characterized genes: cgeB, cgeA, cgeC, cgeD, cgeE,ctpA, and odhA. The deduced products of 25 ORFs were found todisplay significant similarities to proteins in the data banks.We have identified genes involved in detoxification, cell walls,and in the metabolism of biotins, purines, fatty acids, carbohydratesand amino acids. The remaining 22 ORFs showed no similarityto known proteins. Both an attachment site of the SPßprophage and 2 new putative DNA replication terminators wereidentified in this region.     

Cloning and Sequencing of a 36-kb Region of the Bacillus subtilis Genome between the gnt and iol Operons

Within the framework of an international project for the sequencingof the entire Bacillus subtilis genome, a 36-kb chromosome segment,which covers the region between the gnt and iol operons, hasbeen cloned and sequenced. This region (36447 bp) contains 33complete open reading frames (ORFs; genes) including the fourgnt genes and one partial gene. A homology search for the productsof the 33 complete ORFs revealed significant homology to knownproteins in 16 of them such as tetracycline resistance protein(Clostridium perfringens), asparagine synthetase (Arabidopsisthaliana), aldehyde dehydrogenase (Pseudomonas oleovorans),2,5-dichloro-2,5-cyclohexadiene-1,4-diol dehydrogenase (P. paucimobilis),heat shock protein HtpG (Escherichia coli), galactose-protonsymporter (E. coli), auxin-induced protein (common tobacco),glucitol operon repressor (E. coli) and methylmalonate-semialdehydedehydrogenase (P. aeruginosa). Unlike the regions we sequencedso far, this region contained two short sequence multiplications:one was a tandem sequence duplication (409 and 410 bp), andthe other a triplication consisting of two highly conserved118-bp tandem sequences preceded by a less conserved similarsequence (129 bp). The reasons for the presence of these sequencemultiplications in the gnt to iol region were deduced.     

Cloning and Sequencing of a 27.8-kb Nucleotide Sequence of the 79{degrees}-81{degrees}. Region of the Bacillus subtilis Genome Containing the sspE Locus

The nucleotide sequence of a 27830-bp DNA segment in the 79°–81°.region of the Bacillus subtilis genome has been determined.This region contains 29 complete ORFs including the sspE gene,which encodes a small acid-soluble spore protein gamma and locateson the one side terminal of our assigned region. A homologysearch for the products deduced from the 29 ORFs revealed thatnine of them exhibit significant similarity to known proteins,e.g. proteins involved in an iron uptake system, a multidrugresistance protein, a chloramphenicol resistance protein, epoxidehydrolase, adenine glycosylase, and a glucose-1-dehydrogenasehomolog.     

Cloning and Sequencing of a 23-kb Region of the Bacillus subtilis Genome between the iol and hut Operons

Within the framework of an international project for the sequencingof the entire Bacillus subtilis genome, a 23-kb chromosomalsegment, which covers the region between the iol and hut operons,has been cloned and sequenced, creating a 99-kb contig fromthe gnt operon to the wapA locus. This region (23351 bp) contains25 complete open reading frames (ORFs; genes) including deoR,dra, nupC and pdp and two partial ones. The region (5140 bp)containing these four genes, being also sequenced by H. H. Saxildet al., was sequenced by subjecting a long polymerase chainreaction product to random sequencing using phage M13mp19. However,we could detect no conflict, between two independently determinedsequences, which could be attributed to our sequencing method.A homology search for the 24 newly identified gene productsrevealed significant homology to known proteins in 14 of them.It was notable that three proteins, encoded by the successivegenes (yxeMNO), exhibited meaningful homology to the E. coliGlnHPQ products constituting a periplasmic ATP-dependent transportsystem for glutamine.     

Identification of the Coding Region of Saccharomyces cerevisiae Chromosome VI Using the Computer Program GenMark

We searched the nucleotide sequence of budding yeast Saccharomycescerevisiae chromosome VI (270 kb) for candidate coding regions,using the computer program GenMark. One hundred and twenty-nineputative genes were identified, which is almost the same asthe number of ORFs on this chromosome. Nineteen new putativegenes were identified through the GenMark analysis. Most largeORFs were also correctly identified (87% of the predicted putativegenes identified by the GenMark (110 of 127) matched the reportedORFs). The new coding regions were mostly small but they weredistinguished from the more than 2000 ORFs identified by Genetyx.GenMark did not predict 17 ORFs that were over 300 bp long.As these ORFs include known genes, their sequence context maydiffer somewhat from that of typical yeast genes. These analysesrevealed the high potential of GenMark to identify putativegenes from numerous short ORFs and will produce informationon the likelihood of their being actual genes.     

Sequencing of three lambda clones from the genome of alkaliphilic Bacillus sp. strain C-125

The nucleotide sequences of three independent fragments (designated no. 3, 4, and 9; each 15–20 kb in size) of the genome of alkaliphilic Bacillus sp. C-125 cloned in a λ phage vector have been determined. Thirteen putative open reading frames (ORFs) were identified in sequenced fragment no. 3 and 11 ORFs were identified in no. 4. Twenty ORFs were also identified in fragment no. 9. All putative ORFs were analyzed in comparison with the BSORF database and non-redundant protein databases. The functions of 5 ORFs in fragment no. 3 and 3 ORFs in fragment no. 4 were suggested by their significant similarities to known proteins in the database. Among the 20 ORFs in fragment no. 9, the functions of 11 ORFs were similarly suggested. Most of the annotated ORFs in the DNA fragments of the genome of alkaliphilic Bacillus sp. C-125 were conserved in the Bacillus subtilis genome. The organization of ORFs in the genome of strain C-125 was found to differ from the order of genes in the chromosome of B. subtilis, although some gene clusters (ydh, yqi, yer, and yts) were conserved as operon units the same as in B. subtilis. Received: April 17, 1998 / Accepted: June 23, 1998     

Phylogenetic analysis of Bacillus P450 monooxygenases and evaluation of their activity towards steroids

Cytochrome P450 (P450) open reading frames (ORFs) identified in genome sequences of Bacillus species are potential resources for new oxidation biocatalysts. Phylogenetic analysis of 29 Bacillus P450 ORFs revealed that the P450s consist of a limited number of P450 families, CYP102, CYP106, CYP107, CYP109, CYP134, CYP152, and CYP197. Previously, we identified the catalytic activities of three P450s of Bacillus subtilis towards steroids by rapid substrate screening using Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR/MS). Here, we further applied this method to evaluate the activity of Bacillus cereus P450s towards steroids. Five P450 genes were cloned from B. cereus ATCC 10987 based on its genomic sequence and were expressed in Escherichia coli. These P450s were reacted with a mixture of 30 compounds that mainly included steroids, and the reaction mixtures were analyzed using FT-ICR/MS. We found that BCE_2659 (CYP106) catalyzed the monooxygenation of methyltestosterone, progesterone, 11-ketoprogesterone, medroxyprogesterone acetate, and chlormadinone acetate. BCE_2654 (CYP107) monooxygenated testosterone enanthate, and BCE_3250 (CYP109) monooxygenated testosterone and compactin. Based on the phylogenetic relationship and the known substrate specificities including ones identified in this study, we discuss the catalytic potential of Bacillus P450s towards steroids.     

Compilation of All Genes Encoding Two-component Phosphotransfer Signal Transducers in the Genome of Escherichia coli

Bacteria have devised sophisticated His-Asp phosphorelay signalingsystems for eliciting a variety of adaptive responses to theirenvironment, which are generally referred to as the "two-componentregulatory system." The widespread occurrence of the His-Aspphosphorelay signaling in both prokaryotes and eukaryotes impliesthat it is a powerful device for a wide variety of adaptiveresponses of cells to their environment. The two-component signaltransducers contain one or more of three common and characteristicphosphotransfer signaling domains, named the "transmitter, receiver,and histidine-containing phosphotransfer (HPt) domains." Therecently determined entire genomic sequence of Escherichia coliallowed us to compile systematically a complete list of genesencoding such two-component signal transduction proteins. Theresults of such an effort, made in this study, revealed thatat least 62 open reading frames(ORFs) were identified as putativemembers of the two-component signaltransducers in this singlespecies. Among them, 32 were identified as response regulatorand 23 were identified as orthodox sensory kinases. In addition,E. coli has five hybrid sensory kinases. The precise locationof each ORF was mapped on a physical map of the entire E. coligenome. All of these ORFs were then compiled and annotated extensively.     

Sequence Analysis of a 32-kb Region Including the Major Ribosomal Protein Gene Clusters from Alkaliphilic Bacillus sp. Strain C-125

Forty-one open reading frames (ORFs) were identified in a 32-kb DNA fragment of alkaliphilic Bacillus sp. C-125. A similarity search using the BSORF database found 37 ORFs with significant sequence similarity to B. subtilis RNA polymerase subunits, elongation factor G, elongation factor Tu, and ribosomal proteins. Each ORF product showed more than 70% identity to those of B. subtilis. Gene organization in the region of str, S10, spc, and the α cluster was highly conserved among three strains, C-125, B. subtilis, and B. stearothermophilus.     

Nucleotide Sequence and Features of the Bacillus licheniformis gnt Operon

Bacillus licheniformis was able to utilize gluconate as thesole carbon source as efficiently as Bacillus subtilis did.Southern analysis indicated that B. licheniformis likely possessesonly one gnt determinant. The nucleotide sequence (6278 bp)of the B. licheniformis DNA containing the gnt operon was determined,revealing the five complete open reading frames (ORF; genes).The putative product of the first gene, oug, did not show anysignificant homology to known proteins, but those of the secondto fifth genes exhibited striking homology to the gntRKPZ genesof B. subtilis, respectively, indicating that they are the correspondinggnt genes of B. licheniformis. Not only is the organizationof the gnt genes of these two Bacilli highly conserved, butso are the cis regulatory elements of their gnt operon. Sequenceanalysis of the upstream regions of these two gnt operons impliedthat a chromosome rearrangement in B. subtilis might have occurredimmediately upstream of the gnt operon during evolution, causingit to diverge from a common ancestor into B. licheniformis andB. subtilis.     

Exploring the Biosynthesis of Unsaturated Fatty Acids in Bacillus cereus ATCC 14579 and Functional Characterization of Novel Acyl-Lipid Desaturases

At low temperatures, Bacillus cereus synthesizes large amounts of unsaturated fatty acids (UFAs) with double bonds in positions Δ5 and Δ10, as well as Δ5,10 diunsaturated fatty acids. Through sequence homology searches, we identified two open reading frames (ORFs) encoding a putative Δ5 desaturase and a fatty acid acyl-lipid desaturase in the B. cereus ATCC 14579 genome, and these were named BC2983 and BC0400, respectively. Functional characterization of ORFs BC2983 and BC0400 by means of heterologous expression in Bacillus subtilis confirmed that they both encode acyl-lipid desaturases that use phospholipids as the substrates and have Δ5 and Δ10 desaturase activities. Thus, these ORFs were correspondingly named desA (Δ5 desaturase) and desB (Δ10 desaturase). We established that DesA utilizes ferredoxin and flavodoxins (Flds) as electron donors for the desaturation reaction, while DesB preferably employs Flds. In addition, increased amounts of UFAs were found when B. subtilis expressing B. cereus desaturases was subjected to a cold shock treatment, indicating that the activity or the expression of these enzymes is upregulated in response to a decrease in growth temperature. This represents the first work reporting the functional characterization of fatty acid desaturases from B. cereus.     

Hundreds of Out-of-Frame Remodeled Gene Families in the Escherichia coli Pangenome

All genomes include gene families with very limited taxonomic distributions that potentially represent new genes and innovations in protein-coding sequence, raising questions on the origins of such genes. Some of these genes are hypothesized to have formed de novo, from noncoding sequences, and recent work has begun to elucidate the processes by which de novo gene formation can occur. A special case of de novo gene formation, overprinting, describes the origin of new genes from noncoding alternative reading frames of existing open reading frames (ORFs). We argue that additionally, out-of-frame gene fission/fusion events of alternative reading frames of ORFs and out-of-frame lateral gene transfers could contribute to the origin of new gene families. To demonstrate this, we developed an original pattern-search in sequence similarity networks, enhancing the use of these graphs, commonly used to detect in-frame remodeled genes. We applied this approach to gene families in 524 complete genomes of Escherichia coli. We identified 767 gene families whose evolutionary history likely included at least one out-of-frame remodeling event. These genes with out-of-frame components represent ∼2.5% of all genes in the E. coli pangenome, suggesting that alternative reading frames of existing ORFs can contribute to a significant proportion of de novo genes in bacteria.     

A phylogenomic investigation into the origin of metazoa

The evolution of multicellular animals (Metazoa) from theirunicellular ancestors was a key transition that was accompaniedby the emergence and diversification of gene families associatedwith multicellularity. To clarify the timing and order of specificevents in this transition, we conducted expressed sequence tagsurveys on 4 putative protistan relatives of Metazoa includingthe choanoflagellate Monosiga ovata, the ichthyosporeans Sphaeroformaarctica and Amoebidium parasiticum, and the amoeba Capsasporaowczarzaki, and 2 members of Amoebozoa, Acanthamoeba castellaniiand Mastigamoeba balamuthi. We find that homologs of genes involvedin metazoan multicellularity exist in several of these unicellularorganisms, including 1 encoding a membrane-associated guanylatekinase with an inverted arrangement of protein-protein interactiondomains (MAGI) in Capsaspora. In Metazoa, MAGI regulates tightjunctions involved in cell-cell communication. By phylogenomicanalyses of genes encoded in nuclear and mitochondrial genomes,we show that the choanoflagellates are the closest relativesof the Metazoa, followed by the Capsaspora and Ichthyosporealineages, although the branching order between the latter 2groups remains unclear. Understanding the function of "metazoan-specific"proteins we have identified in these protists will clarify theevolutionary steps that led to the emergence of the Metazoa.     

Complete set of ORF clones of Escherichia coli ASKA library (A Complete Set of E. coli K-12 ORF Archive): Unique Resources for Biological Research

Based on the genomic sequence data of Escherichia coli K-12strain, we have constructed a complete set of cloned individualgenes encoding Histidine-tagged proteins with or without GFPfused for functional genomic analysis. Each clone encodes aprotein of predicted ORF attached by Histidines and seven spaceramino acids at the N-terminal end, and five spacer amino acidsand GFP at the C-terminal end. SfiI restriction sites are generatedat both the N- and C-terminal boundaries of ORF upon cloning,which enables easy transfer of ORF to other vector systems bycutting with SfiI. Expression of cloned ORF is under the controlof an IPTG-inducible promoter, which is strictly repressed bylacI^q repressor gene product. The set of cloned ORFs describedhere should provide unique resources for systematic functionalgenomic approaches including (i) construction of DNA microarray,(ii) production and purification of proteins, (iii) analysisof protein localization by monitoring GFP fluorescence and (iv)analysis of protein–protein interaction.