Comparison of the mesophilic cellulosome‐producing Clostridium cellulovorans genome with other cellulosome‐related clostridial genomes

Clostridium cellulovorans, an anaerobic and mesophilic bacterium, degrades native substrates in soft biomass such as corn fibre and rice straw efficiently by producing an extracellular enzyme complex called the cellulosome. Recently, we have reported the whole‐genome sequence of C. cellulovorans comprising 4220 predicted genes in 5.10 Mbp [Y. Tamaru et al., (2010) J. Bacteriol., 192: 901–902]. As a result, the genome size of C. cellulovorans was about 1 Mbp larger than that of other cellulosome‐producing clostridia, mesophilic C. cellulolyticum and thermophilic C. thermocellum. A total of 57 cellulosomal genes were found in the C. cellulovorans genome, and they coded for not only carbohydrate‐degrading enzymes but also a lipase, peptidases and proteinase inhibitors. Interestingly, two novel genes encoding scaffolding proteins were found in the genome. According to KEGG metabolic pathways and their comparison with 11 Clostridial genomes, gene expansion in the C. cellulovorans genome indicated mainly non‐cellulosomal genes encoding hemicellulases and pectin‐degrading enzymes. Thus, by examining genome sequences from multiple Clostridium species, comparative genomics offers new insight into genome evolution and the way natural selection moulds functional DNA sequence evolution. Our analysis, coupled with the genome sequence data, provides a roadmap for constructing enhanced cellulosome‐producing Clostridium strains for industrial applications such as biofuel production.     

Comparative genomic analysis of hyper-ammonia producing Acetoanaerobium sticklandii DSM 519 with purinolytic Gottschalkia acidurici 9a and pathogenic Peptoclostridium difficile 630

Acetoanaerobium sticklandii DSM519 (CST) is a hype-ammonia producing non-pathogenic anaerobe that can use amino acids as important carbon and energy sources through the Stickland reactions. Biochemical aspects of this organism have been extensively studied, but systematic studies addressing its metabolic discrepancy remain scant. In this perspective, we have intensively analyzed its genomic and metabolic characteristics to comprehend the evolutionary conservation of amino acid catabolism by a comparative genomic approach. The whole-genome data indicated that CST has shown a phylogenomic similarity with hyper-ammonia producing, purinolytic, and proteolytic pathogenic Clostridia. CST has shown to common genomic context sharing across the purinolytic Gottschalkia acidurici 9a and pathogenic Peptoclostridium difficile 630. Genome syntenic analysis described that syntenic orthologs might be originated from the recent ancestor at a slow evolution rate and syntenic-out paralogs evolved from either CDF or CAC via α-event and β-event. Collinearity of either gene orders or gene families was adjusted with syntenic out-paralogs across these genomes. The genome-wide metabolic analysis predicted 11 unique putative metabolic subsystems from the CST genome for amino acid catabolism and hydrogen production. The in silico analysis of our study revealed that a characteristic system for amino acid catabolism-directed biofuel synthesis might have slowly evolved and established as a core genomic content of CST.     

Comparative studies on physiology and taxonomy of obligately purinolytic clostridia

Eleven strains of obligately purinolytic clostridia have been studied with respect to their assignment to the three type strains of Clostridium acidiurici, C. cylindrosporum, and C. purinolyticum. DNA/DNA-hybridization proved to be the method of choice for differentiation whereas phenotypic characteristics such as spore morphology, substrate spectra, nutritional requirements, product formation, and sensitivity against various antibiotics did not allow unequivocal identification. All strains depended on selenite for growth.     

Isolation and Characterization of a Native Composite Transposon, Tn14751, Carrying 17.4 Kilobases of Corynebacterium glutamicum Chromosomal DNA

