PGAAS: a prokaryotic genome assembly assistant system

MOTIVATION: In order to accelerate the finishing phase of genome assembly, especially for the whole genome shotgun approach of prokaryotic species, we have developed a software package designated prokaryotic genome assembly assistant system (PGAAS). The approach upon which PGAAS is based is to confirm the order of contigs and fill gaps between contigs through peptide links obtained by searching each contig end with BLASTX against protein databases. RESULTS: We used the contig dataset of the cyanobacterium Synechococcus sp. strain PCC7002 (PCC7002), which was sequenced with six-fold coverage and assembled using the Phrap package. The subject database is the protein database of the cyanobacterium, Synechocystis sp. strain PCC6803 (PCC6803). We found more than 100 non-redundant peptide segments which can link at least 2 contigs. We tested one pair of linked contigs by sequencing and obtained satisfactory result. PGAAS provides a graphic user interface to show the bridge peptides and pier contigs. We integrated Primer3 into our package to design PCR primers at the adjacent ends of the pier contigs. AVAILABILITY: We tested PGAAS on a Linux (Redhat 6.2) PC machine. It is developed with free software (MySQL, PHP and Apache). The whole package is distributed freely and can be downloaded as UNIX compress file: ftp://ftp.cbi.pku.edu.cn/pub/software/unix/pgaas1.0.tar.gz. The package is being continually updated.     

Building a genome engineering toolbox in nonmodel prokaryotic microbes

The realization of a sustainable bioeconomy requires our ability to understand and engineer complex design principles for the development of platform organisms capable of efficient conversion of cheap and sustainable feedstocks (e.g., sunlight, CO₂, and nonfood biomass) into biofuels and bioproducts at sufficient titers and costs. For model microbes, such as Escherichia coli, advances in DNA reading and writing technologies are driving the adoption of new paradigms for engineering biological systems. Unfortunately, microbes with properties of interest for the utilization of cheap and renewable feedstocks, such as photosynthesis, autotrophic growth, and cellulose degradation, have very few, if any, genetic tools for metabolic engineering. Therefore, it is important to develop “design rules” for building a genetic toolbox for novel microbes. Here, we present an overview of our current understanding of these rules for the genetic manipulation of prokaryotic microbes and the available genetic tools to expand our ability to genetically engineer nonmodel systems.     

Principles of prokaryotic genome organization]

The peculiarities of bacterial chromosome organization are discussed, based mainly on the data on Escherichia coli. Highly important for bacterial genome organization is its division into two approx. equal half-genomes undergoing periodically "exchanges" of some kind displayed as continuous inversions including the oriC region of replication initiation. It is believed that short oligonucleotides are comprised in either of half-genomes. The former are predominantly oriented as direct repeats, which ensures the possibility of formation of tandem duplications consisting of identical genes--under conditions when selection for enhancing functions of corresponding genes takes place. Multiple tandem duplications capable of excision of plasmatic gene copies seem to initiate horizontal gene transfer in bacteria. Tandem gene duplications are probably being formed in the process of bacterial genetic recombination as well, when, as a result of non-equal crossing over, gene alleles derived from different strains are united into a tandem.     

Microbial genescapes: a prokaryotic view of the yeast genome

We examine the translated open reading frames (ORFs) of the yeast Saccharomyces cerevisiae, focusing on those that have FASTA matches in phyletically defined sets of completely sequenced genomes. On this basis, we identify archaeal yeast, bacterial yeast, universal yeast, and yeast ORFs that do not have a match in any of nine prokaryote genomes. Similarly, we examine the yeast mitochondrial genome and the subset of the yeast nuclear ORFs identified as being involved in mitochondrial biogenesis. For the yeast ORFs that match one or more ORFs in these prokaryote genomes, we examine the phyletic and functional distributions of these matches as a function of match strength. These results provide genome level insights into the origin of the eukaryotic cell and the origin of mitochondria. More generally, they exemplify how the growing database of prokaryote genome sequences can help us understand eukaryote genomes.     

Predicting prokaryotic ecological niches using genome sequence analysis

Automated DNA sequencing technology is so rapid that analysis has become the rate-limiting step. Hundreds of prokaryotic genome sequences are publicly available, with new genomes uploaded at the rate of approximately 20 per month. As a result, this growing body of genome sequences will include microorganisms not previously identified, isolated, or observed. We hypothesize that evolutionary pressure exerted by an ecological niche selects for a similar genetic repertoire in those prokaryotes that occupy the same niche, and that this is due to both vertical and horizontal transmission. To test this, we have developed a novel method to classify prokaryotes, by calculating their Pfam protein domain distributions and clustering them with all other sequenced prokaryotic species. Clusters of organisms are visualized in two dimensions as 'mountains' on a topological map. When compared to a phylogenetic map constructed using 16S rRNA, this map more accurately clusters prokaryotes according to functional and environmental attributes. We demonstrate the ability of this map, which we term a "niche map", to cluster according to ecological niche both quantitatively and qualitatively, and propose that this method be used to associate uncharacterized prokaryotes with their ecological niche as a means of predicting their functional role directly from their genome sequence.     

