Differentiation of regions with atypical oligonucleotide composition in bacterial genomes

Life Within a Contaminated Niche: Comparative Genomic Analyses of an Integrative Conjugative Element ICEnahCSV86 and Two Genomic Islands From Pseudomonas bharatica CSV86T Suggest Probable Role in Colonization and Adaptation

Comparative genomic and functional analyses revealed the presence of three genomic islands (GIs, >50 Kb size): ICEnahCSV86, Pseudomonas bharatica genomic island-1 (PBGI-1), and PBGI-2 in the preferentially aromatic-degrading soil bacterium, Pseudomonas bharatica CSV86^T. Site-specific genomic integration at or near specific transfer RNAs (tRNAs), near-syntenic structural modules, and phylogenetic relatedness indicated their evolutionary lineage to the type-4 secretion system (T4SS) ICEclc family, thus predicting these elements to be integrative conjugative elements (ICEs). These GIs were found to be present as a single copy in the genome and the encoded phenotypic traits were found to be stable, even in the absence of selection pressure. ICEnahCSV86 harbors naphthalene catabolic (nah-sal) cluster, while PBGI-1 harbors Co-Zn-Cd (czc) efflux genes as cargo modules, whereas PBGI-2 was attributed to as a mixed-function element. The ICEnahCSV86 has been reported to be conjugatively transferred (frequency of 7 × 10^–8/donor cell) to Stenotrophomonas maltophilia CSV89. Genome-wide comparative analyses of aromatic-degrading bacteria revealed nah-sal clusters from several Pseudomonas spp. as part of probable ICEs, syntenic to conjugatively transferable ICEnahCSV86 of strain CSV86^T, suggesting it to be a prototypical element for naphthalene degradation. It was observed that the plasmids harboring nah-sal clusters were phylogenetically incongruent with predicted ICEs, suggesting genetic divergence of naphthalene metabolic clusters in the Pseudomonas population. Gene synteny, divergence estimates, and codon-based Z-test indicated that ICEnahCSV86 is probably derived from PBGI-2, while multiple recombination events masked the ancestral lineage of PBGI-1. Diversifying selection pressure (dN-dS = 2.27–4.31) imposed by aromatics and heavy metals implied the modular exchange-fusion of various cargo clusters through events like recombination, rearrangement, domain reshuffling, and active site optimization, thus allowing the strain to evolve, adapt, and maximize the metabolic efficiency in a contaminated niche. The promoters (Pnah and Psal) of naphthalene cargo modules (nah, sal) on ICEnahCSV86 were proved to be efficient for heterologous protein expression in Escherichia coli. GI-based genomic plasticity expands the metabolic spectrum and versatility of CSV86^T, rendering efficient adaptation to the contaminated niche. Such isolate(s) are of utmost importance for their application in bioremediation and are the probable ideal host(s) for metabolic engineering.

EuGI: a novel resource for studying genomic islands to facilitate horizontal gene transfer detection in eukaryotes

Background

Genomic islands (GIs) are inserts of foreign DNA that have potentially arisen through horizontal gene transfer (HGT). There are evidences that GIs can contribute significantly to the evolution of prokaryotes. The acquisition of GIs through HGT in eukaryotes has, however, been largely unexplored. In this study, the previously developed GI prediction tool, SeqWord Gene Island Sniffer (SWGIS), is modified to predict GIs in eukaryotic chromosomes. Artificial simulations are used to estimate ratios of predicting false positive and false negative GIs by inserting GIs into different test chromosomes and performing the SWGIS v2.0 algorithm. Using SWGIS v2.0, GIs are then identified in 36 fungal, 22 protozoan and 8 invertebrate genomes.

Results

SWGIS v2.0 predicts GIs in large eukaryotic chromosomes based on the atypical nucleotide composition of these regions. Averages for predicting false negative and false positive GIs were 20.1% and 11.01% respectively. A total of 10,550 GIs were identified in 66 eukaryotic species with 5299 of these GIs coding for at least one functional protein. The EuGI web-resource, freely accessible at http://eugi.bi.up.ac.za, was developed that allows browsing the database created from identified GIs and genes within GIs through an interactive and visual interface.

Conclusions

