Software for analysis of bacterial genomes

GenomeExplorer is a program for comparative analysis of regulation in prokaryotic genomes. The program has options for signal search, comparison of gene samples, search for paralogs and orthologs, iterative construction of signal profiles. The program has a convenient graphic interface, allowing for navigation in the annotation window, in the genome map, and in the table of gene similarities. The use of the system clipboard allows one to export the results of analysis into Word and Excel, and to call external programs via the Internet.     

A large-scale analysis of codon usage bias in 4868 bacterial genomes shows association of codon adaptation index with GC content,protein functional domains and bacterial phenotypes

Multiple synonymous codons code for the same amino acid, resulting in the degeneracy of the genetic code and in the preferred used of some codons called codon bias usage (CBU). We performed a large-scale analysis of codon usage bias analysing the distribution of the codon adaptation index (CAI) and the codon relative adaptiveness index (RA) in 4868 bacterial genomes. We found that CAI values differ significantly between protein functional domains and part of the protein outside domains and show how CAI, GC content and preferred usage of polymerase III alpha subunits are related. Additionally, we give evidence of the association between CAI and bacterial phenotypes.     

The Neolithic revolution of bacterial genomes

Current human activities undoubtedly impact natural ecosystems. However, the influence of Homo sapiens on living organisms must have also occurred in the past. Certain genomic characteristics of prokaryotes can be used to study the impact of ancient human activities on microorganisms. By analyzing DNA sequence similarity features of transposable elements, dramatic genomic changes have been identified in bacteria that are associated with large and stable human communities, agriculture and animal domestication: three features unequivocally linked to the Neolithic revolution. It is hypothesized that bacteria specialized in human-associated niches underwent an intense transformation after the social and demographic changes that took place with the first Neolithic settlements. These genomic changes are absent in related species that are not specialized in humans.     

Fold predictions for bacterial genomes

Fold assignments for newly sequenced genomes belong to the most important and interesting applications of the booming field of protein structure prediction. We present a brief survey and a discussion of such assignments completed to date, using as an example several fold assignment projects for proteins from the Escherichia coli genome. This review focuses on steps that are necessary to go beyond the simple assignment projects and into the development of tools extending our understanding of functions of proteins in newly sequenced genomes. This paper also discusses several problems seldom addressed in the literature, such as the problem of domain prediction and complementary predictions (e.g., transmembrane regions and flexible regions) and cross-correlation of predictions from different servers. The influence of sequence and structure database growth on prediction success is also addressed. Finally, we discuss the perspectives of the field in the context of massive sequence and structure determination projects, as well as the development of novel prediction methods.     

Algorithm for large-scale clustering across multiple genomes

Identifying genomic regions that descended from a common ancestor helps us study the gene function and genome evolution. In distantly related genomes, clusters of homologous gene pairs are evidently used in function prediction, operon detection, etc. Currently, there are many kinds of computational methods that have been proposed defining gene clusters to identify gene families and operons. However, most of those algorithms are only available on a data set of small size. We developed an efficient gene clustering algorithm that can be applied on hundreds of genomes at the same time. This approach allows for large-scale study of evolutionary relationships of gene clusters and study of operon formation and destruction. An analysis of proposed algorithms shows that more biological insight can be obtained by analyzing gene clusters across hundreds of genomes, which can help us understand operon occurrences, gene orientations and gene rearrangements.     

A novel system for large-scale gene expression analysis: bacterial colonies array

In the present work, we report the use of bacterial colonies to optimize macroarray technique. The devised system is significantly cheaper than other methods available to detect large-scale differential gene expression. Recombinant Escherichia coli clones containing plasmid-encoded copies of 4,608 individual expressed sequence tag (ESTs) were robotically spotted onto nylon membranes that were incubated for 6 and 12 h to allow the bacteria to grow and, consequently, amplify the cloned ESTs. The membranes were then hybridized with a beta-lactamase gene specific probe from the recombinant plasmid and, subsequently, phosphorimaged to quantify the microbial cells. Variance analysis demonstrated that the spot hybridization signal intensity was similar for 3,954 ESTs (85.8%) after 6 h of bacterial growth. Membranes spotted with bacteria colonies grown for 12 h had 4,017 ESTs (87.2%) with comparable signal intensity but the signal to noise ratio was fivefold higher. Taken together, the results of this study indicate that it is possible to investigate large-scale gene expression using macroarrays based on bacterial colonies grown for 6 h onto membranes.     

Junker: an intergenic explorer for bacterial genomes

In the past few decades, scientists from all over the world have taken a keen interest in novel functional units such as small regulatory RNAs, small open reading frames, pseudogenes, transposons, integrase binding attB/attP sites, repeat elements within the bacterial intergenic regions (IGRs) and in the analysis of those "junk" regions for genomic complexity. Here we have developed a web server, named Junker, to facilitate the in-depth analysis of IGRs for examining their length distribution, four-quadrant plots, GC percentage and repeat details. Upon selection of a particular bacterial genome, the physical genome map is displayed as a multiple loci with options to view any loci of interest in detail. In addition, an IGR statistics module has been created and implemented in the web server to analyze the length distribution of the IGRs and to understand the disordered grouping of IGRs across the genome by generating the four-quadrant plots. The proposed web server is freely available at the URL http://pranag.physics.iisc.ernet.in/junker/.     

CNproScan: Hybrid CNV detection for bacterial genomes