A native composite transposon was isolated from Corynebacterium glutamicum ATCC 14751. This transposon comprises two functional copies of a corynebacterial IS31831-like insertion sequence organized as converging terminal inverted repeats. This novel 20.3-kb element, Tn14751, carries 17.4 kb of C. glutamicum chromosomal DNA containing various genes, including genes involved in purine biosynthesis but not genes related to bacterial warfare, such as genes encoding mediators of antibiotic resistance or extracellular toxins. A derivative of this element carrying a kanamycin resistance cassette, minicomposite Tn14751, transposed into the genome of C. glutamicum at an efficiency of 1.8 × 10² transformants per μg of DNA. Random insertion of the Tn14751 derivative carrying the kanamycin resistance cassette into the chromosome was verified by Southern hybridization. This work paves the way for realization of the concept of minimum genome factories in the search for metabolic engineering via genome-scale directed evolution through a combination of random and directed approaches.     

Complete Genome Sequences of Elephant Endotheliotropic Herpesviruses 1A and 1B Determined Directly from Fatal Cases

A highly lethal hemorrhagic disease associated with infection by elephant endotheliotropic herpesvirus (EEHV) poses a severe threat to Asian elephant husbandry. We have used high-throughput methods to sequence the genomes of the two genotypes that are involved in most fatalities, namely, EEHV1A and EEHV1B (species Elephantid herpesvirus 1, genus Proboscivirus, subfamily Betaherpesvirinae, family Herpesviridae). The sequences were determined from postmortem tissue samples, despite the data containing tiny proportions of viral reads among reads from a host for which the genome sequence was not available. The EEHV1A genome is 180,421 bp in size and consists of a unique sequence (174,601 bp) flanked by a terminal direct repeat (2,910 bp). The genome contains 116 predicted protein-coding genes, of which six are fragmented, and seven paralogous gene families are present. The EEHV1B genome is very similar to that of EEHV1A in structure, size, and gene layout. Half of the EEHV1A genes lack orthologs in other members of subfamily Betaherpesvirinae, such as human cytomegalovirus (genus Cytomegalovirus) and human herpesvirus 6A (genus Roseolovirus). Notable among these are 23 genes encoding type 3 membrane proteins containing seven transmembrane domains (the 7TM family) and seven genes encoding related type 2 membrane proteins (the EE50 family). The EE50 family appears to be under intense evolutionary selection, as it is highly diverged between the two genotypes, exhibits evidence of sequence duplications or deletions, and contains several fragmented genes. The availability of the genome sequences will facilitate future research on the epidemiology, pathogenesis, diagnosis, and treatment of EEHV-associated disease.     

Complete Genomic Structure of the Cultivated Rice Endophyte Azospirillum sp. B510

We determined the nucleotide sequence of the entire genome of a diazotrophic endophyte, Azospirillum sp. B510. Strain B510 is an endophytic bacterium isolated from stems of rice plants (Oryza sativa cv. Nipponbare). The genome of B510 consisted of a single chromosome (3 311 395 bp) and six plasmids, designated as pAB510a (1 455 109 bp), pAB510b (723 779 bp), pAB510c (681 723 bp), pAB510d (628 837 bp), pAB510e (537 299 bp), and pAB510f (261 596 bp). The chromosome bears 2893 potential protein-encoding genes, two sets of rRNA gene clusters (rrns), and 45 tRNA genes representing 37 tRNA species. The genomes of the six plasmids contained a total of 3416 protein-encoding genes, seven sets of rrns, and 34 tRNAs representing 19 tRNA species. Eight genes for plasmid-specific tRNA species are located on either pAB510a or pAB510d. Two out of eight genomic islands are inserted in the plasmids, pAB510b and pAB510e, and one of the islands is inserted into trnfM-CAU in the rrn located on pAB510e. Genes other than the nif gene cluster that are involved in N₂ fixation and are homologues of Bradyrhizobium japonicum USDA110 include fixABCX, fixNOQP, fixHIS, fixG, and fixLJK. Three putative plant hormone-related genes encoding tryptophan 2-monooxytenase (iaaM) and indole-3-acetaldehyde hydrolase (iaaH), which are involved in IAA biosynthesis, and ACC deaminase (acdS), which reduces ethylene levels, were identified. Multiple gene-clusters for tripartite ATP-independent periplasmic-transport systems and a diverse set of malic enzymes were identified, suggesting that B510 utilizes C₄-dicarboxylate during its symbiotic relationship with the host plant.     

