Correction: PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph.

[This corrects the article DOI: 10.1371/journal.pcbi.1007732.].

Tempo and mode of genome evolution in a 50,000-generation experiment.

Adaptation by natural selection depends on the rates, effects and interactions of many mutations, making it difficult to determine what proportion of mutations in an evolving lineage are beneficial. Here we analysed 264 complete genomes from 12 Escherichia coli populations to characterize their dynamics over 50,000 generations. The populations that retained the ancestral mutation rate support a model in which most fixed mutations are beneficial, the fraction of beneficial mutations declines as fitness rises, and neutral mutations accumulate at a constant rate. We also compared these populations to mutation-accumulation lines evolved under a bottlenecking regime that minimizes selection. Nonsynonymous mutations, intergenic mutations, insertions and deletions are overrepresented in the long-term populations, further supporting the inference that most mutations that reached high frequency were favoured by selection. These results illuminate the shifting balance of forces that govern genome evolution in populations adapting to a new environment.

MicroScope-an integrated resource for community expertise of gene functions and comparative analysis of microbial genomic and metabolic data.

The overwhelming list of new bacterial genomes becoming available on a daily basis makes accurate genome annotation an essential step that ultimately determines the relevance of thousands of genomes stored in public databanks. The MicroScope platform (http://www.genoscope.cns.fr/agc/microscope) is an integrative resource that supports systematic and efficient revision of microbial genome annotation, data management and comparative analysis. Starting from the results of our syntactic, functional and relational annotation pipelines, MicroScope provides an integrated environment for the expert annotation and comparative analysis of prokaryotic genomes. It combines tools and graphical interfaces to analyze genomes and to perform the manual curation of gene function in a comparative genomics and metabolic context. In this article, we describe the free-of-charge MicroScope services for the annotation and analysis of microbial (meta)genomes, transcriptomic and re-sequencing data. Then, the functionalities of the platform are presented in a way providing practical guidance and help to the nonspecialists in bioinformatics. Newly integrated analysis tools (i.e. prediction of virulence and resistance genes in bacterial genomes) and original method recently developed (the pan-genome graph representation) are also described. Integrated environments such as MicroScope clearly contribute, through the user community, to help maintaining accurate resources.

PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph.

The use of comparative genomics for functional, evolutionary, and epidemiological studies requires methods to classify gene families in terms of occurrence in a given species. These methods usually lack multivariate statistical models to infer the partitions and the optimal number of classes and don't account for genome organization. We introduce a graph structure to model pangenomes in which nodes represent gene families and edges represent genomic neighborhood. Our method, named PPanGGOLiN, partitions nodes using an Expectation-Maximization algorithm based on multivariate Bernoulli Mixture Model coupled with a Markov Random Field. This approach takes into account the topology of the graph and the presence/absence of genes in pangenomes to classify gene families into persistent, cloud, and one or several shell partitions. By analyzing the partitioned pangenome graphs of isolate genomes from 439 species and metagenome-assembled genomes from 78 species, we demonstrate that our method is effective in estimating the persistent genome. Interestingly, it shows that the shell genome is a key element to understand genome dynamics, presumably because it reflects how genes present at intermediate frequencies drive adaptation of species, and its proportion in genomes is independent of genome size. The graph-based approach proposed by PPanGGOLiN is useful to depict the overall genomic diversity of thousands of strains in a compact structure and provides an effective basis for very large scale comparative genomics. The software is freely available at https://github.com/labgem/PPanGGOLiN.

High-Throughput Sequencing (HTS) of newly synthetized RNAs enables one shot detection and identification of live mycoplasmas and differentiation from inert nucleic acids.

Mycoplasma contamination threatens both the safety of biologics produced in cell substrates as well as the quality of scientific results based on cell-culture observations. Methods currently used to detect contamination of cells include culture, enzymatic activity, immunofluorescence and PCR but suffer from some limitations. High throughput sequencing (HTS) can be used to identify microbes like mycoplasmas in biologics since it enables an unbiased approach to detection without the need to design specific primers to pre-amplify target sequences but it does not enable the confirmation of microbial infection since this could reflect carryover of inert sequences. In order to unambiguously differentiate the presence of live or dead mycoplasmas in biological products, the present method was developed based on metabolic RNA labelling of newly synthetized mycoplasmal RNAs. HTS of labelled RNA detected A549 cell infection with Acholeplasma laidlawii in a manner similar to both PCR and culture and demonstrated that this technique can unambiguously identify bacterial species and differentiates infected cells from cells exposed to a high inoculum of heat-inactivated mycoplasmas. This method therefore combines the advantage of culture (that detects only live microorganisms) with those of molecular tests (rapidity) together with a very broad range of bacterial detection and identification.

The Arginine Deiminase Operon Is Responsible for a Fitness Trade-Off in Extended-Spectrum-ß-Lactamase-Producing Strains of Escherichia coli.

We previously identified an operon involved in an arginine deiminase (ADI) pathway (arc operon) on a CTX-M-producing plasmid from an O102-ST405 strain of Escherichia coli As the ADI pathway was shown to be involved in the virulence of various Gram-positive bacteria, we tested whether the ADI pathway could be involved in the epidemiological success of extended-spectrum-ß-lactamase (ESBL)-producing E. coli strains. We studied two collections of human E. coli isolated in France (n = 493) and England (n = 1,509) and show that the prevalence of the arc operon (i) is higher in ESBL-producing strains (12.1%) than in nonproducers (2.5%), (ii) is higher in CTX-M-producing strains (16%) than in other ESBL producers (3.5%), and (iii) increased over time in ESBL-producing strains from 0% before 2000 to 43.3% in 2011 to 2012. The arc operon, found in strains from various phylogenetic backgrounds, is carried by IncF plasmids (85%) or chromosomes (15%) in regions framed by numerous insertion sequences, indicating multiple arrivals. Competition experiments showed that the arc operon enhances fitness of the strain in vitro in lysogeny broth with arginine. In vivo competition experiments showed that the arc operon is advantageous for the strain in a mouse model of urinary tract infection (UTI), whereas it is a burden in a mouse model of intestinal colonization. In summary, we have identified a trait linked to CTX-M-producing strains that is responsible for a trade-off between two main E. coli lifestyles, UTI and gut commensalism. This trait alone cannot explain the wide spread of ESBLs in E. coli but merits epidemiological surveillance.

MicroScope in 2017: an expanding and evolving integrated resource for community expertise of microbial genomes.

The annotation of genomes from NGS platforms needs to be automated and fully integrated. However, maintaining consistency and accuracy in genome annotation is a challenging problem because millions of protein database entries are not assigned reliable functions. This shortcoming limits the knowledge that can be extracted from genomes and metabolic models. Launched in 2005, the MicroScope platform (http://www.genoscope.cns.fr/agc/microscope) is an integrative resource that supports systematic and efficient revision of microbial genome annotation, data management and comparative analysis. Effective comparative analysis requires a consistent and complete view of biological data, and therefore, support for reviewing the quality of functional annotation is critical. MicroScope allows users to analyze microbial (meta)genomes together with post-genomic experiment results if any (i.e. transcriptomics, re-sequencing of evolved strains, mutant collections, phenotype data). It combines tools and graphical interfaces to analyze genomes and to perform the expert curation of gene functions in a comparative context. Starting with a short overview of the MicroScope system, this paper focuses on some major improvements of the Web interface, mainly for the submission of genomic data and on original tools and pipelines that have been developed and integrated in the platform: computation of pan-genomes and prediction of biosynthetic gene clusters. Today the resource contains data for more than 6000 microbial genomes, and among the 2700 personal accounts (65% of which are now from foreign countries), 14% of the users are performing expert annotations, on at least a weekly basis, contributing to improve the quality of microbial genome annotations.

Expansion of the SOS regulon of Vibrio cholerae through extensive transcriptome analysis and experimental validation.

BACKGROUND: The SOS response is an almost ubiquitous response of cells to genotoxic stresses. The full complement of genes in the SOS regulon for Vibrio species has only been addressed through bioinformatic analyses predicting LexA binding box consensus and in vitro validation. Here, we perform whole transcriptome sequencing from Vibrio cholerae treated with mitomycin C as an SOS inducer to characterize the SOS regulon and other pathways affected by this treatment. RESULTS: Comprehensive transcriptional profiling allowed us to define the full landscape of promoters and transcripts active in V. cholerae. We performed extensive transcription start site (TSS) mapping as well as detection/quantification of the coding and non-coding RNA (ncRNA) repertoire in strain N16961. To improve TSS detection, we developed a new technique to treat RNA extracted from cells grown in various conditions. This allowed for identification of 3078 TSSs with an average 5'UTR of 116 nucleotides, and peak distribution between 16 and 64 nucleotides; as well as 629 ncRNAs. Mitomycin C treatment induced transcription of 737 genes and 28 ncRNAs at least 2 fold, while it repressed 231 genes and 17 ncRNAs. Data analysis revealed that in addition to the core genes known to integrate the SOS regulon, several metabolic pathways were induced. This study allowed for expansion of the Vibrio SOS regulon, as twelve genes (ubiEJB, tatABC, smpA, cep, VC0091, VC1190, VC1369-1370) were found to be co-induced with their adjacent canonical SOS regulon gene(s), through transcriptional read-through. Characterization of UV and mitomycin C susceptibility for mutants of these newly identified SOS regulon genes and other highly induced genes and ncRNAs confirmed their role in DNA damage rescue and protection. CONCLUSIONS: We show that genotoxic stress induces a pervasive transcriptional response, affecting almost 20% of the V. cholerae genes. We also demonstrate that the SOS regulon is larger than previously known, and its syntenic organization is conserved among Vibrio species. Furthermore, this specific co-localization is found in other γ-proteobacteria for genes recN-smpA and rmuC-tatABC, suggesting SOS regulon conservation in this phylum. Finally, we comment on the limitations of widespread NGS approaches for identification of all RNA species in bacteria.

Recruitment of a Lineage-Specific Virulence Regulatory Pathway Promotes Intracellular Infection by a Plant Pathogen Experimentally Evolved into a Legume Symbiont.

Ecological transitions between different lifestyles, such as pathogenicity, mutualism and saprophytism, have been very frequent in the course of microbial evolution, and often driven by horizontal gene transfer. Yet, how genomes achieve the ecological transition initiated by the transfer of complex biological traits remains poorly known. Here, we used experimental evolution, genomics, transcriptomics and high-resolution phenotyping to analyze the evolution of the plant pathogen Ralstonia solanacearum into legume symbionts, following the transfer of a natural plasmid encoding the essential mutualistic genes. We show that a regulatory pathway of the recipient R. solanacearum genome involved in extracellular infection of natural hosts was reused to improve intracellular symbiosis with the Mimosa pudica legume. Optimization of intracellular infection capacity was gained through mutations affecting two components of a new regulatory pathway, the transcriptional regulator efpR and a region upstream from the RSc0965-0967 genes of unknown functions. Adaptive mutations caused the downregulation of efpR and the over-expression of a downstream regulatory module, the three unknown genes RSc3146-3148, two of which encoding proteins likely associated to the membrane. This over-expression led to important metabolic and transcriptomic changes and a drastic qualitative and quantitative improvement of nodule intracellular infection. In addition, these adaptive mutations decreased the virulence of the original pathogen. The complete efpR/RSc3146-3148 pathway could only be identified in the genomes of the pathogenic R. solanacearum species complex. Our findings illustrate how the rewiring of a genetic network regulating virulence allows a radically different type of symbiotic interaction and contributes to ecological transitions and trade-offs.

Genomic study of a novel magnetotactic Alphaproteobacteria uncovers the multiple ancestry of magnetotaxis.

Ecological and evolutionary processes involved in magnetotactic bacteria (MTB) adaptation to their environment have been a matter of debate for many years. Ongoing efforts for their characterization are progressively contributing to understand these processes, including the genetic and molecular mechanisms responsible for biomineralization. Despite numerous culture-independent MTB characterizations, essentially within the Proteobacteria phylum, only few species have been isolated in culture because of their complex growth conditions. Here, we report a newly cultivated magnetotactic, microaerophilic and chemoorganoheterotrophic bacterium isolated from the Mediterranean Sea in Marseille, France: Candidatus Terasakiella magnetica strain PR-1 that belongs to an Alphaproteobacteria genus with no magnetotactic relative. By comparing the morphology and the whole genome shotgun sequence of this MTB with those of closer relatives, we brought further evidence that the apparent vertical ancestry of magnetosome genes suggested by previous studies within Alphaproteobacteria hides a more complex evolutionary history involving horizontal gene transfers and/or duplication events before and after the emergence of Magnetospirillum, Magnetovibrio and Magnetospira genera. A genome-scale comparative genomics analysis identified several additional candidate functions and genes that could be specifically associated to MTB lifestyle in this class of bacteria.

Strand-specific transcriptomes of Enterohemorrhagic Escherichia coli in response to interactions with ground beef microbiota: interactions between microorganisms in raw meat.

BACKGROUND: Enterohemorrhagic Escherichia coli (EHEC) are zoonotic agents associated with outbreaks worldwide. Growth of EHEC strains in ground beef could be inhibited by background microbiota that is present initially at levels greater than that of the pathogen E. coli. However, how the microbiota outcompetes the pathogenic bacteria is unknown. Our objective was to identify metabolic pathways of EHEC that were altered by natural microbiota in order to improve our understanding of the mechanisms controlling the growth and survival of EHECs in ground beef. RESULTS: Based on 16S metagenomics analysis, we identified the microbial community structure in our beef samples which was an essential preliminary for subtractively analyzing the gene expression of the EHEC strains. Then, we applied strand-specific RNA-seq to investigate the effects of this microbiota on the global gene expression of EHEC O2621765 and O157EDL933 strains by comparison with their behavior in beef meat without microbiota. In strain O2621765, the expression of genes connected with nitrate metabolism and nitrite detoxification, DNA repair, iron and nickel acquisition and carbohydrate metabolism, and numerous genes involved in amino acid metabolism were down-regulated. Further, the observed repression of ftsL and murF, involved respectively in building the cytokinetic ring apparatus and in synthesizing the cytoplasmic precursor of cell wall peptidoglycan, might help to explain the microbiota's inhibitory effect on EHECs. For strain O157EDL933, the induced expression of the genes implicated in detoxification and the general stress response and the repressed expression of the peR gene, a gene negatively associated with the virulence phenotype, might be linked to the survival and virulence of O157:H7 in ground beef with microbiota. CONCLUSION: In the present study, we show how RNA-Seq coupled with a 16S metagenomics analysis can be used to identify the effects of a complex microbial community on relevant functions of an individual microbe within it. These findings add to our understanding of the behavior of EHECs in ground beef. By measuring transcriptional responses of EHEC, we could identify putative targets which may be useful to develop new strategies to limit their shedding in ground meat thus reducing the risk of human illnesses.

Experimental evolution of rhizobia may lead to either extra- or intracellular symbiotic adaptation depending on the selection regime.

Experimental evolution is a powerful approach to study the process of adaptation to new environments, including the colonization of eukaryotic hosts. Facultative endosymbionts, including pathogens and mutualists, face changing and spatially structured environments during the symbiotic process, which impose diverse selection pressures. Here, we provide evidence that different selection regimes, involving different times spent in the plant environment, can result in either intra- or extracellular symbiotic adaptations. In previous work, we introduced the symbiotic plasmid of Cupriavidus taiwanensis, the rhizobial symbiont of Mimosa pudica, into the phytopathogen Ralstonia solanacearum and selected three variants able to form root nodules on M. pudica, two (CBM212 and CBM349) being able to rudimentarily infect nodule cells and the third one (CBM356) only capable of extracellular infection of nodules. Each nodulating ancestor was further challenged to evolve using serial ex planta-in planta cycles of either 21 (three short-cycle lineages) or 42 days (three long-cycle lineages). In this study, we compared the phenotype of the 18 final evolved clones. Evolution through short and long cycles resulted in similar adaptive paths on lineages deriving from the two intracellularly infectious ancestors, CBM212 and CBM349. In contrast, only short cycles allowed a stable acquisition of intracellular infection in lineages deriving from the extracellularly infecting ancestor, CBM356. Long cycles, instead, favoured improvement of extracellular infection. Our work highlights the importance of the selection regime in shaping desired traits during host-mediated selection experiments.

Transient hypermutagenesis accelerates the evolution of legume endosymbionts following horizontal gene transfer.

Horizontal gene transfer (HGT) is an important mode of adaptation and diversification of prokaryotes and eukaryotes and a major event underlying the emergence of bacterial pathogens and mutualists. Yet it remains unclear how complex phenotypic traits such as the ability to fix nitrogen with legumes have successfully spread over large phylogenetic distances. Here we show, using experimental evolution coupled with whole genome sequencing, that co-transfer of imuABC error-prone DNA polymerase genes with key symbiotic genes accelerates the evolution of a soil bacterium into a legume symbiont. Following introduction of the symbiotic plasmid of Cupriavidus taiwanensis, the Mimosa symbiont, into pathogenic Ralstonia solanacearum we challenged transconjugants to become Mimosa symbionts through serial plant-bacteria co-cultures. We demonstrate that a mutagenesis imuABC cassette encoded on the C. taiwanensis symbiotic plasmid triggered a transient hypermutability stage in R. solanacearum transconjugants that occurred before the cells entered the plant. The generated burst in genetic diversity accelerated symbiotic adaptation of the recipient genome under plant selection pressure, presumably by improving the exploration of the fitness landscape. Finally, we show that plasmid imuABC cassettes are over-represented in rhizobial lineages harboring symbiotic plasmids. Our findings shed light on a mechanism that may have facilitated the dissemination of symbiotic competency among α- and ß-proteobacteria in natura and provide evidence for the positive role of environment-induced mutagenesis in the acquisition of a complex lifestyle trait. We speculate that co-transfer of complex phenotypic traits with mutagenesis determinants might frequently enhance the ecological success of HGT.

Contrasting effects of historical contingency on phenotypic and genomic trajectories during a two-step evolution experiment with bacteria.

BACKGROUND: The impact of historical contingency, i.e. the past evolutionary history of a population, on further adaptation is mostly unknown at both the phenotypic and genomic levels. We addressed this question using a two-step evolution experiment. First, replicate populations of Escherichia coli were propagated in four different environmental conditions for 1000 generations. Then, all replicate populations were transferred and propagated for further 1000 generations to a single new environment. RESULTS: Using this two-step experimental evolution strategy, we investigated, at both the phenotypic and genomic levels, whether and how adaptation in the initial historical environments impacted evolutionary trajectories in a new environment. We showed that both the growth rate and fitness of the evolved populations obtained after the second step of evolution were contingent upon past evolutionary history. In contrast however, the genes that were modified during the second step of evolution were independent from the previous history of the populations. CONCLUSIONS: Our work suggests that historical contingency affects phenotypic adaptation to a new environment. This was however not reflected at the genomic level implying complex relationships between environmental factors and the genotype-to-phenotype map.

Synonymous Genetic Variation in Natural Isolates of Escherichia coli Does Not Predict Where Synonymous Substitutions Occur in a Long-Term Experiment.

Synonymous genetic differences vary by more than 20-fold among genes in natural isolates of Escherichia coli. One hypothesis to explain this heterogeneity is that genes with high levels of synonymous variation mutate at higher rates than genes with low synonymous variation. If so, then one would expect to observe similar mutational patterns in evolution experiments. In fact, however, the pattern of synonymous substitutions in a long-term evolution experiment with E. coli does not support this hypothesis. In particular, the extent of synonymous variation across genes in that experiment does not reflect the variation observed in natural isolates of E. coli. Instead, gene length alone predicts with high accuracy the prevalence of synonymous changes in the experimental populations. We hypothesize that patterns of synonymous variation in natural E. coli populations are instead caused by differences across genomic regions in their effective population size that, in turn, reflect different histories of recombination, horizontal gene transfer, selection, and population structure.

The revisited genome of Pseudomonas putida KT2440 enlightens its value as a robust metabolic chassis.

By the time the complete genome sequence of the soil bacterium Pseudomonas putida KT2440 was published in 2002 (Nelson et al., ) this bacterium was considered a potential agent for environmental bioremediation of industrial waste and a good colonizer of the rhizosphere. However, neither the annotation tools available at that time nor the scarcely available omics data-let alone metabolic modeling and other nowadays common systems biology approaches-allowed them to anticipate the astonishing capacities that are encoded in the genetic complement of this unique microorganism. In this work we have adopted a suite of state-of-the-art genomic analysis tools to revisit the functional and metabolic information encoded in the chromosomal sequence of strain KT2440. We identified 242 new protein-coding genes and re-annotated the functions of 1548 genes, which are linked to almost 4900 PubMed references. Catabolic pathways for 92 compounds (carbon, nitrogen and phosphorus sources) that could not be accommodated by the previously constructed metabolic models were also predicted. The resulting examination not only accounts for some of the known stress tolerance traits known in P. putida but also recognizes the capacity of this bacterium to perform difficult redox reactions, thereby multiplying its value as a platform microorganism for industrial biotechnology.

RNase J depletion leads to massive changes in mRNA abundance in Helicobacter pylori.

Degradation of RNA as an intermediate message between genes and corresponding proteins is important for rapid attenuation of gene expression and maintenance of cellular homeostasis. This process is controlled by ribonucleases that have different target specificities. In the bacterial pathogen Helicobacter pylori, an exo- and endoribonuclease RNase J is essential for growth. To explore the role of RNase J in H. pylori, we identified its putative targets at a global scale using next generation RNA sequencing. We found that strong depletion for RNase J led to a massive increase in the steady-state levels of non-rRNAs. mRNAs and RNAs antisense to open reading frames were most affected with over 80% increased more than 2-fold. Non-coding RNAs expressed in the intergenic regions were much less affected by RNase J depletion. Northern blotting of selected messenger and non-coding RNAs validated these results. Globally, our data suggest that RNase J of H. pylori is a major RNase involved in degradation of most cellular RNAs.

Mutation rate dynamics in a bacterial population reflect tension between adaptation and genetic load.

Mutations are the ultimate source of heritable variation for evolution. Understanding how mutation rates themselves evolve is thus essential for quantitatively understanding many evolutionary processes. According to theory, mutation rates should be minimized for well-adapted populations living in stable environments, whereas hypermutators may evolve if conditions change. However, the long-term fate of hypermutators is unknown. Using a phylogenomic approach, we found that an adapting Escherichia coli population that first evolved a mutT hypermutator phenotype was later invaded by two independent lineages with mutY mutations that reduced genome-wide mutation rates. Applying neutral theory to synonymous substitutions, we dated the emergence of these mutations and inferred that the mutT mutation increased the point-mutation rate by ∼150-fold, whereas the mutY mutations reduced the rate by ∼40-60%, with a corresponding decrease in the genetic load. Thus, the long-term fate of the hypermutators was governed by the selective advantage arising from a reduced mutation rate as the potential for further adaptation declined.

Comparative genomics of 43 strains of Xanthomonas citri pv. citri reveals the evolutionary events giving rise to pathotypes with different host ranges.

BACKGROUND: The identification of factors involved in the host range definition and evolution is a pivotal challenge in the goal to predict and prevent the emergence of plant bacterial disease. To trace the evolution and find molecular differences between three pathotypes of Xanthomonas citri pv. citri that may explain their distinctive host ranges, 42 strains of X. citri pv. citri and one outgroup strain, Xanthomonas citri pv. bilvae were sequenced and compared. RESULTS: The strains from each pathotype form monophyletic clades, with a short branch shared by the A(w) and A pathotypes. Pathotype-specific recombination was detected in seven regions of the alignment. Using Ancestral Character Estimation, 426 SNPs were mapped to the four branches at the base of the A, A*, A(w) and A/A(w) clades. Several genes containing pathotype-specific nonsynonymous mutations have functions related to pathogenicity. The A pathotype is enriched for SNP-containing genes involved in defense mechanisms, while A* is significantly depleted for genes that are involved in transcription. The pathotypes differ by four gene islands that largely coincide with regions of recombination and include genes with a role in virulence. Both A* and A(w) are missing genes involved in defense mechanisms. In contrast to a recent study, we find that there are an extremely small number of pathotype-specific gene presences and absences. CONCLUSIONS: The three pathotypes of X. citri pv. citri that differ in their host ranges largely show genomic differences related to recombination, horizontal gene transfer and single nucleotide polymorphism. We detail the phylogenetic relationship of the pathotypes and provide a set of candidate genes involved in pathotype-specific evolutionary events that could explain to the differences in host range and pathogenicity between them.

MicroScope--an integrated microbial resource for the curation and comparative analysis of genomic and metabolic data.

MicroScope is an integrated platform dedicated to both the methodical updating of microbial genome annotation and to comparative analysis. The resource provides data from completed and ongoing genome projects (automatic and expert annotations), together with data sources from post-genomic experiments (i.e. transcriptomics, mutant collections) allowing users to perfect and improve the understanding of gene functions. MicroScope (http://www.genoscope.cns.fr/agc/microscope) combines tools and graphical interfaces to analyse genomes and to perform the manual curation of gene annotations in a comparative context. Since its first publication in January 2006, the system (previously named MaGe for Magnifying Genomes) has been continuously extended both in terms of data content and analysis tools. The last update of MicroScope was published in 2009 in the Database journal. Today, the resource contains data for >1600 microbial genomes, of which ∼300 are manually curated and maintained by biologists (1200 personal accounts today). Expert annotations are continuously gathered in the MicroScope database (∼50 000 a year), contributing to the improvement of the quality of microbial genomes annotations. Improved data browsing and searching tools have been added, original tools useful in the context of expert annotation have been developed and integrated and the website has been significantly redesigned to be more user-friendly. Furthermore, in the context of the European project Microme (Framework Program 7 Collaborative Project), MicroScope is becoming a resource providing for the curation and analysis of both genomic and metabolic data. An increasing number of projects are related to the study of environmental bacterial (meta)genomes that are able to metabolize a large variety of chemical compounds that may be of high industrial interest.