SWGIS v2.0 along with the EuGI database, which houses GIs identified in 66 different eukaryotic species, and the EuGI web-resource, provide the first comprehensive resource for studying HGT in eukaryotes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4724-8) contains supplementary material, which is available to authorized users.

SWPhylo - A Novel Tool for Phylogenomic Inferences by Comparison of Oligonucleotide Patterns and Integration of Genome-Based and Gene-Based Phylogenetic Trees

Modern phylogenetic studies may benefit from the analysis of complete genome sequences of various microorganisms. Evolutionary inferences based on genome-scale analysis are believed to be more accurate than the gene-based alternative. However, the computational complexity of current phylogenomic procedures, inappropriateness of standard phylogenetic tools to process genome-wide data, and lack of reliable substitution models which correlates with alignment-free phylogenomic approaches deter microbiologists from using these opportunities. For example, the super-matrix and super-tree approaches of phylogenomics use multiple integrated genomic loci or individual gene-based trees to infer an overall consensus tree. However, these approaches potentially multiply errors of gene annotation and sequence alignment not mentioning the computational complexity and laboriousness of the methods. In this article, we demonstrate that the annotation- and alignment-free comparison of genome-wide tetranucleotide frequencies, termed oligonucleotide usage patterns (OUPs), allowed a fast and reliable inference of phylogenetic trees. These were congruent to the corresponding whole genome super-matrix trees in terms of tree topology when compared with other known approaches including 16S ribosomal RNA and GyrA protein sequence comparison, complete genome-based MAUVE, and CVTree methods. A Web-based program to perform the alignment-free OUP-based phylogenomic inferences was implemented at http://swphylo.bi.up.ac.za/. Applicability of the tool was tested on different taxa from subspecies to intergeneric levels. Distinguishing between closely related taxonomic units may be enforced by providing the program with alignments of marker protein sequences, eg, GyrA.

A Practical Guide for Comparative Genomics of Mobile Genetic Elements in Prokaryotic Genomes

Mobile genetic elements (MGEs) are an important feature of prokaryote genomes but are seldom well annotated and, consequently, are often underestimated. MGEs include transposons (Tn), insertion sequences (ISs), prophages, genomic islands (GEIs), integrons, and integrative and conjugative elements (ICEs). They are intimately involved in genome evolution and promote phenomena such as genomic expansion and rearrangement, emergence of virulence and pathogenicity, and symbiosis. In spite of the annotation bottleneck, there are so far at least 75 different programs and databases dedicated to prokaryotic MGE analysis and annotation, and this number is rapidly growing. Here, we present a practical guide to explore, compare, and visualize prokaryote MGEs using a combination of available software and databases tailored to small scale genome analyses. This protocol can be coupled with expert MGE annotation and exploited for evolutionary and comparative genomic analyses.

Diversity of the cell-wall associated genomic island of the archaeon Haloquadratum walsbyi

BACKGROUND

Haloquadratum walsbyi represents up to 80% of cells in NaCl-saturated brines worldwide, but is notoriously difficult to maintain under laboratory conditions. In order to establish the extent of genetic diversity in a natural population of this microbe, we screened a H. walsbyi enriched metagenomic fosmid library and recovered seven novel version of its cell-wall associated genomic island. The fosmid inserts were sequenced and analysed.

RESULTS

The novel cell-wall associated islands delineated two major clades within H. walsbyi. The islands predominantly contained genes putatively involved in biosynthesis of surface layer, genes encoding cell surface glycoproteins and genes involved in envelope formation. We further found that these genes are maintained in the population and that the diversity of this region arises through homologous recombination but also through the action of mobile genetic elements, including viruses.

CONCLUSIONS

The population of H. walsbyi in the studied saltern brine is composed of numerous clonal lineages that differ in surface structures including the cell wall. This type of variation probably reflects a number of mechanisms that minimize the infection rate of predating viruses.

Pre_GI: a global map of ontological links between horizontally transferred genomic islands in bacterial and archaeal genomes

