Oleg N. Reva

Comparative analysis of the complete genome sequence of the plant growth-promoting bacterium Bacillus amyloliquefaciens FZB42.

Bacillus amyloliquefaciens FZB42 is a Gram-positive, plant-associated bacterium, which stimulates plant growth and produces secondary metabolites that suppress soil-borne plant pathogens. Its 3,918-kb genome, containing an estimated 3,693 protein-coding sequences, lacks extended phage insertions, which occur ubiquitously in the closely related Bacillus subtilis 168 genome. The B. amyloliquefaciens FZB42 genome reveals an unexpected potential to produce secondary metabolites, including the polyketides bacillaene and difficidin. More than 8.5% of the genome is devoted to synthesizing antibiotics and siderophores by pathways not involving ribosomes. Besides five gene clusters, known from B. subtilis to mediate nonribosomal synthesis of secondary metabolites, we identified four giant gene clusters absent in B. subtilis 168. The pks2 gene cluster encodes the components to synthesize the macrolactin core skeleton.

Sequence Analysis of the Mobile Genome Island pKLC102 of Pseudomonas aeruginosa C

The Pseudomonas aeruginosa plasmid pKLC102 coexists as a plasmid and a genome island in clone C strains. Whereas the related plasmid pKLK106 reversibly recombines with P. aeruginosa clone K chromosomes at one of the two tRNA(Lys) genes, pKLC102 is incorporated into the tRNA(Lys) gene only close to the pilA locus. Targeting of the other tRNA(Lys) copy in the chromosome is blocked by a 23,395-bp mosaic of truncated PAO open reading frames, transposons, and pKLC102 homologs. Annotation and phylogenetic analysis of the large 103,532-bp pKLC102 sequence revealed that pKLC102 is a hybrid of plasmid and phage origin. The plasmid lineage conferred oriV and genes for replication, partitioning, and conjugation, including a pil cluster encoding type IV thin sex pili and an 8,524-bp chvB glucan synthetase gene that is known to be a major determinant for host tropism and virulence. The phage lineage conferred integrase, att, and a syntenic set of conserved hypothetical genes also observed in the tRNA(Gly)-associated genome islands of P. aeruginosa clone C chromosomes. In subgroup C isolates from patients with cystic fibrosis, pKLC102 was irreversibly fixed into the chromosome by the insertion of the large 23,061-bp class I transposon TNCP23, which is a composite of plasmid, integron, and IS6100 elements. Intramolecular transposition of a copy of IS6100 led to chromosomal inversions and disruption of plasmid synteny. The case of pKLC102 in P. aeruginosa clone C documents the intraclonal evolution of a genome island from a mobile ancestor via a reversibly integrated state to irreversible incorporation and dissipation in the chromosome.

Diversity of the Abundant pKLC102/PAGI-2 Family of Genomic Islands in Pseudomonas aeruginosa

The known genomic islands of Pseudomonas aeruginosa clone C strains are integrated into tRNA(Lys) (pKLC102) or tRNA(Gly) (PAGI-2 and PAGI-3) genes and differ from their core genomes by distinctive tetranucleotide usage patterns. pKLC102 and the related island PAPI-1 from P. aeruginosa PA14 were spontaneously mobilized from their host chromosomes at frequencies of 10% and 0.3%, making pKLC102 the most mobile genomic island known with a copy number of 30 episomal circular pKLC102 molecules per cell. The incidence of islands of the pKLC102/PAGI-2 type was investigated in 71 unrelated P. aeruginosa strains from diverse habitats and geographic origins. pKLC102- and PAGI-2-like islands were identified in 50 and 31 strains, respectively, and 15 and 10 subtypes were differentiated by hybridization on pKLC102 and PAGI-2 macroarrays. The diversity of PAGI-2-type islands was mainly caused by one large block of strain-specific genes, whereas the diversity of pKLC102-type islands was primarily generated by subtype-specific combination of gene cassettes. Chromosomal loss of PAGI-2 could be documented in sequential P. aeruginosa isolates from individuals with cystic fibrosis. PAGI-2 was present in most tested Cupriavidus metallidurans and Cupriavidus campinensis isolates from polluted environments, demonstrating the spread of PAGI-2 across habitats and species barriers. The pKLC102/PAGI-2 family is prevalent in numerous beta- and gammaproteobacteria and is characterized by high asymmetry of the cDNA strands. This evolutionarily ancient family of genomic islands retained its oligonucleotide signature during horizontal spread within and among taxa.

Genome sequence of B. amyloliquefaciens type strain DSM7(T) reveals differences to plant-associated B. amyloliquefaciens FZB42.

