Interference with histidyl-tRNA synthetase by a CRISPR spacer sequence as a factor in the evolution of Pelobacter carbinolicus
© Aklujkar and Lovley; licensee BioMed Central Ltd. 2010
Received: 28 April 2010
Accepted: 28 July 2010
Published: 28 July 2010
Pelobacter carbinolicus, a bacterium of the family Geobacteraceae, cannot reduce Fe(III) directly or produce electricity like its relatives. How P. carbinolicus evolved is an intriguing problem. The genome of P. carbinolicus contains clustered regularly interspaced short palindromic repeats (CRISPR) separated by unique spacer sequences, which recent studies have shown to produce RNA molecules that interfere with genes containing identical sequences.
CRISPR spacer #1, which matches a sequence within hisS, the histidyl-tRNA synthetase gene of P. carbinolicus, was shown to be expressed. Phylogenetic analysis and genetics demonstrated that a gene paralogous to hisS in the genomes of Geobacteraceae is unlikely to compensate for interference with hisS. Spacer #1 inhibited growth of a transgenic strain of Geobacter sulfurreducens in which the native hisS was replaced with that of P. carbinolicus. The prediction that interference with hisS would result in an attenuated histidyl-tRNA pool insufficient for translation of proteins with multiple closely spaced histidines, predisposing them to mutation and elimination during evolution, was investigated by comparative genomics of P. carbinolicus and related species. Several ancestral genes with high histidine demand have been lost or modified in the P. carbinolicus lineage, providing an explanation for its physiological differences from other Geobacteraceae.
The disappearance of multiheme c-type cytochromes and other genes typical of a metal-respiring ancestor from the P. carbinolicus lineage may be the consequence of spacer #1 interfering with hisS, a condition that can be reproduced in a heterologous host. This is the first successful co-introduction of an active CRISPR spacer and its target in the same cell, the first application of a chimeric CRISPR construct consisting of a spacer from one species in the context of repeats of another species, and the first report of a potential impact of CRISPR on genome-scale evolution by interference with an essential gene.
Clustered regularly interspaced short palindromic repeats (CRISPR), which consist of direct repeats of a short sequence (21-47 bp) separated by nonrepetitive sequences of similar size, have been identified in the genome sequences of almost all archaea and numerous bacteria, with a variable complement of adjacent CRISPR-associated (cas) genes [1–9]. A fraction of the spacer sequences between repeats have been found to match sequences termed "proto-spacers" within genes, from which they may be derived [8, 10, 11], and the fact that many of these genes belong to phage or plasmid entities led to the hypothesis that CRISPR and the Cas proteins may function as an RNA interference-based immune system . The link between specific CRISPR spacers and proto-spacers and phage resistance has been established by mutational analysis in Streptococcus thermophilus [12, 13], and by testing synthetic CRISPR constructs in Escherichia coli . Similarly, resistance of Staphylococcus epidermidis to a conjugative plasmid has been shown to depend on a CRISPR spacer and the corresponding proto-spacer . Expression of CRISPR loci as long transcripts processed into smaller RNA molecules has been observed in several archaea [5, 16–19] and bacteria [14, 20]. A complex of Cas proteins has been shown to carry out this processing in E. coli and to be required for resistance to infection ; a different protein (Cas6) processes CRISPR transcripts in Pyrococcus furiosus . CRISPR-derived RNAs have been shown to form RNA-protein complexes in P. furiosus , which leads to degradation of RNAs containing matching proto-spacers , whereas DNA was shown to be the target of interference by spacer-containing RNAs in S. epidermidis . Although CRISPR are widely regarded as an immunological phenomenon, CRISPR and cas genes have also been implicated in spore development of Myxococcus xanthus [3, 23] and in inhibition of biofilm formation and swarming of Pseudomonas aeruginosa by a lysogenic phage , and there has been speculation that spacers with matches to housekeeping genes represent a novel mechanism of gene regulation .
The Geobacteraceae, a Fe(III)-respiring family of Deltaproteobacteria, are of interest for their role in bioremediation of U(VI)-contaminated environments and their ability to donate electrons directly to graphite electrodes, producing an electrical current [26, 27]. Pelobacter carbinolicus is a member of the Geobacteraceae that grows by fermentation of acetoin and 2,3-butanediol, as well as by indirect Fe(III) respiration with ethanol as the electron donor and acetate as the end product [28, 29]. Unlike its relatives in the genus Geobacter, P. carbinolicus cannot reduce Fe(III) directly in the absence of sulfur or sulfide , or produce electricity . The genome of P. carbinolicus was sequenced for the purpose of comparison to those of Geobacter species, three of which have been extensively curated: Geobacter sulfurreducens , Geobacter metallireducens  and Geobacter bemidjiensis (Aklujkar et al., submitted). This report explores how evolution of the P. carbinolicus genome may have been influenced by a spacer within the CRISPR locus that matches a proto-spacer within histidyl-tRNA synthetase (hisS), resulting in the elimination of ancestral genes containing multiple closely spaced histidines. The interfering nature of the spacer was confirmed by introducing it to a transgenic G. sulfurreducens strain containing the target gene.
Analysis of CRISPR spacers
The CRISPR locus was identified when manual curation of the P. carbinolicus genome revealed a series of suspiciously repetitive predicted genes. The nonredundant nucleotide sequence database was queried with each of the 111 CRISPR spacers of P. carbinolicus using the BLAST algorithm , with the minimum possible word size of 7 bp and without filtering out low-complexity regions of the queries. Alignments with five or fewer mismatches out of 32 bases were considered significant.
Phylogenetic analysis of HisS and HisZ proteins
The sequences of all predicted hisS gene products of the Geobacteraceae, together with HisS and HisZ protein sequences representative of various families of Bacteria and Archaea, were aligned by TCoffee  and trimmed using Mesquite (Maddison, W. P., and Maddison, D. R.. 2006. Mesquite: a modular system for evolutionary analysis. Version 1.12). Phylogenetic trees were constructed using Phylip (Felsenstein, J. 2005. PHYLIP (Phylogeny Inference Package) version 3.6) with 500 bootstrap runs.
Quantitative real-time PCR of reverse-transcribed RNA
Oligonucleotides for QRT-PCR and genetic manipulations.
Primers for QRT-PCR
Sequence (5' to 3')
PCR of hisS
PCR of hisZ
Primers for construction of hisS transgenic strain MA159
Sequence (5' to 3')
PCR on 5'
side of hisS
PCR on 3'
side of hisS
PCR of hisS
Oligonucleotides for construction of chimeric CRISPR expression plasmid pMA35
Sequence (5' to 3')
PCR of lacI-
Recombinant DNA techniques
All restriction enzymes were purchased from New England Biolabs; LA Taq polymerase was from Takara Mirus Bio; plasmids were propagated in E. coli TOP10 cells from Invitrogen; DNA purification kits for plasmids and agarose gel slices were from QIAGEN, and the MasterPure kit for genomic DNA extraction was from EPICENTRE. To construct a transgenic strain of G. sulfurreducens in which the native hisS gene was replaced with hisS from P. carbinolicus, three primer pairs (Table 1) were used to amplify the 5'-side and 3'-side flanking regions of hisS GSU1659 and the coding sequence of hisS Pcar_1041. The two flanking regions were digested with Spe I and ligated; this product and the Pcar_1041 amplicon were separately TOPO-cloned (Invitrogen) and sequenced. Digestion with BspH I (overlapping the start codon) and Spe I (overlapping the stop codon) was used to insert the Pcar_1041 gene between the flanking regions of GSU1659. As a selectable marker, a gentamicin resistance cartridge was amplified from plasmid pCM351  with Xba I site-containing primers (Table 1), maintained as a TOPO clone, and ligated into the Nhe I site between the Spe I site and the 3' flanking region of GSU1659. (A similar construct, in which only the marker was inserted without Pcar_1041, was used in unsuccessful attempts to delete GSU1659.) The entire hisS replacement construct was excised using EcoR I, purified from an agarose gel, and electroporated into the wild type G. sulfurreducens strain DL1 as previously described . An isolated gentamicin-resistant colony was streaked for purity before transfer to liquid. The genotype of this strain, called MA159G, was confirmed by PCR of genomic DNA, with primers MA0334 and MA0335 (which amplify hisS Pcar_1041 but not hisS GSU1659) as well as MA0330 and MA0333 (which give a larger product for MA159G than for DL1, due to the inserted marker). The marker, which had loxP sites on either side, was removed from the chromosome of strain MA159G by introducing the Cre recombinase expression plasmid pCM158  by electroporation and selecting for resistance to kanamycin (Sigma). Four colonies of the resultant strain called MA159 were streaked for purity and their genotypes were confirmed by PCR; amplicons were digested with Pst I to distinguish Pcar_1041 from GSU1659. This strain was electroporated with a plasmid called pRG6 (R. Glaven, personal communication), which is incompatible with pCM158 and selectable with spectinomycin (Sigma); it differs from pRG5  only in that it carries the lacI repressor gene. The chromosomal genotype of this strain was confirmed by PCR and Pst I digestion.
A plasmid vector called pMA36, incompatible with pRG6, was constructed for isopropylthio-β-D-galactopyranoside (IPTG)-inducible expression of a chimeric CRISPR containing spacer #1 from P. carbinolicus between two repeats typical of the CRISPR2 locus of G. sulfurreducens. The lacI repressor gene and taclacUV5 promoter of plasmid pCD341  were amplified by PCR with Pci I and Sph I site-containing primers (Table 1), TOPO-cloned and sequenced, and excised for ligation into plasmid pCM66 , resulting in plasmid pMA36. The chimeric CRISPR, consisting of annealed oligonucleotides MA0429 and MA0430 (Table 1), was ligated between the BamH I and EcoR I sites of pMA36. The sequence of this plasmid, called pMA35-1, was confirmed using the sequencing primer MA36R (Table 1). Serendipitously, two variants were discovered: pMA35-2 in which the chimeric CRISPR had expanded to two copies of spacer #1 with a third copy of the repeat between them, and pMA35-! in which spacer #1 was disrupted by duplication of a single G:C base pair at its exact centre (underlined in Table 1). All three chimeric CRISPR expression plasmids were electroporated into DL1 and MA159(pRG6), in parallel with pMA36 as a control. The genotypes of multiple kanamycin-resistant colonies of each transformation were confirmed by PCR of Pcar_1041 followed by Pst I digestion as well as cloning and sequencing, and by sequencing of plasmids present in the genomic DNA extracts (after transformation into E. coli to improve DNA quality). Another variant of pMA35-1 called pMA35-0 was serendipitously discovered in which spacer #1 had been deleted by recombination of the repeats on either side.
Transformants of G. sulfurreducens were selected on NBAF medium  containing 1.5% Agar Noble (Difco), supplemented with 5 mM cysteine hydrochloride and 0.1% yeast extract in an anaerobic chamber. Growth experiments were carried out with liquid cultures in either NBAF medium supplemented with 1 mM cysteine hydrochloride or FWAFC medium  modified to contain 10 mM acetate and supplemented with 1 mM ferrous ammonium sulfate, in an atmosphere of N2 and CO2 (80%:20%) in rubber-stoppered 26 ml glass tubes at 30°C.
Codon usage was determined using the CodonFrequency algorithm of the Genetics Computer Group Wisconsin Package version 10.3 (Accelrys Inc., San Diego, CA). A script to compute the number of histidines and the distances between them for every protein sequence in a list was written in Perl.
The CRISPR locus of P. carbinolicusincludes a spacer matching its own histidyl-tRNA synthetase
Quantitative detection of CRISPR spacer #1 transcripts
The sense strand spacer #1-containing RNA detected in this experiment may represent the 3' end of a long transcript initiated near the leader sequence, whereas the antisense strand spacer #1-containing RNA may be produced independently from a promoter at the opposite end of the CRISPR. It is also possible that one strand is produced from the other by an unidentified RNA-directed RNA polymerase. If the sense strand spacer #1-containing RNA undergoes processing similarly to the E. coli CRISPR transcript , which belongs to the same phylogenetic cluster of repeat sequences as the P. carbinolicus CRISPR , cleavage within the stem-loops of the repeats flanking spacer #1, followed by 3' end trimming, would release a short RNA with predicted secondary structure (Figure 2b). This or the corresponding antisense strand spacer #1-containing RNA may hybridize to the proto-spacer DNA sequence within the hisS gene (Figure 2c).
Phylogenetic and experimental evidence that interference with hisScannot be compensated
It is surprising that spacer #1 is retained by P. carbinolicus if it interferes with the essential function of histidine activation for protein synthesis. Comparative genome analysis revealed that P. carbinolicus and its close relatives, the Geobacteraceae, possess two full-length hisS-like genes, whereas other bacteria have only one. Interference with Pcar_1041 by spacer #1 might have had negligible effect if Pcar_0202 also produced histidyl-tRNA synthetase activity. However, both phylogenetic and mutational studies suggest that Pcar_1041 is essential, being the only real hisS, as detailed below.
The phylogenetic tree also demonstrates that the hisS gene Pcar_1041 was not acquired laterally; it is clearly an ancestral gene containing a proto-spacer that is not present in its closest relatives.
Spacer #1 inhibits growth of a transgenic G. sulfurreducens strain containing hisS of P. carbinolicus
Despite prior expression of the LacI repressor protein from plasmid pRG6 to prevent premature expression of the chimeric CRISPR, electroporations of MA159(pRG6) with pMA35-1 and two serendipitously obtained variants (pMA35-2 with two copies of spacer #1 and pMA35-! with spacer #1 interrupted by a single base pair insertion as shown in Table 1) were marginally successful (yielding zero to 4 colonies per attempt), whereas electroporation with an equal amount of the empty vector pMA36 produced hundreds of colonies per attempt. Electroporations of the wild type G. sulfurreducens strain DL1, carried out in parallel, yielded hundreds of colonies for all three chimeric CRISPR expression plasmids. These observations suggest that even leaky expression of the chimeric CRISPR containing spacer #1, or its mere presence as DNA, is largely incompatible with the hisS gene of P. carbinolicus containing the matching proto-spacer, which is present in the MA159 host, but not with the native hisS gene of G. sulfurreducens in the DL1 host.
Growth experiments provided further proof that spacer #1 interferes with hisS of P. carbinolicus. The G. sulfurreducens transformants were first checked by PCR, restriction digestion, and sequencing to verify that the hisS transgene and spacer #1 were intact. In one transformant, spontaneous recombination of the repeats on either side of spacer #1, eliminating it from the chimeric CRISPR expression plasmid, resulted in a strain called MA159(pMA35-0) that possessed both hisS of P. carbinolicus and the repeat sequence on the plasmid, but no spacer #1. Compared to this control that grew similarly to MA159, the presence of spacer #1 in the other MA159 transformants (i.e., with hisS of P. carbinolicus) resulted in long lag periods and somewhat reduced growth rates in NBAF medium with fumarate as the electron acceptor (Figure 6a), and very poor growth in FWAFC medium with Fe(III) citrate as the electron acceptor (Figure 6b). This effect was the same with either one copy of spacer #1 in MA159(pMA35-1) or two copies in MA159(pMA35-2), or with a single base pair insertion in spacer #1 in MA159(pMA35-!) - for the sake of clarity, only MA159(pMA35-2) is shown. The only exception was that in one experiment, triplicate cultures of MA159(pMA35-!) grew especially poorly after three transfers in NBAF (not shown). Wild type growth patterns were observed when any of the three plasmids was present in the DL1 host (i.e., with hisS of G. sulfurreducens) - for clarity, only DL1(pMA35-2) is shown (Figure 6). Although expression of spacer #1 is expected to be low in the absence of IPTG, growth inhibition of the MA159 strains was observed, and addition of IPTG had no effect, indicating that expression of chimeric CRISPR RNA was not the limiting factor for inhibition of growth.
Spacer #1 reduces the amount of hisS RNA in transgenic G. sulfurreducens no more than it affects hisZRNA
The P. carbinolicusgenome has fewer genes with numerous or closely spaced histidine codons than closely related genomes
The P. carbinolicusgenome has lost ancestral genes with numerous or closely spaced histidines
Genes of D. acetoxidans, G. sulfurreducens, G. metallireducens and G. bemidjiensis that contain 35 or more histidine codons, or have a histidine demand index above 5.0, were examined in order to identify ancestral genes that have reduced their histidine demand or have been lost specifically in the P. carbinolicus genome (Additional file 1: Table S1). Many genes found in Geobacter species are not necessarily ancestral to P. carbinolicus; they lack homologs in either the unfinished D. acetoxidans genome or the partial genome sequences of a mixture of D. acetoxidans and D. palmitatis (D. R. L. and coworkers, unpublished). Other genes that are present in D. acetoxidans, but not Geobacter species, could have been acquired after divergence from P. carbinolicus. Five gene families actually show increased histidine demand in P. carbinolicus compared to other Geobacteraceae, and in many other cases, a P. carbinolicus gene has lower histidine demand than its orthologs, but is still above the cutoff value of 5.0, or contains a similar number of histidine codons. However, after all these considerations there remain sixteen clearly ancestral gene families with typically high histidine demand that are missing or have reduced histidine demand specifically in P. carbinolicus (Additional file 1: Table S1). The functional annotations of these gene families are briefly described below.
The nuoL-1 gene family, with a high-histidine-demand representative in every Geobacter species and D. acetoxidans, encodes the L subunit of an NADH dehydrogenase I complex that P. carbinolicus has entirely lost. The NuoL protein is thought to have a proton-pumping function , for which the imidazole groups of its clustered histidines are well suited.
The znuA (or zntC) gene encodes the periplasmic protein of an ATP-binding cassette transporter, within which a cluster of histidines is thought to bind zinc with high affinity . The znuA gene of D. acetoxidans has eight histidines in its putative metal-binding cluster, whereas its closest relative, the P. carbinolicus znuA gene (Pcar_3026), has only four.
D. acetoxidans possesses two clearly ancestral genes related to arsenite S-adenosylmethyltransferase (Dace_2134, Dace_3081), with homologs in Geobacter lovleyi and other Deltaproteobacteria that contain numerous histidines. The closest homolog of Dace_3081 is in P. carbinolicus (Pcar_2089), with fewer histidines and a much reduced histidine demand index. P. carbinolicus lacks an ortholog of Dace_2134.
Two ancestral genes encoding polyketide synthase-type enzymes, with representatives in D. acetoxidans, G. lovleyi, Geobacter uraniireducens, Geobacter bemidjiensis, Geobacter sp. M21, Geobacter sp. M18, Geobacter daltonii, and Pelobacter propionicus that contain numerous histidines, have clearly been lost by P. carbinolicus. Both polyketide synthases are composed of unusual domain combinations. Dace_0979 and its orthologs contain just one elongation domain, one acyl carrier protein acylation domain, two to five acyl carrier protein domains, one reductase domain, and no dehydratase or thioesterase domains, suggesting that they build up a long-chain fatty acid with multiple hydroxyl groups, or possibly a storage polymer of a precursor resembling 3-hydroxybutanoate. Dace_1838 and its orthologs contain two elongation domains, one acyl carrier protein acylation domain, three to four dehydratase domains, and no acyl carrier protein or reductase or thioesterase domains, suggesting that their product may be polyhydroxylated, polyunsaturated, or cyclic. The loss of these two enzymes likely means that P. carbinolicus does not make secondary metabolites that are present in most other Geobacteraceae.
Three other high-histidine-demand gene families that P. carbinolicus has lost encode a DUF323 domain-containing methyltransferase (Dace_1886) with homologs in G. uraniireducens, G. bemidjiensis, Geobacter sp. M21, and Geobacter sp. M18, a glycoside hydrolase (GSU2359) with homologs in D. palmitatis, G. bemidjiensis, Geobacter sp. M21, Geobacter sp. M18, G. daltonii, and P. propionicus, and a predicted c-type heme-binding, GAF domain-containing phosphohydrolase (GSU2622) with homologs in D. palmitatis, G. metallireducens, G. uraniireducens, G. bemidjiensis, Geobacter sp. M21, Geobacter sp. M18, and G. daltonii. The specific reactions catalyzed by these enzymes are not known.
The remaining eight ancestral gene families are multiheme c-type cytochromes represented in D. acetoxidans or D. palmitatis and one or both of G. sulfurreducens and G. metallireducens. At least one gene of each family was most probably inherited by P. carbinolicus from a common ancestor of the Geobacteraceae, and then lost. Several other cytochrome families, found in the Geobacter genomes but not the incomplete D. acetoxidans and D. palmitatis genomes, may have been inherited and lost by P. carbinolicus, but the evidence is inconclusive. Nevertheless, it is notable that there are no cytochromes among the eighteen proteins of P. carbinolicus that have either more than 35 histidines or a histidine demand index above 5.0 (Additional file 1: Table S1).
Although CRISPR spacers that match phage/plasmid genes have been shown to confer immunity against infection [12, 15], and it is known that their mode of action is distinct from previously described mechanisms of phage resistance , the significant number of spacers that match host genes have not been investigated . The activity of CRISPR spacers against genes encoding essential housekeeping enzymes has not been demonstrated before, nor have its consequences for genome-scale evolution been examined. The present study establishes that when spacer #1 of the P. carbinolicus CRISPR and its putative target, the hisS gene of the same species, are both present in a transgenic G. sulfurreducens strain, interference occurs that severely affects growth. Very few cells containing the hisS transgene were able to take up the chimeric CRISPR expression plasmids. Transformants in which spacer #1 was present grew more poorly. Attempts to induce expression of spacer #1 with IPTG had no effect, but growth was almost totally inhibited by switching from fumarate to Fe(III) as the electron acceptor, under which condition expression of Cas proteins and the chromosomal CRISPR2 is upregulated . Therefore, the limiting factor for interference with hisS may be the amount of one or more Cas proteins, or the leader sequence, rather than the amount of spacer #1 transcript. Alternatively, growth by respiration of Fe(III) may require protein factors that cannot be expressed properly when histidyl-tRNA synthetase activity is low. Induction of CRISPR expression with IPTG in S. epidermidis has also shown no effect on interference with a target plasmid , but that CRISPR was inactive without a leader sequence in cis, whereas severe interference between spacer #1 and hisS requires neither the native leader sequence of the G. sulfurreducens CRISPR nor that of the P. carbinolicus CRISPR in cis.
The small decrease in the amount of hisS RNA in the presence of spacer #1, comparable to the decrease in the amount of transcript for a control housekeeping gene (hisZ) that accompanies the growth defect, is unlike the extent of decimation that one would expect if spacer #1-containing RNA catalyzed degradation of hisS RNA. It is even less likely that any transgenic strains could be viable if spacer #1-containing RNA targeted hisS DNA for degradation. An alternative possibility is that spacer #1 RNA hybridizes with the hisS gene without marking it for degradation (but perhaps recruiting one or more proteins), and must be displaced by RNA polymerase in order to complete transcription. In this scenario, even leaky expression of spacer #1 would saturate the available targets, and overexpression of spacer #1-containing RNA would not prevent RNA polymerase from displacing the one molecule obstructing transcription of the hisS gene. Consistent with these predictions, induction of spacer #1 expression with IPTG was not required for inhibition of growth, and did not cause stronger inhibition. It will be interesting to examine whether any putative nucleases have mutated in the few viable transformants that carry both spacer #1 and its target. None of the Cmr proteins implicated in RNA-targeted CRISPR function in P. furiosus  has a homolog in G. sulfurreducens or P. carbinolicus, nor does the CRISPR-associated double-stranded DNA/RNA-specific endonuclease of Sulfolobus solfataricus .
Sequencing confirmed that neither hisS nor spacer #1 was mutated in the viable transformants. Therefore, in contrast to earlier studies that showed absolute incompatibility between host spacers and phage proto-spacers, evaded only by mutation of the proto-spacers or proto-spacer-adjacent motifs [12, 13], these results indicate that it is possible to establish a spacer that persistently interferes with an essential housekeeping gene. Both spacer and target can be maintained intact over numerous generations of inhibited growth. Insertion of a single base pair in the middle of the spacer did not reduce its efficacy, indicating that a small bulge in the region of complementarity between spacer and proto-spacer does not necessarily mitigate interference.
The experiments described herein with the chimeric CRISPR construct also show that a spacer from one species can be active in the context of the somewhat different repeats and cas genes of another species. The P. carbinolicus and G. sulfurreducens genomes encode homologs of all five components of the Cascade complex (cse1 CasA Pcar_0957, GSU1385; cse2 CasB Pcar_0958, GSU1386; cse4 CasC Pcar_0959, GSU1387; cas5e CasD Pcar_0960, GSU1388; and cse3 CasE Pcar_0961, GSU1389), which processes CRISPR transcripts into target-active RNAs in E. coli , along with the Cas3 helicase (Pcar_0956, GSU1384) that is required for their activity against targets , the Cas2 endoribonuclease (Pcar_0965, GSU1393) that may degrade them , and the Cas1 endodeoxyribonuclease (Pcar_0964, GSU1392) that may aid in the acquisition of new spacers , and promote DNA/RNA annealing . However, the sequences of most of these Cas proteins are so divergent between the two species (e.g. 28% identity for Cas3) that it is noteworthy that the G. sulfurreducens system seemingly required no context other than its cognate repeat to recognize spacer #1 as a guide and hisS of P. carbinolicus as a target. There is evidence that proto-spacer-adjacent motifs are determinants of target recognition [13, 20]. The proto-spacer within hisS is followed by a CTT motif (Figure 2), typical of sequences captured by CRISPR of phylogenetic cluster 2 , to which the P. carbinolicus CRISPR and G. sulfurreducens CRISPR2 loci both belong . The facile reconstitution of interference between species with similar repeats and proto-spacer-adjacent motifs despite Cas protein divergence is encouraging for future development of CRISPR-based gene silencing technology with synthetic spacers.
Both strands of spacer #1 were detected at similar levels in P. carbinolicus RNA from actively growing cultures, in contrast with previous observations that both strands are present as RNA in Sulfolobus acidocaldarius only during stationary phase , that unequal amounts of RNA were detected with probes for the two strands in P. furiosus  and Thermus thermophilus , and that only transcripts from one strand were detected in E. coli , S. epidermidis , and Xanthomonas oryzae . The presence of RNA representing both strands is consistent with CRISPR spacers on one strand of a cluster matching both the sense and antisense of proto-spacers , and with artificial spacers in both orientations having activity against their targets, although the efficacy may differ by several orders of magnitude .
Together, the evidence of spacer #1 expression and the proof that spacer #1 can inhibit growth in a hisS-dependent manner indicate that interference with hisS by spacer #1 almost certainly occurred during the recent evolutionary history of P. carbinolicus. It is likely that this interference resulted in a growth-limiting shortage of histidyl-tRNA. As expected under this selective pressure, genes with numerous and/or closely spaced histidine codons have evidently been counterselected in the genome of P. carbinolicus. Although loss or mutation of ancestral genes is easiest to surmise from comparative genomics, acquisition of low-histidine-demand genes and failure to acquire new high-histidine-demand genes also contribute to this difference.
Among the ancestral gene losses from the P. carbinolicus genome attributable to interference with hisS by spacer #1, the most obvious are the multiheme c-type cytochromes typical of Geobacter species. The P. carbinolicus genome encodes only 14 c-type cytochromes , whereas the D. acetoxidans genome encodes at least 80 such genes (M. A., unpublished). These cytochromes are not well conserved even between closely related Geobacter species [33, 63], and it has been hypothesized that they have a generic function as capacitors . Therefore, in addition to the loss of eight ancestral multiheme c-type cytochrome families by P. carbinolicus, the non-evolution of new families may also be an effect of interference with hisS.
Of the other eight families of ancestral genes that exhibit loss or reduced histidine demand in P. carbinolicus, two have especially interesting metabolic and physiological implications. The nuoL-1 gene encodes a subunit of NADH dehydrogenase I with a conserved histidine cluster, implicated in proton pumping. Abortive expression of nuoL-1 due to interference with hisS by spacer #1 in an ancestor of P. carbinolicus, leading to a loss of function for the entire NADH dehydrogenase I complex, would have favoured the elimination of all fourteen structural genes of the complex from the genome, which is what has occurred. The reason why this deletion was not lethal is probably that another ancestral NADH dehydrogenase I complex, for which the nuoL-2 gene does not contain a cluster of histidine codons, has been retained by P. carbinolicus. Conceivably, loss of a major respiratory enzyme complex by P. carbinolicus caused it to rely more on laterally acquired fermentative pathways.
Interference with hisS by spacer #1 also offers an explanation for the diminutive histidine cluster in the znuA gene product of P. carbinolicus, which functions to bind zinc in the periplasm for active transport into the cell. Zinc is a cofactor of key metabolic enzymes such as carbonic anhydrase, acetate kinase and phosphotransacetylase [65–67]. If the mutations in the histidine cluster of znuA reduce the efficiency of zinc uptake, it could have an impact on metabolism related to the unexplained inability of P. carbinolicus to oxidize acetate  despite the presence of a complete set of tricarboxylic acid cycle genes in the genome.
This paper reports genetic and comparative genomic evidence that housekeeping genes can be targets of chronic CRISPR interference. Spacer #1 is shown to be transcribed into RNA with comparable amounts of both strands in P. carbinolicus, and to inhibit the growth of a transgenic strain of G. sulfurreducens that relies on hisS of P. carbinolicus, without drastically reducing the level of hisS RNA. The genome of P. carbinolicus exhibits the expected effect of a histidyl-tRNA shortage. Overall, the ancestral genes lost or mutated by the P. carbinolicus genome as part of its tendency towards fewer histidines per gene and lower histidine demand illustrate the de-emphasis of the metal-respiring metabolism that is typical of other Geobacteraceae. More generally, while previous studies have approached CRISPR as a microbial immune system, another important role of the system may be to exert pressure on endogenous and essential genes, resulting in dramatic changes in the genome content and physiology of the host species.
complementary deoxyribonucleic acid
clustered regularly interspaced short palindromic repeats
polymerase chain reaction
quantitative real-time polymerase chain reaction
transfer ribonucleic acid.
We thank S. Haveman, R. DiDonato, and A. Boles for technical assistance, B. Methé, R. Glaven, D. Holmes, L. Villanueva-Alvarez, M. J. Larrahando, J. Butler, and N. Young for guidance and helpful discussions. This research was supported by the Office of Science (Biological and Environmental Research), U.S. Department of Energy (Grant No. DE-FC02-02ER63446).
The Geobacter Project website: http://www.geobacter.org
- Godde JS, Bickerton A: The repetitive DNA elements called CRISPRs and their associated genes: evidence of horizontal transfer among prokaryotes. J Mol Evol. 2006, 62 (6): 718-729. 10.1007/s00239-005-0223-z.View ArticlePubMedGoogle Scholar
- Grissa I, Vergnaud G, Pourcel C: The CRISPRdb database and tools to display CRISPRs and to generate dictionaries of spacers and repeats. BMC Bioinformatics. 2007, 8: 172-10.1186/1471-2105-8-172.PubMed CentralView ArticlePubMedGoogle Scholar
- Haft DH, Selengut J, Mongodin EF, Nelson KE: A guild of 45 CRISPR-associated (Cas) protein families and multiple CRISPR/Cas subtypes exist in prokaryotic genomes. PLoS Comput Biol. 2005, 1 (6): e60-10.1371/journal.pcbi.0010060.PubMed CentralView ArticlePubMedGoogle Scholar
- Jansen R, Embden JD, Gaastra W, Schouls LM: Identification of genes that are associated with DNA repeats in prokaryotes. Mol Microbiol. 2002, 43 (6): 1565-1575. 10.1046/j.1365-2958.2002.02839.x.View ArticlePubMedGoogle Scholar
- Lillestøl RK, Redder P, Garrett RA, Brugger K: A putative viral defence mechanism in archaeal cells. Archaea. 2006, 2 (1): 59-72. 10.1155/2006/542818.PubMed CentralView ArticlePubMedGoogle Scholar
- Makarova K, Grishin N, Shabalina S, Wolf Y, Koonin E: A putative RNA-interference-based immune system in prokaryotes: computational analysis of the predicted enzymatic machinery, functional analogies with eukaryotic RNAi, and hypothetical mechanisms of action. Biology Direct. 2006, 1 (1): 7-10.1186/1745-6150-1-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Makarova KS, Aravind L, Grishin NV, Rogozin IB, Koonin EV: A DNA repair system specific for thermophilic Archaea and bacteria predicted by genomic context analysis. Nucleic Acids Res. 2002, 30 (2): 482-496. 10.1093/nar/30.2.482.PubMed CentralView ArticlePubMedGoogle Scholar
- Mojica FJ, Diez-Villasenor C, Garcia-Martinez J, Soria E: Intervening sequences of regularly spaced prokaryotic repeats derive from foreign genetic elements. J Mol Evol. 2005, 60 (2): 174-182. 10.1007/s00239-004-0046-3.View ArticlePubMedGoogle Scholar
- Mojica FJ, Diez-Villasenor C, Soria E, Juez G: Biological significance of a family of regularly spaced repeats in the genomes of Archaea, Bacteria and mitochondria. Mol Microbiol. 2000, 36 (1): 244-246. 10.1046/j.1365-2958.2000.01838.x.View ArticlePubMedGoogle Scholar
- Bolotin A, Quinquis B, Sorokin A, Ehrlich SD: Clustered regularly interspaced short palindrome repeats (CRISPRs) have spacers of extrachromosomal origin. Microbiology. 2005, 151 (Pt 8): 2551-2561. 10.1099/mic.0.28048-0.View ArticlePubMedGoogle Scholar
- Pourcel C, Salvignol G, Vergnaud G: CRISPR elements in Yersinia pestis acquire new repeats by preferential uptake of bacteriophage DNA, and provide additional tools for evolutionary studies. Microbiology. 2005, 151 (Pt 3): 653-663. 10.1099/mic.0.27437-0.View ArticlePubMedGoogle Scholar
- Barrangou R, Fremaux C, Deveau H, Richards M, Boyaval P, Moineau S, Romero DA, Horvath P: CRISPR provides acquired resistance against viruses in prokaryotes. Science. 2007, 315 (5819): 1709-1712. 10.1126/science.1138140.View ArticlePubMedGoogle Scholar
- Deveau H, Barrangou R, Garneau JE, Labonte J, Fremaux C, Boyaval P, Romero DA, Horvath P, Moineau S: Phage response to CRISPR-encoded resistance in Streptococcus thermophilus. J Bacteriol. 2008, 190 (4): 1390-1400. 10.1128/JB.01412-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Brouns SJ, Jore MM, Lundgren M, Westra ER, Slijkhuis RJ, Snijders AP, Dickman MJ, Makarova KS, Koonin EV, van der Oost J: Small CRISPR RNAs guide antiviral defense in prokaryotes. Science. 2008, 321 (5891): 960-964. 10.1126/science.1159689.View ArticlePubMedGoogle Scholar
- Marraffini LA, Sontheimer EJ: CRISPR interference limits horizontal gene transfer in staphylococci by targeting DNA. Science. 2008, 322 (5909): 1843-1845. 10.1126/science.1165771.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang TH, Bachellerie JP, Rozhdestvensky T, Bortolin ML, Huber H, Drungowski M, Elge T, Brosius J, Huttenhofer A: Identification of 86 candidates for small non-messenger RNAs from the archaeon Archaeoglobus fulgidus. Proc Natl Acad Sci USA. 2002, 99 (11): 7536-7541. 10.1073/pnas.112047299.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang TH, Polacek N, Zywicki M, Huber H, Brugger K, Garrett R, Bachellerie JP, Huttenhofer A: Identification of novel non-coding RNAs as potential antisense regulators in the archaeon Sulfolobus solfataricus. Mol Microbiol. 2005, 55 (2): 469-481. 10.1111/j.1365-2958.2004.04428.x.View ArticlePubMedGoogle Scholar
- Lillestøl RK, Shah SA, Brugger K, Redder P, Phan H, Christiansen J, Garrett RA: CRISPR families of the crenarchaeal genus Sulfolobus: bidirectional transcription and dynamic properties. Mol Microbiol. 2009, 72 (1): 259-272. 10.1111/j.1365-2958.2009.06641.x.View ArticlePubMedGoogle Scholar
- Hale C, Kleppe K, Terns RM, Terns MP: Prokaryotic silencing (psi)RNAs in Pyrococcus furiosus. RNA. 2008, 14 (12): 2572-2579. 10.1261/rna.1246808.PubMed CentralView ArticlePubMedGoogle Scholar
- Semenova E, Nagornykh M, Pyatnitskiy M, Artamonova II, Severinov K: Analysis of CRISPR system function in plant pathogen Xanthomonas oryzae. FEMS Microbiol Lett. 2009, 296 (1): 110-116. 10.1111/j.1574-6968.2009.01626.x.View ArticlePubMedGoogle Scholar
- Carte J, Wang R, Li H, Terns RM, Terns MP: Cas6 is an endoribonuclease that generates guide RNAs for invader defense in prokaryotes. Genes Dev. 2008, 22 (24): 3489-3496. 10.1101/gad.1742908.PubMed CentralView ArticlePubMedGoogle Scholar
- Hale CR, Zhao P, Olson S, Duff MO, Graveley BR, Wells L, Terns RM, Terns MP: RNA-guided RNA cleavage by a CRISPR RNA-Cas protein complex. Cell. 2009, 139 (5): 945-956. 10.1016/j.cell.2009.07.040.PubMed CentralView ArticlePubMedGoogle Scholar
- Viswanathan P, Murphy K, Julien B, Garza AG, Kroos L: Regulation of dev, an operon that includes genes essential for Myxococcus xanthus development and CRISPR-associated genes and repeats. J Bacteriol. 2007, 189 (10): 3738-3750. 10.1128/JB.00187-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Zegans ME, Wagner JC, Cady KC, Murphy DM, Hammond JH, O'Toole GA: Interaction between bacteriophage DMS3 and host CRISPR region inhibits group behaviors of Pseudomonas aeruginosa. J Bacteriol. 2009, 191 (1): 210-219. 10.1128/JB.00797-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Sorek R, Kunin V, Hugenholtz P: CRISPR--a widespread system that provides acquired resistance against phages in bacteria and archaea. Nat Rev Microbiol. 2008, 6 (3): 181-186. 10.1038/nrmicro1793.View ArticlePubMedGoogle Scholar
- Lovley DR: The microbe electric: conversion of organic matter to electricity. Curr Opin Biotechnol. 2008, 19 (6): 564-571. 10.1016/j.copbio.2008.10.005.View ArticlePubMedGoogle Scholar
- Lovley DR, Holmes DE, Nevin KP: Dissimilatory Fe(III) and Mn(IV) reduction. Adv Microb Physiol. 2004, 49: 219-286. 10.1016/S0065-2911(04)49005-5.View ArticlePubMedGoogle Scholar
- Lovley DR, Phillips EJ, Lonergan DJ, Widman PK: Fe(III) and S0 reduction by Pelobacter carbinolicus. Appl Environ Microbiol. 1995, 61 (6): 2132-2138.PubMed CentralPubMedGoogle Scholar
- Schink B: Fermentation of 2,3-butanediol by Pelobacter carbinolicus sp. nov. and Pelobacter propionicus sp. nov., and evidence for propionate formation from C2 compounds. Arch Microbiol. 1984, 137: 33-41. 10.1007/BF00425804.View ArticleGoogle Scholar
- Haveman SA, DiDonato RJ, Villanueva L, Shelobolina ES, Postier BL, Xu B, Liu A, Lovley DR: Genome-wide gene expression patterns and growth requirements suggest that Pelobacter carbinolicus reduces Fe(III) indirectly via sulfide production. Appl Environ Microbiol. 2008, 74 (14): 4277-4284. 10.1128/AEM.02901-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Richter H, Lanthier M, Nevin KP, Lovley DR: Lack of electricity production by Pelobacter carbinolicus indicates that the capacity for Fe(III) oxide reduction does not necessarily confer electron transfer ability to fuel cell anodes. Appl Environ Microbiol. 2007, 73 (16): 5347-5353. 10.1128/AEM.00804-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Methé BA, Nelson KE, Eisen JA, Paulsen IT, Nelson W, Heidelberg JF, Wu D, Wu M, Ward N, Beanan MJ, et al: Genome of Geobacter sulfurreducens: metal reduction in subsurface environments. Science. 2003, 302 (5652): 1967-1969. 10.1126/science.1088727.View ArticlePubMedGoogle Scholar
- Aklujkar M, Krushkal J, DiBartolo G, Lapidus A, Land ML, Lovley DR: The genome sequence of Geobacter metallireducens: features of metabolism, physiology and regulation common and dissimilar to Geobacter sulfurreducens. BMC Microbiol. 2009, 9: 109-10.1186/1471-2180-9-109.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302 (1): 205-217. 10.1006/jmbi.2000.4042.View ArticlePubMedGoogle Scholar
- Haveman SA, Holmes DE, Ding YH, Ward JE, Didonato RJ, Lovley DR: c-type cytochromes in Pelobacter carbinolicus. Appl Environ Microbiol. 2006, 72 (11): 6980-6985. 10.1128/AEM.01128-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Holmes DE, Nevin KP, Lovley DR: In situ expression of nifD in Geobacteraceae in subsurface sediments. Appl Environ Microbiol. 2004, 70 (12): 7251-7259. 10.1128/AEM.70.12.7251-7259.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Holmes DE, Nevin KP, O'Neil RA, Ward JE, Adams LA, Woodard TL, Vrionis HA, Lovley DR: Potential for quantifying expression of the Geobacteraceae citrate synthase gene to assess the activity of Geobacteraceae in the subsurface and on current-harvesting electrodes. Appl Environ Microbiol. 2005, 71 (11): 6870-6877. 10.1128/AEM.71.11.6870-6877.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Coppi MV, Leang C, Sandler SJ, Lovley DR: Development of a genetic system for Geobacter sulfurreducens. Appl Environ Microbiol. 2001, 67 (7): 3180-3187. 10.1128/AEM.67.7.3180-3187.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Marx CJ, Lidstrom ME: Broad-host-range cre-lox system for antibiotic marker recycling in gram-negative bacteria. Biotechniques. 2002, 33 (5): 1062-1067.PubMedGoogle Scholar
- Butler JE, Glaven RH, Esteve-Nunez A, Nunez C, Shelobolina ES, Bond DR, Lovley DR: Genetic characterization of a single bifunctional enzyme for fumarate reduction and succinate oxidation in Geobacter sulfurreducens and engineering of fumarate reduction in Geobacter metallireducens. J Bacteriol. 2006, 188 (2): 450-455. 10.1128/JB.188.2.450-455.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Dehio M, Knorre A, Lanz C, Dehio C: Construction of versatile high-level expression vectors for Bartonella henselae and the use of green fluorescent protein as a new expression marker. Gene. 1998, 215 (2): 223-229. 10.1016/S0378-1119(98)00319-9.View ArticlePubMedGoogle Scholar
- Marx CJ, Lidstrom ME: Development of improved versatile broad-host-range vectors for use in methylotrophs and other Gram-negative bacteria. Microbiology. 2001, 147 (Pt 8): 2065-2075.View ArticlePubMedGoogle Scholar
- Kunin V, Sorek R, Hugenholtz P: Evolutionary conservation of sequence and secondary structures in CRISPR repeats. Genome Biol. 2007, 8 (4): R61-10.1186/gb-2007-8-4-r61.PubMed CentralView ArticlePubMedGoogle Scholar
- Horvath P, Romero DA, Coute-Monvoisin AC, Richards M, Deveau H, Moineau S, Boyaval P, Fremaux C, Barrangou R: Diversity, activity, and evolution of CRISPR loci in Streptococcus thermophilus. J Bacteriol. 2008, 190 (4): 1401-1412. 10.1128/JB.01415-07.PubMed CentralView ArticlePubMedGoogle Scholar
- van der Ploeg JR: Analysis of CRISPR in Streptococcus mutans suggests frequent occurrence of acquired immunity against infection by M102-like bacteriophages. Microbiology. 2009, 155 (Pt 6): 1966-1976. 10.1099/mic.0.027508-0.View ArticlePubMedGoogle Scholar
- Mills S, Griffin C, Coffey A, Meijer WC, Hafkamp B, Ross RP: CRISPR analysis of bacteriophage-insensitive mutants (BIMs) of industrial Streptococcus thermophilus - implications for starter design. J Appl Microbiol. 2009Google Scholar
- Cui Y, Li Y, Gorge O, Platonov ME, Yan Y, Guo Z, Pourcel C, Dentovskaya SV, Balakhonov SV, Wang X, et al: Insight into microevolution of Yersinia pestis by clustered regularly interspaced short palindromic repeats. PLoS One. 2008, 3 (7): e2652-10.1371/journal.pone.0002652.PubMed CentralView ArticlePubMedGoogle Scholar
- Salzberg SL, Sommer DD, Schatz MC, Phillippy AM, Rabinowicz PD, Tsuge S, Furutani A, Ochiai H, Delcher AL, Kelley D, et al: Genome sequence and rapid evolution of the rice pathogen Xanthomonas oryzae pv. oryzae PXO99A. BMC Genomics. 2008, 9: 204-10.1186/1471-2164-9-204.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersson AF, Banfield JF: Virus population dynamics and acquired virus resistance in natural microbial communities. Science. 2008, 320 (5879): 1047-1050. 10.1126/science.1157358.View ArticlePubMedGoogle Scholar
- Tyson GW, Banfield JF: Rapidly evolving CRISPRs implicated in acquired resistance of microorganisms to viruses. Environ Microbiol. 2008, 10 (1): 200-207.PubMedGoogle Scholar
- Sissler M, Delorme C, Bond J, Ehrlich SD, Renault P, Francklyn C: An aminoacyl-tRNA synthetase paralog with a catalytic role in histidine biosynthesis. Proc Natl Acad Sci USA. 1999, 96 (16): 8985-8990. 10.1073/pnas.96.16.8985.PubMed CentralView ArticlePubMedGoogle Scholar
- Bond JP, Francklyn C: Proteobacterial histidine-biosynthetic pathways are paraphyletic. J Mol Evol. 2000, 50 (4): 339-347.PubMedGoogle Scholar
- Holt PJ, Morgan DJ, Sazanov LA: The location of NuoL and NuoM subunits in the membrane domain of the Escherichia coli complex I: implications for the mechanism of proton pumping. J Biol Chem. 2003, 278 (44): 43114-43120. 10.1074/jbc.M308247200.View ArticlePubMedGoogle Scholar
- Patzer SI, Hantke K: The ZnuABC high-affinity zinc uptake system and its regulator Zur in Escherichia coli. Mol Microbiol. 1998, 28 (6): 1199-1210. 10.1046/j.1365-2958.1998.00883.x.View ArticlePubMedGoogle Scholar
- Methé BA, Webster J, Nevin K, Butler J, Lovley DR: DNA microarray analysis of nitrogen fixation and Fe(III) reduction in Geobacter sulfurreducens. Appl Environ Microbiol. 2005, 71 (5): 2530-2538. 10.1128/AEM.71.5.2530-2538.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Han D, Krauss G: Characterization of the endonuclease SSO2001 from Sulfolobus solfataricus P2. FEBS Lett. 2009, 583 (4): 771-776. 10.1016/j.febslet.2009.01.024.View ArticlePubMedGoogle Scholar
- Beloglazova N, Brown G, Zimmerman MD, Proudfoot M, Makarova KS, Kudritska M, Kochinyan S, Wang S, Chruszcz M, Minor W, et al: A novel family of sequence-specific endoribonucleases associated with the clustered regularly interspaced short palindromic repeats. J Biol Chem. 2008, 283 (29): 20361-20371. 10.1074/jbc.M803225200.PubMed CentralView ArticlePubMedGoogle Scholar
- Wiedenheft B, Zhou K, Jinek M, Coyle SM, Ma W, Doudna JA: Structural basis for DNase activity of a conserved protein implicated in CRISPR-mediated genome defense. Structure. 2009, 17 (6): 904-912. 10.1016/j.str.2009.03.019.View ArticlePubMedGoogle Scholar
- Han D, Lehmann K, Krauss G: SSO1450--a CAS1 protein from Sulfolobus solfataricus P2 with high affinity for RNA and DNA. FEBS Lett. 2009, 583 (12): 1928-1932. 10.1016/j.febslet.2009.04.047.View ArticlePubMedGoogle Scholar
- Mojica FJ, Diez-Villasenor C, Garcia-Martinez J, Almendros C: Short motif sequences determine the targets of the prokaryotic CRISPR defence system. Microbiology. 2009, 155 (Pt 3): 733-740. 10.1099/mic.0.023960-0.View ArticlePubMedGoogle Scholar
- Agari Y, Sakamoto K, Tamakoshi M, Oshima T, Kuramitsu S, Shinkai A: Transcription profile of Thermus thermophilus CRISPR systems after phage infection. J Mol Biol. 2010, 395 (2): 270-281. 10.1016/j.jmb.2009.10.057.View ArticlePubMedGoogle Scholar
- Butler JE, Young ND, Lovley DR: Evolution of electron transfer out of the cell: comparative genomics of six Geobacter genomes. BMC Genomics. 2010, 11: 40-10.1186/1471-2164-11-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Esteve-Nunez A, Sosnik J, Visconti P, Lovley DR: Fluorescent properties of c-type cytochromes reveal their potential role as an extracytoplasmic electron sink in Geobacter sulfurreducens. Environ Microbiol. 2008, 10 (2): 497-505. 10.1111/j.1462-2920.2007.01470.x.View ArticlePubMedGoogle Scholar
- Katayama A, Tsujii A, Wada A, Nishino T, Ishihama A: Systematic search for zinc-binding proteins in Escherichia coli. Eur J Biochem. 2002, 269 (9): 2403-2413. 10.1046/j.1432-1033.2002.02900.x.View ArticlePubMedGoogle Scholar
- Matsuyama A, Yamamoto-Otake H, Hewitt J, MacGillivray RT, Nakano E: Nucleotide sequence of the phosphotransacetylase gene of Escherichia coli strain K12. Biochim Biophys Acta. 1994, 1219 (2): 559-562.View ArticlePubMedGoogle Scholar
- Smith KS, Ferry JG: Prokaryotic carbonic anhydrases. FEMS Microbiol Rev. 2000, 24 (4): 335-366. 10.1111/j.1574-6976.2000.tb00546.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.