Molecular cloning with bifunctional plasmid vectors in Bacillus subtilis

Summary Cloning in Escherichia coli and Bacillus subtilis was carried out using the bifunctional plasmid pDH5060. B. subtilis chromosomal DNA and pDH5060 DNA were digested with either BamHI or SalI, then annealed, ligated, and transformed into E. coli SK2267. Transformants containing sequences ligated into the BamHI or SalI sites in the Tc^r gene of pDH5060 were selected directly using a modification of the fusaric acid technique. The BamHI and SalI clone banks contain about 250 and 140 B. subtilis fragments, respectively, with an average insert size of 8–9 Kbp in the BamHI and 4–5 Kbp in the SalI bank. The inserts ranged in size from 0.3 Kbp to greater than 20 Kbp. The vector used here therefore accepts inserts which are significantly larger than previously reported for other B. subtilis cloning systems. All individual cloned B. subtilis sequences examined were stably propagated in E. coli SK2267. Eight of eighteen B. subtilis auxotrophic markers tested (aroG, gltA, glyB, ilvA, metC, purA, pyrD, and thrA) were transformed to prototrophy with BamHI or SalI clone bank DNA. All or part of the hybrid plasmid DNA recombined at the sites of homology in the chromosome of these Rec⁺ recipients. Loss of sequences from hybrid plasmids was not prevented in a r ^- m ^- recE4 recipient strain of B. subtilis. Although the recE4 background prevented recombination between homologous chromosomal DNA, a variety of cloned fragments were shown to be unstable and undergo deletions of both insert and plasmid sequences. In addition, B. subtilis sequences propagated in E. coli transformed B. subtilis recE4 recipients with a 500-1,000-fold reduced efficiency.     

Molecular analysis of some genes from plasmid p19 of the Bacillus subtilis 19 soil strain involved in conjugation

Two fragments of conjugative plasmid p19 (95 kb) from the Bacillus subtilis 19 soil strain were cloned and sequenced; these fragments carry genes, products of which are indispensable for the conjugative transfer. One of the fragments 4518 bp in size carries five open reading frames and their fragments (ORF1–ORF5). The protein corresponding to ORF4 is homologous to proteins from the family VirD4. Inactivation of ORF4 and ORF1 by insertional mutagenesis caused a three-to-fivefold decrease in the frequency of plasmid p19 conjugative transfer. Another 2932-bp fragment of p19 was shown to possess a rep region homologous to the rep region of plasmid pBS72 from the B. subtilis 72 soil strain and a novel ORF (ORF6); the protein corresponding to this ORF contains the HTH motif typical for DNA-binding proteins.     

Genome of Bacillus subtilis Bacteriophage SPP1: Structure and Nucleotide Sequence of pac, the Origin of DNA Packaging

The DNA of Bacillus subtilis bacteriophage SPP1 is terminally redundant and partially circularly permuted. To explain these parameters, we followed the Streisinger-Botstein models of phage maturation and assumed that packaging of SPP1 DNA begins at a unique genomic site (“pac”) and proceeds sequentially from there. We describe the sequence of about 1,000 nucleotides surrounding pac. This together with size determinations of small, pac-terminated restriction fragments has revealed heterogeneity of the natural pac ends of SPP1 DNA. Such ends fell in each DNA strand into a region of five to seven nucleotides. However, within this range more than 50% of all molecules terminated with defined cytosines on both strands, generating a 3′ protruding terminus. The nucleotide sequence of the DNA segment surrounding pac did not reveal any features which would distinguish this region.     

Metabolic Engineering of Carotenoid Biosynthesis in Escherichia coli by Ordered Gene Assembly in Bacillus subtilis

We attempted to optimize the production of zeaxanthin in Escherichia coli by reordering five biosynthetic genes in the natural carotenoid cluster of Pantoea ananatis. Newly designed operons for zeaxanthin production were constructed by the ordered gene assembly in Bacillus subtilis (OGAB) method, which can assemble multiple genes in one step using an intrinsic B. subtilis plasmid transformation system. The highest level of production of zeaxanthin in E. coli (820 μg/g [dry weight]) was observed in the transformant with a plasmid in which the gene order corresponds to the order of the zeaxanthin metabolic pathway (crtE-crtB-crtI-crtY-crtZ), among a series of plasmids with circularly permuted gene orders. Although two of five operons using intrinsic zeaxanthin promoters failed to assemble in B. subtilis, the full set of operons was obtained by repressing operon expression during OGAB assembly with a p_R promoter-cI repressor system. This result suggests that repressing the expression of foreign genes in B. subtilis is important for their assembly by the OGAB method. For all tested operons, the abundance of mRNA decreased monotonically with the increasing distance of the gene from the promoter in E. coli, and this may influence the yield of zeaxanthin. Our results suggest that rearrangement of biosynthetic genes in the order of the metabolic pathway by the OGAB method could be a useful approach for metabolic engineering.     