Genome Wide Analysis of the Potato Soft Rot Pathogen Pectobacterium carotovorum Strain ICMP 5702 to Predict Novel Insights into Its Genetic Features

Pectobacterium carotovorum subsp. carotovorum (Pcc) is a gram-negative, broad host range bacterial pathogen which causes soft rot disease in potatoes as well as other vegetables worldwide. While Pectobacterium infection relies on the production of major cell wall degrading enzymes, other virulence factors and the mechanism of genetic adaptation of this pathogen is not yet clear. In the present study, we have performed an in-depth genome-wide characterization of Pcc strain ICMP5702 isolated from potato and compared it with other pathogenic bacteria from the Pectobacterium genus to identify key virulent determinants. The draft genome of Pcc ICMP5702 contains 4,774,457 bp with a G + C content of 51.90% and 4,520 open reading frames. Genome annotation revealed prominent genes encoding key virulence factors such as plant cell wall degrading enzymes, flagella-based motility, phage proteins, cell membrane structures, and secretion systems. Whereas, a majority of determinants were conserved among the Pectobacterium strains, few notable genes encoding AvrE-family type III secretion system effectors, pectate lyase and metalloprotease in addition to the CRISPR-Cas based adaptive immune system were uniquely represented. Overall, the information generated through this study will contribute to decipher the mechanism of infection and adaptive immunity in Pcc.     

Complete chloroplast genome sequences of Dipterygium glaucum and Cleome chrysantha and other Cleomaceae Species,comparative analysis and phylogenetic relationships

This current study presents, for the first time, the complete chloroplast genome of two Cleomaceae species: Dipterygium glaucum and Cleome chrysantha in order to evaluate the evolutionary relationship. The cp genome is 158,576 bp in length with 35.74% GC content in D. glaucum and 158,111 bp with 35.96% GC in C. chrysantha. Inverted repeats IR 26,209 bp, 26,251 bp each, LSC of 87,738 bp, 87,184 bp and SSC of 18,420 bp, 18,425 bp respectively. There are 136 genes in the genome, which includes 80 protein coding genes, 31 tRNA genes and four rRNA genes were observed in both chloroplast genomes. 117 genes are unique while the remaining 19 genes are duplicated in IR regions. The analysis of repeats shows that the cp genome includes all types of repeats with more frequent occurrences of palindromic; Also, this analysis indicates that the total number of simple sequence repeats (SSR) were 323 in D. glaucum, and 313 in C. chrysantha, of which the majority of the SSRs in these plastid genomes were mononucleotide repeats A/T which are located in the intergenic spacer. Moreover, the comparative analysis of the four cp sequences revealed four hotspot genes (atpF, rpoC2, rps19, and ycf1), these variable regions could be used as molecular makers for the species authentication as well as resources for inferring phylogenetic relationships of the species. All the relationships in the phylogenetic tree are with high support, this indicate that the complete chloroplast genome is a useful data for inferring phylogenetic relationship within the Cleomaceae and other families. The simple sequence repeats identified will be useful for identification, genetic diversity, and other evolutionary studies of the species. This study reported the first cp genome of the genus Dipterygium and Cleome. The finding of this study will be beneficial for biological disciplines such as evolutionary and genetic diversity studies of the species within the core Cleomaceae.     

Genetic Analysis of Two Bile Salt Hydrolase Activities in Lactobacillus acidophilus NCFM