The Predicted Genomic Islands database (Pre_GI) is a comprehensive repository of prokaryotic genomic islands (islands, GIs) freely accessible at http://pregi.bi.up.ac.za/index.php. Pre_GI, Version 2015, catalogues 26 744 islands identified in 2407 bacterial/archaeal chromosomes and plasmids. It provides an easy-to-use interface which allows users the ability to query against the database with a variety of fields, parameters and associations. Pre_GI is constructed to be a web-resource for the analysis of ontological roads between islands and cartographic analysis of the global fluxes of mobile genetic elements through bacterial and archaeal taxonomic borders. Comparison of newly identified islands against Pre_GI presents an alternative avenue to identify their ontology, origin and relative time of acquisition. Pre_GI aims to aid research on horizontal transfer events and materials through providing data and tools for holistic investigation of migration of genes through ecological niches and taxonomic boundaries.

The genome of Pseudomonas fluorescens strain R124 demonstrates phenotypic adaptation to the mineral environment

Microbial adaptation to environmental conditions is a complex process, including acquisition of positive traits through horizontal gene transfer or the modification of existing genes through duplication and/or mutation. In this study, we examined the adaptation of a Pseudomonas fluorescens isolate (R124) from the nutrient-limited mineral environment of a silica cave in comparison with P. fluorescens isolates from surface soil and the rhizosphere. Examination of metal homeostasis gene pathways demonstrated a high degree of conservation, suggesting that such systems remain functionally similar across chemical environments. The examination of genomic islands unique to our strain revealed the presence of genes involved in carbohydrate metabolism, aromatic carbon metabolism, and carbon turnover, confirmed through phenotypic assays, suggesting the acquisition of potentially novel mechanisms for energy metabolism in this strain. We also identified a twitching motility phenotype active at low-nutrient concentrations that may allow alternative exploratory mechanisms for this organism in a geochemical environment. Two sets of candidate twitching motility genes are present within the genome, one on the chromosome and one on a plasmid; however, a plasmid knockout identified the functional gene as being present on the chromosome. This work highlights the plasticity of the Pseudomonas genome, allowing the acquisition of novel nutrient-scavenging pathways across diverse geochemical environments while maintaining a core of functional stress response genes.

Genome sequence of the plant growth promoting strain Bacillus amyloliquefaciens subsp. plantarum B9601-Y2 and expression of mersacidin and other secondary metabolites

The plant-associated Bacillus amyloliquefaciens subsp. plantarum strain B9601-Y2, isolated from wheat rhizosphere, is a powerful plant growth-promoting rhizobacterium. Its relative large genome size of 4.24Mbp, exceeding that of other representatives of the B. amyloliquefaciens subsp. plantarum taxon, is mainly due to the presence of 18 DNA-islands containing remnants of phages, a unique restriction modification system, a gene cluster for mersacidin synthesis, and an orphan gene cluster devoted to non-ribosomal synthesis of an unidentified peptide. Like other members of the taxon, the Y2 genome contains giant gene clusters for non-ribosomal synthesis of the polyketides macrolactin, difficidin, and bacillaene, the antifungal lipopeptides bacillomycin D, and fengycin, the siderophore bacillibactin, and the dipeptide bacilysin. A gene cluster encoding enzymes for a degradative pathway with 2-keto-3-deoxygluconate and 2-keto-3-deoxy-phosphogluconate as intermediates was explored by genome mining and found as being a unique feature for representatives of the plantarum subspecies. A survey of the Y2 genome against other B. amyloliquefaciens genomes revealed 130 genes only occurring in subsp. plantarum but not in subsp. amyloliquefaciens. Notably, the surfactin gene cluster is not functional due to a large deletion removing parts of the Srf synthetases B and C. Expression of polyketides, lipopeptides, mersacidin, and of the growth hormone indole-3-acetic acid in Y2 was demonstrated by matrix-assisted laser desorption ionization-time of flight mass spectroscopy and high-performance liquid chromatography, respectively.

Mainstreams of horizontal gene exchange in enterobacteria: consideration of the outbreak of enterohemorrhagic E. coli O104:H4 in Germany in 2011

BACKGROUND

Escherichia coli O104:H4 caused a severe outbreak in Europe in 2011. The strain TY-2482 sequenced from this outbreak allowed the discovery of its closest relatives but failed to resolve ways in which it originated and evolved. On account of the previous statement, may we expect similar upcoming outbreaks to occur recurrently or spontaneously in the future? The inability to answer these questions shows limitations of the current comparative and evolutionary genomics methods.

