Cloning, sequencing, and expression of RecA proteins from three distantly related thermophilic eubacteria

Genomewide comparison and novel ncRNAs of Aquificales

Background

The Aquificales are a diverse group of thermophilic bacteria that thrive in terrestrial and marine hydrothermal environments. They can be divided into the families Aquificaceae, Desulfurobacteriaceae and Hydrogenothermaceae. Although eleven fully sequenced and assembled genomes are available, only little is known about this taxonomic order in terms of RNA metabolism.

Results

In this work, we compare the available genomes, extend their protein annotation, identify regulatory sequences, annotate non-coding RNAs (ncRNAs) of known function, predict novel ncRNA candidates, show idiosyncrasies of the genetic decoding machinery, present two different types of transfer-messenger RNAs and variations of the CRISPR systems. Furthermore, we performed a phylogenetic analysis of the Aquificales based on entire genome sequences, and extended this by a classification among all bacteria using 16S rRNA sequences and a set of orthologous proteins.

Combining several in silico features (e.g. conserved and stable secondary structures, GC-content, comparison based on multiple genome alignments) with an in vivo dRNA-seq transcriptome analysis of Aquifex aeolicus, we predict roughly 100 novel ncRNA candidates in this bacterium.

Conclusions

We have here re-analyzed the Aquificales, a group of bacteria thriving in extreme environments, sharing the feature of a small, compact genome with a reduced number of protein and ncRNA genes. We present several classical ncRNAs and riboswitch candidates. By combining in silico analysis with dRNA-seq data of A. aeolicus we predict nearly 100 novel ncRNA candidates.

The hierarchical system of the 'Alphaproteobacteria': description of Hyphomonadaceae fam. nov., Xanthobacteraceae fam. nov. and Erythrobacteraceae fam. nov

Phylogenetic analysis of the class 'Alphaproteobacteria', including physiologically diverse species, was conducted by using small-subunit rRNA gene sequences. The 16S rRNA gene sequences of 261 species in the class 'Alphaproteobacteria' were obtained from GenBank/EMBL/DDBJ for constructing a phylogenetic tree by using maximum-likelihood analysis. In the resulting tree, members of the class 'Alphaproteobacteria' were subdivided into five major clusters, which were compared with the taxonomic outline of Bergey's Manual of Systematic Biology and the arb tree. Based on this phylogenetic tree, three novel families are proposed: Hyphomonadaceae fam. nov. to accommodate the bacterial genera Hyphomonas, Hirschia, Maricaulis and Oceanicaulis, Xanthobacteraceae fam. nov. to include the genera Xanthobacter, Azorhizobium, Ancylobacter, Labrys and Starkeya, and Erythrobacteraceae fam. nov. to accommodate the genera Erythrobacter, Porphyrobacter and Erythromicrobium. The phylogenetic tree of 16S rRNA gene sequences established in this study may provide a sound basis for future taxonomic reconstruction of the class 'Alphaproteobacteria'.

Temperature-dependent hypermutational phenotype in recA mutants of Thermus thermophilus HB27

The recA gene from Thermus thermophilus HB27 was cloned and engineered to obtain insertion (recA::kat) and deletion (deltarecA) derivatives. Transcription of recA in this extreme thermophile was induced by mitomycin C, leading to the synthesis of a monocistronic mRNA. This DNA damage-mediated induction was dependent on the integrity of recA. In addition to UV sensitivity, the recA mutants of T. thermophilus showed severe pleiotropic defects, ranging from irregular nucleoid condensation and segregation to a dramatic reduction in viability during culture. An increase in the frequency of both carotenoidless and auxotrophic mutants within surviving cells of the deltarecA strain indicated a high mutation rate. As RecA is not required for plasmid transformation, we have used the alpha-lacZ gene fragment and the ampicillin resistance gene from Escherichia coli as passenger reporters to confirm such high mutation rates. Our data support the idea that the absence of RecA results in a hypermutational phenotype in T. thermophilus. Furthermore, a direct relationship is deduced between the growth temperature and mutation rate, which finally has a deleterious effect on cell survival in the absence of RecA.

Intein-mediated affinity-fusion purification of the Escherichia coli RecA protein

