- Open Access
Comparison of single-molecule sequencing and hybrid approaches for finishing the genome of Clostridium autoethanogenum and analysis of CRISPR systems in industrial relevant Clostridia
Biotechnology for Biofuels volume 7, Article number: 40 (2014)
Clostridium autoethanogenum strain JA1-1 (DSM 10061) is an acetogen capable of fermenting CO, CO2 and H2 (e.g. from syngas or waste gases) into biofuel ethanol and commodity chemicals such as 2,3-butanediol. A draft genome sequence consisting of 100 contigs has been published.
A closed, high-quality genome sequence for C. autoethanogenum DSM10061 was generated using only the latest single-molecule DNA sequencing technology and without the need for manual finishing. It is assigned to the most complex genome classification based upon genome features such as repeats, prophage, nine copies of the rRNA gene operons. It has a low G + C content of 31.1%. Illumina, 454, Illumina/454 hybrid assemblies were generated and then compared to the draft and PacBio assemblies using summary statistics, CGAL, QUAST and REAPR bioinformatics tools and comparative genomic approaches. Assemblies based upon shorter read DNA technologies were confounded by the large number repeats and their size, which in the case of the rRNA gene operons were ~5 kb. CRISPR (Clustered Regularly Interspaced Short Paloindromic Repeats) systems among biotechnologically relevant Clostridia were classified and related to plasmid content and prophages. Potential associations between plasmid content and CRISPR systems may have implications for historical industrial scale Acetone-Butanol-Ethanol (ABE) fermentation failures and future large scale bacterial fermentations. While C. autoethanogenum contains an active CRISPR system, no such system is present in the closely related Clostridium ljungdahlii DSM 13528. A common prophage inserted into the Arg-tRNA shared between the strains suggests a common ancestor. However, C. ljungdahlii contains several additional putative prophages and it has more than double the amount of prophage DNA compared to C. autoethanogenum. Other differences include important metabolic genes for central metabolism (as an additional hydrogenase and the absence of a phophoenolpyruvate synthase) and substrate utilization pathway (mannose and aromatics utilization) that might explain phenotypic differences between C. autoethanogenum and C. ljungdahlii.
Single molecule sequencing will be increasingly used to produce finished microbial genomes. The complete genome will facilitate comparative genomics and functional genomics and support future comparisons between Clostridia and studies that examine the evolution of plasmids, bacteriophage and CRISPR systems.
The development of next-generation DNA sequencing technologies since the first human genome sequence was completed has led to remarkable increases in sequencing efficiency on the order of approximately 100,000-fold . Costs have dropped dramatically and computational methods have advanced along with sequencing technology, leading to large increases in DNA sequencing output and in the number of available genome sequences [2, 3]. A variety of assembly algorithms and methods for quality evaluation have been developed [3–13]. However, the majority of sequenced genomes are incomplete due to technical difficulties, time, and expense leading to an increasing disparity between the number of finished and draft genomes in databases [1–3, 5, 14].
The PacBio sequencing system  is the only long-read, single-molecule sequencer available at present and the performance of the PacBio RS system was compared to two short-read sequencing platforms also released in 2011 . The original RS system with C1 chemistry generated mean read-lengths in the range of 1,500 bp and yielded approximately 100 Mb of sequence data per run, and reads in this range were useful in generating improved scaffolds for de novo assemblies. However, the original system was not optimal for de novo assembly applications  and hybrid assembly approaches have been developed to overcome limitations in short-read technologies and higher error rates associated with third-generation technology [17, 18].
Repetitive stretches of DNA are abundant and are one of the main technical challenges that hinder accurate sequencing and genome assembly efforts . In the case of bacteria, the rRNA gene operon is often the largest region of repetitive sequence and range in size between 5 and 7 kb . Last year, the longest PacBio RS reads were reported as being approximately 14 kb and these longer reads are useful in resolving repeats during genome assemblies . The PacBio RS II system was released last year and it produces more and longer reads. In a recent study, the longest read before correction was 15,634 bp and the genomes of six bacteria were sequenced and assembled using single-molecule sequencing based on C2 chemistry . Koren et al.  suggested that the majority of bacterial genomes could be assembled into finished-grade quality, that is, without gaps, and with data derived from a single PacBio sequencing library per sample . The combination of the longer reads, depth of coverage and random nature of sequencing errors facilitates de novo assemblies for microbial isolates [15, 20, 21]. The advantages of single-molecule sequencing have been discussed . To date, relatively few genomes sequences have been determined exclusively using single-molecule technology and only a handful represent finished genomes [14, 20, 21, 23–25].
In this study, a finished genome sequence for Clostridium autoethanogenum strain JA1-1 (DSM 10061) was generated using the latest PacBio RS II instrument. This represents one of the first de novo genomes finished into a single contiguous sequence using RS II data alone (that is, without addition of other next-generation sequence data or manual finishing steps). To offer insights into this technology, the PacBio assembly was compared to assemblies based on 454 GS FLX Titanium and Illumina MiSeq data and an earlier draft genome sequence of 100 contigs for this strain obtained from 454 GS FLX Titanium and Ion Torrent data .
C. autoethanogenum is an anaerobic, Gram-positive, mesophilic, acetogenic bacterium isolated using carbon monoxide (CO) . Other substrates include the greenhouse gas CO2 plus H2, pyruvate, xylose, arabinose, fructose, rhamnose, and L-glutamate. There is significant biotechnological interest in this organism as well as other acetogenic bacteria for their ability to use gases containing CO, H2 and CO2 as the sole source of carbon and energy for the production of fuel and chemicals at scale. The ability to use these gases in fermentative processes enables acetogens to potentially provide a route to more sustainable fuel and chemical production from a range of feedstocks including biomass and municipal solid waste-derived syngas, reformed biogas and industrial waste gases derived, for example, from steel production facilities [28–33].
Results and discussion
Sequencing output and assembly statistics for C. autoethanogenum DSM 10061
Sequencing statistics show that for each platform a large number of raw reads were attained that resulted in high degrees of genome coverage (Table 1). Raw Illumina data were trimmed and filtered before assembly, but in the case of the 454 and PacBio assemblers raw instrument output files were used. Bruno-Barcena et al. used a combination of 454 GS FLX Titanium and Ion Torrent Personal Genome Machine (PGM) data to generate a genome reported as 4.5 Mb for C. autoethanogenum DSM 10061 . The number of 454 reads (452,052) and genome coverage (39×) from the earlier study was similar to this one (Table 1), although addition of the PGM reads resulted in 905,738 raw reads being used to generate the preliminary assembly by Newbler (version 2.6). The record [GenBank: ASZX00000000.1] for strain DSM 10061 draft genome is reported as 4,323,309 bp.
In this study, Newbler (version 2.8) was used to assemble new 454 paired-end reads from a 3-kb insert length library (Table 1) into a draft genome sequence that consisted of 32 contigs (Table 2). The lower number of contigs (32 versus 100) from the new 454-only assembly compared to the draft version  is likely due to differences in library types (paired-end versus shotgun) and software versions. Assembly of Illumina-only data was conducted using the SPAdes , Velvet , Abyss  and the CLC Genomics Workbench (CLC Bio) assemblers and the best results were obtained by the Velvet assembler (Table 2). Previously, we have assembled genome sequences for a range of bacteria using a combination of 454 and Illumina technologies, whereby initial Illumina consensus sequences were shredded into 1.5-kbp overlapped fake reads and assembled together with the 454 data [37–42]. The best genome assembly obtained for strain DSM 10061 using second generation sequencing technologies employed such a hybrid approach, which is reflected in the lowest number of contigs, the largest single contig and highest N50 value (Table 2). Preliminary studies using the Clostridium ljungdahlii DSM 13528 genome as a reference and a PCR/Sanger sequencing strategy showed contigs could be joined by such an approach (Additional file 1). As manual finishing is time consuming the potential of PacBio data to generate finished microbial genome sequences was assessed.
Remarkably, one PacBio library preparation and two single molecule real-time sequencing (SMRT) cells produced sufficient sequence such that it could be assembled into one contiguous DNA fragment that represented the DSM 10061 genome. The PacBio genome assembly is a similar size to the other assemblies (Tables 1 and 2) and genome completeness was confirmed by sequence wrap-around. This is one of the first de novo sequenced genomes we are aware of that has been closed without manual finishing or additional data, despite the complexity of the C. autoethanogenum genome.
A comparison of the 454/Illumina hybrid assembly to the PacBio assembly showed there were small regions of overlap in the hybrid assembly that weakly joined contigs, and were supported by PCR and Sanger data, but there was insufficient support for the Newbler software to join them (Additional file 1A). PCR and Sanger data joined small gaps between contigs (for example, see Additional file 1B) in line with predictions using C. ljungdahlii DSM 13528 as a reference but in other examples much larger products were obtained compared to the predicted PCR product sizes (Additional file 1C). Other challenges involved using a related but different species, or strain from manual finishing included instances of software not being able to design PCR primers, not obtaining PCR products, and instances of obtaining multiple PCR products of different sizes and/or DNA smears.
Assembly quality assessments and comparisons
The complexity of the C. autoethanogenum DSM 10061 genome sequence was assessed and it is classified as a class III genome, according to previously described criteria for repeat sequence content and type . Class III genomes are defined as containing repeats that can include rRNA gene operons, many mid-scale repeats, such as insertion sequences and simple sequence repeats, and large phage-mediated repeats, duplications, or large tandem arrays that are considerably larger than the rRNA gene operon.
PacBio sequencing technology has a high error rate, which has been reported as being approximately 18% . Due to the random nature of the error , it is however, possible to get a highly accurate consensus sequence when there is high coverage [14, 20, 21]. For genomes such as C. autoethanogenu m with extreme guanine and cytosine (G + C) contents (31.1 mol% G + C content) and long homonucleotide stretches this provides an advantage over other sequencing technologies.
Beyond simple metrics, such as contig number, N50 and largest contig size, several bioinformatics approaches have been developed to assess assembly quality. The computing genome assembly likelihoods (CGAL) method is one recent approach that assesses uniformity of read coverage for assemblies and also evaluates the read errors, library insert size distribution and the degree of unassembled data . At present, CGAL is only able to utilize Illumina reads for its assembly assessment and using Illumina reads it ranked the assemblies in the order of best to worst as Illumina only, Illumina/454 hybrid, 454, published draft, to PacBio, respectively (Additional file 2). The CGAL likelihood principle is based on the possibility that a read is produced from every single location in the assembly. Regions of repetitive DNA were to be sequenced by longer reads, which were at times not resolved by the Illumina reads (Figure 1) and this may have contributed to the lower CGAL scores for assemblies that contained longer reads and no Illumina data. QUAST , which used the PacBio assembly as the reference, ranked the Illumina/454 hybrid, 454, published draft, and Illumina only assemblies in the order of best to worst, respectively and details are provided (Additional file 3).
The tool, recognising errors in assemblies using paired reads (REAPR) for genome assembly evaluation  detected no collapsed repeats in the PacBio assembly and five in the hybrid assembly and four in each of the other assemblies (Additional file 4). The fragment coverage distribution (FCD) error detected by REAPR in PacBio assembly was at location 3872494 to 3873407 (913 bp). This region contains an rRNA gene operon and had very low Illumina coverage (40× as compared to the average of 127×). Hence, REAPR reported an error (based on Illumina reads only). Even 454 coverage was low in this region (19× as compared to average of 46×). However, there was 108× PacBio reads covering this (913 bp) region and for the first 392 bp there was also high-quality Sanger sequence support indicating it is unlikely that there is an issue for the PacBio assembly in this region. The hybrid and PacBio assemblies contained the fewest warnings (83 and 96, respectively), followed by the Illumina assembly (182) and then published draft assembly contained the most (190).
A multiple genome alignment was conducted by aligning contigs from the different assemblies to the PacBio reference assembly to identify conserved regions and to evaluate gaps in the different DSM 10061 assemblies. Regions with no or partial 454 or Illumina contig coverage predominantly contained predicted rRNA gene operons and other duplicated genes (Figure 1 and Table 3). While the draft genome sequence for strain DSM 10061 predicts one copy of the 16S rRNA gene , nine rRNA clusters were predicted using the DSM 10061 PacBio assembly, which is the same number of rRNA operons as in the closely related C. ljungdahlii DSM 13528 . Based on findings in this study and earlier ones [1, 3, 14], the large number of DSM 10061 rRNA clusters and their repetitive nature confounded assembly of the shorter reads.
The latest PacBio RS II SMRT cells are designed to select for larger read-lengths when long insert libraries (10 to 20 kb) are being prepared, however, preferential loading of smaller fragments can still occur and this limits sequence output. In this study, smaller fragments were removed from the PacBio library by size exclusion leading to longer read-lengths and greater amounts of sequence data than otherwise might have been attained. The long reads produced by the new PacBio RS II system, combined with sequence depth meant that the principal regions of complexity could be resolved using one library preparation and two SMRT cells to generate a complete genome sequence. The application of long, single-molecule sequencing data will lead to a greater number of finished genomes and quality improvements in microbial genome databases , however the application of the newest version of this technology requires more evaluation before its full potential can be assessed for complex genomes.
General features of the C. autoethanogenum genome, its metabolism and comparison to C. ljungdahlii
The finished genome of C. autoethanogenum DSM 10061 consists of one chromosome of 4,352,205 bp in size with a G + C content of 31.1 mol% and consists of 89 RNA genes (Additional file 5). Of the 4,161 genes predicted for this strain, 4,042 are protein-coding genes (CDSs) and 18 are pseudogenes. The distribution of genes into COG functional categories is presented (Additional file 6). The previously published draft DSM 10061 genome annotation included 4,135 predicted coding sequences  and the related finished C. ljungdahlii DSM 13528 genome which is 277,860 bp larger in size contained 4,184 protein coding genes . Predicted gene content differences reflect the use of different gene-calling algorithms, that draft sequences can split genes in two, and genotypic differences. The methodology, accuracy, and specificity of the Prodigal gene prediction algorithm used in this study has been described previously .
Phenotypic and metabolic differences have been reported for C. autoethanogenum and C. ljungdahlii[27, 44–47]. The two are indistinguishable at the 16S rRNA gene level  and have high scores for similarity based on in silico average nucleotide identity comparisons across the genomes (0.9977 ANIb) . To evaluate potential coding sequence differences between the two organisms, OrthoMCL , a genome-scale algorithm for grouping orthologous protein sequences, was used to compare all the C. autoethanogenum proteins to those in C. ljungdahlii and for the reciprocal evaluation (Additional file 7). A general description for all OrthoMCL proteins is provided (Additional file 7: Table S1). Putative paralogs were identified (1_taxa tab in Additional file 7) along with putative orthologs (2_taxa file tab in Additional file 7). Proteins without orthologs or paralogs were identified using the default settings (C. autoethanogenum unique or C. ljungdahlii unique tabs in Additional file 7). This analysis revealed that over 10% of the proteome is unique to each bacterium when comparing C. autoethanogenum (427 proteins out of 4,134) and C. ljungdahlii (447 out of 4,198). The 427 proteins with unique genes to DSM 10061 (as listed by OrthoMCL) were searched against the entire C. ljungdahlii proteome using BLASTP and an e-value similarity criteria of 1e-5 to identify proteins with truly unique function and no homolog, which reduced the number of dissimilar or unique proteins to 221 (BLASTP analysis tab in Additional file 7). From the proteins identified as unique to each bacterium, the majority were proteins with hypothetical functions or proteins related to particular phage, transposon or CRISPR sequences, but proteins with key functions in the metabolism were also identified that could explain different phenotypes. These differences are discussed below.
The Wood-Ljungdahl pathway (Figure 2) plays a key role in the acetogenic metabolism by allowing the formation of acetyl-CoA from CO or CO2, and thus, is essential for autotrophic growth. Under heterotrophic growth conditions it permits utilization of produced CO2 and reducing equivalents generated during glycolysis to form an additional molecule of acetyl-CoA . The genes encoding for the enzymes of the Wood-Ljungdahl pathway are co-localized in one large cluster (CAETHG_ 1606–1621). The same organization is also found in other acetogens such as C. ljungdahlii, C. ragsdalei or C. difficile, but significant differences at the sequence level are present as described earlier [28, 51]. This cluster also includes the genes for the bifunctional carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS) enzyme complex, the key enzyme in the Wood-Ljungdahl pathway. As in C. ljungdahlii, two additional monofunctional carbon monoxide dehydrogenases (CAETHG_3005 and CAETHG_3899) are encoded in the genome of C. autoethanogenum that may also be involved in utilization of CO and CO2. Although CO can be both a carbon and an energy source for the bacteria, CO2 can only be used as a carbon source. Additional energy can be generated from hydrogen, via hydrogenase enzymes. The genome of C. autoethanogenum encodes for six hydrogenases, one (NiFe) hydrogenase and five (FeFe) hydrogenases. Interestingly, C. ljungdahlii only has five hydrogenases, lacking one of the iron-only hydrogenases that are present in C. autoethanogenum. The genes for this unique (FeFe) hydrogenase are in an operon with two genes for NuoF-like oxidoreductases (CAETHG_1575-78). The presence of an additional hydrogenase enzyme complex could represent a significant advantage for C. autoethanogenum during autotrophic growth on CO, CO2 and H2 containing gases. Preliminary RNA-Seq experiments show that this cluster is highly expressed under such conditions, underlining the importance of this enzyme (Additional file 8). Of the other C. autoethanogenum hydrogenases, a second (FeFe) hydrogenase gene cluster was also found to be highly expressed. This nicotinamide adenine dinucleotide phosphate-oxidase (NADPH)-specific electron-bifurcating Hyt hydrogenase was recently characterized and found to form a functional complex with a formate dehydrogenase . Formate dehydrogenase activates CO2 to formate in the Wood-Ljungdahl pathway and additional formate dehydrogenase genes are present in the C. autoethanogenum genome (Figure 2). C. autoethanogenum also has a predicted formate transporter (CAETHG_1601) that is not present in C. ljungdahlii.
During autotrophic growth, all biomass and products must be derived of acetyl-CoA from the Wood-Ljungdahl pathway (Figure 2). Fatty acid biosynthesis starts directly from acetyl-CoA, whereas production of nucleic acids, amino acids, vitamins, cofactors and secondary metabolites proceed via pyruvate and gluconeogenesis or the TCA cycle. The C. autoethanogenum genome encodes for two pyruvate:ferredoxin oxidoreductases (PFOR) that catalyze the conversion of acetyl-CoA into pyruvate. C. autoethanogenum has a pyruvate, phosphate dikinase (PPDK), but interestingly no phophoenolpyruvate synthase (PPSA) and C. ljungdahlii contains two such enzymes (CLJU_c14340 and CLJU_c38600). The rest of gluconeogenesis is similar in both organisms. C. autoethanogenum has an incomplete TCA cycle to succinate and 3-oxogluterate (Figure 2). C. autoethanogenum products such as 2,3-butanediol and lactate are also derived from pyruvate , whereas ethanol is produced from acetyl-CoA via acetaldehyde, either directly via bifunctional aldehyde/alcohol dehydrogenases or via acetate using phosphotransacetylase and acetate kinase and an aldehyde:ferredoxin oxidoreductase (Figure 2). Several additional alcohol dehydrogenases are present in the genome of C. autoethanogenum.
Heterotrophic growth on a range of other products such as a range of C5 and C6 sugars has been described for C. autoethanogenum. A PTS system and other respective genes could be identified in the genome (Figure 2). In contrast to C. ljungdahlii, some extra genes involved in mannose metabolism are present in C. autoethanogenum as well as genes for aromatic compound degradation. C. autoethanogenum also has an additional predicted nitrate reductase (CAETHG_0085) and both organisms differ in some of their transport systems.
Other differences between C. autoethanogenum and C. ljungdahlii include variations in the sporulation program, with several unique predicted sporulation proteins and regulators present in C. autoethanogenum strain DSM 10061, and different defense systems, such as restriction/methylation systems and a CRISPR system that is present in C. autoethanogenum but not in C. ljungdahlii. Insertion sequence (IS) elements, are usually unique to a strain, and one is found in C. autoethanogenum between 4,345,780 and 4,347,448 bp that is 100% identical to one in C. ljungdahlii. C. autoethanogenum DSM 10061 was enriched from rabbit feces in Belgium  and C. ljungdahlii DSM 13528 was isolated in the US from chicken yard waste . Despite the geographical separation of the isolates, the overall degree of similarity between C. ljungdahlii and C. autoethanogenum suggests a common ancestor.
C. autoethanogenum CRISPR system
CRISPR are prokaryotic DNA loci that carry the memory of past bacterial infections of phages and plasmids to provide immunity against mobile genetic elements [53, 54]. In the last decade, several studies have unraveled CRISPR defense molecular details and mechanisms of action [53, 55, 56]. Briefly, CRISPR loci are composed of arrays of 24 to 47 bp partially palindromic, highly conserved repeats separated by variable spacers specific to the infecting DNA. CRISPR-associated (cas) genes are involved in spacer acquisition, expression and interference to phage or plasmid. cas gene operons are classified into three types and several subtypes, and can target either DNA or RNA, or both . CRISPR and cas gene operons are proposed to be transferred between distinctly related strains by horizontal gene transfer and/or by transposons , and the latter can be identified by the presence of insertion elements and transposase/mutase in its vicinity. Thus, CRISPR appear to be dynamic heritable defense systems in bacteria against plasmids and phages that are ever fast-evolving and play important roles in the co-evolution of both bacteria and phages.
The genome of C. autoethanogenum is found to contain eight cas genes of Type-I B, all predicted to be in one operon on the antisense strand with a predicted transcription terminator at the end of the cas2 gene, and it is flanked by three CRISPR arrays (Table 4, Additional files 9 and 10) with a total of 93 30-bp-repeats (consensus 5′-GTTGAACCTCAACATGAGATGTATTTAAAT-3′) and 90 spacers of 35 to 38 bp. All three CRISPR arrays are preceded by a 177-bp-leader sequence, which is essential for array transcription and a fragment of the leader sequence is co-transcribed with the array [58, 59]. The three putative C. autoethanogenum CRISPR arrays leader sequences share 82 to 91% sequence similarity between them (Additional file 9). Interestingly, 10 kb downstream of the three CRISPR arrays, an incomplete leader sequence of 65 bp was found that has high sequence identity to the other leaders in close proximity (100 bp) to an imperfect CRISPR repeat (5′-GTTGAACCTtAACATGAGATGTAaaggtAa-3′). In addition to the three CRISPR arrays flanking the cas genes, a putative extra CRISPR array was identified in the genome, consisting of three 55-bp-repeats and two 16-bp-spacer (Additional file 10).
Expression of cas genes and CRISPR arrays along with their leader sequence were studied by Reverse Transcriptase PCR (RT-PCR) and RNA-Seq during logarithmic growth under autotrophic conditions. PCR amplification of fragments of expected sizes were observed only with cDNA template and not with RNA, showing the absence of genomic DNA contamination in RNA preparations (Additional file 9). All eight predicted cas genes appear to be co-expressed and from a single operon. Expression of spacers distal to the leader sequences in all three arrays was also assessed. Based on sequence similarity between the three leader sequences, a common reverse primer was designed to align to the conserved region in the leader sequences of all three arrays and forward primers aligning specifically to spacers proximal to the leader sequence in each array (Additional file 9). Both the leader proximal and distal spacers of array 2 were found to be expressed, whereas expression was detected only for spacers proximal to the leader sequence in array 3 and no expression was detected from array 1 or the identified extra leader. Preliminary RNA-Seq data showed expression of all CRISPR RNAs (crRNAs), with different abundances. A few transcripts corresponding to the leader region of the three CRISPR arrays were also detected. Similar to the cas gene operon, the three CRISPR arrays including their leader regions were transcribed from the antisense strand. The three CRISPR arrays were found to be constitutively transcribed and processed into crRNAs of varying lengths (Additional file 8). However, the processed crRNAs appeared to have a well-defined eight-nucleotide 5′ handle, 5′-ATTTAAAT-3′, originating from the repeat region followed by the spacer sequence (Additional file 8). The 3′ end of these processed crRNAs had varying tags (Additional file 8). Based on these findings, a scheme for CRISPR processing in C. autoethanogenum is proposed (Additional file 8). Similar processing of crRNA was observed across the three samples collected at different time points.
CRISPR spacer sequences in C. autoethanogenum were analyzed to identify potential target DNA sequence. A BLAST search did not result in high identity hits within the National Center for Biotechnology Information (NCBI) database or against its own genome. A comparison of regions of DNA from putative C. autoethanogenum processing crRNAs from all three arrays identified the sequence 5′-ATTTAAAT-3′ (Additional file 9), which is similar to sequences from Clostridium thermocellum, Methanococcus maripaludis, Escherichia coli and Pyrococcus furiosus, which also have type-IB CRISPR systems. In these organisms the processing of crRNA is mediated by the cas6 gene, which is also found in C. autoethanogenum. Unlike in C. thermocellum, M. maripaludis, Sulfolobus acidocaldarius and P. furiosus, C. autoethanogenum crRNAs were transcribed only from the antisense strand and no anti-crRNA transcripts originating from the complementary strand were detected.
Identification and classification of CRISPR systems in industrial relevant clostridia
The presence of a CRISPR system in C. autoethanogenum compared to C. ljungdahlii could provide an advantage in industrial fermentations. The C. autoethanogenum CRISPR system was compared to those from other industrial relevant Clostridia strains to better understand their characteristics and their potential physiological and applied roles. In particular, the Clostridial ABE fermentation process has a history of phage infections . CRISPR systems from 14 Clostridium species were examined for the first time including those used in ABE fermentation processes: C. acetobutylicum, C. beijerinckii, C. saccharobutylicum and C. saccharoperbutylicum, cellulose degrading C. thermocellum, C. cellulolyticum, C. cellulovorans, and C. phytofermentans, and the acetogens C. autoethanogenum, C. ljungdahlii and C. carboxidivorans. CRISPR elements were identified only in 8 of the 14 Clostridium species analyzed by PILER  and CRISPRdb . All of the loci were found on chromosomes and none on any plasmids or megaplasmids (Table 4). From the ABE fermentation-Clostridia examined, only C. saccharobutylicum DSM13864 has a CRISPR system, but not several strains of the more commonly used C. acetobutylicum, C. beijerinckii and C. saccharobutylicum. This may be one of the reasons why the ABE fermentation process was historically found to be prone to phage infections . From the three acetogenic strains investigated only C. autoethanogenum had a CRISPR system, whereas all analyzed cellulolytic Clostridia, but C. phytofermentans contain CRISPR systems.
In all Clostridium species that harbor CRISPR arrays, cas genes were identified. C. cellulolyticum and C. thermocellum had two and four different cas operons, respectively (Additional files 10 and 11). These cas operons were classified based on a recently proposed classification system  and their target molecule(s) inferred. A phylogenetic analysis of cas1 genes was performed and compared to the 16S rRNA phylogeny (Additional file 12). In C. cellulolyticum, arrays 1 and 2 are associated with the Type I-C cas system and Adb with the Type II cas system. The two arrays are separated by a transposase and mutase genes (Additional file 11). C. cellulolyticum has two different sets of cas genes, both of which appear to target DNA. The C. cellulovorans cas operon could not be classified according by these criteria, nor could the target of its cas genes be inferred. C. thermocellum appears to have a Type III cas genes system (Additional file 11). The Type III cas system contains more than one type of cas gene operon belonging to either Type I or II or the repeat-associated mysterious proteins (RAMP) module operon and are predicted to target both DNA and RNA . The arrays 3 and 4 in C. thermocellum are associated with the Type I-B cas system  and arrays 5 and 6 to a cas system similar to Type I-B but interrupted by insertion of multiple other genes that separate cas1, cas2, cas4 genes (possibly involved in spacer acquisition) from cas3, 5, 7 and 8b (predicted to be involved in DNA interference) (Additional file 11). The array 1 is not associated with any cas gene cassette and is flanked by a integrase and mutase genes at the 3′ end. Apart from these two cas systems, C. thermocellum also has a RAMP module gene cassette that may be involved in RNA interference. This cassette is not associated with any array and could be acting in trans. As in a few lactic acid bacteria , integrase and mutase genes were frequently found near the cas gene cassette, particularly flanking the Type 1 cas gene cassette found in C. cellulovorans, C. thermocellum and C. autoethanogenum, suggesting possible horizontal gene transfer. These genes were also found next to array-1 in C. thermocellum.
The C. autoethanogenum CRISPR repeat DNA was not found in any of the other Clostridium species included in this study. A search for organisms with repeats similar to C. autoethanogenum in the CRISPRdb database  resulted in Clostridium novyi, Eubacterium limosum, along with a few Clostridium botulinum substrains. A comparison of the repeat sequences showed very high sequence similarity (Additional file 12). The cas genes operon in C. autoethanogenum, C. novyi and E. limosum were all of Type-I B, whereas the C. botulinum substrains had different sets of cas genes. The cas gene operon architecture, the arrangement of arrays on the chromosome and the presence of two hypothetical genes separating arrays 2 and 3 in C. autoethanogenum and C. novyi are strikingly alike, suggesting a common lineage of these two CRISPR-cas systems. This observation was further strengthened by the phylogenetic classification placing C. autoethanogenum cas1 gene together with cas1 genes from C. novyi and E. limosum and apart from the other Clostridium species (Additional file 12). Even though the repeat and the cas genes operon in C. autoethanogenum, C. novyi and E. limosum are largely identical, no similarity was found between the spacers.
Comparison of strains with/without the CRISPR system to plasmid and prophage content
Correlation between the presence of CRISPR and the occurrence of prophages or plasmids has been reported . To assess the 14 Clostridium species for this correlation their genome sequences were analyzed for presence of potential prophage regions. In C. autoethanogenum, four putative prophages were identified, an incomplete prophage similar to a Singapore grouper iridovirus, an intact prophage similar to a Geobacillus E2 virus and two intact prophages inserted into tRNAs (Additional file 10). One prophage was identified a Trp-tRNA and the other in an Arg-tRNA. The latter is in almost identical form also present in C. ljungdahlii, suggesting a shared lineage. Prophage regions were detected in all species irrespective of the presence of CRISPR modules (Table 4 and Additional file 10). Although there seems to be no general trend and it cannot be determined whether a prophage infection occurred before or after a CRISPR system was acquired, in a few cases bacteria that lacked CRISPR systems appeared to have more abundant prophage sequences.
When looking for plasmid content, only one out of seven strains containing CRISPR systems was found to contain a plasmid. Likewise, only one out of five plasmid-carrying strains contained a CRISPR system. CRISPR-mediated immunity has been shown experimentally to block conjugative plasmid acquisition , although the role of CRISPR in driving plasmid and phage evolution for industrially relevant Clostridia and other microorganisms remains to be fully elucidated.
A comparative genomic analysis revealed short-read technologies were unable to overcome C. autoethanogenum DSM 10061 repeat regions largely associated with nine copies of the rRNA gene operons. A previous study suggested that long single-molecule reads are sufficient to assemble most known microbial genomes based on a bioinformatics analysis of 2,267 complete genomes for bacteria and archaea and sequencing results for six bacteria . The genome sequence of C. autoethanogenum DSM 10061 is classified as within the most complex class of bacterial genomes and a complete genome sequence was generated for it using long single-molecule reads and without the need for manual finishing. The relatively low cost to generate the PacBio data (approximately US$1,500) and the outcome of this study support the assertion this technology will be valuable in future studies where a complete genome sequence is important and for complex genomes that contain large repeat elements.
Clostridia are known for their substrate and metabolic flexibility, which makes them attractive biocatalysts for biofuel and biorefinery applications . Acetogenic Clostridia, such as C. autoethanogenum, are of interest due to their ability to ferment abundant syngas or waste gases to useful products . The C. autoethanogenum genome sequence will facilitate strain development for biofuels and biochemicals production and comparative genomics in the future. A comparison between C. autoethanogenum and C. ljungdahlii identified distinct differences, notably the presence of a CRISPR system, an additional C. autoethanogenum hydrogenase, and several differences in central metabolism, although the two bacteria likely descend from a common ancestor. Comparative genomic analysis and characterization of CRISPR, plasmid content and prophage among Clostridia with biotechnological interest was performed. Notably, the classic ABE fermentation strains C. acetobutylicum and C. beijeinckii are reported to be prone to bacteriophage infections  and all lack a CRISPR system and only one of the analyzed 14 strains contain both a plasmid and a CRISPR system. From the acetogenic Clostridium strains sequenced to date, only C. autoethanogenum possesses a CRISPR system. Further consideration of Clostridia CRISPR systems may be informative for bioprocess development strategies and for ecological studies.
DNA sequence data generation
C. autoethanogenum strain JA1-1 was obtained from the Deutsche Sammlung von Mikroorganismen und Zellkulturen (DSMZ) culture collection (DSM 10061). C. autoethanogenum strain JA1-1 was cultured in PETC medium as described . Single colony was purified and 16S rDNA sequence confirmed before genomic DNA was prepared. High molecular weight genomic DNA was prepared as described earlier , quantified with a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and quality was assessed with Agilent Bioanalyzer (Agilent, Santa Clara, CA, USA).
Pyrosequencing was conducted using the Roche 454 GS FLX System (Roche 454 Life Sciences, Branford, CT, USA) with the method of paired-end DNA library preparation and average insert sizes in the 3-kb range and Titanium chemistry, according to the manufacturer’s instructions as described previously [38, 80]. Sequence data were also generated using a MiSeq instrument (Illumina, San Diego, CA, USA)  and a paired-end approach with an approximate insert library size of 500 bp and read lengths of 151 bp, as described previously  and according to the manufacturer’s instructions. DNA for PacBio sequencing was sheared with G-tubes (Covaris, Inc., Woburn, MA, USA), targeting 20-kb fragments. PacBio libraries were prepared with the DNA Template Prep Kit 2.0 (Pacific Biosciences, Menlo Park, CA, USA) and library fragments above 4 kb were isolated using the Blue Pippin system (Sage Science, Inc., Beverly, MA, USA). The average PacBio library insert size (including adapters) was approximately 19 kb and samples were sequenced using Magbead loading, C2 chemistry, Polymerase version P4, and software version 2.02. Raw next-generation sequence data available through the NCBI SRA database [SRX352885; SRX352888; SRP030033]. PCR and Sanger sequencing were conducted using standard approaches as described previously  and primer sequences are described (Additional file 1).
Sequence data trimming, filtering, annotation and assembly
The CLC Genomics Workbench (version 6.0.2) (CLC bio, Cambridge, MA, USA) was used to trim and filter Illumina reads for quality sequence data and the subsequent Illumina assembly. The Newbler application (version 2.8) in the 454 GS FLX software package (Roche 454 Life Sciences) was used to assemble reads generated from the GS FLX instrument and in combination with reads from the Illumina instrument, as described previously . The consensus Illumina sequences were processed before inputting into the Newbler assembler by generating 1.5-kb overlapping fake reads using the fb_dice.pl script, which is part of the FragBlast module (http://www.clarkfrancis.com/codes/fb_dice.pl). The PacBio reads were assembled through SMRTanalysis v 2.0 (Pacific Biosciences) using the HGAP protocol . The DSM 10061 PacBio assembly was annotated using the Prodigal gene-calling algorithm  and deposited in the NCBI database [GenBank: CP006763].
Assessment of genome assembly quality
The in silico evaluation of genome assemblies was performed using CGAL (version 0.9.6) , REAPR (version 1.0.16) , QUAST (version 2.2)  and Circos . The genomic repeats were identified using Nucmer ; genome complexity was determined based on count and length of the repeats as suggested earlier . Gaps in the 454/Illumina hybrid and published draft assemblies were determined by performing multiple genome alignment through Mauve (version 2.3.1)  with PacBio assembly used as reference genome. The order of contigs in 454/Illumina hybrid assembly and alignment of Sanger sequences was determined using Genious software (version 6.1.5) (Biomatters, Auckland, New Zealand).
Analysis, classification and comparison of CRISPR, plasmid, and prophage content in C. autoethanogenum and other fuel-producing Clostridia
The genome of C. autoethanogenum (NC_022592) and genome sequences of C. acetobutylicum ATCC824 (NC_003030), DSM1731 (NC_015687) and EA2018 (NC_017295), C. beijerinckii NCIMB8052 (NC_009617), C. saccharobutylicum (NC_022571), C. saccharoperbutylacetonicum (NC_020291), C. cellulolyticum H10 (NC_011898), C. cellulovorans 743B (NC_014393), C. thermocellum ATCC27405 (NC_009012) and DSM1313 (NC_017304), C. phytofermentans Isdg (NC_010001), C. ljungdahlii DSM13528 (NC_014328) and C. carboxidivorans (ACVI01000000; ADEK01000000) were retrieved from NCBI Genbank. The genome sequences of all these organisms were analyzed for CRISPR repeats using the PILER algorithm  and CRISPRdb . The sequence of plasmids found in C. acetobutylicum ATCC824 (NC_001988), DSM1731 (NC_015686 and NC_015688) and EA2018 (CP002119), C. saccharoperbutylacetonicum (NC_020292), and C. carboxidivorans (NC_014565) were also analyzed for the presence of CRISPR loci. The repeat sequences detected by PILER and CRISPRdb were combined and manually compared at sequence level. Degenerated CRISPR repeats with several mismatches were not taken into account. The genome sequences of these species were also analyzed for prophage regions using PHAST , Phage_Finder  and the outputs program manually analyzed. Multiple sequence alignment of repeats and their sequence logos were generated using CLUSTALW and Weblogo. Phylogenetic analyses based on 16S rRNA and cas1 genes were made using Geneious software. The phylogenetic tree was constructed using the neighbor-joining method with 100 bootstrap steps.
RT-PCR was performed to study the expression and operon structure of cas genes and the expression CRISPR arrays. Briefly, RNA was isolated (RNeasy Mini Kit, Qiagen, Valencia, CA, USA) from 20 mL C. autoethanogenum culture growing in serum bottles at an optical density (OD)600 of 0.2 in PETC media and steel mill waste gas (composition: 42% CO, 36% N2, 20% CO2, and 2% H2; collected from a New Zealand Steel site in Glenbrook, New Zealand) as the sole energy and carbon source. cDNA was synthesized using 500 ng of DNaseI (Ambion Inc., Austin, TX, USA)-treated RNA, Superscript III reverse transcriptase and random primers (Life Technologies, Grand Island, NY, USA). PCR was set with 30 ng cDNA- and DNaseI-treated RNA (control) as templates and iproof DNA polymerase (Biorad, Hercules, CA, USA). The primers used in this study are listed (Additional file 9).
RNA-Seq was performed from C. autoethanogenum growing in continuous culture in a 1.5-L continuous-stirred tank reactor (CSTR) with steel mill waste gas (composition: 42% CO, 36% N2, 20% CO2, and 2% H2; collected from a New Zealand Steel site in Glenbrook, New Zealand) as the sole energy and carbon source as described previously . A 20-ml sample was centrifuged at 4,000 × rpm for 10 minutes at 4°C. The supernatant was discarded and the pellet was stabilized adding 5 ml of RNAlater® (Ambion Inc). Total RNA was isolated from the cell pellet using RiboPureTM-Bacteria Kit (Ambion Inc.) according to the manufacturer’s standard protocol. DNA was removed using the TURBO DNA-free kit (Ambion Inc.) and RNA quality was assessed using a 2100 bioanalyzer (Agilent Technologies). RNA concentration was determined with a nanodrop 2000 (Thermo Fischer Scientific, Waltham, MA, USA). Ribodepletion was conducted using MICROBEExpressTM kit (Ambion Inc.). cDNA libraries were prepared and sequenced by standard procedures using SOLiD2 sequencing technology. Output from the SOLID run was processed in LifeScope v2.5.1., as specified by the manufacturer (http://downloads.lifetechnologies.com/Analysis_Software/GS/LifeScope/v2.5.1/LifeScope-v2.5.1_4476538_AUG.pdf). Processed reads were mapped to a reference assembly, and the resulting BAM files were imported, displayed, and manually inspected in the Geneious genome browser, v7.0.3 (Biomatters).
Steven D. Brown: corresponding author for DNA sequencing and comparative genomics
Michael Köpke: corresponding author for Clostridium autoethanogenum
computing genome assembly likelihoods
carbon monoxide dehydrogenase/acetyl-CoA synthase
Clusters of Orthologous Groups
clustered regularly interspaced short palindromic repeats
fragment coverage distribution
nicotinamide adenine dinucleotide phosphate-oxidase
Personal Genome Machine
repeat-associated mysterious proteins
reverse transcriptase PCR
single molecule real-time sequencing
tricarboxcylic acid cycle.
Treangen TJ, Salzberg SL: Repetitive DNA and next-generation sequencing: computational challenges and solutions. Nat Rev Genet 2011, 13: 36-46.
Chain PS, Grafham DV, Fulton RS, Fitzgerald MG, Hostetler J, Muzny D, Ali J, Birren B, Bruce DC, Buhay C, Cole JR, Ding Y, Dugan S, Field D, Garrity GM, Gibbs R, Graves T, Han CS, Harrison SH, Highlander S, Hugenholtz P, Khouri HM, Kodira CD, Kolker E, Kyrpides NC, Lang D, Lapidus A, Malfatti SA, Markowitz V, Metha T, et al.: Genomics. Genome project standards in a new era of sequencing. Science 2009, 326: 236-237. 10.1126/science.1180614
Nagarajan N, Pop M: Sequence assembly demystified. Nat Rev Genet 2013, 14: 157-167. 10.1038/nrg3367
Salzberg SL, Phillippy AM, Zimin A, Puiu D, Magoc T, Koren S, Treangen TJ, Schatz MC, Delcher AL, Roberts M, Marcais G, Pop M, Yorke JA: GAGE: A critical evaluation of genome assemblies and assembly algorithms. Genome Res 2012, 22: 557-567. 10.1101/gr.131383.111
Mavromatis K, Land ML, Brettin TS, Quest DJ, Copeland A, Clum A, Goodwin L, Woyke T, Lapidus A, Klenk HP, Cottingham RW, Kyrpides NC: The fast changing landscape of sequencing technologies and their impact on microbial genome assemblies and annotation. PLoS ONE 2012, 7: e48837. 10.1371/journal.pone.0048837
Kisand V, Lettieri T: Genome sequencing of bacteria: sequencing, de novo assembly and rapid analysis using open source tools. BMC Genomics 2013, 14: 211. 10.1186/1471-2164-14-211
Magoc T, Pabinger S, Canzar S, Liu X, Su Q, Puiu D, Tallon LJ, Salzberg SL: GAGE-B: an evaluation of genome assemblers for bacterial organisms. Bioinformatics 2013, 29: 1718-1725. 10.1093/bioinformatics/btt273
Medini D, Serruto D, Parkhill J, Relman DA, Donati C, Moxon R, Falkow S, Rappuoli R: Microbiology in the post-genomic era. Nat Rev Micro 2008, 6: 419-430.
Brown SD, Gilmour CC, Kucken AM, Wall JD, Elias DA, Brandt CC, Podar M, Chertkov O, Held B, Bruce DC, Detter JC, Tapia R, Han CS, Goodwin LA, Cheng J-F, Pitluck S, Woyke T, Mikhailova N, Ivanova NN, Han J, Lucas S, Lapidus AL, Land ML, Hauser LJ, Palumbo AV: Genome sequence of the mercury-methylating strain Desulfovibrio desulfurican s ND132. J Bacteriol 2078–2079, 2011: 193.
Vezzi F, Narzisi G, Mishra B: Feature-by-feature – evaluating de novo sequence assembly. PLoS One 2012, 7: e31002. 10.1371/journal.pone.0031002
Hunt M, Kikuchi T, Sanders M, Newbold C, Berriman M, Otto T: REAPR: a universal tool for genome assembly evaluation. Genome Biol 2013, 14: R47. 10.1186/gb-2013-14-5-r47
Gurevich A, Saveliev V, Vyahhi N, Tesler G: QUAST: quality assessment tool for genome assemblies. Bioinformatics 2013, 29: 1072-1075. 10.1093/bioinformatics/btt086
Rahman A, Pachter L: CGAL: computing genome assembly likelihoods. Genome Biol 2013, 14: R8. 10.1186/gb-2013-14-1-r8
Koren S, Harhay G, Smith T, Bono J, Harhay D, Mcvey S, Radune D, Bergman N, Phillippy A: Reducing assembly complexity of microbial genomes with single-molecule sequencing. Genome Biol 2013, 14: R101. 10.1186/gb-2013-14-9-r101
Eid J, Fehr A, Gray J, Luong K, Lyle J, Otto G, Peluso P, Rank D, Baybayan P, Bettman B, Bibillo A, Bjornson K, Chaudhuri B, Christians F, Cicero R, Clark S, Dalal R, Dewinter A, Dixon J, Foquet M, Gaertner A, Hardenbol P, Heiner C, Hester K, Holden D, Kearns G, Kong X, Kuse R, Lacroix Y, Lin S: Real-time DNA sequencing from single polymerase molecules. Science 2009, 323: 133-138. 10.1126/science.1162986
Quail M, Smith M, Coupland P, Otto T, Harris S, Connor T, Bertoni A, Swerdlow H, Gu Y: A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers. BMC Genomics 2012, 13: 341. 10.1186/1471-2164-13-341
Bashir A, Klammer AA, Robins WP, Chin CS, Webster D, Paxinos E, Hsu D, Ashby M, Wang S, Peluso P, Sebra R, Sorenson J, Bullard J, Yen J, Valdovino M, Mollova E, Luong K, Lin S, LaMay B, Joshi A, Rowe L, Frace M, Tarr CL, Turnsek M, Davis BM, Kasarskis A, Mekalanos JJ, Waldor MK, Schadt EE: A hybrid approach for the automated finishing of bacterial genomes. Nat Biotechnol 2012, 30: 701-707. 10.1038/nbt.2288
English AC, Richards S, Han Y, Wang M, Vee V, Qu J, Qin X, Muzny DM, Reid JG, Worley KC, Gibbs RA: Mind the gap: upgrading genomes with Pacific Biosciences RS long-read sequencing technology. PLoS One 2012, 7: e47768. 10.1371/journal.pone.0047768
Treangen TJ, Abraham A-L, Touchon M, Rocha EPC: Genesis, effects and fates of repeats in prokaryotic genomes. FEMS Microbiol Rev 2009, 33: 539-571. 10.1111/j.1574-6976.2009.00169.x
Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, Copeland A, Huddleston J, Eichler EE, Turner SW, Korlach J: Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods 2013, 10: 563-569. 10.1038/nmeth.2474
Koren S, Schatz M, Walenz B, Martin J, Howard J, Ganapathy G, Wang Z, Rasko D, McCombie W, Jarvis E, Phillippy A: Hybrid error correction and de novo assembly of single-molecule sequencing reads. Nat Biotechnol 2012, 30: 693-700. 10.1038/nbt.2280
Roberts R, Carneiro M, Schatz M: The advantages of SMRT sequencing. Genome Biol 2013, 14: 405. 10.1186/gb-2013-14-6-405
Hoefler BC, Konganti K, Straight PD: De novo assembly of the Streptomyces sp. strain Mg1 genome using PacBio single-molecule sequencing. Genome Announc 2013, 1: 1:e00535-00513.
Rasko DA, Webster DR, Sahl JW, Bashir A, Boisen N, Scheutz F, Paxinos EE, Sebra R, Chin C-S, Iliopoulos D, Klammer A, Peluso P, Lee L, Kislyuk AO, Bullard J, Kasarskis A, Wang S, Eid J, Rank D, Redman JC, Steyert SR, Frimodt-Møller J, Struve C, Petersen AM, Krogfelt KA, Nataro JP, Schadt EE, Waldor MK: Origins of the E. coli strain causing an outbreak of hemolytic–uremic syndrome in Germany. N Eng J Med 2011, 365: 709-717. 10.1056/NEJMoa1106920
Powers JG, Weigman VJ, Shu J, Pufky JM, Cox D, Hurban P: Efficient and accurate whole genome assembly and methylome profiling of E. coli . BMC Genomics 2013, 14: 675. 10.1186/1471-2164-14-675
Bruno-Barcena JM, Chinn MS, Grunden AM: Genome sequence of the autotrophic acetogen Clostridium autoethanogenum JA1-1 strain DSM 10061, a producer of ethanol from carbon monoxide. Genome Announc 2013, 1: e00628-00613.
Abrini J, Naveau H, Nyns EJ: Clostridium autoethanogenum , sp. nov., an anaerobic bacterium that produces ethanol from carbon monoxide. Arch Microbiol 1994, 161: 345-351. 10.1007/BF00303591
Köpke M, Held C, Hujer S, Liesegang H, Wiezer A, Wollherr A, Ehrenreich A, Liebl W, Gottschalk G, Dürre P: Clostridium ljungdahlii represents a microbial production platform based on syngas. Proc Natl Acad Sci USA 2010, 107: 13087-13092. 10.1073/pnas.1004716107
Köpke M, Mihalcea C, Bromley JC, Simpson SD: Fermentative production of ethanol from carbon monoxide. Curr Opin Biotechnol 2011, 22: 320-325. 10.1016/j.copbio.2011.01.005
Köpke M, Mihalcea C, Liew F-M, Tizard JH, Ali MS, Conolly JJ, Al-Sinawi B, Simpson SD: 2,3-butanediol production by acetogenic bacteria, an alternative route to chemical synthesis, using industrial waste gas. Appl Environ Microbiol 2011, 77: 5467-5475. 10.1128/AEM.00355-11
Tirado-Acevedo O, Chinn MS, Grunden AM: Production of Biofuels from Synthesis Gas Using Microbial Catalysts. Advances in Applied Microbiology. In Production of Biofuels from Synthesis Gas Using Microbial Catalysts, Volume 70 Edited by: Laskin AI, Sariaslani S, Gadd GM. 2010, 57-92.
Mohammadi M, Najafpour GD, Younesi H, Lahijani P, Uzir MH, Mohamed AR: Bioconversion of synthesis gas to second generation biofuels: A review. Renew Sustain Energ Rev 2011, 15: 4255-4273. 10.1016/j.rser.2011.07.124
Munasinghe PC, Khanal SK: Biomass-derived syngas fermentation into biofuels: Opportunities and challenges. Biores Technol 2010, 101: 5013-5022. 10.1016/j.biortech.2009.12.098
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko S, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA: SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 2012, 19: 455-477. 10.1089/cmb.2012.0021
Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 2008, 18: 821-829. 10.1101/gr.074492.107
Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol İ: ABySS: A parallel assembler for short read sequence data. Genome Res 2009, 19: 1117-1123. 10.1101/gr.089532.108
Brown SD, Klingeman DM, Lu TY, Johnson CM, Utturkar SM, Land ML, Schadt CW, Doktycz MJ, Pelletier DA: Draft genome sequence of Rhizobium sp. strain PDO1–076, a bacterium isolated from Populus deltoides . J Bacteriol 2012, 194: 2383-2384. 10.1128/JB.00198-12
Brown SD, Lamed R, Morag E, Borovok I, Shoham Y, Klingeman DM, Johnson CM, Yang Z, Land ML, Utturkar SM, Keller M, Bayer EA: Draft genome sequences for Clostridium thermocellum wild-type strain YS and derived cellulose adhesion-defective mutant strain AD2. J Bacteriol 2012, 194: 3290-3291. 10.1128/JB.00473-12
Brown SD, Podar M, Klingeman DM, Johnson CM, Yang ZK, Utturkar SM, Land ML, Mosher JJ, Hurt RA Jr, Phelps TJ, Palumbo AV, Arkin AP, Hazen TC, Elias DA: Draft genome sequences for two metal-reducing Pelosinus fermentans strains isolated from a Cr(VI)-contaminated site and for type strain R7. J Bacteriol 2012, 194: 5147-5148. 10.1128/JB.01174-12
Elkins JG, Lochner A, Hamilton-Brehm SD, Davenport KW, Podar M, Brown SD, Land ML, Hauser LJ, Klingeman DM, Raman B, Goodwin LA, Tapia R, Meincke LJ, Detter JC, Bruce DC, Han CS, Palumbo AV, Cottingham RW, Keller M, Graham DE: Complete genome sequence of the cellulolytic thermophile Caldicellulosiruptor obsidiansis OB47T. J Bacteriol 2010, 192: 6099-6100. 10.1128/JB.00950-10
Utturkar SM, Bollmann A, Brzoska RM, Klingeman DM, Epstein SE, Palumbo AV, Brown SD: Draft Genome Sequence for Caulobacter sp. Strain OR37, a Bacterium Tolerant to Heavy Metals. Genome Announc 2013, 1: e00322-00313.
Utturkar SM, Bollmann A, Brzoska RM, Klingeman DM, Epstein SE, Palumbo AV, Brown SD: Draft Genome Sequence for Ralstonia sp. Strain OR214, a Bacterium with Potential for Bioremediation. Genome Announc 2013, 1: e00321-00313.
Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ: Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 2010, 11: 119. 10.1186/1471-2105-11-119
Cotter JL, Chinn MS, Grunden AM: Influence of process parameters on growth of Clostridium ljungdahlii and Clostridium autoethanogenum on synthesis gas. Enzym Microb Tech 2009, 44: 281-288. 10.1016/j.enzmictec.2008.11.002
Cotter JL, Chinn MS, Grunden AM: Ethanol and acetate production by Clostridium ljungdahlii and Clostridium autoethanogenum using resting cells. Bioproc Biosyst Eng 2009, 32: 369-380. 10.1007/s00449-008-0256-y
Tirado-Acevedo O, Cotter JL, Chinn MS, Grunden AM: Influence of carbon source preadaptation on Clostridium ljungdahlii growth and product formation. J Bioprocess Biotechniq 2011, S2: 001.
Tanner RS, Miller LM, Yang D: Clostridium ljungdahlii sp. nov., an acetogenic species in Clostridial rRNA homology group I. Int J Syst Bacteriol 1993, 43: 232-236. 10.1099/00207713-43-2-232
Stackebrandt E, Kramer I, Swiderski J, Hippe H: Phylogenetic basis for a taxonomic dissection of the genus Clostridium . FEMS Immunol Med Microbiol 1999, 24: 253-258. 10.1111/j.1574-695X.1999.tb01291.x
Chen F, Mackey AJ, Stoeckert CJ Jr, Roos DS: OrthoMCL-DB: querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res 2006, 34: D363-D368. 10.1093/nar/gkj123
Wood HG: Life with CO or CO2 and H2 as a source of carbon and energy. FASEB J 1991, 5: 156-163.
Köpke M, Straub M, Dürre P: Clostridium difficile is an autotrophic bacterial pathogen. PLoS One 2013, 8: e62157. 10.1371/journal.pone.0062157
Wang S, Huang H, Kahnt J, Mueller AP, Köpke M, Thauer RK: NADP-specific electron-bifurcating [FeFe]-hydrogenase in a functional complex with formate dehydrogenase in Clostridium autoethanogenum grown on CO. J Bacteriol 2013, 195: 4373-4386. 10.1128/JB.00678-13
Bhaya D, Davison M, Barrangou R: CRISPR-Cas systems in Bacteria and Archaea: Versatile small RNAs for adaptive defense and regulation. Ann Rev Genet 2011, 45: 273-297. 10.1146/annurev-genet-110410-132430
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: 181-186. 10.1038/nrmicro1793
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 Comp Biol 2005, 1: 474-483.
Makarova KS, Haft DH, Barrangou R, Brouns SJJ, Charpentier E, Horvath P, Moineau S, Mojica FJM, Wolf YI, Yakunin AF, van der Oost J, Koonin EV: Evolution and classification of the CRISPR-Cas systems. Nat Rev Microbiol 2011, 9: 467-477. 10.1038/nrmicro2577
Horvath P, Coute-Monvoisin AC, Romero DA, Boyaval P, Fremaux C, Barrangou R: Comparative analysis of CRISPR loci in lactic acid bacteria genomes. Int J Food Microbiol 2009, 131: 62-70. 10.1016/j.ijfoodmicro.2008.05.030
Lillestøl RK, Redder P, Garrett RA, Brügger K: A putative viral defence mechanism in archaeal cells. Archaea 2006, 2: 59-72. 10.1155/2006/542818
Pougach K, Semenova E, Bogdanova E, Datsenko KA, Djordjevic M, Wanner BL, Severinov K: Transcription, processing and function of CRISPR cassettes in Escherichia coli . Mol Microbiol 2010, 77: 1367-1379. 10.1111/j.1365-2958.2010.07265.x
Nölling J, Breton G, Omelchenko MV, Makarova KS, Zeng Q, Gibson R, Lee HM, Dubois J, Qiu D, Hitti J, Wolf YI, Tatusov RL, Sabathe F, Doucette-Stamm L, Soucaille P, Daly MJ, Bennett GN, Koonin EV, Smith DR, GTC Sequencing Center Production Finishing Bioinformatics Teams: Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum . J Bacteriol 2001, 183: 4823-4838. 10.1128/JB.183.16.4823-4838.2001
Hu S, Zheng H, Gu Y, Zhao J, Zhang W, Yang Y, Wang S, Zhao G, Yang S, Jiang W: Comparative genomic and transcriptomic analysis revealed genetic characteristics related to solvent formation and xylose utilization in Clostridium acetobutylicum EA 2018. BMC Genomics 2011, 12: 93. 10.1186/1471-2164-12-93
Bao G, Wang R, Zhu Y, Dong H, Mao S, Zhang Y, Chen Z, Li Y, Ma Y: Complete genome sequence of Clostridium acetobutylicum DSM 1731, a solvent-producing strain with multireplicon genome architecture. J Bacteriol 2011, 193: 5007-5008. 10.1128/JB.05596-11
Poehlein A, Hartwich K, Krabben P, Ehrenreich A, Liebl W, Dürre P, Gottschalk G, Daniel R: Complete genome sequence of the solvent producer Clostridium saccharobutylicum NCP262 (DSM 13864). Genome Announc 2013, 1: e00997-13-6. 10.1128/genomeA.00997-13
del Cerro C, Felpeto-Santero C, Rojas A, Tortajada M, Ramón D, García JL: Genome sequence of the butanol hyperproducer Clostridium saccharoperbutylacetonicum N1-4. Genome Announc 2013, 1: e0007013.
Tamaru Y, Miyake H, Kuroda K, Nakanishi A, Kawade Y, Yamamoto K, Uemura M, Fujita Y, Doi RH, Ueda M: Genome sequence of the cellulosome-producing mesophilic organism Clostridium cellulovorans 743B. J Bacteriol 2010, 192: 901-902. 10.1128/JB.01450-09
Wilson CM, Rodriguez M, Johnson CM, Martin SL, Chu T-M, Wolfinger RD, Hauser LJ, Land ML, Klingeman DM, Syed M, Ragauskas AJ, Tschaplinski TJ, Mielenz JR, Brown SD: Global transcriptome analysis of Clostridium thermocellum ATCC 27405 during growth on dilute acid pretreated Populus and switchgrass. Biotechnol Biofuels 2013, 6: 179. 10.1186/1754-6834-6-179
Feinberg L, Foden J, Barrett T, Davenport KW, Bruce D, Detter C, Tapia R, Han C, Lapidus A, Lucas S, Cheng J-F, Pitluck S, Woyke T, Ivanova N, Mikhailova N, Land M, Hauser L, Argyros DA, Goodwin L, Hogsett D, Caiazza N: Complete genome sequence of the cellulolytic thermophile Clostridium thermocellum DSM1313. J Bacteriol 2011, 193: 2906-2907. 10.1128/JB.00322-11
Paul D, Austin FW, Arick T, Bridges SM, Burgess SC, Dandass YS, Lawrence ML: Genome sequence of the solvent-producing bacterium Clostridium carboxidivorans strain P7T. J Bacteriol 2010, 192: 5554-5555. 10.1128/JB.00877-10
Richter H, Zoephel J, Schermuly J, Maticzka D, Backofen R, Randau L: Characterization of CRISPR RNA processing in Clostridium thermocellum and Methanococcus maripaludis . Nucleic Acids Res 2012, 40: 9887-9896. 10.1093/nar/gks737
Brouns SJJ, Jore MM, Lundgren M, Westra ER, Slijkhuis RJH, Snijders APL, Dickman MJ, Makarova KS, Koonin EV, van der Oost J: Small CRISPR RNAs guide antiviral defense in prokaryotes. Science 2008, 321: 960-964. 10.1126/science.1159689
Carte J, Wang RY, Li H, Terns RM, Terns MP: Cas6 is an endoribonuclease that generates guide RNAs for invader defense in prokaryotes. Genes Dev 2008, 22: 3489-3496. 10.1101/gad.1742908
Lillestol 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: 259-272. 10.1111/j.1365-2958.2009.06641.x
Hale CR, Majumdar S, Elmore J, Pfister N, Compton M, Olson S, Resch AM, Glover CVC, Graveley BR, Terns RM, Terns MP: Essential features and rational design of CRISPR RNAs that function with the Cas RAMP module complex to cleave RNAs. Mol Cell 2012, 45: 292-302. 10.1016/j.molcel.2011.10.023
Jones DT, Shirley M, Wu XY, Keis S: Bacteriophage infections in the industrial acetone butanol (AB) fermentation process. J Mol Microbiol Biotechnol 2000, 2: 21-26.
Edgar R: PILER-CR: Fast and accurate identification of CRISPR repeats. BMC Bioinformatics 2007, 8: 18. 10.1186/1471-2105-8-18
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
Nozawa T, Furukawa N, Aikawa C, Watanabe T, Haobam B, Kurokawa K, Maruyama F, Nakagawa I: CRISPR Inhibition of Prophage Acquisition in Streptococcus pyogenes . PLoS One 2011, 6(5):e19543. 10.1371/journal.pone.0019543
Marraffini LA, Sontheimer EJ: CRISPR Interference Limits Horizontal Gene Transfer in Staphylococci by Targeting DNA. Science 1843–1845, 2008: 322.
Tracy BP, Jones SW, Fast AG, Indurthi DC, Papoutsakis ET: Clostridia: the importance of their exceptional substrate and metabolite diversity for biofuel and biorefinery applications. Curr Opin Biotechnol 2012, 23: 364-381. 10.1016/j.copbio.2011.10.008
Mardis ER: Next-Generation DNA sequencing methods. Annu Rev Genomics Hum Genet 2008, 9: 387-402. 10.1146/annurev.genom.9.081307.164359
Brown SD, Hurt RA, Gilmour CC, Elias DA: Draft genome sequences for three mercury-methylating, sulfate-reducing bacteria. Genome Announc 2013, 1: e00618-00613.
Brown SD, Guss A, Karpinets T, Parks J, Smolin N, Yang S, Land ML, Klingeman DM, Bhandiwad A, Rodriguez JM, Raman B, Shao X, Mielenz JR, Smith JC MK, Lynd LR: Mutant alcohol dehydrogenase leads to improved ethanol tolerance in Clostridium thermocellum . Proc Natl Acad Sci USA 2011, 108: 13752-13757. 10.1073/pnas.1102444108
Krzywinski M, Schein J, Birol İ, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: An information aesthetic for comparative genomics. Genome Res 2009, 19: 1639-1645. 10.1101/gr.092759.109
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol 2004, 5: R12. 10.1186/gb-2004-5-2-r12
Darling AE, Mau B, Perna NT: Progressive Mauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 2010, 5: e11147. 10.1371/journal.pone.0011147
Zhou Y, Liang YJ, Lynch KH, Dennis JJ, Wishart DS: PHAST: a fast phage search tool. Nucleic Acids Res 2011, 39: W347-W352. 10.1093/nar/gkr485
Fouts DE: Phage_Finder: Automated identification and classification of prophage regions in complete bacterial genome sequences. Nucleic Acids Res 2006, 34: 5839-5851. 10.1093/nar/gkl732
This study was funded by LanzaTech, NZ and leveraged resources supported by the BioEnergy Science Center (BESC). BESC is a US Department of Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research (BER) in the DOE Office of Science. The University of Maryland Institute for Genome Sciences generated the PacBio data and Dawn Klingeman (ORNL) generated Illumina data. The RNA sequencing was performed by Jon K Magnuson, David E Culley and James R Collet at Battelle Memorial Institute Pacific Northwest Division. SU is supported by the BER Genomic Science Program as part of the Plant Microbe Interfaces Scientific Focus Area. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the DOE under Contract DE-AC05-00OR22725.
LanzaTech NZ, Ltd has a commercial interest in Clostridium autoethanogenum.
SDB, SU, WM, ML, AD and MK carried out bioinformatics analysis. SDB, SN, SU, WM and MK conceived and designed the study and prepared the manuscript. SN carried out the work around the CRISPR system and RT-PCR. SU performed de novo and hybrid assemblies and evaluations, OrthoMCL analysis, and ring figure creation, and assisted with PCR. SDT prepared genomic DNA. SS carried out the RNA-Seq experiment. All authors read and approved the final manuscript.
Steven D Brown, Shilpa Nagaraju contributed equally to this work.
Electronic supplementary material
Examples of preliminary PCR and Sanger sequencing studies to close DSM 10061 genome compared to PacBio assembly.
Additional file 1: Small regions of overlap in the hybrid assembly weakly joined contigs, and were supported by PCR and Sanger data, but had insufficient support for the Newbler assembly to join contigs (A), PCR and Sanger data joined small gaps between contigs in line with predictions using C. ljungdahlii DSM 13528 as a reference (B), and in other examples much larger products were obtained compared to the predicted PCR product sizes (C). (PNG 948 KB)
Additional file 2: CGAL scores for C. autoethanogenum DSM 10061 assemblies. (DOCX 15 KB)
Additional file 3: QUAST analysis of C. autoethanogenum DSM 10061 assemblies. (DOCX 16 KB)
Additional file 4: REAPR analysis of C. autoethanogenum DSM 10061 assemblies. (DOCX 120 KB)
Additional file 5: General genome statistics for DSM 10061 PacBio assembly. (DOCX 17 KB)
Additional file 6: Number of genes associated with COG functional categories for DSM 10061 PacBio assembly. (DOCX 16 KB)
Additional file 7: C. autoethanogenum from C. ljungdahlii orthologs. The OrthoMCL algorithm  was used for ortholog analysis. The 1_taxa tab contains putative paralogs, 2_taxa file contains putative orthologs, Unique tabs contains Unique proteins using default setting, and Table 1 file is general descriptor for all proteins. The 427C. autoethanogenum from the unique tab were compared to the C. ljungdahlii proteome using BLASTP and 221 proteins were considered dissimilar or unique to DSM 10061 based on e-value scores of ≤1e-5. (XLSX 485 KB)
RNA-Seq data for Hydrogenase operon CAETHG_1575-78 and clustered regularly interspaced short paloindromic repeats-associated (CRISPR-cas) array of
Additional file 8: C. autoethanogenum . Mapped RNA-Seq reads for (FeFe) hydrogenase operon CAETHG_1575-78 and CRISPR-cas system of C. autoethanogenum. Processing of crRNA in C. autoethanogenum. (DOCX 326 KB)
Additional file 9: C. autoethanogenum clustered regularly interspaced short paloindromic repeats-associated (CRISPR- cas ) system. Organization and expression of the C. autoethanogenum CRISPR system as determined by RT-PCR and primer list for RT-PCR. (PDF 242 KB)
clustered regularly interspaced short paloindromic repeats (CRISPR) arrays and prophage regions in analyzed genomes.
Additional file 10: Overview of position and sequences of identified CRISPR arrays and prophages in industrial relevant Clostridia. (XLSX 13 MB)
Graphical representation of clustered regularly interspaced short paloindromic repeats-associated (CRISPR-
Additional file 11: cas ) loci in Clostridium species.(PDF 344 KB)
Phylogenetic classification of industrial relevant Clostridia based on 16S rRNA and clustered regularly interspaced short paloindromic repeats (CRISPR) systems.
Additional file 12: Classification of the CRISPR-cas system of C. autoethanogenum and phylogenetic classification based on 16S rRNA and cas1 genes. (PDF 511 KB)
About this article
Cite this article
Brown, S.D., Nagaraju, S., Utturkar, S. et al. Comparison of single-molecule sequencing and hybrid approaches for finishing the genome of Clostridium autoethanogenum and analysis of CRISPR systems in industrial relevant Clostridia. Biotechnol Biofuels 7, 40 (2014). https://doi.org/10.1186/1754-6834-7-40
- Clostridium Species
- Gene Operon
- Hybrid Assembly
- CRISPR Locus
- CRISPR System