PRINCIPAL FINDINGS

The study revealed oscillations of gene exchange in enterobacteria, which originated from marine γ-Proteobacteria. These mobile genetic elements have become recombination hotspots and effective 'vehicles' ensuring a wide distribution of successful combinations of fitness and virulence genes among enterobacteria. Two remarkable peculiarities of the strain TY-2482 and its relatives were observed: i) retaining the genetic primitiveness by these strains as they somehow avoided the main fluxes of horizontal gene transfer which effectively penetrated other enetrobacteria; ii) acquisition of antibiotic resistance genes in a plasmid genomic island of β-Proteobacteria origin which ontologically is unrelated to the predominant genomic islands of enterobacteria.

CONCLUSIONS

Oscillations of horizontal gene exchange activity were reported which result from a counterbalance between the acquired resistance of bacteria towards existing mobile vectors and the generation of new vectors in the environmental microflora. We hypothesized that TY-2482 may originate from a genetically primitive lineage of E. coli that has evolved in confined geographical areas and brought by human migration or cattle trade onto an intersection of several independent streams of horizontal gene exchange. Development of a system for monitoring the new and most active gene exchange events was proposed.

Pseudomonas aeruginosa Genomic Structure and Diversity

The Pseudomonas aeruginosa genome (G + C content 65–67%, size 5.5–7 Mbp) is made up of a single circular chromosome and a variable number of plasmids. Sequencing of complete genomes or blocks of the accessory genome has revealed that the genome encodes a large repertoire of transporters, transcriptional regulators, and two-component regulatory systems which reflects its metabolic diversity to utilize a broad range of nutrients. The conserved core component of the genome is largely collinear among P. aeruginosa strains and exhibits an interclonal sequence diversity of 0.5–0.7%. Only a few loci of the core genome are subject to diversifying selection. Genome diversity is mainly caused by accessory DNA elements located in 79 regions of genome plasticity that are scattered around the genome and show an anomalous usage of mono- to tetradecanucleotides. Genomic islands of the pKLC102/PAGI-2 family that integrate into tRNA^Lys or tRNA^Gly genes represent hotspots of inter- and intraclonal genomic diversity. The individual islands differ in their repertoire of metabolic genes that make a large contribution to the pangenome. In order to unravel intraclonal diversity of P. aeruginosa, the genomes of two members of the PA14 clonal complex from diverse habitats and geographic origin were compared. The genome sequences differed by less than 0.01% from each other. One hundred ninety-eight of the 231 single nucleotide substitutions (SNPs) were non-randomly distributed in the genome. Non-synonymous SNPs were mainly found in an integrated Pf1-like phage and in genes involved in transcriptional regulation, membrane and extracellular constituents, transport, and secretion. In summary, P. aeruginosa is endowed with a highly conserved core genome of low sequence diversity and a highly variable accessory genome that communicates with other pseudomonads and genera via horizontal gene transfer.

The IncP-1 plasmid backbone adapts to different host bacterial species and evolves through homologous recombination

Plasmids are important members of the bacterial mobile gene pool, and are among the most important contributors to horizontal gene transfer between bacteria. They typically harbour a wide spectrum of host beneficial traits, such as antibiotic resistance, inserted into their backbones. Although these inserted elements have drawn considerable interest, evolutionary information about the plasmid backbones, which encode plasmid related traits, is sparse. Here we analyse 25 complete backbone genomes from the broad-host-range IncP-1 plasmid family. Phylogenetic analysis reveals seven clades, in which two plasmids that we isolated from a marine biofilm represent a novel clade. We also found that homologous recombination is a prominent feature of the plasmid backbone evolution. Analysis of genomic signatures indicates that the plasmids have adapted to different host bacterial species. Globally circulating IncP-1 plasmids hence contain mosaic structures of segments derived from several parental plasmids that have evolved in, and adapted to, different, phylogenetically very distant host bacterial species.

Plasmids are present in many bacteria and are often transferred between different species causing horizontal gene transfer. By comparing the sequences of 25 plasmid DNA backbones, the authors show that homologous recombination is prevalent in plasmids and that the plasmids have adapted to persist in different host bacteria.