In Vivo Random Mutagenesis of Bacillus subtilis by Use of TnYLB-1, a mariner-Based Transposon

This report describes the construction and characterization of a mariner-based transposon system designed to be used in Bacillus subtilis, but potentially applicable to other gram-positive bacteria. Two pUC19-derived plasmids were created that contain the mariner-Himar1 transposase gene, modified for expression in B. subtilis, under the control of either σ^A- or σ^B-dependent promoters. Both plasmids also contain a transposable element (TnYLB-1) consisting of a Kan^r cassette bracketed by the Himar1-recognized inverse terminal repeats, as well as the temperature-sensitive replicon and Erm^r gene of pE194ts. TnYLB-1 transposes into the B. subtilis chromosome with high frequency (10⁻²) from either plasmid. Southern hybridization analyses of 15 transposants and sequence analyses of the insertion sites of 10 of these are consistent with random transposition, requiring only a “TA” dinucleotide as the essential target in the recipient DNA. Two hundred transposants screened for sporulation proficiency and auxotrophy yielded five Spo⁻ clones, three with insertions in known sporulation genes (kinA, spoVT, and yqfD) and two in genes (ybaN and yubB) with unknown functions. Two auxotrophic mutants were identified among the 200 transposants, one with an insertion in lysA and another in a gene (yjzB) whose function is unknown.     

Combinatorial Engineering of 1-Deoxy-D-Xylulose 5-Phosphate Pathway Using Cross-Lapping In Vitro Assembly (CLIVA) Method

The ability to assemble multiple fragments of DNA into a plasmid in a single step is invaluable to studies in metabolic engineering and synthetic biology. Using phosphorothioate chemistry for high efficiency and site specific cleavage of sequences, a novel ligase independent cloning method (cross-lapping in vitro assembly, CLIVA) was systematically and rationally optimized in E. coli. A series of 16 constructs combinatorially expressing genes encoding enzymes in the 1-deoxy-D-xylulose 5-phosphate (DXP) pathway were assembled using multiple DNA modules. A plasmid (21.6 kb) containing 16 pathway genes, was successfully assembled from 7 modules with high efficiency (2.0 x 10³ cfu/ µg input DNA) within 2 days. Overexpressions of these constructs revealed the unanticipated inhibitory effects of certain combinations of genes on the production of amorphadiene. Interestingly, the inhibitory effects were correlated to the increase in the accumulation of intracellular methylerythritol cyclodiphosphate (MEC), an intermediate metabolite in the DXP pathway. The overexpression of the iron sulfur cluster operon was found to modestly increase the production of amorphadiene. This study demonstrated the utility of CLIVA in the assembly of multiple fragments of DNA into a plasmid which enabled the rapid exploration of biological pathways.     

Homology of pCS1 Plasmid Sequences with Chromosomal DNA in Clavibacter michiganense subsp. sepedonicum: Evidence for the Presence of a Repeated Sequence and Plasmid Integration

Restriction fragments of pCS1, a 50.6-kilobase (kb) plasmid present in many strains of Clavibacter michiganense subsp. sepedonicum (“Corynebacterium sepedonicum”), have been cloned in an M13mp11 phage vector. Radiolabeled forms of these cloned fragments have been used as Southern hybridization probes for the presence of plasmid sequences in chromosomal DNA of this organism. These studies have shown that all tested strains lacking the covalently closed circular form of pCS1 contain the plasmid in integrated form. In each case the site of integration exists on a single plasmid restriction fragment with a size of 5.1 kb. Southern hybridizations with these probes have also revealed the existence of a major repeated sequence in C. michiganense subsp. sepedonicum. Hybridizations of chromosomal DNA with deletion subclones of a 2.9-kb plasmid fragment containing the repeated sequence indicate that the size of the repeated sequence is approximately 1.3 kb. One of the copies of the repeated sequence is on the plasmid fragment containing the site of integration.     

