GC-rich Genomes Research Articles

Quantitative frame analysis and the annotation of GC-rich (and other) prokaryotic genomes. An application to Anaeromyxobacter dehalogenans

Motivation: Graphical representations of contrasts in GC usage among codon frame positions (frame analysis) provide evidence of genes missing from the annotations of prokaryotic genomes of high GC content but the qualitative approach of visual frame analysis prevents its applicability on a genomic scale.Results: We developed two quantitative methods for the identification and statistical characterization in sequence regions of three-base periodicity (hits) associated with open reading frame structures. The methods were implemented in the N-Profile Analysis Computational Tool (NPACT), which highlights in graphical representations inconsistencies between newly identified ORFs and pre-existing annotations of coding-regions. We applied the NPACT procedures to two recently annotated strains of the deltaproteobacterium Anaeromyxobacter dehalogenans, identifying in both genomes numerous conserved ORFs not included in the published annotation of coding regions.Availability and implementation: NPACT is available as a web-based service and for download at http://genome.ufl.edu/npact.Contact: lucianob@ufl.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Coupling Between Protein Level Selection and Codon Usage Optimization in the Evolution of Bacteria and Archaea

ABSTRACTThe relationship between the selection affecting codon usage and selection on protein sequences of orthologous genes in diverse groups of bacteria and archaea was examined by using the Alignable Tight Genome Clusters database of prokaryote genomes. The codon usage bias is generally low, with 57.5% of the gene-specific optimal codon frequencies (Fopt) being below 0.55. This apparent weak selection on codon usage contrasts with the strong purifying selection on amino acid sequences, with 65.8% of the gene-specific dN/dS ratios being below 0.1. For most of the genomes compared, a limited but statistically significant negative correlation between Fopt and dN/dS was observed, which is indicative of a link between selection on protein sequence and selection on codon usage. The strength of the coupling between the protein level selection and codon usage bias showed a strong positive correlation with the genomic GC content. Combined with previous observations on the selection for GC-rich codons in bacteria and archaea with GC-rich genomes, these findings suggest that selection for translational fine-tuning could be an important factor in microbial evolution that drives the evolution of genome GC content away from mutational equilibrium. This type of selection is particularly pronounced in slowly evolving, “high-status” genes. A significantly stronger link between the two aspects of selection is observed in free-living bacteria than in parasitic bacteria and in genes encoding metabolic enzymes and transporters than in informational genes. These differences might reflect the special importance of translational fine-tuning for the adaptability of gene expression to environmental changes. The results of this work establish the coupling between protein level selection and selection for translational optimization as a distinct and potentially important factor in microbial evolution.

Comprehensive analysis of DNA polymerase III α subunits and their homologs in bacterial genomes

The analysis of ∼2000 bacterial genomes revealed that they all, without a single exception, encode one or more DNA polymerase III α-subunit (PolIIIα) homologs. Classified into C-family of DNA polymerases they come in two major forms, PolC and DnaE, related by ancient duplication. While PolC represents an evolutionary compact group, DnaE can be further subdivided into at least three groups (DnaE1-3). We performed an extensive analysis of various sequence, structure and surface properties of all four polymerase groups. Our analysis suggests a specific evolutionary pathway leading to PolC and DnaE from the last common ancestor and reveals important differences between extant polymerase groups. Among them, DnaE1 and PolC show the highest conservation of the analyzed properties. DnaE3 polymerases apparently represent an ‘impaired’ version of DnaE1. Nonessential DnaE2 polymerases, typical for oxygen-using bacteria with large GC-rich genomes, have a number of features in common with DnaE3 polymerases. The analysis of polymerase distribution in genomes revealed three major combinations: DnaE1 either alone or accompanied by one or more DnaE2s, PolC + DnaE3 and PolC + DnaE1. The first two combinations are present in Escherichia coli and Bacillus subtilis, respectively. The third one (PolC + DnaE1), found in Clostridia, represents a novel, so far experimentally uncharacterized, set.

Amino Acid Usage Is Asymmetrically Biased in AT- and GC-Rich Microbial Genomes