The RecA protein of Escherichia coli plays important roles in homologous recombination, recombinational DNA repair, and SOS induction. Because its functions are conserved among the phylogenetic kingdoms, RecA investigations have provided a paradigm for understanding these biological processes. The RecA protein has been overproduced in E. coli and purified using a variety of purification schemes requiring multiple, time-intensive steps. The purification schemes share a dependence on appropriate RecA structure and/or function at one or more steps. In this report, we used a modified protein splicing element (intein) and a chitin-binding domain, fused to the C-terminus of RecA, to facilitate a one-step affinity purification of RecA protein without modification of the native protein sequence. Following the single chromatographic step, RecA protein that is greater than 95% physical purity at a concentration of greater than microM was obtained. The protein displays in vitro activities that are identical to those of protein isolated using classical procedures. The purification strategy described here promises to yield mutant RecA proteins in sufficient quantity for rigorous biophysical characterization without dependence on intrinsic RecA function.

The ruv proteins of Thermotoga maritima: branch migration and resolution of Holliday junctions

In homologous recombination in bacteria, the RuvAB Holliday junction-specific helicase catalyzes Holliday junction branch migration, and the RuvC Holliday junction resolvase catalyzes formation of spliced or patched structures. RuvAB and RuvC from the hyperthermophile Thermotoga maritima were expressed in Escherichia coli and purified to homogeneity. An inverted repeat sequence with unique termini was produced by PCR, restriction endonuclease cleavage, and head-to-tail ligation. A second inverted repeat sequence was derived by amplification of a second template containing a three-nucleotide insertion. Reassociation products from a mixture of these two sequences were homoduplex linear molecules and heteroduplex heat-stable Holliday junctions, which acted as substrates for both T. maritima RuvAB and RuvC. The T. maritima RuvAB helicase catalyzed energy-dependent Holliday junction branch migration at 70 degrees C, leading to heteroduplex linear duplex molecules with two three-nucleotide loops. Either ATP or ATP gamma S hydrolysis served as the energy source. T. maritima RuvC resolved Holliday junctions at 70 degrees C. Remarkably, the cleavage site was identical to the preferred cleavage site for E. coli RuvC [(A/T)TT(downward arrow)(G/C)]. The conservation of function and the ease of purification of wild-type and mutant thermophilic proteins argues for the use of T. maritima proteins for additional biochemical and structural studies.

Molecular analysis of the Deinococcus radiodurans recA locus and identification of a mutation site in a DNA repair-deficient mutant, rec30

Deinococcus radiodurans strain rec30, which is a DNA damage repair-deficient mutant, has been estimated to be defective in the deinococcal recA gene. To identify the mutation site of strain rec30 and obtain information about the region flanking the gene, a 4.4-kb fragment carrying the wild-type recA gene was sequenced. It was revealed that the recA locus forms a polycistronic operon with the preceding cistrons (orf105a and orf105b). Predicted amino acid sequences of orf105a and orf105b showed substantial similarity to the competence-damage inducible protein (cinA gene product) from Streptococcus pneumoniae and the 2'-5' RNA ligase from Escherichia coli, respectively. By analyzing polymerase chain reaction (PCR) fragments derived from the genomic DNA of strain rec30, the mutation site in the strain was identified as a single G:C to A:T transition which causes an amino acid substitution at position 224 (Gly to Ser) of the deinococcal RecA protein. Furthermore, we succeeded in expressing both the wild-type and mutant recA genes of D. radiodurans in E. coli without any obvious toxicity or death. The gamma-ray resistance of an E. coli recA1 strain was fully restored by the expression of the wild-type recA gene of D. radiodurans that was cloned in an E. coli vector plasmid. This result is consistent with evidence that RecA proteins from many bacterial species can functionally complement E. coli recA mutants. In contrast with the wild-type gene, the mutant recA gene derived from strain rec30 did not complement E. coli recA1, suggesting that the mutant RecA protein lacks functional activity for recombinational repair.

Biochemical properties of a high fidelity DNA ligase from Thermus species AK16D

NAD+-dependent DNA ligases from thermophilic bacteria Thermus species are highly homologous with amino acid sequence identities ranging from 85 to 98%. Thermus species AK16D ligase, the most divergent of the seven Thermus isolates collected worldwide, was cloned, expressed in Escherichia coli and purified to homogeneity. This Thermus ligase is similar to Thermus thermophilus HB8 ligase with respect to pH, salt, NAD+, divalent cation profiles and steady-state kinetics.However, the former is more discriminative toward T/G mismatches at the 3'-side of the ligation junction, as judged by the ratios of initial ligation rates of matched and mismatched substrates. The two wild-type Thermus ligases and a Tth ligase mutant (K294R) demonstrate 1-2 orders of magnitude higher fidelity than viral T4 DNA ligase. Both Thermus ligases are active with either the metal cofactor Mg2+, Mn2+or Ca2+but not with Co2+, Ni2+, Cu2+or Zn2+. While the nick closure step with Ca2+becomes rate-limiting which results in the accumulation of DNA-adenylate intermediate, Ni2+only supports intermediate formation to a limited extent. Both Thermus ligases exhibit enhanced mismatch ligation when Mn2+is substituted for Mg2+, but the Tsp. AK16D ligase remains more specific toward perfectly matched substrate.

