Impaired removal of H3K4 methylation affects cell fate determination and gene transcription.

Methylation of histone 3 lysine 4 (H3K4) is largely associated with promoters and enhancers of actively transcribed genes and is finely regulated during development by the action of histone methyltransferases and demethylases. H3K4me3 demethylases of the KDM5 family have been previously implicated in development, but how the regulation of H3K4me3 level controls developmental processes is not fully established. Here, we show that the H3K4 demethylase RBR-2, the unique member of the KDM5 family in C. elegans, acts cell-autonomously and in a catalytic-dependent manner to control vulva precursor cells fate acquisition, by promoting the LIN-12/Notch pathway. Using genome-wide approaches, we show that RBR-2 reduces the H3K4me3 level at transcription start sites (TSSs) and in regions upstream of the TSSs, and acts both as a transcription repressor and activator. Analysis of the lin-11 genetic locus, a direct RBR-2 target gene required for vulva precursor cell fate acquisition, shows that RBR-2 controls the epigenetic signature of the lin-11 vulva-specific enhancer and lin-11 expression, providing in vivo evidence that RBR-2 can positively regulate transcription and cell fate acquisition by controlling enhancer activity.

Loss of PRDM11 promotes MYC-driven lymphomagenesis.

The PR-domain (PRDM) family of genes encodes transcriptional regulators, several of which are deregulated in cancer. By using a functional screening approach, we sought to identify novel tumor suppressors among the PRDMs. Here we demonstrate oncogenic collaboration between depletion of the previously uncharacterized PR-domain family member Prdm11 and overexpression of MYC. Overexpression of PRDM11 inhibits proliferation and induces apoptosis. Prdm11 knockout mice are viable, and loss of Prdm11 accelerates MYC-driven lymphomagenesis in the Eµ-Myc mouse model. Moreover, we show that patients with PRDM11-deficient diffuse large B-cell lymphomas (DLBCLs) have poorer overall survival and belong to the nongerminal center B-cell-like subtype. Mechanistically, genome-wide mapping of PRDM11 binding sites coupled with transcriptome sequencing in human DLBCL cells evidenced that PRDM11 associates with transcriptional start sites of target genes and regulates important oncogenes such as FOS and JUN. Hence, we characterize PRDM11 as a putative novel tumor suppressor that controls the expression of key oncogenes, and we add new mechanistic insight into B-cell lymphomagenesis.

H3K23me2 is a new heterochromatic mark in Caenorhabditis elegans.

Genome-wide analyses in Caenorhabditis elegans show that post-translational modifications (PTMs) of histones are evolutionary conserved and distributed along functionally distinct genomic domains. However, a global profile of PTMs and their co-occurrence on the same histone tail has not been described in this organism. We used mass spectrometry based middle-down proteomics to analyze histone H3 N-terminal tails from C. elegans embryos for the presence, the relative abundance and the potential cross-talk of co-existing PTMs. This analysis highlighted that the lysine 23 of histone H3 (H3K23) is extensively modified by methylation and that tri-methylated H3K9 (H3K9me3) is exclusively detected on histone tails with di-methylated H3K23 (H3K23me2). Chromatin immunoprecipitation approaches revealed a positive correlation between H3K23me2 and repressive marks. By immunofluorescence analyses, H3K23me2 appears differentially regulated in germ and somatic cells, in part by the action of the histone demethylase JMJD-1.2. H3K23me2 is enriched in heterochromatic regions, localizing in H3K9me3 and heterochromatin protein like-1 (HPL-1)-positive foci. Biochemical analyses indicated that HPL-1 binds to H3K23me2 and interacts with a conserved CoREST repressive complex. Thus, our study suggests that H3K23me2 defines repressive domains and contributes to organizing the genome in distinct heterochromatic regions during embryogenesis.

Genomic variation in Salmonella enterica core genes for epidemiological typing.

BACKGROUND: Technological advances in high throughput genome sequencing are making whole genome sequencing (WGS) available as a routine tool for bacterial typing. Standardized procedures for identification of relevant genes and of variation are needed to enable comparison between studies and over time. The core genes--the genes that are conserved in all (or most) members of a genus or species--are potentially good candidates for investigating genomic variation in phylogeny and epidemiology. RESULTS: We identify a set of 2,882 core genes clusters based on 73 publicly available Salmonella enterica genomes and evaluate their value as typing targets, comparing whole genome typing and traditional methods such as 16S and MLST. A consensus tree based on variation of core genes gives much better resolution than 16S and MLST; the pan-genome family tree is similar to the consensus tree, but with higher confidence. The core genes can be divided into two categories: a few highly variable genes and a larger set of conserved core genes, with low variance. For the most variable core genes, the variance in amino acid sequences is higher than for the corresponding nucleotide sequences, suggesting that there is a positive selection towards mutations leading to amino acid changes. CONCLUSIONS: Genomic variation within the core genome is useful for investigating molecular evolution and providing candidate genes for bacterial genome typing. Identification of genes with different degrees of variation is important especially in trend analysis.