Two genes, bshA and bshB, encoding bile salt hydrolase enzymes (EC 3.5.1.24) were identified in the genome sequence of Lactobacillus acidophilus NCFM. Targeted inactivation of these genes via chromosomal insertion of an integration vector demonstrated different substrate specificities for these two enzymes.     

Genome Sequence of the Cellulosome-Producing Mesophilic Organism Clostridium cellulovorans 743B

Clostridium cellulovorans 743B was isolated from a wood chip pile and is an anaerobic and mesophilic spore-forming bacterium. This organism degrades native substrates in soft biomass such as corn fiber and rice straw efficiently by producing an extracellular enzyme complex called the cellulosome. Here we report the genome sequence of C. cellulovorans 743B.The biotechnological potential of polysaccharolytic enzymes has resulted in the isolation and characterization of a large number of anaerobic, Gram-positive, spore-forming bacteria, the majority of which have been allocated to the genus Clostridium. Among clostridia, the cellulosomes produced by Clostridium species are particularly designed for efficient degradation of plant cell wall polysaccharides. The component parts of the multicomponent complex are integrated by virtue of a unique family of integrating modules, the cohesins and the dockerins, whose distribution and specificity dictate the overall cellulosome architecture (3). The cellulosome system in Clostridium cellulovorans 743B (ATCC 35296) has been studied extensively for the last 20 years and has resulted in providing basic information about mesophilic cellulosomes. This organism was isolated from a wood chip pile and is an anaerobic spore-forming bacterium whose optimal growth temperature is 37°C (9). It has the ability to utilize cellulose, xylan, pectin, cellobiose, glucose, fructose, galactose, and mannose as carbon sources for growth. Its fermentation products include H₂, CO₂, acetate, butyrate, formate, lactate, and ethanol. When it is grown in the presence of cellulose, electron micrographs have shown that large protuberances are present on its cell surface (4), while small or no protuberances are evident when cells are grown in the presence of glucose or cellobiose (5).We sequenced a total length of 101,749,598 bp and analyzed 381,514 reads by Genome Sequencer FLX 454./Roche sequencing (8) (GS-FLX version) to highly oversample the genome (20× coverage) and generated 123,892 paired-end sequence tags, to enable the assembly of all tags using the GS De Novo Assembler version 1.1.03.24 (Roche Diagnostics) and the Genome Analyzer II and sequencing kit 36-Cycle Run (Illumina). Finally, we assembled 30 scaffolds (sets of 601 ordered and oriented contigs; total length of 5,123,527 bp) to generate approximately 5.1 Mbp of nearly contiguous Clostridium botulinum E3 strain Alaska E43 (accession no. {"type":"entrez-nucleotide","attrs":{"text":"NC_010723","term_id":"188587536","term_text":"NC_010723"}}NC_010723) complete genome sequence. We analyzed a number of predicted genes included in the C. cellulovorans genome using CRITICA (version 1.05b) (2) and Glimmer 2 (version 2.10) (6) to find regions in proteins with known functions. We annotated and classified according to Gene Ontology (GO) (1). In silico Molecular Cloning Genomic Edition ver. 3.0.26 software (In Silico Biology, Co., Ltd., Japan) was used for individual genomic analysis.The C. cellulovorans 743B (ATCC 35296) genome consists of 5,123,527 bp. A total of 4,220 polypeptide-encoding open reading frames (ORFs) were identified using CRITICA, while 4,297 ORFs were identified using Glimmer 2. The number of ORFs identical between CRITICA and Glimmer 2 was 2,773. Sixty-three tRNAs and 33 anticodons were also identified using tRNAscan-SE (7). In comparison of the genome sizes among cellulosomal clostridia such as Clostridium cellulolyticum H10 (4.07 Mbp) (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"CP001348","term_id":"219997787","term_text":"CP001348"}}CP001348) and Clostridium thermocellum ATCC 27405 (3.84 Mbp) (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"CP000568","term_id":"125712750","term_text":"CP000568"}}CP000568), the C. cellulovorans genome was over 1 Mbp larger than the other genomes. Moreover, the number of predicted genes (4,220 by CRITICA) in the C. cellulovorans genome was the largest among them. On the other hand, the G+C content in C. cellulovorans was 31.1%, similar to that (30.9%) in Clostridium acetobutylicum ATCC 824 (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"AE001437","term_id":"25168256","term_text":"AE001437"}}AE001437), while the G+C contents in C. cellulolyticum and C. thermocellum were 37.7% and 39.0%, respectively.A protein BLAST search against the database of clusters of orthologous groups (COGs) of proteins indicated that 4,171 genes were annotated by 4,220 predicted coding sequences using CRITICA while 4,098 genes were identified by 4,297 predicted coding sequences using Glimmer 2. On the other hand, a protein BLAST search against the NCBI nr database indicated that 4,184 genes were annotated by 4,220 predicted coding sequences using CRITICA while 4,071 genes were identified by 4,297 predicted coding sequences using Glimmer 2. Interestingly, 57 cellulosomal genes were found in the C. cellulovorans genome and coded for not only carbohydrate-active enzymes but also lipases, peptidases, and proteinase inhibitors. Moreover, two novel genes encoding a scaffolding protein were found in the genome. Thus, by examining genome sequences from multiple Clostridium species, comparative genomics offers new insight into genome evolution and the way in which natural selection molds functional DNA sequence evolution. Our analysis, coupled with the genome sequence data, will provide a road map for constructing enhanced cellulosome-producing Clostridium strains for industrial applications such as biofuel production.     

