Intrinsic and extrinsic approaches for detecting genes in a bacterial genome.

The unannotated regions of the Escherichia coli genome DNA sequence from the EcoSeq6 database, totaling 1,278 'intergenic' sequences of the combined length of 359,279 basepairs, were analyzed using computer-assisted methods with the aim of identifying putative unknown genes. The proposed strategy for finding new genes includes two key elements: i) prediction of expressed open reading frames (ORFs) using the GeneMark method based on Markov chain models for coding and non-coding regions of Escherichia coli DNA, and ii) search for protein sequence similarities using programs based on the BLAST algorithm and programs for motif identification. A total of 354 putative expressed ORFs were predicted by GeneMark. Using the BLASTX and TBLASTN programs, it was shown that 208 ORFs located in the unannotated regions of the E. coli chromosome are significantly similar to other protein sequences. Identification of 182 ORFs as probable genes was supported by GeneMark and BLAST, comprising 51.4% of the GeneMark 'hits' and 87.5% of the BLAST 'hits'. 73 putative new genes, comprising 20.6% of the GeneMark predictions, belong to ancient conserved protein families that include both eubacterial and eukaryotic members. This value is close to the overall proportion of highly conserved sequences among eubacterial proteins, indicating that the majority of the putative expressed ORFs that are predicted by GeneMark, but have no significant BLAST hits, nevertheless are likely to be real genes. The majority of the putative genes identified by BLAST search have been described since the release of the EcoSeq6 database, but about 70 genes have not been detected so far. Among these new identifications are genes encoding proteins with a variety of predicted functions including dehydrogenases, kinases, several other metabolic enzymes, ATPases, rRNA methyltransferases, membrane proteins, and different types of regulatory proteins.     

Local de novo assembly of RAD paired-end contigs using short sequencing reads

Despite the power of massively parallel sequencing platforms, a drawback is the short length of the sequence reads produced. We demonstrate that short reads can be locally assembled into longer contigs using paired-end sequencing of restriction-site associated DNA (RAD-PE) fragments. We use this RAD-PE contig approach to identify single nucleotide polymorphisms (SNPs) and determine haplotype structure in threespine stickleback and to sequence E. coli and stickleback genomic DNA with overlapping contigs of several hundred nucleotides. We also demonstrate that adding a circularization step allows the local assembly of contigs up to 5 kilobases (kb) in length. The ease of assembly and accuracy of the individual contigs produced from each RAD site sequence suggests RAD-PE sequencing is a useful way to convert genome-wide short reads into individually-assembled sequences hundreds or thousands of nucleotides long.     

Frameshift mutagenesis by eucaryotic DNA polymerases in vitro

The frequency and specificity of frameshift errors produced during a single round of in vitro DNA synthesis by DNA polymerases-alpha, -beta, and -gamma (pol-alpha, -beta, and -gamma, respectively) have been determined. DNA polymerase-beta is the least accurate enzyme, producing frameshift errors at an average frequency of one error for each 1,000-3,000 nucleotides polymerized, a frequency similar to its average base substitution accuracy. DNA polymerase-alpha is approximately 10-fold more accurate, producing frameshifts at an average frequency of one error for every 10,000-30,000 nucleotides polymerized, a frequency which is about 2- to 6-fold lower than the average pol-alpha base substitution accuracy. DNA polymerase-gamma is highly accurate, producing on the average less than one frameshift error for every 200,000-400,000 nucleotides polymerized. This represents a more than 10-fold higher fidelity than for base substitutions. Among the collection of sequenced frameshifts produced by DNA polymerases-alpha and beta, both common features and distinct specificities are apparent. These specificities suggest a major role for eucaryotic DNA polymerases in modulating frameshift fidelity. Possible mechanisms for production of frameshifts are discussed in relation to the observed biases. One of these models has been experimentally supported using site-directed mutagenesis to change the primary DNA sequence of the template. Alteration of a pol-beta frameshift hotspot sequence TTTT to CTCT reduced the frequency of pol-beta-dependent minus-one-base errors at this site by more than 30-fold, suggesting that more than 97% of the errors at the TTTT run involve a slippage mechanism.     

Strand asymmetry of +1 frameshift mutagenesis at a homopolymeric run by DNA polymerase III holoenzyme of Escherichia coli

