Set Of Orthologous Genes Research Articles

Computational identification of functional introns: high positional conservation of introns that harbor RNA genes

An appreciable fraction of introns is thought to have some function, but there is no obvious way to predict which specific intron is likely to be functional. We hypothesize that functional introns experience a different selection regime than non-functional ones and will therefore show distinct evolutionary histories. In particular, we expect functional introns to be more resistant to loss, and that this would be reflected in high conservation of their position with respect to the coding sequence. To test this hypothesis, we focused on introns whose function comes about from microRNAs and snoRNAs that are embedded within their sequence. We built a data set of orthologous genes across 28 eukaryotic species, reconstructed the evolutionary histories of their introns and compared functional introns with the rest of the introns. We found that, indeed, the position of microRNA- and snoRNA-bearing introns is significantly more conserved. In addition, we found that both families of RNA genes settled within introns early during metazoan evolution. We identified several easily computable intronic properties that can be used to detect functional introns in general, thereby suggesting a new strategy to pinpoint non-coding cellular functions.

Comparative genome analysis and gene finding in Candida species using CGOB.

The Candida Gene Order Browser (CGOB) was developed as a tool to visualize and analyze synteny relationships in multiple Candida species, and to provide an accurate, manually curated set of orthologous Candida genes for evolutionary analyses. Here, we describe major improvements to CGOB. The underlying structure of the database has been changed significantly. Genomic features are now based directly on genome annotations rather than on protein sequences, which allows non-protein features such as centromere locations in Candida albicans and tRNA genes in all species to be included. The data set has been expanded to 13 species, including genomes of pathogens (C. albicans, C. parapsilosis, C. tropicalis, and C. orthopsilosis), and those of xylose-degrading species with important biotechnological applications (C. tenuis, Scheffersomyces stipitis, and Spathaspora passalidarum). Updated annotations of C. parapsilosis, C. dubliniensis, and Debaryomyces hansenii have been incorporated. We discovered more than 1,500 previously unannotated genes among the 13 genomes, ranging in size from 29 to 3,850 amino acids. Poorly conserved and rapidly evolving genes were also identified. Re-analysis of the mating type loci of the xylose degraders suggests that C. tenuis is heterothallic, whereas both Spa. passalidarum and S. stipitis are homothallic. As well as hosting the browser, the CGOB website (http://cgob.ucd.ie) gives direct access to all the underlying genome annotations, sequences, and curated orthology data.

Tracking recent adaptive evolution in microbial species using TimeZone

An important goal of the analysis of sequenced genomes of microbial pathogens is to improve the therapy of infectious diseases. In this context, a major challenge is to detect genomic-level evolutionary changes that increase microbial virulence. TimeZone, a genome analysis software package, is designed to detect footprints of positive selection for functionally adaptive point mutations. The uniqueness of TimeZone lies in its ability to predict recent adaptive mutations that are overlooked by conventional microevolutionary tools. This protocol describes the use of TimeZone to analyze adaptive footprints in either individual genes or in sets of genomes. Three major workflows are described: (i) extraction of orthologous gene sets from multiple genomes; (ii) alignment and phylogenetic analysis of genes; and (iii) identification of candidate genes under positive selection for point mutations, taking into account the effect of recombination events. This software package can be downloaded free from http://sourceforge.net/projects/timezone1/. In the case, for example, of the analysis of 14 Escherichia coli genomes, the protocol described here can be completed in ∼32 h.

Universal Pacemaker of Genome Evolution

A fundamental observation of comparative genomics is that the distribution of evolution rates across the complete sets of orthologous genes in pairs of related genomes remains virtually unchanged throughout the evolution of life, from bacteria to mammals. The most straightforward explanation for the conservation of this distribution appears to be that the relative evolution rates of all genes remain nearly constant, or in other words, that evolutionary rates of different genes are strongly correlated within each evolving genome. This correlation could be explained by a model that we denoted Universal PaceMaker (UPM) of genome evolution. The UPM model posits that the rate of evolution changes synchronously across genome-wide sets of genes in all evolving lineages. Alternatively, however, the correlation between the evolutionary rates of genes could be a simple consequence of molecular clock (MC). We sought to differentiate between the MC and UPM models by fitting thousands of phylogenetic trees for bacterial and archaeal genes to supertrees that reflect the dominant trend of vertical descent in the evolution of archaea and bacteria and that were constrained according to the two models. The goodness of fit for the UPM model was better than the fit for the MC model, with overwhelming statistical significance, although similarly to the MC, the UPM is strongly overdispersed. Thus, the results of this analysis reveal a universal, genome-wide pacemaker of evolution that could have been in operation throughout the history of life.

Comparative Analysis of Syntenic Genes in Grass Genomes Reveals Accelerated Rates of Gene Structure and Coding Sequence Evolution in Polyploid Wheat

Cycles of whole-genome duplication (WGD) and diploidization are hallmarks of eukaryotic genome evolution and speciation. Polyploid wheat (Triticum aestivum) has had a massive increase in genome size largely due to recent WGDs. How these processes may impact the dynamics of gene evolution was studied by comparing the patterns of gene structure changes, alternative splicing (AS), and codon substitution rates among wheat and model grass genomes. In orthologous gene sets, significantly more acquired and lost exonic sequences were detected in wheat than in model grasses. In wheat, 35% of these gene structure rearrangements resulted in frame-shift mutations and premature termination codons. An increased codon mutation rate in the wheat lineage compared with Brachypodium distachyon was found for 17% of orthologs. The discovery of premature termination codons in 38% of expressed genes was consistent with ongoing pseudogenization of the wheat genome. The rates of AS within the individual wheat subgenomes (21%-25%) were similar to diploid plants. However, we uncovered a high level of AS pattern divergence between the duplicated homeologous copies of genes. Our results are consistent with the accelerated accumulation of AS isoforms, nonsynonymous mutations, and gene structure rearrangements in the wheat lineage, likely due to genetic redundancy created by WGDs. Whereas these processes mostly contribute to the degeneration of a duplicated genome and its diploidization, they have the potential to facilitate the origin of new functional variations, which, upon selection in the evolutionary lineage, may play an important role in the origin of novel traits.

What Can Animal Models Teach Us About the Potential Safety and Efficacy of Regulatory Macrophage Therapy in Patients?

Cell-based immunosuppressive therapy is very much in its infancy, but with the inception of The ONE Study, a coordinated effort to bring such treatments to clinical trial is being made; therefore, we must ask, how great are the benefits of regulatory cell therapy in transplantation expected to be? Adoptive transfer of transplant tolerance with isolated populations of regulatory cell types is reliably achieved in homologous animal models, but whether such observations can be extrapolated to patients hinges on proving the equivalence between the animal and human cells in question. Our work centers on human regulatory macrophages (M reg) as a means of attenuating responses against alloantigen. The mouse proves a convenient model to test the combined effect of M reg treatment and conventional immunosuppression, but is only clinically relevant if mouse and human M regs are functional counterparts. M regs from both species arise from equivalent populations of monocytes under directly analogous culture conditions and share many morphological features and express unique markers, distinguishing them from other macrophages. Mouse and human M regs profoundly inhibit allogeneic T cell proliferation in co-culture, through deletion and phagocytosis of activated T cells, and a direct suppressive activity. Yet, differences between mouse and human M regs exist, and weighing these differences against the likenesses is necessarily subjective, so gene expression profiling by microarray was adopted as an unbiased approach to evaluating their true similarity. In co-clustering analyses, M regs from the two species grouped together when the entire set of orthologous genes was considered, so identifying human M regs as the closest counterpart of mouse M regs.[Figure]In a mouse heart transplant model, infusion of donor-strain M regs prior to transplantation significantly prolonged allograft survival in unconditioned, non-immunosuppressed, fully-allogeneic recipients, whereas recipient and third party-derived M regs were not effective. Co-treatment with M regs and rapamycin for 10-days post-transplant enhanced the effect of M regs compared to M reg or rapamycin treatment alone. After injection, M regs from CD45.2+ B6 mice or CD45.2+ BALB/c mice were detectable in the lung, blood, liver and spleen of CD45.1+ B6 recipients for up to two weeks. Donor-derived human M regs were administered to two renal transplant recipients 7-days prior to transplantation. Both patients were minimised to low-dose tacrolimus monotherapy within 24 weeks of transplantation and subsequently maintained excellent graft function. To assess their fate after infusion, M regs were labelled with Indium111-oxine for tracking in serial whole-body SPECT studies, which showed most M regs remained viable and trafficked from the pulmonary vasculature via the blood to liver, spleen and bone marrow.

DISCO-SCA and Properly Applied GSVD as Swinging Methods to Find Common and Distinctive Processes

BackgroundIn systems biology it is common to obtain for the same set of biological entities information from multiple sources. Examples include expression data for the same set of orthologous genes screened in different organisms and data on the same set of culture samples obtained with different high-throughput techniques. A major challenge is to find the important biological processes underlying the data and to disentangle therein processes common to all data sources and processes distinctive for a specific source. Recently, two promising simultaneous data integration methods have been proposed to attain this goal, namely generalized singular value decomposition (GSVD) and simultaneous component analysis with rotation to common and distinctive components (DISCO-SCA).ResultsBoth theoretical analyses and applications to biologically relevant data show that: (1) straightforward applications of GSVD yield unsatisfactory results, (2) DISCO-SCA performs well, (3) provided proper pre-processing and algorithmic adaptations, GSVD reaches a performance level similar to that of DISCO-SCA, and (4) DISCO-SCA is directly generalizable to more than two data sources. The biological relevance of DISCO-SCA is illustrated with two applications. First, in a setting of comparative genomics, it is shown that DISCO-SCA recovers a common theme of cell cycle progression and a yeast-specific response to pheromones. The biological annotation was obtained by applying Gene Set Enrichment Analysis in an appropriate way. Second, in an application of DISCO-SCA to metabolomics data for Escherichia coli obtained with two different chemical analysis platforms, it is illustrated that the metabolites involved in some of the biological processes underlying the data are detected by one of the two platforms only; therefore, platforms for microbial metabolomics should be tailored to the biological question.ConclusionsBoth DISCO-SCA and properly applied GSVD are promising integrative methods for finding common and distinctive processes in multisource data. Open source code for both methods is provided.

Genes Associated with Honey Bee Behavioral Maturation Affect Clock-Dependent and -Independent Aspects of Daily Rhythmic Activity in Fruit Flies

BackgroundIn the honey bee, the age-related and socially regulated transition of workers from in-hive task performance (e.g., caring for young) to foraging (provisioning the hive) is associated with changes in many behaviors including the 24-hour pattern of rhythmic activity. We have previously shown that the hive-bee to forager transition is associated with extensive changes in brain gene expression. In this study, we test the possible function of a subset of these genes in daily rhythmic activity pattern using neural-targeted RNA interference (RNAi) of an orthologous gene set in Drosophila melanogaster.Principal FindingsOf 10 genes tested, knockdown of six affected some aspect of locomotor activity under a 12 h∶12 h light:dark regime (LD). Inos affected anticipatory activity preceding lights-off, suggesting a possible clock-dependent function. BM-40-SPARC, U2af50 and fax affected peak activity at dawn without affecting anticipation or overall inactivity (proportion of 15-min intervals without activity), suggesting that these effects may depend on the day-night light cycle. CAH1 affected overall inactivity. The remaining gene, abl, affected peak activity levels but was not clearly time-of-day-specific. No gene tested affected length of period or strength of rhythmicity in constant dark (DD), suggesting that these genes do not act in the core clock.SignificanceTaking advantage of Drosophila molecular genetic tools, our study provides an important step in understanding the large set of gene expression changes that occur in the honey bee transition from hive bee to forager. We show that orthologs of many of these genes influence locomotor activity in Drosophila, possibly through both clock-dependent and -independent pathways. Our results support the importance of both circadian clock and direct environmental stimuli (apart from entrainment) in shaping the bee’s 24-hour pattern of activity. Our study also outlines a new approach to dissecting complex behavior in a social animal.

Development of conserved ortholog set markers linked to the restorer gene Rfp1 in rye

Restoration of male fertility is a prerequisite for hybrid rye breeding and currently the most straightforward approach to minimize ergot infection in hybrid rye varieties. Molecular markers are important tools for the efficient introgression and management of restorer genes like Rfp1 originating from unadapted genetic resources. Furthermore, closely linked markers flanking Rfp1 are indispensible for identifying and selecting individuals with haplotypes showing recombination between Rfp1 and other gene(s) that reside in close proximity and have a negative influence on yield. We identified orthologous gene sets in rice, Brachypodium, and Sorghum and used these gene models as templates to establish conserved ortholog set (COS) markers for the restorer gene Rfp1 on the long arm of rye chromosome 4R. The novel co-dominant markers delimit Rfp1 within a 0.7-cM interval and allow prediction of Rfp1 genotypes with a precision not feasible before. The COS markers enabled an alignment of the improved genetic map of rye chromosome 4R with wheat and barley maps and allowed identification of regions orthologous to Rfp1 in wheat and barley on the short arms of chromosomes 6D and 6H, respectively. Results obtained in this study revealed that micro-collinearity around the Rfp1 locus in rye is affected by rearrangements relative to other grass genomes. The impact of the novel COS markers for practical hybrid rye breeding is discussed.

Development of gene-based markers and construction of an integrated linkage map in eggplant by using Solanum orthologous (SOL) gene sets

We constructed an integrated DNA marker linkage map of eggplant (Solanum melongena L.) using DNA marker segregation data sets obtained from two independent intraspecific F(2) populations. The linkage map consisted of 12 linkage groups and encompassed 1,285.5 cM in total. We mapped 952 DNA markers, including 313 genomic SSR markers developed by random sequencing of simple sequence repeat (SSR)-enriched genomic libraries, and 623 single-nucleotide polymorphisms (SNP) and insertion/deletion polymorphisms (InDels) found in eggplant-expressed sequence tags (ESTs) and related genomic sequences [introns and untranslated regions (UTRs)]. Because of their co-dominant inheritance and their highly polymorphic and multi-allelic nature, the SSR markers may be more versatile than the SNP and InDel markers for map-based genetic analysis of any traits of interest using segregating populations derived from any intraspecific crosses of practical breeding materials. However, we found that the distribution of microsatellites in the genome was biased to some extent, and therefore a considerable part of the eggplant genome was first detected when gene-derived SNP and InDel markers were mapped. Of the 623 SNP and InDel markers mapped onto the eggplant integrated map, 469 were derived from eggplant unigenes contained within Solanum orthologous (SOL) gene sets (i.e., sets of orthologous unigenes from eggplant, tomato, and potato). Out of the 469 markers, 326 could also be mapped onto the tomato map. These common markers will be informative landmarks for the transfer of tomato's more saturated genomic information to eggplant and will also provide comparative information on the genome organization of the two solanaceous species. The data are available from the DNA marker database of vegetables, VegMarks (http://vegmarks.nivot.affrc.go.jp).

Proteome adaptation to high temperatures in the ectothermic hydrothermal vent Pompeii worm.

Taking advantage of the massive genome sequencing effort made on thermophilic prokaryotes, thermal adaptation has been extensively studied by analysing amino acid replacements and codon usage in these unicellular organisms. In most cases, adaptation to thermophily is associated with greater residue hydrophobicity and more charged residues. Both of these characteristics are positively correlated with the optimal growth temperature of prokaryotes. In contrast, little information has been collected on the molecular ‘adaptive’ strategy of thermophilic eukaryotes. The Pompeii worm A. pompejana, whose transcriptome has recently been sequenced, is currently considered as the most thermotolerant eukaryote on Earth, withstanding the greatest thermal and chemical ranges known. We investigated the amino-acid composition bias of ribosomal proteins in the Pompeii worm when compared to other lophotrochozoans and checked for putative adaptive changes during the course of evolution using codon-based Maximum likelihood analyses. We then provided a comparative analysis of codon usage and amino-acid replacements from a greater set of orthologous genes between the Pompeii worm and Paralvinella grasslei, one of its closest relatives living in a much cooler habitat. Analyses reveal that both species display the same high GC-biased codon usage and amino-acid patterns favoring both positively-charged residues and protein hydrophobicity. These patterns may be indicative of an ancestral adaptation to the deep sea and/or thermophily. In addition, the Pompeii worm displays a set of amino-acid change patterns that may explain its greater thermotolerance, with a significant increase in Tyr, Lys and Ala against Val, Met and Gly. Present results indicate that, together with a high content in charged residues, greater proportion of smaller aliphatic residues, and especially alanine, may be a different path for metazoans to face relatively ‘high’ temperatures and thus a novelty in thermophilic metazoans.

Genome-Wide Analysis of Syntenic Gene Deletion in the Grasses

The grasses, Poaceae, are one of the largest and most successful angiosperm families. Like many radiations of flowering plants, the divergence of the major grass lineages was preceded by a whole-genome duplication (WGD), although these events are not rare for flowering plants. By combining identification of syntenic gene blocks with measures of gene pair divergence and different frequencies of ancient gene loss, we have separated the two subgenomes present in modern grasses. Reciprocal loss of duplicated genes or genomic regions has been hypothesized to reproductively isolate populations and, thus, speciation. However, in contrast to previous studies in yeast and teleost fishes, we found very little evidence of reciprocal loss of homeologous genes between the grasses, suggesting that post-WGD gene loss may not be the cause of the grass radiation. The sets of homeologous and orthologous genes and predicted locations of deleted genes identified in this study, as well as links to the CoGe comparative genomics web platform for analyzing pan-grass syntenic regions, are provided along with this paper as a resource for the grass genetics community.

PGAT: a multistrain analysis resource for microbial genomes

Motivation: The Prokaryotic-genome Analysis Tool (PGAT) is a web-based database application for comparing gene content and sequence across multiple microbial genomes facilitating the discovery of genetic differences that may explain observed phenotypes. PGAT supports database queries to identify genes that are present or absent in user-selected genomes, comparison of sequence polymorphisms in sets of orthologous genes, multigenome display of regions surrounding a query gene, comparison of the distribution of genes in metabolic pathways and manual community annotation.Availability and Implementation:The PGAT website may be accessed at http://nwrce.org/pgat.Contact: mbrittna@uw.edu

Phylogeny-wide analysis of social amoeba genomes highlights ancient origins for complex intercellular communication.

Dictyostelium discoideum (DD), an extensively studied model organism for cell and developmental biology, belongs to the most derived group 4 of social amoebas, a clade of altruistic multicellular organisms. To understand genome evolution over long time periods and the genetic basis of social evolution, we sequenced the genomes of Dictyostelium fasciculatum (DF) and Polysphondylium pallidum (PP), which represent the early diverging groups 1 and 2, respectively. In contrast to DD, PP and DF have conventional telomere organization and strongly reduced numbers of transposable elements. The number of protein-coding genes is similar between species, but only half of them comprise an identifiable set of orthologous genes. In general, genes involved in primary metabolism, cytoskeletal functions and signal transduction are conserved, while genes involved in secondary metabolism, export, and signal perception underwent large differential gene family expansions. This most likely signifies involvement of the conserved set in core cell and developmental mechanisms, and of the diverged set in niche- and species-specific adaptations for defense and food, mate, and kin selection. Phylogenetic dating using a concatenated data set and extensive loss of synteny indicate that DF, PP, and DD split from their last common ancestor at least 0.6 billion years ago.

MultiMSOAR 2.0: An Accurate Tool to Identify Ortholog Groups among Multiple Genomes

The identification of orthologous genes shared by multiple genomes plays an important role in evolutionary studies and gene functional analyses. Based on a recently developed accurate tool, called MSOAR 2.0, for ortholog assignment between a pair of closely related genomes based on genome rearrangement, we present a new system MultiMSOAR 2.0, to identify ortholog groups among multiple genomes in this paper. In the system, we construct gene families for all the genomes using sequence similarity search and clustering, run MSOAR 2.0 for all pairs of genomes to obtain the pairwise orthology relationship, and partition each gene family into a set of disjoint sets of orthologous genes (called super ortholog groups or SOGs) such that each SOG contains at most one gene from each genome. For each such SOG, we label the leaves of the species tree using 1 or 0 to indicate if the SOG contains a gene from the corresponding species or not. The resulting tree is called a tree of ortholog groups (or TOGs). We then label the internal nodes of each TOG based on the parsimony principle and some biological constraints. Ortholog groups are finally identified from each fully labeled TOG. In comparison with a popular tool MultiParanoid on simulated data, MultiMSOAR 2.0 shows significantly higher prediction accuracy. It also outperforms MultiParanoid, the Roundup multi-ortholog repository and the Ensembl ortholog database in real data experiments using gene symbols as a validation tool. In addition to ortholog group identification, MultiMSOAR 2.0 also provides information about gene births, duplications and losses in evolution, which may be of independent biological interest. Our experiments on simulated data demonstrate that MultiMSOAR 2.0 is able to infer these evolutionary events much more accurately than a well-known software tool Notung. The software MultiMSOAR 2.0 is available to the public for free.

Computational methods for Gene Orthology inference

Accurate inference of orthologous genes is a pre-requisite for most comparative genomics studies, and is also important for functional annotation of new genomes. Identification of orthologous gene sets typically involves phylogenetic tree analysis, heuristic algorithms based on sequence conservation, synteny analysis, or some combination of these approaches. The most direct tree-based methods typically rely on the comparison of an individual gene tree with a species tree. Once the two trees are accurately constructed, orthologs are straightforwardly identified by the definition of orthology as those homologs that are related by speciation, rather than gene duplication, at their most recent point of origin. Although ideal for the purpose of orthology identification in principle, phylogenetic trees are computationally expensive to construct for large numbers of genes and genomes, and they often contain errors, especially at large evolutionary distances. Moreover, in many organisms, in particular prokaryotes and viruses, evolution does not appear to have followed a simple 'tree-like' mode, which makes conventional tree reconciliation inapplicable. Other, heuristic methods identify probable orthologs as the closest homologous pairs or groups of genes in a set of organisms. These approaches are faster and easier to automate than tree-based methods, with efficient implementations provided by graph-theoretical algorithms enabling comparisons of thousands of genomes. Comparisons of these two approaches show that, despite conceptual differences, they produce similar sets of orthologs, especially at short evolutionary distances. Synteny also can aid in identification of orthologs. Often, tree-based, sequence similarity- and synteny-based approaches can be combined into flexible hybrid methods.

Comparative genomics of the social amoebae Dictyostelium discoideum and Dictyostelium purpureum

BackgroundThe social amoebae (Dictyostelia) are a diverse group of Amoebozoa that achieve multicellularity by aggregation and undergo morphogenesis into fruiting bodies with terminally differentiated spores and stalk cells. There are four groups of dictyostelids, with the most derived being a group that contains the model species Dictyostelium discoideum.ResultsWe have produced a draft genome sequence of another group dictyostelid, Dictyostelium purpureum, and compare it to the D. discoideum genome. The assembly (8.41 × coverage) comprises 799 scaffolds totaling 33.0 Mb, comparable to the D. discoideum genome size. Sequence comparisons suggest that these two dictyostelids shared a common ancestor approximately 400 million years ago. In spite of this divergence, most orthologs reside in small clusters of conserved synteny. Comparative analyses revealed a core set of orthologous genes that illuminate dictyostelid physiology, as well as differences in gene family content. Interesting patterns of gene conservation and divergence are also evident, suggesting function differences; some protein families, such as the histidine kinases, have undergone little functional change, whereas others, such as the polyketide synthases, have undergone extensive diversification. The abundant amino acid homopolymers encoded in both genomes are generally not found in homologous positions within proteins, so they are unlikely to derive from ancestral DNA triplet repeats. Genes involved in the social stage evolved more rapidly than others, consistent with either relaxed selection or accelerated evolution due to social conflict.ConclusionsThe findings from this new genome sequence and comparative analysis shed light on the biology and evolution of the Dictyostelia.

Community transcriptomics reveals universal patterns of protein sequence conservation in natural microbial communities

BackgroundCombined metagenomic and metatranscriptomic datasets make it possible to study the molecular evolution of diverse microbial species recovered from their native habitats. The link between gene expression level and sequence conservation was examined using shotgun pyrosequencing of microbial community DNA and RNA from diverse marine environments, and from forest soil.ResultsAcross all samples, expressed genes with transcripts in the RNA sample were significantly more conserved than non-expressed gene sets relative to best matches in reference databases. This discrepancy, observed for many diverse individual genomes and across entire communities, coincided with a shift in amino acid usage between these gene fractions. Expressed genes trended toward GC-enriched amino acids, consistent with a hypothesis of higher levels of functional constraint in this gene pool. Highly expressed genes were significantly more likely to fall within an orthologous gene set shared between closely related taxa (core genes). However, non-core genes, when expressed above the level of detection, were, on average, significantly more highly expressed than core genes based on transcript abundance normalized to gene abundance. Finally, expressed genes showed broad similarities in function across samples, being relatively enriched in genes of energy metabolism and underrepresented by genes of cell growth.ConclusionsThese patterns support the hypothesis, predicated on studies of model organisms, that gene expression level is a primary correlate of evolutionary rate across diverse microbial taxa from natural environments. Despite their complexity, meta-omic datasets can reveal broad evolutionary patterns across taxonomically, functionally, and environmentally diverse communities.

ODB: a database for operon organizations, 2011 update

ODB (Operon DataBase) aims to collect data of all known and conserved operons in completely sequenced genomes. Three newly updated features of this database have been added as follows: (i) Data from included operons were updated. The genome-wide analysis of transcription and transcriptional units has become popular recently and ODB successfully integrates these high-throughput operon data, including genome-wide transcriptional units of five prokaryotes and two eukaryotes. The current version of our database contains information from about 10 000 known operons in more than 50 genomes, and more than 400 000 conserved operons obtained from more than 1000 bacterial genomes. (ii) ODB proposes the idea of reference operons as a new operon prediction tool. A reference operon, a set of possible orthologous genes that organize operons, is defined by clustering all known operons. A large number of known operons, including the recently added genome-wide analysis of operons, allowed us to define more reliable reference operons. (iii) ODB also provides new graphical interfaces. One is for comparative analyses of operon structures in multiple genomes. The other is for visualization of possible operons in multiple genomes obtained from the reference operons. The 2011 updated version of ODB is now available at http://operondb.jp/.

Comparative genomics of Lactobacillus

SummaryThe genus Lactobacillus includes a diverse group of bacteria consisting of many species that are associated with fermentations of plants, meat or milk. In addition, various lactobacilli are natural inhabitants of the intestinal tract of humans and other animals. Finally, several Lactobacillus strains are marketed as probiotics as their consumption can confer a health benefit to host. Presently, 154 Lactobacillus species are known and a growing fraction of these are subject to draft genome sequencing. However, complete genome sequences are needed to provide a platform for detailed genomic comparisons. Therefore, we selected a total of 20 genomes of various Lactobacillus strains for which complete genomic sequences have been reported. These genomes had sizes varying from 1.8 to 3.3 Mb and other characteristic features, such as G+C content that ranged from 33% to 51%. The Lactobacillus pan genome was found to consist of approximately 14 000 protein‐encoding genes while all 20 genomes shared a total of 383 sets of orthologous genes that defined the Lactobacillus core genome (LCG). Based on advanced phylogeny of the proteins encoded by this LCG, we grouped the 20 strains into three main groups and defined core group genes present in all genomes of a single group, signature group genes shared in all genomes of one group but absent in all other Lactobacillus genomes, and Group‐specific ORFans present in core group genes of one group and absent in all other complete genomes. The latter are of specific value in defining the different groups of genomes. The study provides a platform for present individual comparisons as well as future analysis of new Lactobacillus genomes.