Conserved Plasmid Hydrogen-Uptake (hup)-Specific Sequences within Hup+Rhizobium leguminosarum Strains

Thirteen Rhizobium leguminosarum strains previously reported as H₂-uptake hydrogenase positive (Hup⁺) or negative (Hup⁻) were analyzed for the presence and conservation of DNA sequences homologous to cloned Bradyrhizobium japonicum hup-specific DNA from cosmid pHU1 (M. A. Cantrell, R. A. Haugland, and H. J. Evans, Proc. Natl. Acad. Sci. USA 80:181-185, 1983). The Hup phenotype of these strains was reexamined by determining hydrogenase activity induced in bacteroids from pea nodules. Five strains, including H₂ oxidation-ATP synthesis-coupled and -uncoupled strains, induced significant rates of H₂-uptake hydrogenase activity and contained DNA sequences homologous to three probe DNA fragments (5.9-kilobase [kb] HindIII, 2.9-kb EcoRI, and 5.0-kb EcoRI) from pHU1. The pattern of genomic DNA HindIII and EcoRI fragments with significant homology to each of the three probes was identical in all five strains regardless of the H₂-dependent ATP generation trait. The restriction fragments containing the homology totalled about 22 kb of DNA common to the five strains. In all instances the putative hup sequences were located on a plasmid that also contained nif genes. The molecular sizes of the identified hup-sym plasmids ranged between 184 and 212 megadaltons. No common DNA sequences homologous to B. japonicum hup DNA were found in genomic DNA from any of the eight remaining strains showing no significant hydrogenase activity in pea bacteroids. These results suggest that the identified DNA region contains genes essential for hydrogenase activity in R. leguminosarum and that its organization is highly conserved within Hup⁺ strains in this symbiotic species.     

Construction of a recombinant plasmid composed of B. subtilis leucine genes and a B. subtilis (natto) plasmid: its use as cloning vehicle in B. subtilis 168.

Summary Recombinant plasmids composed of Bacillus subtilis 168 leucine genes and a B. subtilis (natto) plasmid have been constructed in a recombination deficient (recE4) mutant of Bacillus subtilis 168. The process involved EcoRI fragmentation and ligation of a B. subtilis (natto) plasmid and a composite plasmid RSF2124-B · leu in which B. subtilis 168 leucine genes are linked to the R-factor RSF2124. A constructed plasmid (pLS102) was found to be composed of an EcoRI fragment derived from the vector plasmid and two tandemly repeated EcoRI fragments carrying the leucine genes. A derivative plasmid (pLS101 or pLS103) consisting of one molecule each of the EcoRI fragments was obtained by in vivo intramolecular recombination between the repeated leucine gene fragments in pLS102. pLS103 was cleaved once with BamNI, SmaI and HpaI. Insertion of foreign DNA (Escherichia coli plasmid pBR322) into the BamNI site inactivated leuA but not the leuC function which thus can serve as selective marker if the plasmid is used as vector in molecular cloning. The penicillin resistance carried in pBR322 was not functionally expressed in B. subtilis cells. By partial digestion of pLS103 with HindIII followed by ligation with T4-induced ligase, pLS107 was obtained which contained only one EcoRI site. However, insertion of exogenous DNA (pBR322) into this EcoRI site inactivated both leuA and leuC functions.     

Whole-Genome Sequencing and Comparative Genome Analysis of Bacillus subtilis Strains Isolated from Non-Salted Fermented Soybean Foods

Bacillus subtilis is the main component in the fermentation of soybeans. To investigate the genetics of the soybean-fermenting B. subtilis strains and its relationship with the productivity of extracellular poly-γ-glutamic acid (γPGA), we sequenced the whole genome of eight B. subtilis stains isolated from non-salted fermented soybean foods in Southeast Asia. Assembled nucleotide sequences were compared with those of a natto (fermented soybean food) starter strain B. subtilis BEST195 and the laboratory standard strain B. subtilis 168 that is incapable of γPGA production. Detected variants were investigated in terms of insertion sequences, biotin synthesis, production of subtilisin NAT, and regulatory genes for γPGA synthesis, which were related to fermentation process. Comparing genome sequences, we found that the strains that produce γPGA have a deletion in a protein that constitutes the flagellar basal body, and this deletion was not found in the non-producing strains. We further identified diversity in variants of the bio operon, which is responsible for the biotin auxotrophism of the natto starter strains. Phylogenetic analysis using multilocus sequencing typing revealed that the B. subtilis strains isolated from the non-salted fermented soybeans were not clustered together, while the natto-fermenting strains were tightly clustered; this analysis also suggested that the strain isolated from “Tua Nao” of Thailand traces a different evolutionary process from other strains.     