Concept of a bacterial species and the evolution of the prokaryotic genome

The present concepts of evolution and species delineation in prokaryotes are considered. Recently a considerable extension of knowledge on the processes of microevolution of medically significant bacteria was noted alongside with the importance of horizontal and lateral transfer of genes. The phylophenetic concept of species was considered in detail. The inclusion of the ecological criterion into a phylophenetic concept of a species is supposed to facilitate the development of more adequate notion on the evolution of bacteria, the improvement of species delineation in prokaryotes, their classification and nomenclature.     

Prophage Finder: a prophage loci prediction tool for prokaryotic genome sequences

Prophage loci often remain under-annotated or even unrecognized in prokaryotic genome sequencing projects. A PHP application, Prophage Finder, has been developed and implemented to predict prophage loci, based upon clusters of phage-related gene products encoded within DNA sequences. This application provides results detailing several facets of these clusters to facilitate rapid prediction and analysis of prophage sequences. Prophage Finder was tested using previously annotated prokaryotic genomic sequences with manually curated prophage loci as benchmarks. Additional analyses from Prophage Finder searches of several draft prokaryotic genome sequences are available through the Web site (http://bioinformatics.uwp.edu/~phage/DOEResults.php) to illustrate the potential of this application.     

Large-scale prokaryotic gene prediction and comparison to genome annotation

MOTIVATION: Prokaryotic genomes are sequenced and annotated at an increasing rate. The methods of annotation vary between sequencing groups. It makes genome comparison difficult and may lead to propagation of errors when questionable assignments are adapted from one genome to another. Genome comparison either on a large or small scale would be facilitated by using a single standard for annotation, which incorporates a transparency of why an open reading frame (ORF) is considered to be a gene. RESULTS: A total of 143 prokaryotic genomes were scored with an updated version of the prokaryotic genefinder EasyGene. Comparison of the GenBank and RefSeq annotations with the EasyGene predictions reveals that in some genomes up to approximately 60% of the genes may have been annotated with a wrong start codon, especially in the GC-rich genomes. The fractional difference between annotated and predicted confirms that too many short genes are annotated in numerous organisms. Furthermore, genes might be missing in the annotation of some of the genomes. We predict 41 of 143 genomes to be over-annotated by >5%, meaning that too many ORFs are annotated as genes. We also predict that 12 of 143 genomes are under-annotated. These results are based on the difference between the number of annotated genes not found by EasyGene and the number of predicted genes that are not annotated in GenBank. We argue that the average performance of our standardized and fully automated method is slightly better than the annotation.     

MetaGene: prokaryotic gene finding from environmental genome shotgun sequences

Exhaustive gene identification is a fundamental goal in all metagenomics projects. However, most metagenomic sequences are unassembled anonymous fragments, and conventional gene-finding methods cannot be applied. We have developed a prokaryotic gene-finding program, MetaGene, which utilizes di-codon frequencies estimated by the GC content of a given sequence with other various measures. MetaGene can predict a whole range of prokaryotic genes based on the anonymous genomic sequences of a few hundred bases, with a sensitivity of 95% and a specificity of 90% for artificial shotgun sequences (700 bp fragments from 12 species). MetaGene has two sets of codon frequency interpolations, one for bacteria and one for archaea, and automatically selects the proper set for a given sequence using the domain classification method we propose. The domain classification works properly, correctly assigning domain information to more than 90% of the artificial shotgun sequences. Applied to the Sargasso Sea dataset, MetaGene predicted almost all of the annotated genes and a notable number of novel genes. MetaGene can be applied to wide variety of metagenomic projects and expands the utility of metagenomics.     

Complex prokaryotic genome structure: rapid evolution of chromosome II

Although many bacteria with two chromosomes have been sequenced, the roles of such complex genome structuring are still unclear. To uncover levels of chromosome I (CI) and chromosome II (CII) sequence divergence, Mauve 2.2.0 was used to align the CI- and CII-specific sequences of bacteria with complex genome structuring in two sets of comparisons: the first set was conducted among the CI and CII of bacterial strains of the same species, while the second set was conducted among the CI and CII of species in Alphaproteobacteria that possess two chromosomes. The analyses revealed a rapid evolution of CII-specific DNA sequences compared with CI-specific sequences in a majority of organisms. In addition, levels of protein divergence between CI-specific and CII-specific genes were determined using phylogenetic analyses and confirmed the DNA alignment findings. Analysis of synonymous and nonsynonymous substitutions revealed that the structural and functional constraints on CI and CII genes are not significantly different. Also, horizontal gene transfer estimates in selected organisms demonstrated that CII in many species has acquired higher levels of horizontally transferred segments than CI. In summary, rapid evolution of CII may perform particular roles for organisms such as aiding in adapting to specialized niches.     

PeerGAD: a peer-review-based and community-centric web application for viewing and annotating prokaryotic genome sequences

PeerGAD is a web-based database-driven application that allows community-wide peer-reviewed annotation of prokaryotic genome sequences. The application was developed to support the annotation of the Pseudomonas syringae pv. tomato strain DC3000 genome sequence and is easily portable to other genome sequence annotation projects. PeerGAD incorporates several innovative design and operation features and accepts annotations pertaining to gene naming, role classification, gene translation and annotation derivation. The annotator tool in PeerGAD is built around a genome browser that offers users the ability to search and navigate the genome sequence. Because the application encourages annotation of the genome sequence directly by researchers and relies on peer review, it circumvents the need for an annotation curator while providing added value to the annotation data. Support for the Gene Ontology vocabulary, a structured and controlled vocabulary used in classification of gene roles, is emphasized throughout the system. Here we present the underlying concepts integral to the functionality of PeerGAD.     

The elusive short gene--an ensemble method for recognition for prokaryotic genome

Glucagon-like peptide-1 (GLP-1) and leptin are anorectic hormones produced in the small intestine and white adipose tissue, respectively. Investigating how these hormones act together as an integrated anorectic signal is important to elucidate a mechanism to maintain energy balance. In the present study, coadministration of subthreshold GLP-1 and leptin dramatically reduced feeding in rats. Although coadministration of GLP-1 with leptin did not enhance leptin signal transduction in the hypothalamus, it significantly decreased phosphorylation of AMP-activated protein kinase (AMPK). In addition, coadministration of GLP-1 with leptin significantly increased proopiomelanocortin (POMC) mRNA levels. Considering that α-melanocortin stimulating hormone (α-MSH) is derived from POMC and functions through the melanocortin-4-receptor (MC4-R) as a key molecule involved in feeding reduction, the interaction of GLP-1 and leptin on feeding reduction may be mediated through the α-MSH/MC4-R system. As expected, the interaction of GLP-1 and leptin was abolished by intracerebroventricular preadministration of the MC4-R antagonists agouti-related peptide and SHU9119. Taken together, GLP-1 and leptin cooperatively reduce feeding at least in part via inhibition of AMPK following binding of α-MSH to MC4-R.     

Evolution of prokaryotic gene order: genome rearrangements in closely related species

Conservation of gene order in prokaryotes has become important in predicting protein function because, over the evolutionary timescale, genomes are shuffled so that local gene-order conservation reflects the functional constraints within the protein. Here, we compare closely related genomes to identify the rate with which gene order is disrupted and to infer the genes involved in the genome rearrangement.     

原核基因组的进化方向和癌／正常细胞的熵产生

本文评述了理论生物学的两个基础实验。一是用离子束辐照加速大肠杆菌进化的方法研究原核基因组的进化方向,证明进化中的缺失偏好性和编码信息量扩增律不矛盾。二是提出用导热法测量细胞的熵产生,比较研究了癌细胞和正常细胞的熵产生和外加电场的关系,证明在一定强度的电场作用下正常细胞的熵产生可以明显超过癌细胞,从而实现改变两类细胞间熵流方向的目的。     

Energy beyond food: foraging theory informs time spent in thermals by a large soaring bird

Current understanding of how animals search for and exploit food resources is based on microeconomic models. Although widely used to examine feeding, such constructs should inform other energy-harvesting situations where theoretical assumptions are met. In fact, some animals extract non-food forms of energy from the environment, such as birds that soar in updraughts. This study examined whether the gains in potential energy (altitude) followed efficiency-maximising predictions in the world's heaviest soaring bird, the Andean condor (Vultur gryphus). Animal-attached technology was used to record condor flight paths in three-dimensions. Tracks showed that time spent in patchy thermals was broadly consistent with a strategy to maximise the rate of potential energy gain. However, the rate of climb just prior to leaving a thermal increased with thermal strength and exit altitude. This suggests higher rates of energetic gain may not be advantageous where the resulting gain in altitude would lead to a reduction in the ability to search the ground for food. Consequently, soaring behaviour appeared to be modulated by the need to reconcile differing potential energy and food energy distributions. We suggest that foraging constructs may provide insight into the exploitation of non-food energy forms, and that non-food energy distributions may be more important in informing patterns of movement and residency over a range of scales than previously considered.     

Towards a prokaryotic genomic taxonomy

One of the most interesting developments in the field of modern-day microbiology is the ever increasing number of whole-genome sequences that is publicly available. There is an increasing interest in the use of these genome sequences to assess evolutionary relationships among microbial taxa, as it is anticipated that much additional taxonomic information can be extracted from these sequences. In a first part of the present review, mechanisms that are responsible for the evolution of genomes will be discussed. Subsequently, we will give an overview of approaches that are presently available to assess the taxonomic relationships between prokaryotic species based on complete genome sequences, followed by a brief discussion of the potential implications of these novel approaches for bacterial taxonomy in general and our thinking about the bacterial species concept in particular.