We have recently shown that single-base frameshifts were predominant among mutations induced within the rpsL target sequence upon oriC plasmid DNA replication in vitro. We found that the occurrence of +1 frameshifts at a run of 6 residues of dA/dT could be increased proportionally by increasing the concentration of dATP present in the in vitro replication. Using single-stranded circular DNA containing either the coding sequence of the rpsL gene or its complementary sequence, the +1 frameshift mutagenesis by DNA polymerase III holoenzyme of Escherichia coli was extensively examined. A(6) --> A(7) frameshifts occurred 30 to 90 times more frequently during DNA synthesis with the noncoding sequence (dT tract) template than with the coding sequence (dA tract). Excess dATP enhanced the occurrence of +1 frameshifts during DNA synthesis with the dT tract template, but no other dNTPs showed such an effect. In the presence of 0.1 mM dATP, the A(6) --> A(7) mutagenesis with the dT tract template was not inhibited by 1.5 mM dCTP, which is complementary to the residue immediately upstream of the dT tract. These results strongly suggested that the A(6) --> A(7) frameshift mutagenesis possesses an asymmetric strand nature and that slippage errors leading to the +1 frameshift are made during chain elongation within the tract rather than by misincorporation of nucleotides opposite residues next to the tract.     

De novo assembly of the Pseudomonas syringae pv. syringae B728a genome using Illumina/Solexa short sequence reads

Illumina's Genome Analyzer generates ultra-short sequence reads, typically 36 nucleotides in length, and is primarily intended for resequencing. We tested the potential of this technology for de novo sequence assembly on the 6 Mbp genome of Pseudomonas syringae pv. syringae B728a with several freely available assembly software packages. Using an unpaired data set, velvet assembled >96% of the genome into contigs with an N50 length of 8289 nucleotides and an error rate of 0.33%. edena generated smaller contigs (N50 was 4192 nucleotides) and comparable error rates. ssake and vcake yielded shorter contigs with very high error rates. Assembly of paired-end sequence data carrying 400 bp inserts produced longer contigs (N50 up to 15 628 nucleotides), but with increased error rates (0.5%). Contig length and error rate were very sensitive to the choice of parameter values. Noncoding RNA genes were poorly resolved in de novo assemblies, while >90% of the protein-coding genes were assembled with 100% accuracy over their full length. This study demonstrates that, in practice, de novo assembly of 36-nucleotide reads can generate reasonably accurate assemblies from about 40 × deep sequence data sets. These draft assemblies are useful for exploring an organism's proteomic potential, at a very economic low cost.     

Roles of RecA protein in spontaneous mutagenesis in Escherichia coli

To verify the extent of contribution of spontaneous DNA lesions to spontaneous mutagenesis, we have developed a new genetic system to examine simultaneously both forward mutations and recombination events occurring within about 600 base pairs of a transgenic rpsL target sequence located on Escherichia coli chromosome. In a wild-type strain, the recombination events were occurring at a frequency comparable to that of point mutations within the rpsL sequence. When the cells were UV-irradiated, the recombination events were induced much more sharply than point mutations. In a recA null mutant, no recombination event was observed. These data suggest that the blockage of DNA replication, probably caused by spontaneous DNA lesions, occurs often in normally growing E. coli cells and is mainly processed by cellular functions requiring the RecA protein. However, the recA mutant strain showed elevated frequencies of single-base frameshifts and large deletions, implying a novel mutator action of this strain. A similar mutator action of the recA mutant was also observed with a plasmid-based rpsL mutation assay. Therefore, if the recombinogenic problems in DNA replication are not properly processed by the RecA function, these would be a potential source for mutagenesis leading to single-base frameshift and large deletion in E. coli. Furthermore, the single-base frameshifts induced in the recA-deficient cells appeared to be efficiently suppressed by the mutS-dependent mismatch repair system. Thus, it seems likely that the single-base frameshifts are derived from slippage errors that are not directly caused by DNA lesions but made indirectly during some kind of error-prone DNA synthesis in the recA mutant cells.     

Cloning of the argF gene encoding the ornithine carbamoyltransferase from Corynebacterium glutamicum.

The argF gene encoding ornithine carbamoyl-transferase (OTCase; EC2.1.3.3) has been cloned from Corynebacterium glutamicum by transforming the Escherichia coli arginine auxotroph with the genomic DNA library. The cloned DNA also complements the E. coli argG mutant, suggesting a clustered organization of the genes in the genome. We have determined the DNA sequence of the minimal fragment complementing the E. coli argF mutant. The coding region of the cloned gene is 957 nucleotides long with a deduced molecular mass of about 35 kDa polypeptide. The enzyme activity and size of the expressed protein in the E. coli auxotroph carrying the argF gene revealed that the cloned gene indeed codes for OTCase. Analysis of the amino acid sequence of the predicted protein revealed a strong similarity to the corresponding protein of other bacteria.     

The sequence of the 6S RNA gene of Pseudomonas aeruginosa.