Discovering copy number variation (CNV) in bacteria is not in the spotlight compared to the attention focused on CNV detection in eukaryotes. However, challenges arising from bacterial drug resistance bring further interest to the topic of CNV and its role in drug resistance. General CNV detection methods do not consider bacteria's features and there is space to improve detection accuracy. Here, we present a CNV detection method called CNproScan focused on bacterial genomes. CNproScan implements a hybrid approach and other bacteria-focused features and depends only on NGS data. We benchmarked our method and compared it to the previously published methods and we can resolve to achieve a higher detection rate together with providing other beneficial features, such as CNV classification. Compared with other methods, CNproScan can detect much shorter CNV events.     

The most deviated codon position in AT-rich bacterial genomes: a function related analysis

We have performed systematic study on more than 120 archaeal and bacterial genomes. Based on the index proposed in the current paper, clear patterns are observed showing the relation between the base compositional deviation at three codon positions and the genomic GC content. For AT-rich genomes, the Most Deviated Codon Position (MDCP) is the 1st codon position, while for GC-rich genomes, MDCP appears at the 2nd or 3rd codon position alternatively. According to MDCP, the CDSs of a genome can be classified into two types: typical and atypical. In AT-rich genomes the typical represent the majority and account for about 3/4 of all the CDSs. Based on the functional classification of COG database, the two types of CDSs are examined. An apparent bias of distribution is observed that the CDSs with the function of 'information processing' are more likely to present in typical type.     

The generation of variation in bacterial genomes

In the context of a general overview of molecular mechanisms of microbial evolution, several genetic systems known to either promote or restrain the generation of genetic variations are discussed. Particular attention is given to functions involved in DNA rearrangements and DNA acquisition. Sporadic actions by a variety of such systems influencing genetic stability in either way result in a level of genetic plasticity which is tolerable to the overall wealth of microbial populations but which allows for evolutionary change needed for a steady adaptation to variable selective forces. Although these evolutionarily relevant biological functions are encoded by the genome of each individual, their actions are exerted to some degree randomly in rare individuals and are therefore seemingly nondeterministic and become manifest at the population level.     

Conformational analysis of invariant peptide sequences in bacterial genomes

The functional significance of evolutionarily conserved motifs/patterns of short regions in proteins is well documented. Although a large number of sequences are conserved, only a small fraction of these are invariant across several organisms. Here, we have examined the structural features of the functionally important peptide sequences, which have been found invariant across diverse bacterial genera. Ramachandran angles (phi,psi) have been used to analyze the conformation, folding patterns and geometrical location (buried/exposed) of these invariant peptides in different crystal structures harboring these sequences. The analysis indicates that the peptides preferred a single conformation in different protein structures, with the exception of only a few longer peptides that exhibited some conformational variability. In addition, it is noticed that the variability of conformation occurs mainly due to flipping of peptide units about the virtual C(alpha)...C(alpha) bond. However, for a given invariant peptide, the folding patterns are found to be similar in almost all the cases. Over and above, such peptides are found to be buried in the protein core. Thus, we can safely conclude that these invariant peptides are structurally important for the proteins, since they acquire unique structures across different proteins and can act as structural determinants (SD) of the proteins. The location of these SD peptides on the protein chain indicated that most of them are clustered towards the N-terminal and middle region of the protein with the C-terminal region exhibiting low preference. Another feature that emerges out of this study is that some of these SD peptides can also play the roles of "fold boundaries" or "hinge nucleus" in the protein structure. The study indicates that these SD peptides may act as chain-reversal signatures, guiding the proteins to adopt appropriate folds. In some cases the invariant signature peptides may also act as folding nuclei (FN) of the proteins.     

Mining bacterial genomes for antimicrobial targets

The elucidation of whole-genome sequences is expected to have a revolutionary impact on the discovery of novel medicines. With the availability of complete genome sequences of more than 30 different species, the field of antimicrobial drug discovery has the opportunity to access a remarkable diversity of genomic information. In this review, I summarize how microbial genomics has changed strategies of drug discovery by applying bioinformatics, novel genetic approaches and genomics-based technologies, including analysis of gene expression using DNA microarrays.     

The evolution of domain-content in bacterial genomes

Background  

Across all sequenced bacterial genomes, the number of domains n _c in different functional categories c scales as a power-law in the total number of domains n, i.e. , with exponents α _c that vary across functional categories. Here we investigate the implications of these scaling laws for the evolution of domain-content in bacterial genomes and derive the simplest evolutionary model consistent with these scaling laws.     

Evolution of bacterial genomes

This review examines evolution of bacterial genomes with an emphasis on RNA based life, the transition to functional DNA and small evolving genomes (possibly plasmids) that led to larger, functional bacterial genomes.     

Statistical analysis of complete bacterial genomes: Avoidance of palindromes and restriction-modification systems

Avoidance of 4-, 5-, and 6-letter palindromes is observed in many prokaryotic genomes. A large fraction of such palindromes is formed by restriction sites of the species itself or a closely related species. One possible reason for that is the horizontal transfer of genes encoding restriction-modification systems. In organisms isolated from the action of such systems (e.g., in Mycoplasma), palindromes are not avoided. The general tendencies in preferences and avoidance of palindromes were studied for 33 available prokaryotic genomes. The results obtained provide additional insight into the relationships within and between taxonomic groups.