Estimating variation within the genes and inferring the phylogeny of 186 sequenced diverse Escherichia coli genomes.

BACKGROUND: Escherichia coli exists in commensal and pathogenic forms. By measuring the variation of individual genes across more than a hundred sequenced genomes, gene variation can be studied in detail, including the number of mutations found for any given gene. This knowledge will be useful for creating better phylogenies, for determination of molecular clocks and for improved typing techniques. RESULTS: We find 3,051 gene clusters/families present in at least 95% of the genomes and 1,702 gene clusters present in 100% of the genomes. The former 'soft core' of about 3,000 gene families is perhaps more biologically relevant, especially considering that many of these genome sequences are draft quality. The E. coli pan-genome for this set of isolates contains 16,373 gene clusters.A core-gene tree, based on alignment and a pan-genome tree based on gene presence/absence, maps the relatedness of the 186 sequenced E. coli genomes. The core-gene tree displays high confidence and divides the E. coli strains into the observed MLST type clades and also separates defined phylotypes. CONCLUSION: The results of comparing a large and diverse E. coli dataset support the theory that reliable and good resolution phylogenies can be inferred from the core-genome. The results further suggest that the resolution at the isolate level may, subsequently be improved by targeting more variable genes. The use of whole genome sequencing will make it possible to eliminate, or at least reduce, the need for several typing steps used in traditional epidemiology.

snpTree--a web-server to identify and construct SNP trees from whole genome sequence data.

BACKGROUND: The advances and decreasing economical cost of whole genome sequencing (WGS), will soon make this technology available for routine infectious disease epidemiology. In epidemiological studies, outbreak isolates have very little diversity and require extensive genomic analysis to differentiate and classify isolates. One of the successfully and broadly used methods is analysis of single nucletide polymorphisms (SNPs). Currently, there are different tools and methods to identify SNPs including various options and cut-off values. Furthermore, all current methods require bioinformatic skills. Thus, we lack a standard and simple automatic tool to determine SNPs and construct phylogenetic tree from WGS data. RESULTS: Here we introduce snpTree, a server for online-automatic SNPs analysis. This tool is composed of different SNPs analysis suites, perl and python scripts. snpTree can identify SNPs and construct phylogenetic trees from WGS as well as from assembled genomes or contigs. WGS data in fastq format are aligned to reference genomes by BWA while contigs in fasta format are processed by Nucmer. SNPs are concatenated based on position on reference genome and a tree is constructed from concatenated SNPs using FastTree and a perl script. The online server was implemented by HTML, Java and python script.The server was evaluated using four published bacterial WGS data sets (V. cholerae, S. aureus CC398, S. Typhimurium and M. tuberculosis). The evaluation results for the first three cases was consistent and concordant for both raw reads and assembled genomes. In the latter case the original publication involved extensive filtering of SNPs, which could not be repeated using snpTree. CONCLUSIONS: The snpTree server is an easy to use option for rapid standardised and automatic SNP analysis in epidemiological studies also for users with limited bioinformatic experience. The web server is freely accessible at http://www.cbs.dtu.dk/services/snpTree-1.0/.

Multilocus sequence typing of total-genome-sequenced bacteria.

Accurate strain identification is essential for anyone working with bacteria. For many species, multilocus sequence typing (MLST) is considered the "gold standard" of typing, but it is traditionally performed in an expensive and time-consuming manner. As the costs of whole-genome sequencing (WGS) continue to decline, it becomes increasingly available to scientists and routine diagnostic laboratories. Currently, the cost is below that of traditional MLST. The new challenges will be how to extract the relevant information from the large amount of data so as to allow for comparison over time and between laboratories. Ideally, this information should also allow for comparison to historical data. We developed a Web-based method for MLST of 66 bacterial species based on WGS data. As input, the method uses short sequence reads from four sequencing platforms or preassembled genomes. Updates from the MLST databases are downloaded monthly, and the best-matching MLST alleles of the specified MLST scheme are found using a BLAST-based ranking method. The sequence type is then determined by the combination of alleles identified. The method was tested on preassembled genomes from 336 isolates covering 56 MLST schemes, on short sequence reads from 387 isolates covering 10 schemes, and on a small test set of short sequence reads from 29 isolates for which the sequence type had been determined by traditional methods. The method presented here enables investigators to determine the sequence types of their isolates on the basis of WGS data. This method is publicly available at www.cbs.dtu.dk/services/MLST.