IntroductionGenomic base composition ranges from less than 25% AT to more than 85% AT in prokaryotes. Since only a small fraction of prokaryotic genomes is not protein coding even a minor change in genomic base composition will induce profound protein changes. We examined how amino acid and codon frequencies were distributed in over 2000 microbial genomes and how these distributions were affected by base compositional changes. In addition, we wanted to know how genome-wide amino acid usage was biased in the different genomes and how changes to base composition and mutations affected this bias. To carry this out, we used a Generalized Additive Mixed-effects Model (GAMM) to explore non-linear associations and strong data dependences in closely related microbes; principal component analysis (PCA) was used to examine genomic amino acid- and codon frequencies, while the concept of relative entropy was used to analyze genomic mutation rates.ResultsWe found that genomic amino acid frequencies carried a stronger phylogenetic signal than codon frequencies, but that this signal was weak compared to that of genomic %AT. Further, in contrast to codon usage bias (CUB), amino acid usage bias (AAUB) was differently distributed in AT- and GC-rich genomes in the sense that AT-rich genomes did not prefer specific amino acids over others to the same extent as GC-rich genomes. AAUB was also associated with relative entropy; genomes with low AAUB contained more random mutations as a consequence of relaxed purifying selection than genomes with higher AAUB.ConclusionGenomic base composition has a substantial effect on both amino acid- and codon frequencies in bacterial genomes. While phylogeny influenced amino acid usage more in GC-rich genomes, AT-content was driving amino acid usage in AT-rich genomes. We found the GAMM model to be an excellent tool to analyze the genomic data used in this study.

Phyletic distribution and conservation of the bacterial transcription termination factor Rho

Transcription termination factor Rho is a ring-shaped, ATP-dependent molecular motor that targets hundreds of transcription units in Escherichia coli. Interest in Rho was renewed recently on the realization that this essential factor is involved in multiple interactions and cellular processes that protect the E. coli genome and regulate its expression on a global scale. Yet it is currently unknown if (and how) Rho-dependent mechanisms are conserved throughout the bacterial kingdom. Here, we mined public databases to assess the distribution, expression and structural conservation of Rho across bacterial phyla. We found that rho is present in more than 90 % of sequenced bacterial genomes, although Cyanobacteria, Mollicutes and a fraction of Firmicutes are totally devoid of rho. Genomes lacking rho tend to be small and AT-rich and often belong to species with parasitic/symbiotic lifestyles (such as Mollicutes). By contrast, large GC-rich genomes, such as those of Actinobacteria, often contain rho duplicates and/or encode Rho proteins that bear insertion domains of unknown function(s). Notwithstanding, most Rho sequences across taxa contain canonical RNA-binding and ATP hydrolysis signature motifs, a feature suggestive of largely conserved mechanism(s) of action. Mutations that impair binding of bicyclomycin are present in ~5 % of rho sequences, implying that species from diverse ecosystems have developed resistance against this natural antibiotic. Altogether, these findings assert that Rho function is widespread among bacteria and suggest that it plays a particularly relevant role in the expression of complex genomes and/or bacterial adaptation to changing environments.

An efficient rRNA removal method for RNA sequencing in GC-rich bacteria

BackgroundNext generation sequencing (NGS) technologies have revolutionized gene expression studies and functional genomics analysis. However, further improvement of RNA sequencing protocols is still desirable, in order to reduce NGS costs and to increase its accuracy. In bacteria, a major problem in RNA sequencing is the abundance of ribosomal RNA (rRNA), which accounts for 95-98% of total RNA and can therefore hinder sufficient coverage of mRNA, the main focus of transcriptomic studies. Thus, efficient removal of rRNA is necessary to achieve optimal coverage, good detection sensitivity and reliable results. An additional challenge is presented by microorganisms with GC-rich genomes, in which rRNA removal is less efficient.ResultsIn this work, we tested two commercial kits for rRNA removal, either alone or in combination, on Burkholderia thailandensis. This bacterium, chosen as representative of the important Burkholderia genus, which includes both pathogenic and environmental bacteria, has a rather large (6.72 Mb) and GC-rich (67.7%) genome. Each enriched mRNA sample was sequenced through paired-end Illumina GAIIx run in duplicate, yielding between 10 and 40 million reads. We show that combined treatment with both kits allows an mRNA enrichment of more than 238-fold, enabling the sequencing of almost all (more than 90%) B. thailandensis transcripts from less than 10 million reads, without introducing any bias in mRNA relative abundance, thus preserving differential expression profile.ConclusionsThe mRNA enrichment protocol presented in this work leads to an increase in detection sensitivity up to 770% compared to total RNA; such increased sensitivity allows for a corresponding reduction in the number of sequencing reads necessary for the complete analysis of whole transcriptome expression profiling. Thus we can conclude that the MICROBExpress/Ovation combined rRNA removal method could be suitable for RNA sequencing of whole transcriptomes of microorganisms with high GC content and complex genomes enabling at the same time an important scaling down of sequencing costs.