The complete genome sequence of Bacillus amyloliquefaciens type strain DSM7(T) is presented. A comparative analysis between the genome sequences of the plant associated strain FZB42 (Chen et al., 2007) with the genome of B. amyloliquefaciens DSM7(T) revealed obvious differences in the variable part of the genomes, whilst the core genomes were found to be very similar. The strains FZB42 and DSM7(T) have in common 3345 genes (CDS) in their core genomes; whilst 547 and 344 CDS were found to be unique in DSM7(T) and FZB42, respectively. The core genome shared by both strains exhibited 97.89% identity on amino acid level. The number of genes representing the core genome of the strains FZB42, DSM7(T), and Bacillus subtilis DSM10(T) was calculated as being 3098 and their identity was 92.25%. The 3,980,199 bp genome of DSM7(T) contains numerous genomic islands (GI) detected by different methods. Many of them were located in vicinity of tRNA, glnA, and glmS gene copies. In contrast to FZB42, but similar to B. subtilis DSM10(T), the GI were enriched in prophage sequences and often harbored transposases, integrases and recombinases. Compared to FZB42, B. amyloliquefaciens DSM7(T) possessed a reduced potential to non-ribosomally synthesize secondary metabolites with antibacterial and/or antifungal action. B. amyloliquefaciens DSM7(T) did not produce the polyketides difficidin and macrolactin and was impaired in its ability to produce lipopeptides other than surfactin. Differences established within the variable part of the genomes, justify our proposal to discriminate the plant-associated ecotype represented by FZB42 from the group of type strain related B. amyloliquefaciens soil bacteria.

Genome sequence and functional genomic analysis of the oil-degrading bacterium Oleispira antarctica.

Ubiquitous bacteria from the genus Oleispira drive oil degradation in the largest environment on Earth, the cold and deep sea. Here we report the genome sequence of Oleispira antarctica and show that compared with Alcanivorax borkumensis—the paradigm of mesophilic hydrocarbonoclastic bacteria—O. antarctica has a larger genome that has witnessed massive gene-transfer events. We identify an array of alkane monooxygenases, osmoprotectants, siderophores and micronutrient-scavenging pathways. We also show that at low temperatures, the main protein-folding machine Cpn60 functions as a single heptameric barrel that uses larger proteins as substrates compared with the classical double-barrel structure observed at higher temperatures. With 11 protein crystal structures, we further report the largest set of structures from one psychrotolerant organism. The most common structural feature is an increased content of surface-exposed negatively charged residues compared to their mesophilic counterparts. Our findings are relevant in the context of microbial cold-adaptation mechanisms and the development of strategies for oil-spill mitigation in cold environments.

Diversity of Glycosyl Hydrolases from Cellulose-Depleting Communities Enriched from Casts of Two Earthworm Species†

ABSTRACT The guts and casts of earthworms contain microbial assemblages that process large amounts of organic polymeric substrates from plant litter and soil; however, the enzymatic potential of these microbial communities remains largely unexplored. In the present work, we retrieved carbohydrate-modifying enzymes through the activity screening of metagenomic fosmid libraries from cellulose-depleting microbial communities established with the fresh casts of two earthworm species, Aporrectodea caliginosa and Lumbricus terrestris, as inocula. Eight glycosyl hydrolases (GHs) from the A. caliginosa-derived community were multidomain endo-β-glucanases, β-glucosidases, β-cellobiohydrolases, β-galactosidase, and β-xylosidases of known GH families. In contrast, two GHs derived from the L. terrestris microbiome had no similarity to any known GHs and represented two novel families of β-galactosidases/α-arabinopyranosidases. Members of these families were annotated in public databases as conserved hypothetical proteins, with one being structurally related to isomerases/dehydratases. This study provides insight into their biochemistry, domain structures, and active-site architecture. The two communities were similar in bacterial composition but significantly different with regard to their eukaryotic inhabitants. Further sequence analysis of fosmids and plasmids bearing the GH-encoding genes, along with oligonucleotide usage pattern analysis, suggested that those apparently originated from Gammaproteobacteria (pseudomonads and Cellvibrio-like organisms), Betaproteobacteria (Comamonadaceae), and Alphaproteobacteria (Rhizobiales).

Functional Metagenomics Unveils a Multifunctional Glycosyl Hydrolase from the Family 43 Catalysing the Breakdown of Plant Polymers in the Calf Rumen