The Salmonella enterica pan-genome.

Salmonella enterica is divided into four subspecies containing a large number of different serovars, several of which are important zoonotic pathogens and some show a high degree of host specificity or host preference. We compare 45 sequenced S. enterica genomes that are publicly available (22 complete and 23 draft genome sequences). Of these, 35 were found to be of sufficiently good quality to allow a detailed analysis, along with two Escherichia coli strains (K-12 substr. DH10B and the avian pathogenic E. coli (APEC O1) strain). All genomes were subjected to standardized gene finding, and the core and pan-genome of Salmonella were estimated to be around 2,800 and 10,000 gene families, respectively. The constructed pan-genomic dendrograms suggest that gene content is often, but not uniformly correlated to serotype. Any given Salmonella strain has a large stable core, whilst there is an abundance of accessory genes, including the Salmonella pathogenicity islands (SPIs), transposable elements, phages, and plasmid DNA. We visualize conservation in the genomes in relation to chromosomal location and DNA structural features and find that variation in gene content is localized in a selection of variable genomic regions or islands. These include the SPIs but also encompass phage insertion sites and transposable elements. The islands were typically well conserved in several, but not all, isolates--a difference which may have implications in, e.g., host specificity.

Proteome analysis of pod and seed development in the model legume Lotus japonicus.

Legume pods serve important functions during seed development and are themselves sources of food and feed. Compared to seeds, the metabolism and development of pods are not well-defined. The present characterization of pods from the model legume Lotus japonicus, together with the detailed analyses of the pod and seed proteomes in five developmental stages, paves the way for comparative pathway analysis and provides new metabolic information. Proteins were analyzed by two-dimensional gel electrophoresis and tandem-mass spectrometry. These analyses lead to the identification of 604 pod proteins and 965 seed proteins, including 263 proteins distinguishing the pod. The complete data set is publicly available at http://www.cbs.dtu.dk/cgi-bin/lotus/db.cgi , where spots in a reference map are linked to experimental data, such as matched peptides, quantification values, and gene accessions. Identified pod proteins represented enzymes from 85 different metabolic pathways, including storage globulins and a late embryogenesis abundant protein. In contrast to seed maturation, pod maturation was associated with decreasing total protein content, especially proteins involved in protein biosynthesis and photosynthesis. Proteins detected only in pods included three enzymes participating in the urea cycle and four in nitrogen and amino group metabolism, highlighting the importance of nitrogen metabolism during pod development. Additionally, five legume seed proteins previously unassigned in the glutamate metabolism pathway were identified.

Comparative profiling of the transcriptional response to iron restriction in six serotypes of Actinobacillus pleuropneumoniae with different virulence potential.

BACKGROUND: Comparative analysis of gene expression among serotypes within a species can provide valuable information on important differences between related genomes. For the pig lung pathogen Actinobacillus pleuropneumoniae, 15 serotypes with a considerable variation in virulence potential and immunogenicity have been identified. This serotypic diversity can only partly be explained by amount of capsule and differences in the RTX toxin genes in their genomes. Iron acquisition in vivo is an important bacterial function and in pathogenic bacteria, iron-limitation is often a signal for the induction of virulence genes. We used a pan-genomic microarray to study the transcriptional response to iron restriction in vitro in six serotypes of A. pleuropneumoniae (1, 2, 3, 5b, 6, and 7), representing at least two levels of virulence. RESULTS: In total, 45 genes were significantly (p < 0.0001) up-regulated and 67 genes significantly down-regulated in response to iron limitation. Not previously observed in A. pleuropneumoniae was the up-regulation of a putative cirA-like siderophore in all six serotypes. Three genes, recently described in A. pleuropneumoniae as possibly coding for haemoglobin-haptoglobin binding proteins, displayed significant serotype related up-regulation to iron limitation. For all three genes, the expression appeared at its lowest in serotype 3, which is generally considered one of the least virulent serotypes of A. pleuropneumoniae. The three genes share homology with the hmbR haemoglobin receptor of Neisseria meningitidis, a possible virulence factor which contributes to bacterial survival in rats. CONCLUSIONS: By comparative analysis of gene expression among 6 different serotypes of A. pleuropneumoniae we identified a common set of presumably essential core genes, involved in iron regulation. The results support and expand previous observations concerning the identification of new potential iron acquisition systems in A. pleuropneumoniae, showing that this bacterium has evolved several strategies for scavenging the limited iron resources of the host. The combined effect of iron-depletion and serotype proved to be modest, indicating that serotypes of both moderate and high virulence at least in vitro are reacting almost identical to iron restriction. One notable exception, however, is the haemoglobin-haptoglobin binding protein cluster which merits further investigation.