Genome Sequence of the Chemolithoautotrophic Nitrite-Oxidizing Bacterium Nitrobacter winogradskyi Nb-255

The alphaproteobacterium Nitrobacter winogradskyi (ATCC 25391) is a gram-negative facultative chemolithoautotroph capable of extracting energy from the oxidation of nitrite to nitrate. Sequencing and analysis of its genome revealed a single circular chromosome of 3,402,093 bp encoding 3,143 predicted proteins. There were extensive similarities to genes in two alphaproteobacteria, Bradyrhizobium japonicum USDA110 (1,300 genes) and Rhodopseudomonas palustris CGA009 CG (815 genes). Genes encoding pathways for known modes of chemolithotrophic and chemoorganotrophic growth were identified. Genes encoding multiple enzymes involved in anapleurotic reactions centered on C₂ to C₄ metabolism, including a glyoxylate bypass, were annotated. The inability of N. winogradskyi to grow on C₆ molecules is consistent with the genome sequence, which lacks genes for complete Embden-Meyerhof and Entner-Doudoroff pathways, and active uptake of sugars. Two gene copies of the nitrite oxidoreductase, type I ribulose-1,5-bisphosphate carboxylase/oxygenase, cytochrome c oxidase, and gene homologs encoding an aerobic-type carbon monoxide dehydrogenase were present. Similarity of nitrite oxidoreductases to respiratory nitrate reductases was confirmed. Approximately 10% of the N. winogradskyi genome codes for genes involved in transport and secretion, including the presence of transporters for various organic-nitrogen molecules. The N. winogradskyi genome provides new insight into the phylogenetic identity and physiological capabilities of nitrite-oxidizing bacteria. The genome will serve as a model to study the cellular and molecular processes that control nitrite oxidation and its interaction with other nitrogen-cycling processes.     

Isolation and characterization of a bile acid inducible 7α-dehydroxylating operon in Clostridium hylemonae TN271