Characterization of thermostable RecA protein and analysis of its interaction with single-stranded DNA

Thermostable RecA protein (ttRecA) from Thermus thermophilus HB8 showed strand exchange activity at 65 degrees C but not at 37 degrees C, although nucleoprotein complex was observed at both temperatures. ttRecA showed single-stranded DNA (ssDNA)-dependent ATPase activity, and its activity was maximal at 65 degrees C. The kinetic parameters, K(m) and kcat, for adenosine triphosphate (ATP) hydrolysis with poly(dT) were 1.4 mM and 0.60 s-1 at 65 degrees C, and 0.34 mM and 0.28 s-1 at 37 degrees C, respectively. Substrate cooperativity was observed at both temperatures, and the Hill coefficient was about 2. At 65 degrees C, all tested ssDNAs were able to stimulate the ATPase activity. The order of ATPase stimulation was: poly(dC) > poly(dT) > M13 ssDNA > poly(dA). Double-stranded DNAs (dsDNA), poly(dT).poly(dA) and M13 dsDNA, were unable to activate the enzyme at 65 degrees C. At 37 degrees C, however, not only dsDNAs but also poly(dA) and M13 ssDNA showed poor stimulating ability. At 25 degrees C, poly(dA) and M13 ssDNA gave circular dichroism (CD) peaks at around 192 nm, which reflect a particular structure of DNA. The conformation was changed by an upshift of temperature or binding to Escherichia coli RecA protein (ecRecA), but not to ttRecA. The dissociation constant between ecRecA and poly(dA) was estimated to be 44 microM at 25 degrees C by the change in the CD. These observations suggest that the capability to modify the conformation of ssDNA may be different between ttRecA and ecRecA. The specific structure of ssDNA was altered by heat or binding of ecRecA. After this alteration, ttRecA and ecRecA can express their activities at each physiological temperature.

Protein phylogenies and signature sequences: A reappraisal of evolutionary relationships among archaebacteria, eubacteria, and eukaryotes

The presence of shared conserved insertion or deletions (indels) in protein sequences is a special type of signature sequence that shows considerable promise for phylogenetic inference. An alternative model of microbial evolution based on the use of indels of conserved proteins and the morphological features of prokaryotic organisms is proposed. In this model, extant archaebacteria and gram-positive bacteria, which have a simple, single-layered cell wall structure, are termed monoderm prokaryotes. They are believed to be descended from the most primitive organisms. Evidence from indels supports the view that the archaebacteria probably evolved from gram-positive bacteria, and I suggest that this evolution occurred in response to antibiotic selection pressures. Evidence is presented that diderm prokaryotes (i.e., gram-negative bacteria), which have a bilayered cell wall, are derived from monoderm prokaryotes. Signature sequences in different proteins provide a means to define a number of different taxa within prokaryotes (namely, low G+C and high G+C gram-positive, Deinococcus-Thermus, cyanobacteria, chlamydia-cytophaga related, and two different groups of Proteobacteria) and to indicate how they evolved from a common ancestor. Based on phylogenetic information from indels in different protein sequences, it is hypothesized that all eukaryotes, including amitochondriate and aplastidic organisms, received major gene contributions from both an archaebacterium and a gram-negative eubacterium. In this model, the ancestral eukaryotic cell is a chimera that resulted from a unique fusion event between the two separate groups of prokaryotes followed by integration of their genomes.

The complete genome of the hyperthermophilic bacterium Aquifex aeolicus

Aquifex aeolicus was one of the earliest diverging, and is one of the most thermophilic, bacteria known. It can grow on hydrogen, oxygen, carbon dioxide, and mineral salts. The complex metabolic machinery needed for A. aeolicus to function as a chemolithoautotroph (an organism which uses an inorganic carbon source for biosynthesis and an inorganic chemical energy source) is encoded within a genome that is only one-third the size of the E. coli genome. Metabolic flexibility seems to be reduced as a result of the limited genome size. The use of oxygen (albeit at very low concentrations) as an electron acceptor is allowed by the presence of a complex respiratory apparatus. Although this organism grows at 95 degrees C, the extreme thermal limit of the Bacteria, only a few specific indications of thermophily are apparent from the genome. Here we describe the complete genome sequence of 1,551,335 base pairs of this evolutionarily and physiologically interesting organism.