A virulent parent with probiotic progeny: comparative genomics of Escherichia coli strains CFT073, Nissle 1917 and ABU 83972.

Escherichia coli is a highly versatile species encompassing a diverse spectrum of strains, i.e. from highly virulent isolates causing serious infectious diseases to commensals and probiotic strains. Although much is known about bacterial pathogenicity in E. coli, the understanding of which genetic determinants differentiates a virulent from an avirulent strain still remains limited. In this study we designed a new comparative genomic hybridization microarray based on 31 sequenced E. coli strains and used it to compare two E. coli strains used as prophylactic agents (i.e. Nissle 1917 and 83972) with the highly virulent uropathogen CFT073. Only relatively minor genetic variations were found between the isolates, suggesting that the three strains may have originated from the same virulent ancestral parent. Interestingly, Nissle 1917 (a gut commensal strain) was more similar to CFT073 with respect to genotype and phenotype than 83972 (an asymptomatic bacteriuria strain). The study indicates that genetic variations (e.g. mutations) and expression differences, rather than genomic content per se, contribute to the divergence in disease-causing ability between these strains. This has implications for the use of virulence factors in epidemiological research, and emphasizes the need for more comparative genomic studies of closely related strains to compare their virulence potential.

A pathway-specific microarray analysis highlights the complex and co-ordinated transcriptional networks of the developing grain of field-grown barley.

The aim of the study was to describe the molecular and biochemical interactions associated with amino acid biosynthesis and storage protein accumulation in the developing grains of field-grown barley. Our strategy was to analyse the transcription of genes associated with the biosynthesis of storage products during the development of field-grown barley grains using a grain-specific microarray assembled in our laboratory. To identify co-regulated genes, a distance matrix was constructed which enabled the identification of three clusters corresponding to early, middle, and late grain development. The gene expression pattern associated with the clusters was investigated using pathway-specific analysis with specific reference to the temporal expression levels of a range of genes involved mainly in the photosynthesis process, amino acid and storage protein metabolism. It is concluded that the grain-specific microarray is a reliable and cost-effective tool for monitoring temporal changes in the transcriptome of the major metabolic pathways in the barley grain. Moreover, it was sensitive enough to monitor differences in the gene expression profiles of different homologues from the storage protein families. The study described here should provide a strong complement to existing knowledge assisting further understanding of grain development and thereby provide a foundation for plant breeding towards storage proteins with improved nutritional quality.

Antisense-mediated suppression of C-hordein biosynthesis in the barley grain results in correlated changes in the transcriptome, protein profile, and amino acid composition.

Antisense- or RNAi-mediated suppression of the biosynthesis of nutritionally inferior storage proteins is a promising strategy for improving the amino acid profile of seeds. However, the potential pleiotropic effects of this on interconnected pathways and the agronomic quality traits need to be addressed. In the current study, a transcriptomic analysis of an antisense C-hordein line of barley was performed, using a grain-specific cDNA array. The C-hordein antisense line is characterized by marked changes in storage protein and amino acid profiles, while the seed weight is within the normal range and no external morphological irregularities were observed. The results of the transcriptome analysis showed excellent correlation with data on changes in the relative proportions of storage proteins and amino acid composition. The antisense line had a lower C-hordein level and down-regulated transcript encoding C-hordein. The production of the S-rich B/gamma- and D-hordeins was increased and significantly higher steady-state expression levels of the corresponding genes were observed. The increased synthesis of S-rich hordeins appeared to increase the demand for sulphur and the S-rich amino acids (cysteine and methionine), resulting in an up-regulation of key genes in the appropriate biosynthetic pathways. This study demonstrated the utility of the grain-specific cDNA microarray analysis to detect perturbations induced by antisense suppression of plant processes.

GenePublisher: Automated analysis of DNA microarray data.