Nuclear DNA content variation among perennial taxa of the genus Cyanus (Asteraceae) in Central Europe and adjacent areas

The genome size of 265 plants and the GC content of 126 plants from 63 populations of the Cyanus triumfetti and Cyanus montanus groups, collected across the Carpathians, Pannonia, Bohemian Massif, and Western and Dinaric Alps were determined by PI and DAPI flow cytometry. Variation of the nuclear DNA content among homoploid species, and intraspecific and interpopulation variation were confirmed in simultaneous analyses. The 2C-value at the diploid level (the C. triumfetti group) varied from 2.53 for Cyanus dominii subsp. sokolensis to 3.06 pg for C. triumfetti s.s. (1.21-fold range). At the tetraploid level (the C. montanus group), the 2C-value varied from 5.19 for Cyanus mollis to 5.84 pg for C. montanus (1.13-fold range). High intraspecific and interpopulation variation in the amount of nuclear DNA in the C. triumfetti group correlates with the extensive morphological variation found in this group. Significant between-species differences in genome size indicate that this attribute may be used as a supportive taxonomic marker for both of the groups studied. The GC content varied by 2.93 %, from 39.46 % for “Cyanus axillaris” to 40.61 % for Cyanus adscendens; this character is of no value for taxonomic purposes. Genome size of the studied populations is significantly higher in southern parts of the distribution area and at higher elevations. Plants with smaller genomes tend to occur in dry areas at low altitudes with high diurnal and annual temperature oscillations. The GC content of the populations studied is significantly correlated with longitude, increasing from east to west; and plants with GC-rich genomes are concentrated in the coldest areas with low minimum temperatures.

Unique Cost Dynamics Elucidate the Role of Frameshifting Errors in Promoting Translational Robustness

There is now considerable evidence supporting the view that codon usage is frequently under selection for translational accuracy. There are, however, multiple forms of inaccuracy (missense, premature termination, and frameshifting errors) and pinpointing a particular error process behind apparently adaptive mRNA anatomy is rarely straightforward. Understanding differences in the fitness costs associated with different types of translational error can help us devise critical tests that can implicate one error process to the exclusion of others. To this end, we present a model that captures distinct features of frameshifting cost and apply this to 641 prokaryotic genomes. We demonstrate that, although it is commonly assumed that the ribosome encounters an off-frame stop codon soon after the frameshift and costs of mis-elongation are therefore limited, genomes with high GC content typically incur much larger per-error costs. We go on to derive the prediction, unique to frameshifting errors, that differences in translational robustness between the 5′ and 3′ ends of genes should be less pronounced in genomes with higher GC content. This prediction we show to be correct. Surprisingly, this does not mean that GC-rich organisms necessarily carry a greater fitness burden as a consequence of accidental frameshifting. Indeed, increased per-error costs are often more than counterbalanced by lower predicted error rates owing to more diverse anticodon repertoires in GC-rich genomes. We therefore propose that selection on tRNA repertoires may operate to reduce frameshifting errors.

Identifying translation initiation sites in prokaryotes using support vector machine