Insights into thermal resistance of proteins from the intrinsic stability of their alpha-helices

To investigate the role of alpha helices in protein thermostability, we compared energy characteristics of alpha helices from thermophilic and mesophilic proteins belonging to four protein families of known three-dimensional structure, for at least one member of each family. The changes in intrinsic free energy of alpha-helix formation were estimated using the statistical mechanical theory for describing helix/coil transitions in peptide helices [Munoz, V., Serrano, L. Nature Struc. Biol. 1:399-409, 1994; Munoz, V., Serrano, L. J. Mol. Biol. 245:275-296, 1995; Munoz, V., Serrano, L. J. Mol. Biol. 245:297-308, 1995]. Based on known sequences of mesophilic and thermophilic RecA proteins we found that (1) a high stability of alpha helices is necessary but is not a sufficient condition for thermostability of RecA proteins, (2) the total helix stability, rather than that of individual helices, is the factor determining protein thermostability, and (3) two facets of intrahelical interactions, the intrinsic helical propensities of amino acids and the side chain-side chain interactions, are the main contributors to protein thermostability. Similar analysis applied to families of L-lactate dehydrogenases, seryl-tRNA synthetases, and aspartate amino transferases led us to conclude that an enhanced total stability of alpha helices is a general feature of many thermophilic proteins. The magnitude of the observed decrease in intrinsic free energy on alpha-helix formation of several thermoresistant proteins was found to be sufficient to explain the experimentally determined increase of their thermostability. Free energies of intrahelical interactions of different RecA proteins calculated at three temperatures that are thought to be close to its normal environmental conditions were found to be approximately equal. This indicates that certain flexibility of RecA protein structure is an essential factor for protein function. All RecA proteins analyzed fell into three temperature-dependent classes of similar alpha-helix stability (delta G(int) = 45.0 +/- 2.0 kcal/mol). These classes were consistent with the natural origin of the proteins. Based on the sequences of protein alpha helices with optimized arrangement of stabilizing interactions, a natural reserve of RecA protein thermoresistance was estimated to be sufficient for conformational stability of the protein at nearly 200 degrees C.

The recA gene from Streptomyces rimosus R6: sequence and expression in Escherichia coli

The recA gene from Streptomyces rimusus encodes a 376-amino acids polypeptide (M(r) 39,702) that is one of the largest bacterial RecA proteins observed. Detailed analyses of the Streptomyces RecA proteins showed that all possess an additional and unique C-terminal, rich in lysines and alanines, which can form an additional terminal alpha helix. Expression of the S. rimosus RecA protein in Escherichia coli FR333 (delta recA306) was demonstrated using antibodies raised against E. coli RecA protein; expression was possible only from the S. rimosus promoter. A Streptomyces-E. coli-like promoter sequence (TTGACA-18bp-TCTTAT) was found in the A+ T-rich region 135-165 base pairs upstream from the initiation codon and was related to Bacillus subtilis DNA damage-inducible promoters.

Cloning, sequencing, and expression of ruvB and characterization of RuvB proteins from two distantly related thermophilic eubacteria

The ruvB genes of the highly divergent thermophilic eubacteria Thermus thermophilus and Thermotoga maritima were cloned, sequenced, and expressed in Escherichia coli. Both thermostable RuvB proteins were purified to homogeneity. Like E. coli RuvB protein, both purified thermostable RuvB proteins showed strong double-stranded DNA-dependent ATPase activity at their temperature optima (> or = 70 degrees C). In the absence of ATP, T. thermophilus RuvB protein bound to linear double-stranded DNA with a preference for the ends. Addition of ATP or gamma-S-ATP destabilized the T. thermophilus RuvB-DNA complexes. Both thermostable RuvB proteins displayed helicase activity on supercoiled DNA. Expression of thermostable T. thermophilus RuvB protein in the E. coli ruvB recG mutant strain N3395 partially complemented the UV-sensitive phenotype, suggesting that T. thermophilus RuvB protein has a function similar to that of E. coli RuvB in vivo.

The recA gene of Borrelia burgdorferi

The nucleotide sequence of the Borrelia burgdorferi (Bb) Sh-2-82 recA gene has been determined using PCR-based approaches without the construction of a genomic library. The gene should encode a protein of 365 amino acids which is highly homologous to other known RecA proteins. It represents a new homolog from a distinct phylogenetic branch of eubacteria. Although, previous reports concluded that recA is absent from Bb, the identification presented here conclusively shows its presence and reaffirms the ubiquity of RecA in prokaryotes.