Diversity of glycosyl hydrolases from cellulose-depleting communities enriched from casts of two earthworm species

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-beta-glucanases, beta-glucosidases, beta-cellobiohydrolases, beta-galactosidase, and beta-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 beta-galactosidases/alpha-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).

Analysis of intra-genomic GC content homogeneity within prokaryotes

BACKGROUND

Bacterial genomes possess varying GC content (total guanines (Gs) and cytosines (Cs) per total of the four bases within the genome) but within a given genome, GC content can vary locally along the chromosome, with some regions significantly more or less GC rich than on average. We have examined how the GC content varies within microbial genomes to assess whether this property can be associated with certain biological functions related to the organism's environment and phylogeny. We utilize a new quantity GCVAR, the intra-genomic GC content variability with respect to the average GC content of the total genome. A low GCVAR indicates intra-genomic GC homogeneity and high GCVAR heterogeneity.

RESULTS

The regression analyses indicated that GCVAR was significantly associated with domain (i.e. archaea or bacteria), phylum, and oxygen requirement. GCVAR was significantly higher among anaerobes than both aerobic and facultative microbes. Although an association has previously been found between mean genomic GC content and oxygen requirement, our analysis suggests that no such association exits when phylogenetic bias is accounted for. A significant association between GCVAR and mean GC content was also found but appears to be non-linear and varies greatly among phyla.

CONCLUSIONS

Our findings show that GCVAR is linked with oxygen requirement, while mean genomic GC content is not. We therefore suggest that GCVAR should be used as a complement to mean GC content.

Comparative genomics of green sulfur bacteria

Eleven completely sequenced Chlorobi genomes were compared in oligonucleotide usage, gene contents, and synteny. The green sulfur bacteria (GSB) are equipped with a core genome that sustains their anoxygenic phototrophic lifestyle by photosynthesis, sulfur oxidation, and CO(2) fixation. Whole-genome gene family and single gene sequence comparisons yielded similar phylogenetic trees of the sequenced chromosomes indicating a concerted vertical evolution of large gene sets. Chromosomal synteny of genes is not preserved in the phylum Chlorobi. The accessory genome is characterized by anomalous oligonucleotide usage and endows the strains with individual features for transport, secretion, cell wall, extracellular constituents, and a few elements of the biosynthetic apparatus. Giant genes are a peculiar feature of the genera Chlorobium and Prosthecochloris. The predicted proteins have a huge molecular weight of 10(6), and are probably instrumental for the bacteria to generate their own intimate (micro)environment.

Differential decay of parent-of-origin-specific genomic sharing in cystic fibrosis-affected sib pairs maps a paternally imprinted locus to 7q34

Cystic fibrosis (CF) is a monogenic disease characterized by a high variability of disease severity and outcome that points to the role of environmental factors and modulating genes that shape the course of this multiorgan disease. We genotyped families of cystic fibrosis sib pairs homozygous for F508del-CFTR who represent extreme clinical phenotypes at informative microsatellite markers spanning a 38 Mb region between CFTR and 7qtel. Recombination events on both parental chromosomes were compared between siblings with concordant clinical phenotypes and siblings with discordant clinical phenotypes. Monitoring parent-of-origin-specific decay of genomic sharing delineated a 2.9-Mb segment on 7q34 in which excess of recombination on paternal chromosomes in discordant pairs was observed compared with phenotypically concordant sibs. This 2.9-Mb core candidate region was enriched in imprinting-related elements such as predicted CCCTC-binding factor consensus sites and CpG islands dense in repetitive elements. Moreover, allele frequencies at a microsatellite marker within the core candidate region differed significantly comparing mildly and severely affected cystic fibrosis sib pairs. The identification of this paternally imprinted locus on 7q34 as a modulator of cystic fibrosis disease severity shows that imprinted elements can be identified by straightforward fine mapping of break points in sib pairs with informative contrasting phenotypes.

Community-wide analysis of microbial genome sequence signatures

BACKGROUND

Analyses of DNA sequences from cultivated microorganisms have revealed genome-wide, taxa-specific nucleotide compositional characteristics, referred to as genome signatures. These signatures have far-reaching implications for understanding genome evolution and potential application in classification of metagenomic sequence fragments. However, little is known regarding the distribution of genome signatures in natural microbial communities or the extent to which environmental factors shape them.