Motivation: Gene identification in genomes has been a fundamental and long-standing task in bioinformatics and computational biology. Many computational methods have been developed to predict genes in prokaryote genomes by identifying translation initiation site (TIS) in transcript data. However, the pseudo-TISs at the genome level make these methods suffer from a high number of false positive predictions. In addition, most of the existing tools use an unsupervised learning framework, whose predictive accuracy may depend on the choice of specific organism. Results: In this paper, we present a supervised learning method, support vector machine (SVM), to identify translation initiation site at the genome level. The features are extracted from the sequence data by modeling the sequence segment around predicted TISs as a position specific weight matrix (PSWM). We train the parameters of our SVM through well constructed positive and negative TIS datasets. Then we apply the method to recognize translation initiation sites in E. coli, B. subtilis, and validate our method on two GC-rich bacteria genomes: Pseudomonas aeruginosa and Burkholderia pseudomallei K96243. We show that translation initiation sites can be recognized accurately at the genome level by our method, irrespective of their GC content. Furthermore, we compare our method with four existing methods and demonstrate that our method outperform these methods by obtaining better performance in all the four organisms.

Extreme diversity in noncalcifying haptophytes explains a major pigment paradox in open oceans

The current paradigm holds that cyanobacteria, which evolved oxygenic photosynthesis more than 2 billion years ago, are still the major light harvesters driving primary productivity in open oceans. Here we show that tiny unicellular eukaryotes belonging to the photosynthetic lineage of the Haptophyta are dramatically diverse and ecologically dominant in the planktonic photic realm. The use of Haptophyta-specific primers and PCR conditions adapted for GC-rich genomes circumvented biases inherent in classical genetic approaches to exploring environmental eukaryotic biodiversity and led to the discovery of hundreds of unique haptophyte taxa in 5 clone libraries from subpolar and subtropical oceanic waters. Phylogenetic analyses suggest that this diversity emerged in Paleozoic oceans, thrived and diversified in the permanently oxygenated Mesozoic Panthalassa, and currently comprises thousands of ribotypic species, belonging primarily to low-abundance and ancient lineages of the "rare biosphere." This extreme biodiversity coincides with the pervasive presence in the photic zone of the world ocean of 19'-hexanoyloxyfucoxanthin (19-Hex), an accessory photosynthetic pigment found exclusively in chloroplasts of haptophyte origin. Our new estimates of depth-integrated relative abundance of 19-Hex indicate that haptophytes dominate the chlorophyll a-normalized phytoplankton standing stock in modern oceans. Their ecologic and evolutionary success, arguably based on mixotrophy, may have significantly impacted the oceanic carbon pump. These results add to the growing evidence that the evolution of complex microbial eukaryotic cells is a critical force in the functioning of the biosphere.

The Consequences of Base Pair Composition Biases for Regulatory Network Organization in Prokaryotes

Given the dramatic variation in guanine-cytosine (GC) content observed in prokaryotes, from approximately 20% to approximately 75% GC, one wonders if these extreme biases in base pair composition affect the evolution of transcription factor-binding sites (BS). This letter shows that, along the wide range of GC content variation in bacteria, bacterial BS keep a high frequency of AT bases, roughly independently of the background (BG) base pair composition of intergenic regions. As a result, the equilibrium base pair frequencies of BS depart the most from those of BS DNA in GC-rich genomes. This not only implies a higher specificity but also a higher coding barrier for BS in GC-rich genomes. In accordance, we observe that the average percentage of divergently transcribed regions increases with the GC content of the genome, suggesting the use of a more efficient coding strategy.

Genome information management and integrated data analysis with HaloLex

HaloLex is a software system for the central management, integration, curation, and web-based visualization of genomic and other -omics data for any given microorganism. The system has been employed for the manual curation of three haloarchaeal genomes, namely Halobacterium salinarum (strain R1), Natronomonas pharaonis, and Haloquadratum walsbyi. HaloLex, in particular, enables the integrated analysis of genome-wide proteomic results with the underlying genomic data. This has proven indispensable to generate reliable gene predictions for GC-rich genomes, which, due to their characteristically low abundance of stop codons, are known to be hard targets for standard gene finders, especially concerning start codon assignment. The proteomic identification of more than 600 N-terminal peptides has greatly increased the reliability of the start codon assignment for Halobacterium salinarum. Application of homology-based methods to the published genome of Haloarcula marismortui allowed to detect 47 previously unidentified genes (a problem that is particularly serious for short protein sequences) and to correct more than 300 start codon misassignments.