Genome engineering using a synthetic gene circuit in Bacillus subtilis

Genome engineering without leaving foreign DNA behind requires an efficient counter-selectable marker system. Here, we developed a genome engineering method in Bacillus subtilis using a synthetic gene circuit as a counter-selectable marker system. The system contained two repressible promoters (B. subtilis xylA (P_xyl) and spac (P_spac)) and two repressor genes (lacI and xylR). P_xyl-lacI was integrated into the B. subtilis genome with a target gene containing a desired mutation. The xylR and P_spac–chloramphenicol resistant genes (cat) were located on a helper plasmid. In the presence of xylose, repression of XylR by xylose induced LacI expression, the LacIs repressed the P_spac promoter and the cells become chloramphenicol sensitive. Thus, to survive in the presence of chloramphenicol, the cell must delete P_xyl-lacI by recombination between the wild-type and mutated target genes. The recombination leads to mutation of the target gene. The remaining helper plasmid was removed easily under the chloramphenicol absent condition. In this study, we showed base insertion, deletion and point mutation of the B. subtilis genome without leaving any foreign DNA behind. Additionally, we successfully deleted a 2-kb gene (amyE) and a 38-kb operon (ppsABCDE). This method will be useful to construct designer Bacillus strains for various industrial applications.     

Crystal structure and functional insights into uracil-DNA glycosylase inhibition by phage ?29 DNA mimic protein p56

Uracil-DNA glycosylase (UDG) is a key repair enzyme responsible for removing uracil residues from DNA. Interestingly, UDG is the only enzyme known to be inhibited by two different DNA mimic proteins: p56 encoded by the Bacillus subtilis phage ϕ29 and the well-characterized protein Ugi encoded by the B. subtilis phage PBS1/PBS2. Atomic-resolution crystal structures of the B. subtilis UDG both free and in complex with p56, combined with site-directed mutagenesis analysis, allowed us to identify the key amino acid residues required for enzyme activity, DNA binding and complex formation. An important requirement for complex formation is the recognition carried out by p56 of the protruding Phe191 residue from B. subtilis UDG, whose side-chain is inserted into the DNA minor groove to replace the flipped-out uracil. A comparative analysis of both p56 and Ugi inhibitors enabled us to identify their common and distinctive features. Thereby, our results provide an insight into how two DNA mimic proteins with different structural and biochemical properties are able to specifically block the DNA-binding domain of the same enzyme.     

Evolutionary Capture of Viral and Plasmid DNA by Yeast Nuclear Chromosomes