From the gram-negative eubacterium Pseudomonas aeruginosa we have isolated a stable 6S RNA, approximately 180 nucleotides in length. The RNA was partially sequenced and identified by comparison with the known Escherichia coli 6S RNA sequence. Southern hybridizations revealed a single copy gene coding for the 6S RNA. DNA from other prokaryotes, i.e. E. coli, Thermus thermophilus, Bacillus subtilis, Bacillus stearothermophilus and Halobacterium maris mortui, did not give detectable hybridization signals. The 6S RNA gene was cloned in E. coli and its complete primary structure was determined. Although the 6S RNA sequences from P. aeruginosa and E. coli share only a 60.4% homology, we are able to propose a common secondary structural model.     

Role of proofreading and mismatch repair in maintaining the stability of nucleotide repeats in DNA.

The role of the proofreading exonuclease in maintaining the stability of multiply repeated units in DNA was studied in Escherichia coli. Reversion of plasmids in which the beta-galactosidase alpha complementing sequence was moved +2 out of frame by inserts containing (CA)14, (CA)5, (CA)2 or (TA)6 or +1 by creating a run of 8 C was compared in mutS and mutSdnaQ strains. Proofreading corrects at least half of the frameshift errors for all the plasmids and at least 99% of the errors in the (CA)2 plasmid. The (CA)2 plasmid reverts mostly by +1 frameshifts in the restriction sites flanking the insert. With the (CA)14, (TA)6, (CA)5 and 8C plasmids, reversion is mainly by loss of a repeat unit. The data support the hypothesis that the dnaQgene product recognizes frameshifts close to the DNA growing point. Frameshifts distal to the growing point are mainly corrected by mismatch repair.We speculate that mismatches in mononucleotide repeats are susceptible to proofreading because they can either migrate to a point where they are recognized by the exonuclease or, alternatively, because single nucleotide distortions are more readily detected than dinucleotides.     

High frequencies of short frameshifts in poly-CA/TG tandem repeats borne by bacteriophage M13 in Escherichia coli K-12.

Slipped-strand mispairing (SSM) may play an major role in repetitive DNA sequence evolution by generating large numbers of short frameshift mutations within simple tandem repeats. Here we examine the frequency and size spectrum of frameshifts generated within poly-CA/TG sequences inserted into bacteriophage M13 in Escherichia coli hosts. The frequency of detectable frameshifts within a 40 bp tract of poly-CA/TG is greater than one percent and increases more than linearly with length, being lower by a factor of four in a 22 bp target sequence. The frequency increases more than 13-fold in mutL and mutS host cells, suggesting that a high proportion of frameshift events are normally repaired by methyl-directed mismatch repair. Of the 87 sequenced frameshifts in this study, 96% result from deletion or insertion of only or two 2 bp repeat units. The most frequent events are 2 bp deletions, 2 bp insertions, and 4 bp deletions, the relative frequencies of these events being about 18:6:1.     

A molecular characterization of spontaneous frameshift mutagenesis within the trpA gene of Escherichia coli

Spontaneous frameshift mutations are an important source of genetic variation in all species and cause a large number of genetic disorders in humans. To enhance our understanding of the molecular mechanisms of frameshift mutagenesis, 583 spontaneous Trp+ revertants of two trpA frameshift alleles in Escherichia coli were isolated and DNA sequenced. In order to measure the contribution of methyl-directed mismatch repair to frameshift production, mutational spectra were constructed for both mismatch repair-proficient and repair-defective strains. The molecular origins of practically all of the frameshifts analyzed could be explained by one of six simple models based upon misalignment of the template or nascent DNA strands with or without misincorporation of primer nucleotides during DNA replication. Most frameshifts occurred within mononucleotide runs as has been shown often in previous studies but the location of the 76 frameshift sites was usually outside of runs. Mismatch repair generally was most effective in preventing the occurrence of frameshifts within runs but there was much variation from site to site. Most frameshift sites outside of runs appear to be refractory to mismatch repair although the small number of occurrences at most of these sites make firm conclusions impossible. There was a dense pattern of reversion sites within the trpA DNA region where reversion events could occur, suggesting that, in general, most DNA sequences are capable of undergoing spontaneous mutational events during replication that can lead to small deletions and insertions. Many of these errors are likely to occur at low frequencies and be tolerated as events too costly to prevent or repair. These studies also revealed an unpredicted flexibility in the primary amino acid sequence of the trpA product, the alpha subunit of tryptophan synthase.     

Polypeptide chain elongation factor 1 alpha (EF-1 alpha) from yeast: nucleotide sequence of one of the two genes for EF-1 alpha from Saccharomyces cerevisiae.