Inversing the natural hydrogen bonding rule to selectively amplify GC-rich ADAR-edited RNAs

DNA complementarity is expressed by way of three hydrogen bonds for a G:C base pair and two for A:T. As a result, careful control of the denaturation temperature of PCR allows selective amplification of AT-rich alleles. Yet for the same reason, the converse is not possible, selective amplification of GC-rich alleles. Inosine (I) hydrogen bonds to cytosine by two hydrogen bonds while diaminopurine (D) forms three hydrogen bonds with thymine. By substituting dATP by dDTP and dGTP by dITP in a PCR reaction, DNA is obtained in which the natural hydrogen bonding rule is inversed. When PCR is performed at limiting denaturation temperatures, it is possible to recover GC-rich viral genomes and inverted Alu elements embedded in cellular mRNAs resulting from editing by dsRNA dependent host cell adenosine deaminases. The editing of Alu elements in cellular mRNAs was strongly enhanced by type I interferon induction indicating a novel link mRNA metabolism and innate immunity.

Investigations of Oligonucleotide Usage Variance Within and Between Prokaryotes

Oligonucleotide usage in archaeal and bacterial genomes can be linked to a number of properties, including codon usage (trinucleotides), DNA base-stacking energy (dinucleotides), and DNA structural conformation (di- to tetranucleotides). We wanted to assess the statistical information potential of different DNA ‘word-sizes’ and explore how oligonucleotide frequencies differ in coding and non-coding regions. In addition, we used oligonucleotide frequencies to investigate DNA composition and how DNA sequence patterns change within and between prokaryotic organisms. Among the results found was that prokaryotic chromosomes can be described by hexanucleotide frequencies, suggesting that prokaryotic DNA is predominantly short range correlated, i.e., information in prokaryotic genomes is encoded in short oligonucleotides. Oligonucleotide usage varied more within AT-rich and host-associated genomes than in GC-rich and free-living genomes, and this variation was mainly located in non-coding regions. Bias (selectional pressure) in tetranucleotide usage correlated with GC content, and coding regions were more biased than non-coding regions. Non-coding regions were also found to be approximately 5.5% more AT-rich than coding regions, on average, in the 402 chromosomes examined. Pronounced DNA compositional differences were found both within and between AT-rich and GC-rich genomes. GC-rich genomes were more similar and biased in terms of tetranucleotide usage in non-coding regions than AT-rich genomes. The differences found between AT-rich and GC-rich genomes may possibly be attributed to lifestyle, since tetranucleotide usage within host-associated bacteria was, on average, more dissimilar and less biased than free-living archaea and bacteria.

Large-Scale Identification of N-Terminal Peptides in the Halophilic Archaea Halobacterium salinarum and Natronomonas pharaonis

Characterization of protein N-terminal peptides supports the quality assessment of data derived from genomic sequences (e.g., the correct assignment of start codons) and hints to in vivo N-terminal modifications such as N-terminal acetylation and removal of the initiator methionine. The current work represents the first large-scale identification of N-terminal peptides from prokaryotes, of the two halophilic euryarchaeota Halobacterium salinarum and Natronomonas pharaonis. Two methods were used that specifically allow the characterization of protein N-terminal peptides: combined fractional diagonal chromatography (COFRADIC) and strong cation exchange chromatography (SCX), both known to enrich for N-terminally blocked peptides. In addition to these specific methods, N-terminal peptide identifications were extracted from our previous genome-wide proteomic data. Combining all data, 606 N-terminal peptides from Hbt. salinarum and 328 from Nmn. pharaonis were reliably identified. These results constitute the largest available dataset holding identified and characterized protein N-termini for prokaryotes (archaea and bacteria). They allowed the validation/improvement of start codon assignments as automatic gene finders tend to misassign start codons for GC-rich genomes. In addition, the dataset allowed unravelling N-terminal protein maturation in archaea, showing that 60% of the proteins undergo methionine cleavage and that-in contrast to current knowledge-Nalpha-acetylation is common in the archaeal domain of life with 13-18% of the proteins being Nalpha-acetylated. The protein sets described in this paper are available by FTP and might be used as reference sets to test the performance of new gene finders.