RESULTS

We analyzed metagenomic sequence data from two acidophilic biofilm communities, including composite genomes reconstructed for nine archaea, three bacteria, and numerous associated viruses, as well as thousands of unassigned fragments from strain variants and low-abundance organisms. Genome signatures, in the form of tetranucleotide frequencies analyzed by emergent self-organizing maps, segregated sequences from all known populations sharing < 50 to 60% average amino acid identity and revealed previously unknown genomic clusters corresponding to low-abundance organisms and a putative plasmid. Signatures were pervasive genome-wide. Clusters were resolved because intra-genome differences resulting from translational selection or protein adaptation to the intracellular (pH approximately 5) versus extracellular (pH approximately 1) environment were small relative to inter-genome differences. We found that these genome signatures stem from multiple influences but are primarily manifested through codon composition, which we propose is the result of genome-specific mutational biases.

CONCLUSIONS

An important conclusion is that shared environmental pressures and interactions among coevolving organisms do not obscure genome signatures in acid mine drainage communities. Thus, genome signatures can be used to assign sequence fragments to populations, an essential prerequisite if metagenomics is to provide ecological and biochemical insights into the functioning of microbial communities.

Visualization of Pseudomonas genomic structure by abundant 8-14mer oligonucleotides

Under- and over-represented mono- to hexanucleotides are signatures of bacterial genomes, but the compositional biases of octa- to tetradecanucleotides have not yet been explored. Thirteen completely sequenced genomes of the Pseudomonas genus were searched for highly overrepresented 8-14mers. Between 59-989 overrepresented 8-14mers were found to exceed the applied threshold value. All genomic data sets of the 13 strains showed a consistent pattern, with individual oligomers clustering in either non-coding or coding regions. Non-coding oligonucleotides were typically part of longer repeats. Coding oligonucleotides were evenly distributed in the core genome, preferred one reading frame and matched with the local tetranucleotide usage patterns. Genomic islands were recognized by the depletion of overrepresented oligonucleotides. Several mainly coding 8-14mers occurred in genomes on average every 10 000 bp or less. Such frequently occurring 8-14mers could become useful markers for species identification. In the future of next-generation ultra-high throughput DNA sequencing, the composition of bacterial metagenomes may be quantified by scanning the primary sequence reads for these 8-14mer markers.

The SeqWord Genome Browser: an online tool for the identification and visualization of atypical regions of bacterial genomes through oligonucleotide usage

BACKGROUND

Data mining in large DNA sequences is a major challenge in microbial genomics and bioinformatics. Oligonucleotide usage (OU) patterns provide a wealth of information for large scale sequence analysis and visualization. The purpose of this research was to make OU statistical analysis available as a novel web-based tool for functional genomics and annotation. The tool is also available as a downloadable package.

RESULTS

The SeqWord Genome Browser (SWGB) was developed to visualize the natural compositional variation of DNA sequences. The applet is also used for identification of divergent genomic regions both in annotated sequences of bacterial chromosomes, plasmids, phages and viruses, and in raw DNA sequences prior to annotation by comparing local and global OU patterns. The applet allows fast and reliable identification of clusters of horizontally transferred genomic islands, large multi-domain genes and genes for ribosomal RNA. Within the majority of genomic fragments (also termed genomic core sequence), regions enriched with housekeeping genes, ribosomal proteins and the regions rich in pseudogenes or genetic vestiges may be contrasted.

CONCLUSION

The SWGB applet presents a range of comprehensive OU statistical parameters calculated for a range of bacterial species, plasmids and phages. It is available on the Internet at http://www.bi.up.ac.za/SeqWord/mhhapplet.php.

Utility of computational methods to identify the apoptosis machinery in unicellular eukaryotes

Apoptosis is the phenotypic result of an active, regulated process of self-destruction. Following various cellular insults, apoptosis has been demonstrated in numerous unicellular eukaryotes, but very little is known about the genes and proteins that initiate and execute this process in this group of organisms. A bioinformatic approach presents an array of powerful methods to direct investigators in the identification of the apoptosis machinery in protozoans. In this review, we discuss some of the available computational methods and illustrate how they may be applied using the identification of a Plasmodium falciparum metacaspase gene as an example.