GenePublisher, a system for automatic analysis of data from DNA microarray experiments, has been implemented with a web interface at http://www.cbs.dtu.dk/services/GenePublisher. Raw data are uploaded to the server together with a specification of the data. The server performs normalization, statistical analysis and visualization of the data. The results are run against databases of signal transduction pathways, metabolic pathways and promoter sequences in order to extract more information. The results of the entire analysis are summarized in report form and returned to the user.

Genomics of an emerging clone of Salmonella serovar Typhimurium ST313 from Nigeria and the Democratic Republic of Congo.

INTRODUCTION: Salmonella enterica serovar Typhimurium ST313 is an invasive and phylogenetically distinct lineage present in sub-Saharan Africa. We report the presence of S. Typhimurium ST313 from patients in the Democratic Republic of Congo and Nigeria. METHODOLOGY: Eighteen S. Typhimurium ST313 isolates were characterized by antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), and multilocus sequence typing (MLST). Additionally, six of the isolates were characterized by whole genome sequence typing (WGST). The presence of a putative virulence determinant was examined in 177 Salmonella isolates belonging to 57 different serovars. RESULTS: All S. Typhimurium ST313 isolates harbored resistant genes encoded by blaTEM1b, catA1, strA/B, sul1, and dfrA1. Additionally, aac(6')1aa gene was detected. Phylogenetic analyses revealed close genetic relationships among Congolese and Nigerian isolates from both blood and stool. Comparative genomic analyses identified a putative virulence fragment (ST313-TD) unique to S. Typhimurium ST313 and S. Dublin. CONCLUSION: We showed in a limited number of isolates that S. Typhimurium ST313 is a prevalent sequence-type causing gastrointestinal diseases and septicemia in patients from Nigeria and DRC. We found three distinct phylogenetic clusters based on the origin of isolation suggesting some spatial evolution. Comparative genomics showed an interesting putative virulence fragment (ST313-TD) unique to S. Typhimurium ST313 and invasive S. Dublin.

Transcriptional portrait of Actinobacillus pleuropneumoniae during acute disease--potential strategies for survival and persistence in the host.

BACKGROUND: Gene expression profiles of bacteria in their natural hosts can provide novel insight into the host-pathogen interactions and molecular determinants of bacterial infections. In the present study, the transcriptional profile of the porcine lung pathogen Actinobacillus pleuropneumoniae was monitored during the acute phase of infection in its natural host. METHODOLOGY/PRINCIPAL FINDINGS: Bacterial expression profiles of A. pleuropneumoniae isolated from lung lesions of 25 infected pigs were compared in samples taken 6, 12, 24 and 48 hours post experimental challenge. Within 6 hours, focal, fibrino hemorrhagic lesions could be observed in the pig lungs, indicating that A. pleuropneumoniae had managed to establish itself successfully in the host. We identified 237 differentially regulated genes likely to encode functions required by the bacteria for colonization and survival in the host. This group was dominated by genes involved in various aspects of energy metabolism, especially anaerobic respiration and carbohydrate metabolism. Remodeling of the bacterial envelope and modifications of posttranslational processing of proteins also appeared to be of importance during early infection. The results suggested that A. pleuropneumoniae is using various strategies to increase its fitness, such as applying Na+ pumps as an alternative way of gaining energy. Furthermore, the transcriptional data provided potential clues as to how A. pleuropneumoniae is able to circumvent host immune factors and survive within the hostile environment of host macrophages. This persistence within macrophages may be related to urease activity, mobilization of various stress responses and active evasion of the host defenses by cell surface sialylation. CONCLUSIONS/SIGNIFICANCE: The data presented here highlight the importance of metabolic adjustments to host conditions as virulence factors of infecting microorganisms and help to provide insight into the mechanisms behind the efficient colonization and persistence of A. pleuropneumoniae during acute disease.

Genomic characterization of Campylobacter jejuni strain M1.

Campylobacter jejuni strain M1 (laboratory designation 99/308) is a rarely documented case of direct transmission of C. jejuni from chicken to a person, resulting in enteritis. We have sequenced the genome of C. jejuni strain M1, and compared this to 12 other C. jejuni sequenced genomes currently publicly available. Compared to these, M1 is closest to strain 81116. Based on the 13 genome sequences, we have identified the C. jejuni pan-genome, as well as the core genome, the auxiliary genes, and genes unique between strains M1 and 81116. The pan-genome contains 2,427 gene families, whilst the core genome comprised 1,295 gene families, or about two-thirds of the gene content of the average of the sequenced C. jejuni genomes. Various comparison and visualization tools were applied to the 13 C. jejuni genome sequences, including a species pan- and core genome plot, a BLAST Matrix and a BLAST Atlas. Trees based on 16S rRNA sequences and on the total gene families in each genome are presented. The findings are discussed in the background of the proven virulence potential of M1.

The proteome of seed development in the model legume Lotus japonicus.

We have characterized the development of seeds in the model legume Lotus japonicus. Like soybean (Glycine max) and pea (Pisum sativum), Lotus develops straight seed pods and each pod contains approximately 20 seeds that reach maturity within 40 days. Histological sections show the characteristic three developmental phases of legume seeds and the presence of embryo, endosperm, and seed coat in desiccated seeds. Furthermore, protein, oil, starch, phytic acid, and ash contents were determined, and this indicates that the composition of mature Lotus seed is more similar to soybean than to pea. In a first attempt to determine the seed proteome, both a two-dimensional polyacrylamide gel electrophoresis approach and a gel-based liquid chromatography-mass spectrometry approach were used. Globulins were analyzed by two-dimensional polyacrylamide gel electrophoresis, and five legumins, LLP1 to LLP5, and two convicilins, LCP1 and LCP2, were identified by matrix-assisted laser desorption ionization quadrupole/time-of-flight mass spectrometry. For two distinct developmental phases, seed filling and desiccation, a gel-based liquid chromatography-mass spectrometry approach was used, and 665 and 181 unique proteins corresponding to gene accession numbers were identified for the two phases, respectively. All of the proteome data, including the experimental data and mass spectrometry spectra peaks, were collected in a database that is available to the scientific community via a Web interface (http://www.cbs.dtu.dk/cgi-bin/lotus/db.cgi). This database establishes the basis for relating physiology, biochemistry, and regulation of seed development in Lotus. Together with a new Web interface (http://bioinfoserver.rsbs.anu.edu.au/utils/PathExpress4legumes/) collecting all protein identifications for Lotus, Medicago, and soybean seed proteomes, this database is a valuable resource for comparative seed proteomics and pathway analysis within and beyond the legume family.

Motif decomposition of the phosphotyrosine proteome reveals a new N-terminal binding motif for SHIP2.

Advances in mass spectrometry-based proteomics have yielded a substantial mapping of the tyrosine phosphoproteome and thus provided an important step toward a systematic analysis of intracellular signaling networks in higher eukaryotes. In this study we decomposed an uncharacterized proteomics data set of 481 unique phosphotyrosine (Tyr(P)) peptides by sequence similarity to known ligands of the Src homology 2 (SH2) and the phosphotyrosine binding (PTB) domains. From 20 clusters we extracted 16 known and four new interaction motifs. Using quantitative mass spectrometry we pulled down Tyr(P)-specific binding partners for peptides corresponding to the extracted motifs. We confirmed numerous previously known interaction motifs and found 15 new interactions mediated by phosphosites not previously known to bind SH2 or PTB. Remarkably, a novel hydrophobic N-terminal motif ((L/V/I)(L/V/I)pY) was identified and validated as a binding motif for the SH2 domain-containing inositol phosphatase SHIP2. Our decomposition of the in vivo Tyr(P) proteome furthermore suggests that two-thirds of the Tyr(P) sites mediate interaction, whereas the remaining third govern processes such as enzyme activation and nucleic acid binding.

An environmental signature for 323 microbial genomes based on codon adaptation indices.

BACKGROUND: Codon adaptation indices (CAIs) represent an evolutionary strategy to modulate gene expression and have widely been used to predict potentially highly expressed genes within microbial genomes. Here, we evaluate and compare two very different methods for estimating CAI values, one corresponding to translational codon usage bias and the second obtained mathematically by searching for the most dominant codon bias. RESULTS: The level of correlation between these two CAI methods is a simple and intuitive measure of the degree of translational bias in an organism, and from this we confirm that fast replicating bacteria are more likely to have a dominant translational codon usage bias than are slow replicating bacteria, and that this translational codon usage bias may be used for prediction of highly expressed genes. By analyzing more than 300 bacterial genomes, as well as five fungal genomes, we show that codon usage preference provides an environmental signature by which it is possible to group bacteria according to their lifestyle, for instance soil bacteria and soil symbionts, spore formers, enteric bacteria, aquatic bacteria, and intercellular and extracellular pathogens. CONCLUSION: The results and the approach described here may be used to acquire new knowledge regarding species lifestyle and to elucidate relationships between organisms that are far apart evolutionarily.