Messenger RNA for yeast cytosolic polypeptide chain elongation factor 1 alpha (EF-1 alpha) was partially purified from Saccharomyces cerevisiae. Double-stranded complementary DNA (cDNA) was synthesized and cloned in Escherichia coli with pBR327 as a vector. Recombinant plasmid carrying yEF-1 alpha cDNA was identified by cross-hybridization with the E. coli tufB gene and the yeast mitochondrial EF-Tu gene (tufM) under non-stringent conditions. A yeast gene library was then screened with the EF-1 alpha cDNA and several clones containing the chromosomal gene for EF-1 alpha were isolated. Restriction analysis of DNA fragments of these clones as well as the Southern hybridization of yeast genomic DNA with labelled EF-1 alpha cDNA indicated that there are two EF-1 alpha genes in S. cerevisiae. The nucleotide sequence of one of the two EF-1 alpha genes (designated as EF1 alpha A) was established together with its 5'- and 3'-flanking sequences. The sequence contained 1374 nucleotides coding for a protein of 458 amino acids with a calculated mol. wt. of 50 300. The derived amino acid sequence showed homologies of 31% and 32% with yeast mitochondrial EF-Tu and E. coli EF-Tu, respectively.     

A common structural feature in promoter sequences of E. coli.

We have searched promoter regions of E. coli, structural genes of the same organism, and computer-generated random sequence DNA for the occurrence of common structural features. This is done by converting the base sequence to a series of numbers representing the sequence of helix twist angles and examining these numerical sequences statistically. Common structural features are shared by the promoter regions with a much higher frequency than are found in structural genes or in random sequences. These structures appear to be scattered randomly throughout the promoters, both in terms of the number of such structures per promoter and in terms of location within each promoter. One particular structure consisting of five successive helix twist angles is reported, along with a list of 60 different hexanucleotide sequences that share this structure. The locations of these structural elements in 61 E. coli promoters are also tabulated.     

Frameshifts and deletions during in vitro translesion synthesis past Pt-DNA adducts by DNA polymerases beta and eta

DNA polymerases beta (pol beta ) and eta (pol eta ) are the only two eukaryotic polymerases known to efficiently bypass cisplatin and oxaliplatin adducts in vitro. Frameshift errors are an important aspect of mutagenesis. We have compared the types of frameshifts that occur during translesion synthesis past cisplatin and oxaliplatin adducts in vitro by pol beta and pol eta on a template containing multiple runs of nucleotides flanking a single platinum-GG adduct. Translesion synthesis past platinum adducts by pol beta resulted in approximately 50% replication products containing single-base deletions. For both adducts the majority of -1 frameshifts occurred in a TTT sequence 3-5 bp upstream of the DNA lesion. For pol eta, all of the bypass products for both cisplatin and oxaliplatin adducts contained -1 frameshifts in the upstream TTT sequence and most of the products of replication on oxaliplatin-damaged templates had multiple replication errors, both frameshifts and misinsertions. In addition, on platinated templates both polymerases generated replication products 4-8 bp shorter than the full-length products. The majority of short cisplatin-induced products contained an internal deletion which included the adduct. In contrast, the majority of oxaliplatin-induced short products contained a 3' terminal deletion. The implications of these in vitro results for in vivo mutagenesis are discussed.     

Escherichia coli mutator (Delta)polA is defective in base mismatch correction: the nature of in vivo DNA replication errors

We constructed a set of Escherichia coli strains containing deletions in genes encoding three SOS polymerases, and defective in MutS and DNA polymerase I (PolI) mismatch repair, and estimated the rate and specificity of spontaneous endogenous tonB(+)-->tonB- mutations. The rate and specificity of mutations in strains proficient or deficient in three SOS polymerases was compared and found that there was no contribution of SOS polymerases to the chromosomal tonB mutations. MutS-deficient strains displayed elevated spontaneous mutation rates, consisting of dominantly minus frameshifts and transitions. Minus frameshifts are dominated by warm spots at run-bases. Among 57 transitions (both G:C-->A:T and A:T-->G:C), 35 occurred at two hotspot sites. PolI-deficient strains possessed an increased rate of deletions and frameshifts, because of a deficiency in postreplicative deletion and frameshift mismatch corrections. Frameshifts in PolI-deficient strains occurred within the entire tonB gene at non-run and run sequences. MutS and PolI double deficiency indicated a synergistic increase in the rate of deletions, frameshifts and transitions. In this case, mutS-specific hotspots for frameshifts and transitions disappeared. The results suggested that, unlike the case previously known pertaining to postreplicative MutS mismatch repair for frameshifts and transitions and PolI mismatch repair for frameshifts and deletions, PolI can recognize and correct transition mismatches. Possible mechanisms for distinct MutS and PolI mismatch repair are discussed. A strain containing deficiencies in three SOS polymerases, MutS mismatch repair and PolI mismatch repair was also constructed. The spectrum of spontaneous mutations in this strain is considered to represent the spectrum of in vivo DNA polymerase III replication errors. The mutation rate of this strain was 219x10(-8), about a 100-fold increase relative to the wild-type strain. Uncorrected polymerase III replication errors were predominantly frameshifts and base substitutions followed by deletions.