Large-scale prokaryotic gene prediction and comparison to genome annotation

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. 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.

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 1stcodon 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.

Comparative Analysis of the Base Composition and Codon Usages in Fourteen Mycobacteriophage Genomes

To study the possible codon usage and base composition variation in the bacteriophages, fourteen mycobacteriophages were used as a model system here and both the parameters in all these phages and their plating bacteria, M. smegmatis had been determined and compared. As all the organisms are GC-rich, the GC contents at third codon positions were found in fact higher than the second codon positions as well as the first + second codon positions in all the organisms indicating that directional mutational pressure is strongly operative at the synonymous third codon positions. Nc plot indicates that codon usage variation in all these organisms are governed by the forces other than compositional constraints. Correspondence analysis suggests that: (i) there are codon usage variation among the genes and genomes of the fourteen mycobacteriophages and M. smegmatis, i.e, codon usage patterns in the mycobacteriophages is phage-specific but not the M. smegmatis-specific; (ii) synonymous codon usage patterns of Barnyard, Che8, Che9d, and Omega are more similar than the rest mycobacteriophages and M. smegmatis; (iii) codon usage bias in the mycobacteriophages are mainly determined by mutational pressure; and (iv) the genes of comparatively GC rich genomes are more biased than the GC poor genomes. Translational selection in determining the codon usage variation in highly expressed genes can be invoked from the predominant occurrences of C ending codons in the highly expressed genes. Cluster analysis based on codon usage data also shows that there are two distinct branches for the fourteen mycobacteriophages and there is codon usage variation even among the phages of each branch.

GeneLook: a novel ab initio gene identification system suitable for automated annotation of prokaryotic sequences

With the rapid increases in the amounts of sequence data for prokaryotic genomes, it has become important to develop systems for automated and accurate genome annotation. We present herein a novel ab initio gene identification system, GeneLook, that predicts protein-coding open reading frames (ORFs) with high sensitivity and specificity with no prior knowledge of the sequence composition. The system predicts protein-coding ORFs in two stages, seed ORF selection and main prediction. In the selection of reliable seed ORFs containing at least 200 codons, GeneLook predicts translation start sites and operon structures through searches for ribosome-binding sites and a novel operon prediction algorithm. The codon and nucleotide frequencies of seed ORFs are then used to determine values for two new coding-potential parameters for identification of protein-coding ORFs of at least 34 codons and for another parameter that improves the prediction accuracy for GC-rich genomes. In the main prediction, GeneLook uses these parameters to identify the most likely genes of a given minimal length. We assessed the performance of GeneLook with two indices, sensitivity and specificity that are defined as true positives (TP)/(TP+false negatives) and TP/(TP+false positives), respectively. This system predicted protein-coding ORFs for Escherichia coli and Bacillus subtilis with sensitivities of 96.5% and 96.2%, respectively, and specificities of 96.9% and 96.1%, respectively. The system also identified 94.1% of annotated genes of the Pseudomonas aeruginosa genome, which is GC-rich, with high specificity (97.2%). Furthermore, GeneLook identified protein-coding ORFs with high accuracy from a wide variety of prokaryotic genomes.

FrameD: A flexible program for quality check and gene prediction in prokaryotic genomes and noisy matured eukaryotic sequences.

We describe FrameD, a program that predicts coding regions in prokaryotic and matured eukaryotic sequences. Initially targeted at gene prediction in bacterial GC rich genomes, the gene model used in FrameD also allows to predict genes in the presence of frameshifts and partially undetermined sequences which makes it also very suitable for gene prediction and frameshift correction in unfinished sequences such as EST and EST cluster sequences. Like recent eukaryotic gene prediction programs, FrameD also includes the ability to take into account protein similarity information both in its prediction and its graphical output. Its performances are evaluated on different bacterial genomes. The web site (http://genopole.toulouse.inra.fr/bioinfo/FrameD/FD) allows direct prediction, sequence correction and translation and the ability to learn new models for new organisms.