A 10-kb region of the nuclear genome of the yeast Vanderwaltozyma polyspora contains an unusual cluster of five pseudogenes homologous to five different genes from yeast killer viruses, killer plasmids, the 2μm plasmid, and a Penicillium virus. By further database searches, we show that this phenomenon is not unique to V. polyspora but that about 40% of the sequenced genomes of Saccharomycotina species contain integrated copies of genes from DNA plasmids or RNA viruses. We propose the name NUPAVs (nuclear sequences of plasmid and viral origin) for these objects, by analogy to NUMTs (nuclear copies of mitochondrial DNA) and NUPTs (nuclear copies of plastid DNA, in plants) of organellar origin. Although most of the NUPAVs are pseudogenes, one intact and active gene that was formed in this way is the KHS1 chromosomal killer locus of Saccharomyces cerevisiae. We show that KHS1 is a NUPAV related to M2 killer virus double-stranded RNA. Many NUPAVs are located beside tRNA genes, and some contain sequences from a mixture of different extrachromosomal sources. We propose that NUPAVs are sequences that were captured by the nuclear genome during the repair of double-strand breaks that occurred during evolution and that some of their properties may be explained by repeated breakage at fragile chromosomal sites.It is well known that the nuclear genomes of most eukaryotes contain integrated fragments of organellar DNA called NUMTs (nuclear copies of mitochondrial DNA) and NUPTs (nuclear copies of plastid DNA, in plants) (26, 29, 44, 45, 57). These fragments are usually pseudogenes, although some NUMTs and NUPTs have become incorporated into functional nuclear genes (38). The NUMTs present in the nuclear genomes of Saccharomycotina yeast species were recently analyzed by Sacerdot et al. (48).In addition to their mitochondrial genomes, yeast species contain a variety of other extranuclear DNA and RNA elements, including viruses and plasmids. These extrachromosomal elements are usually considered to be autonomous entities that do not interact with nuclear DNA. When our laboratory sequenced the genome of the yeast Vanderwaltozyma polyspora (synonym: Kluyveromyces polysporus) (49), we were therefore surprised to find the genomic region we describe here, which contains integrated fragments of several plasmid- and virus-like sequences. We propose that this region was formed by the capture of plasmid and viral sequences by the same mechanism that captures mitochondrial DNA to form NUMTs (43, 65). In a literature search, we could find only one previous report of a similar finding: Utatsu et al. (59) reported the sequences of two regions of nuclear DNA from Zygosaccharomyces rouxii that were highly similar to parts of the 2μm-like plasmid pSR1 from that species, but rearranged.Before describing the V. polyspora region, and similar regions found in other species, we will first briefly introduce the extrachromosomal RNA and DNA entities that are known to exist in yeasts. Extrachromosomal nucleic acids are relatively uncommon in yeasts: a broad survey of 1,800 strains from 600 species by Fukuhara (14) found that 196 strains (11%) contained some sort of extrachromosomal entity. Among these, 105 strains had a double-stranded RNA (dsRNA), 28 had a linear dsDNA plasmid, and 53 had a circular DNA plasmid of the 2μm family. These elements typically also have a patchy distribution within a species, being found in some individuals or strains but not in others. For instance, Nakayashiki et al. (37) surveyed 70 “wild” strains of Saccharomyces (mostly S. cerevisiae) for the presence of five extrachromosomal elements (2μm DNA plasmid, L-A and L-BC helper RNA viruses, and W and T RNA entities) and found each element to be present in between 1 and 38 of the strains, with 1 strain even containing all five elements simultaneously.     

Complete genome sequence of the alkaliphilic bacterium Bacillus halodurans and genomic sequence comparison with Bacillus subtilis

The 4 202 353 bp genome of the alkaliphilic bacterium Bacillus halodurans C-125 contains 4066 predicted protein coding sequences (CDSs), 2141 (52.7%) of which have functional assignments, 1182 (29%) of which are conserved CDSs with unknown function and 743 (18.3%) of which have no match to any protein database. Among the total CDSs, 8.8% match sequences of proteins found only in Bacillus subtilis and 66.7% are widely conserved in comparison with the proteins of various organisms, including B.subtilis. The B.halodurans genome contains 112 transposase genes, indicating that transposases have played an important evolutionary role in horizontal gene transfer and also in internal genetic rearrangement in the genome. Strain C-125 lacks some of the necessary genes for competence, such as comS, srfA and rapC, supporting the fact that competence has not been demonstrated experimentally in C-125. There is no paralog of tupA, encoding teichuronopeptide, which contributes to alkaliphily, in the C-125 genome and an ortholog of tupA cannot be found in the B.subtilis genome. Out of 11 σ factors which belong to the extracytoplasmic function family, 10 are unique to B.halodurans, suggesting that they may have a role in the special mechanism of adaptation to an alkaline environment.     

Transposon Insertion Reveals pRM, a Plasmid of Rickettsia monacensis

Until the recent discovery of pRF in Rickettsia felis, the obligate intracellular bacteria of the genus Rickettsia (Rickettsiales: Rickettsiaceae) were thought not to possess plasmids. We describe pRM, a plasmid from Rickettsia monacensis, which was detected by pulsed-field gel electrophoresis and Southern blot analyses of DNA from two independent R. monacensis populations transformed by transposon-mediated insertion of coupled green fluorescent protein and chloramphenicol acetyltransferase marker genes into pRM. Two-dimensional electrophoresis showed that pRM was present in rickettsial cells as circular and linear isomers. The 23,486-nucleotide (31.8% G/C) pRM plasmid was cloned from the transformant populations by chloramphenicol marker rescue of restriction enzyme-digested transformant DNA fragments and PCR using primers derived from sequences of overlapping restriction fragments. The plasmid was sequenced. Based on BLAST searches of the GenBank database, pRM contained 23 predicted genes or pseudogenes and was remarkably similar to the larger pRF plasmid. Two of the 23 genes were unique to pRM and pRF among sequenced rickettsial genomes, and 4 of the genes shared by pRM and pRF were otherwise found only on chromosomes of R. felis or the ancestral group rickettsiae R. bellii and R. canadensis. We obtained pulsed-field gel electrophoresis and Southern blot evidence for a plasmid in R. amblyommii isolate WB-8-2 that contained genes conserved between pRM and pRF. The pRM plasmid may provide a basis for the development of a rickettsial transformation vector.     