Reliability and applications of statistical methods based on oligonucleotide frequencies in bacterial and archaeal genomes

Background

The increasing number of sequenced prokaryotic genomes contains a wealth of genomic data that needs to be effectively analysed. A set of statistical tools exists for such analysis, but their strengths and weaknesses have not been fully explored. The statistical methods we are concerned with here are mainly used to examine similarities between archaeal and bacterial DNA from different genomes. These methods compare observed genomic frequencies of fixed-sized oligonucleotides with expected values, which can be determined by genomic nucleotide content, smaller oligonucleotide frequencies, or be based on specific statistical distributions. Advantages with these statistical methods include measurements of phylogenetic relationship with relatively small pieces of DNA sampled from almost anywhere within genomes, detection of foreign/conserved DNA, and homology searches. Our aim was to explore the reliability and best suited applications for some popular methods, which include relative oligonucleotide frequencies (ROF), di- to hexanucleotide zero'th order Markov methods (ZOM) and 2.order Markov chain Method (MCM). Tests were performed on distant homology searches with large DNA sequences, detection of foreign/conserved DNA, and plasmid-host similarity comparisons. Additionally, the reliability of the methods was tested by comparing both real and random genomic DNA.

Results

Our findings show that the optimal method is context dependent. ROFs were best suited for distant homology searches, whilst the hexanucleotide ZOM and MCM measures were more reliable measures in terms of phylogeny. The dinucleotide ZOM method produced high correlation values when used to compare real genomes to an artificially constructed random genome with similar %GC, and should therefore be used with care. The tetranucleotide ZOM measure was a good measure to detect horizontally transferred regions, and when used to compare the phylogenetic relationships between plasmids and hosts, significant correlation (R² = 0.4) was found with genomic GC content and intra-chromosomal homogeneity.

Conclusion

The statistical methods examined are fast, easy to implement, and powerful for a number of different applications involving genomic sequence comparisons. However, none of the measures examined were superior in all tests, and therefore the choice of the statistical method should depend on the task at hand.

Global features of the Alcanivorax borkumensis SK2 genome

The global feature of the completely sequenced Alcanivorax borkumensis SK2 type strain chromosome is its symmetry and homogeneity. The origin and terminus of replication are located opposite to each other in the chromosome and are discerned with high signal to noise ratios by maximal oligonucleotide usage biases on the leading and lagging strand. Genomic DNA structure is rather uniform throughout the chromosome with respect to intrinsic curvature, position preference or base stacking energy. The orthologs and paralogs of A. borkumensis genes with the highest sequence homology were found in most cases among gamma-Proteobacteria, with Acinetobacter and P. aeruginosa as closest relatives. A. borkumensis shares a similar oligonucleotide usage and promoter structure with the Pseudomonadales. A comparatively low number of only 18 genome islands with atypical oligonucleotide usage was detected in the A. borkumensis chromosome. The gene clusters that confer the assimilation of aliphatic hydrocarbons, are localized in two genome islands which were probably acquired from an ancestor of the Yersinia lineage, whereas the alk genes of Pseudomonas putida still exhibit the typical Alcanivorax oligonucleotide signature indicating a complex evolution of this major hydrocarbonoclastic trait.

Circuit assemblages derived from net dinucleotide values provide a succinct identity for the HIV-1 genome and each of its genes

Dinucleotide composition has been recognized as a species-specific characteristic of organisms for more than 20 years. Lang (2000, Bioinformatics, 16, 212-221), found that in Monilinia rRNA a species-specific identity is conserved when dinucleotide counts are compressed into net dinucleotide counts (e.g., 50AC + 20CA = 30nAC) and clusters of net dinucleotides of equal value (e.g., 30nAC + 30nCT + 30nTA = 30ACTA) which were called circuits. This study evaluates circuit assemblages (CAs)--the collection of all net dinucleotide circuits derived from a sequence--in a diverse set of 110 HIV-1 genomes. The circuit composition, which is often based on <or= 15% of the total dinucleotides of a sequence, uniquely characterizes each gene and genome, although the pairwise similarity of the sequences is as low as 70%. Variations in net dinucleotide distributions are associated with structural and functional features of the genome and its proteins. Circuit values of the env signal sequence are different between subtypes that have remained localized and those that have become pandemic. CAs of complete genomes of HIV-1 are similar to other retro-transcribing viruses, and distinct from viroids and single- and double-stranded DNA and RNA viruses. CAs provide a succinct, quantitative, and species-specific description of DNA composition that is consistent with the results of traditional analytic methods at multiple levels of genome organization.