Microbial communities from cow rumen are known for their ability to degrade diverse plant polymers at high rates. In this work, we identified 15 hydrolases through an activity-centred metagenome analysis of a fibre-adherent microbial community from dairy cow rumen. Among them, 7 glycosyl hydrolases (GHs) and 1 feruloyl esterase were successfully cloned, expressed, purified and characterised. The most striking result was a protein of GH family 43 (GHF43), hereinafter designated as R_09-02, which had characteristics very distinct from the other proteins in this family with mono-functional β-xylosidase, α-xylanase, α-L-arabinase and α-L-arabinofuranosidase activities. R_09-02 is the first multifunctional enzyme to exhibit β-1,4 xylosidase, α-1,5 arabinofur(pyr)anosidase, β-1,4 lactase, α-1,6 raffinase, α-1,6 stachyase, β-galactosidase and α-1,4 glucosidase activities. The R_09-02 protein appears to originate from the chromosome of a member of Clostridia, a class of phylum Firmicutes, members of which are highly abundant in ruminal environment. The evolution of R_09-02 is suggested to be driven from the xylose- and arabinose-specific activities, typical for GHF43 members, toward a broader specificity to the glucose- and galactose-containing components of lignocellulose. The apparent capability of enzymes from the GHF43 family to utilise xylose-, arabinose-, glucose- and galactose-containing oligosaccharides has thus far been neglected by, or could not be predicted from, genome and metagenome sequencing data analyses. Taking into account the abundance of GHF43-encoding gene sequences in the rumen (up to 7% of all GH-genes) and the multifunctional phenotype herein described, our findings suggest that the ecological role of this GH family in the digestion of ligno-cellulosic matter should be significantly reconsidered.

Think big--giant genes in bacteria.

Long genes should be rare in archaea and eubacteria because of the demanding costs of time and resources for protein production. The search in 580 sequenced prokaryotic genomes, however, revealed 0.2% of all genes to be longer than 5 kb (absolute number: 3732 genes). Eighty giant bacterial genes of more than 20 kb in length were identified in 47 taxa that belong to the phyla Thermotogae (1), Chlorobi (3), Planctomycetes (1), Cyanobacteria (2), Firmicutes (7), Actinobacteria (9), Proteobacteria (23) or Euryarchaeota (1) (number of taxa in brackets). Giant genes are strain-specific, differ in their tetranucleotide usage from the bulk genome and occur preferentially in non-pathogenic environmental bacteria. The two longest bacterial genes known to date were detected in the green sulfur bacterium Chlorobium chlorochromatii CaD3 encoding proteins of 36 806 and 20 647 amino acids, being surpassed in length only by the human titin coding sequence. More than 90% of bacterial giant genes either encode a surface protein or a polyketide/non-ribosomal peptide synthetase. Most surface proteins are acidic, threonine-rich, lack cystein and harbour multiple amino acid repeats. Giant proteins increase bacterial fitness by the production of either weapons towards or shields against animate competitors or hostile environments.

Global features of sequences of bacterial chromosomes, plasmids and phages revealed by analysis of oligonucleotide usage patterns.

BackgroundOligonucleotide frequencies were shown to be conserved signatures for bacterial genomes, however, the underlying constraints have yet not been resolved in detail. In this paper we analyzed oligonucleotide usage (OU) biases in a comprehensive collection of 155 completely sequenced bacterial chromosomes, 316 plasmids and 104 phages.ResultsTwo global features were analyzed: pattern skew (PS) and variance of OU deviations normalized by mononucleotide content of the sequence (OUV). OUV reflects the strength of OU biases and taxonomic signals. PS denotes asymmetry of OU in direct and reverse DNA strands. A trend towards minimal PS was observed for almost all complete sequences of bacterial chromosomes and plasmids, however, PS was substantially higher in separate genomic loci and several types of plasmids and phages characterized by long stretches of non-coding DNA and/or asymmetric gene distribution on the two DNA strands. Five of the 155 bacterial chromosomes have anomalously high PS, of which the chromosomes of Xylella fastidiosa 9a5c and Prochlorococcus marinus MIT9313 exhibit extreme PS values suggesting an intermediate unstable state of these two genomes.ConclusionsStrand symmetry as indicated by minimal PS is a universally conserved feature of complete bacterial genomes that results from the matching mutual compensation of local OU biases on both replichors while OUV is more a taxon specific feature. Local events such as inversions or the incorporation of genome islands are balanced by global changes in genome organization to minimize PS that may represent one of the leading evolutionary forces driving bacterial genome diversification.

Differentiation of regions with atypical oligonucleotide composition in bacterial genomes

BackgroundComplete sequencing of bacterial genomes has become a common technique of present day microbiology. Thereafter, data mining in the complete sequence is an essential step. New in silico methods are needed that rapidly identify the major features of genome organization and facilitate the prediction of the functional class of ORFs. We tested the usefulness of local oligonucleotide usage (OU) patterns to recognize and differentiate types of atypical oligonucleotide composition in DNA sequences of bacterial genomes.ResultsA total of 163 bacterial genomes of eubacteria and archaea published in the NCBI database were analyzed. Local OU patterns exhibit substantial intrachromosomal variation in bacteria. Loci with alternative OU patterns were parts of horizontally acquired gene islands or ancient regions such as genes for ribosomal proteins and RNAs. OU statistical parameters, such as local pattern deviation (D), pattern skew (PS) and OU variance (OUV) enabled the detection and visualization of gene islands of different functional classes.ConclusionA set of approaches has been designed for the statistical analysis of nucleotide sequences of bacterial genomes. These methods are useful for the visualization and differentiation of regions with atypical oligonucleotide composition prior to or accompanying gene annotation.