Rapid metabolic pathway assembly and modification using serine integrase site-specific recombination

Synthetic biology requires effective methods to assemble DNA parts into devices and to modify these devices once made. Here we demonstrate a convenient rapid procedure for DNA fragment assembly using site-specific recombination by ϕC31 integrase. Using six orthogonal attP/attB recombination site pairs with different overlap sequences, we can assemble up to five DNA fragments in a defined order and insert them into a plasmid vector in a single recombination reaction. ϕC31 integrase-mediated assembly is highly efficient, allowing production of large libraries suitable for combinatorial gene assembly strategies. The resultant assemblies contain arrays of DNA cassettes separated by recombination sites, which can be used to manipulate the assembly by further recombination. We illustrate the utility of these procedures to (i) assemble functional metabolic pathways containing three, four or five genes; (ii) optimize productivity of two model metabolic pathways by combinatorial assembly with randomization of gene order or ribosome binding site strength; and (iii) modify an assembled metabolic pathway by gene replacement or addition.     

Degradation of nanoRNA is performed by multiple redundant RNases in Bacillus subtilis

Escherichia coli possesses only one essential oligoribonuclease (Orn), an enzyme that can degrade oligoribonucleotides of five residues and shorter in length (nanoRNA). Firmicutes including Bacillus subtilis do not have an Orn homolog. We had previously identified YtqI (NrnA) as functional analog of Orn in B. subtilis. Screening a genomic library from B. subtilis for genes that can complement a conditional orn mutant, we identify here YngD (NrnB) as a second nanoRNase in B. subtilis. Like NrnA, NrnB is a member of the DHH/DHHA1 protein family of phosphoesterases. NrnB degrades nanoRNA 5-mers in vitro similarily to Orn. Low expression levels of NrnB are sufficient for orn complementation. YhaM, a known RNase present in B. subtilis, degrades nanoRNA efficiently in vitro but requires high levels of expression for only partial complementation of the orn^– strain. A triple mutant (nrnA^–, nrnB^–, yhaM^–) in B. subtilis is viable and shows almost no impairment in growth. Lastly, RNase J1 seems also to have some 5′-to-3′ exoribonuclease activity on nanoRNA and thus can potentially finish degradation of RNA. We conclude that, unlike in E. coli, degradation of nanoRNA is performed in a redundant fashion in B. subtilis.     

Thermoadaptation-Directed Enzyme Evolution in an Error-Prone Thermophile Derived from Geobacillus kaustophilus HTA426

Thermostability is an important property of enzymes utilized for practical applications because it allows long-term storage and use as catalysts. In this study, we constructed an error-prone strain of the thermophile Geobacillus kaustophilus HTA426 and investigated thermoadaptation-directed enzyme evolution using the strain. A mutation frequency assay using the antibiotics rifampin and streptomycin revealed that G. kaustophilus had substantially higher mutability than Escherichia coli and Bacillus subtilis. The predominant mutations in G. kaustophilus were A · T→G · C and C · G→T · A transitions, implying that the high mutability of G. kaustophilus was attributable in part to high-temperature-associated DNA damage during growth. Among the genes that may be involved in DNA repair in G. kaustophilus, deletions of the mutSL, mutY, ung, and mfd genes markedly enhanced mutability. These genes were subsequently deleted to construct an error-prone thermophile that showed much higher (700- to 9,000-fold) mutability than the parent strain. The error-prone strain was auxotrophic for uracil owing to the fact that the strain was deficient in the intrinsic pyrF gene. Although the strain harboring Bacillus subtilis pyrF was also essentially auxotrophic, cells became prototrophic after 2 days of culture under uracil starvation, generating B. subtilis PyrF variants with an enhanced half-denaturation temperature of >10°C. These data suggest that this error-prone strain is a promising host for thermoadaptation-directed evolution to generate thermostable variants from thermolabile enzymes.