Secrets of soil survival revealed by the genome sequence of Arthrobacter aurescens TC1

Arthrobacter sp. strains are among the most frequently isolated, indigenous, aerobic bacterial genera found in soils. Member of the genus are metabolically and ecologically diverse and have the ability to survive in environmentally harsh conditions for extended periods of time. The genome of Arthrobacter aurescens strain TC1, which was originally isolated from soil at an atrazine spill site, is composed of a single 4,597,686 basepair (bp) circular chromosome and two circular plasmids, pTC1 and pTC2, which are 408,237 bp and 300,725 bp, respectively. Over 66% of the 4,702 open reading frames (ORFs) present in the TC1 genome could be assigned a putative function, and 13.2% (623 genes) appear to be unique to this bacterium, suggesting niche specialization. The genome of TC1 is most similar to that of Tropheryma, Leifsonia, Streptomyces, and Corynebacterium glutamicum, and analyses suggest that A. aurescens TC1 has expanded its metabolic abilities by relying on the duplication of catabolic genes and by funneling metabolic intermediates generated by plasmid-borne genes to chromosomally encoded pathways. The data presented here suggest that Arthrobacter's environmental prevalence may be due to its ability to survive under stressful conditions induced by starvation, ionizing radiation, oxygen radicals, and toxic chemicals.

Soil systems contain the greatest diversity of microorganisms on earth, with 5,000–10,000 species of microorganism per gram of soil. Arthrobacter sp. strains have a primitive life cycle and are among the most frequently isolated, indigenous soil bacteria, found in common and deep subsurface soils, arctic ice, and environments contaminated with industrial chemicals and radioactive materials. To better understand how these bacteria survive in environmentally harsh conditions, the authors used a structural genomics approach to identify genes involved in soil survival of Arthrobacter aurescens strain TC1, a bacterium originally isolated for its ability to degrade the herbicide atrazine. They found that the genome of this bacterium comprises a single circular chromosome and two plasmids that encode for a large number proteins involved in stress responses due to starvation, desiccation, oxygen radicals, and toxic chemicals. A. aurescens' metabolic versatility is in part due to the presence of duplicated catabolic genes and its ability to funnel plasmid-derived intermediates into chromosomally encoded pathways. Arthrobacter's array of genes that allow for survival in stressful conditions and its ability to produce a temperature-tolerant “cyst”-like resting cell render this soil microorganism able to survive and prosper in a variety of environmental conditions.

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.

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.

Widespread occurrence and genomic context of unusually small polyketide synthase genes in microbial consortia associated with marine sponges

Numerous marine sponges harbor enormous amounts of as-yet-uncultivated bacteria in their tissues. There is increasing evidence that these symbionts play an important role in the synthesis of protective metabolites, many of which are of great pharmacological interest. In this study, genes for the biosynthesis of polyketides, one of the most important classes of bioactive natural products, were systematically investigated in 20 demosponge species from different oceans. Unexpectedly, the sponge metagenomes were dominated by a ubiquitously present, evolutionarily distinct, and highly sponge-specific group of polyketide synthases (PKSs). Open reading frames resembling animal fatty acid genes were found on three corresponding DNA regions isolated from the metagenomes of Theonella swinhoei and Aplysina aerophoba. Their architecture suggests that methyl-branched fatty acids are the metabolic product. According to a phylogenetic analysis of housekeeping genes, at least one of the PKSs belongs to a bacterium of the Deinococcus-Thermus phylum. The results provide new insights into the chemistry of sponge symbionts and allow inference of a detailed phylogeny of the diverse functional PKS types present in sponge metagenomes. Based on these qualitative and quantitative data, we propose a significantly simplified strategy for the targeted isolation of biomedically relevant PKS genes from complex sponge-symbiont associations.