Primary bile acids are synthesized from cholesterol in the liver, conjugated to either glycine or taurine and secreted into bile. Bile salts undergo enterohepatic circulation several times each day. During this process, they are biotransformed into a variety of metabolites by gut bacteria. The major biotransformation is the 7α-dehydroxylation of cholic acid and chenodeoxycholic acid yielding deoxycholic acid and lithocholic acid, respectively. 7α-Dehydroxylation is a multi-step pathway. The genes encoding enzymes in this pathway have been identified in two species of “high” activity strains of clostridia. Here, we report the isolation and characterization of a bile acid inducible (bai) operon in Clostridium hylemonae, a “low” activity 7α-dehydroxylating strain. The gene organization and sequence of the baiBCDEFGHI operon was highly conserved between C. hylemonae and “high” activity strains. Surprisingly, the baiA gene was missing from the bai operon of C. hylemonae. The baiA gene was isolated using PCR and degenerate oligonucleotide primers. The mRNA start site for the large bai operon was determined and shown to be only 11 bp from the initiation codon of the first gene. It was also discovered that allocholic acid (5α) induced the bai operon and stimulated the conversion of [24-¹⁴C] cholic acid to [24-¹⁴C] allodeoxycholic acid in cultures of C. scindens and C. hylemonae allodeoxycholic acid. Finally, it was discovered that the addition of testosterone to the growth medium markedly increased 7α-dehydroxylation of cholic acid in Clostridium scindens and C. hylemonae. We hypothesize that testosterone may be a gratuitous inducer of genes involved in the reductive arm of the bile acid 7α-dehydroxylation pathway.     

Development of a genome editing technique using the CRISPR/Cas9 system in the industrial filamentous fungus <Emphasis Type="Italic">Aspergillus oryzae</Emphasis>

Objectives

To develop a genome editing method using the CRISPR/Cas9 system in Aspergillus oryzae, the industrial filamentous fungus used in Japanese traditional fermentation and for the production of enzymes and heterologous proteins.

Results

To develop the CRISPR/Cas9 system as a genome editing technique for A. oryzae, we constructed plasmids expressing the gene encoding Cas9 nuclease and single guide RNAs for the mutagenesis of target genes. We introduced these into an A. oryzae strain and obtained transformants containing mutations within each target gene that exhibited expected phenotypes. The mutational rates ranged from 10 to 20 %, and 1 bp deletions or insertions were the most commonly induced mutations.

Conclusions

We developed a functional and versatile genome editing method using the CRISPR/Cas9 system in A. oryzae. This technique will contribute to the use of efficient targeted mutagenesis in many A. oryzae industrial strains.

     

Complete genome sequence of the industrial bacterium Bacillus licheniformis and comparisons with closely related Bacillus species

Background

Bacillus licheniformis is a Gram-positive, spore-forming soil bacterium that is used in the biotechnology industry to manufacture enzymes, antibiotics, biochemicals and consumer products. This species is closely related to the well studied model organism Bacillus subtilis, and produces an assortment of extracellular enzymes that may contribute to nutrient cycling in nature.

Results

We determined the complete nucleotide sequence of the B. licheniformis ATCC 14580 genome which comprises a circular chromosome of 4,222,336 base-pairs (bp) containing 4,208 predicted protein-coding genes with an average size of 873 bp, seven rRNA operons, and 72 tRNA genes. The B. licheniformis chromosome contains large regions that are colinear with the genomes of B. subtilis and Bacillus halodurans, and approximately 80% of the predicted B. licheniformis coding sequences have B. subtilis orthologs.

Conclusions

Despite the unmistakable organizational similarities between the B. licheniformis and B. subtilis genomes, there are notable differences in the numbers and locations of prophages, transposable elements and a number of extracellular enzymes and secondary metabolic pathway operons that distinguish these species. Differences include a region of more than 80 kilobases (kb) that comprises a cluster of polyketide synthase genes and a second operon of 38 kb encoding plipastatin synthase enzymes that are absent in the B. licheniformis genome. The availability of a completed genome sequence for B. licheniformis should facilitate the design and construction of improved industrial strains and allow for comparative genomics and evolutionary studies within this group of Bacillaceae.