The comprehensive microbial resource

Looking for the mechanism of arsenate respiration of Fusibacter sp. strain 3D3, independent of ArrAB

The literature has reported the isolation of arsenate-dependent growing microorganisms which lack a canonical homolog for respiratory arsenate reductase, ArrAB. We recently isolated an arsenate-dependent growing bacterium from volcanic arsenic-bearing environments in Northern Chile, Fusibacter sp. strain 3D3 (Fas) and studied the arsenic metabolism in this Gram-positive isolate. Features of Fas deduced from genome analysis and comparative analysis with other arsenate-reducing microorganisms revealed the lack of ArrAB coding genes and the occurrence of two arsC genes encoding for putative cytoplasmic arsenate reductases named ArsC-1 and ArsC-2. Interestingly, ArsC-1 and ArsC-2 belong to the thioredoxin-coupled family (because of the redox-active disulfide protein used as reductant), but they conferred differential arsenate resistance to the E. coli WC3110 ΔarsC strain. PCR experiments confirmed the absence of arrAB genes and results obtained using uncouplers revealed that Fas growth is linked to the proton gradient. In addition, Fas harbors ferredoxin-NAD⁺ oxidoreductase (Rnf) and electron transfer flavoprotein (etf) coding genes. These are key molecular markers of a recently discovered flavin-based electron bifurcation mechanism involved in energy conservation, mainly in anaerobic metabolisms regulated by the cellular redox state and mostly associated with cytoplasmic enzyme complexes. At least three electron-bifurcating flavoenzyme complexes were evidenced in Fas, some of them shared in conserved genomic regions by other members of the Fusibacter genus. These physiological and genomic findings permit us to hypothesize the existence of an uncharacterized arsenate-dependent growth metabolism regulated by the cellular redox state in the Fusibacter genus.

Validation of a random Vibrio parahaemolyticus genomic library by selection of quinolone resistance in a heterologous host

Vibrio parahaemolyticus is a shellfish-borne pathogen that is a highly prevalent causative agent of inflammatory gastroenteritis in humans. Genomic libraries have proven useful for the identification of novel gene functions in many bacterial species. In this study we prepared a library containing 40 kb fragments of randomly sheared V. parahaemolyticus genomic DNA and introduced this into Escherichia coli HB101 using a commercially available low copy cosmid system. In order to estimate coverage and suitability of the library and potentially identify novel antimicrobial resistance determinants, we screened for the acquisition of resistance to the fluoroquinolone norfloxacin - a phenotype exhibited by V. parahaemolyticus but not the heterologous E. coli host. Upon selection on solid medium containing norfloxacin, 0.52% of the library population was resistant, consistent with the selection of a single resistance locus. End-sequencing identified six distinct insert fragments. All clones displayed fourfold increased norfloxacin MIC compared with E. coli HB101 carrying an empty vector. The common locus contained within resistant clones included qnr, a previously described quinolone resistance gene. These results indicate that the library was unbiased, of sufficient coverage and that heterologous expression was possible. While we hope that this library proves useful for identifying the genetic determinants of complex phenotypes such as those related to virulence, not all norfloxacin resistance genes were detected in our screen. As such, we discuss the benefits and limitations of this approach for identifying the genetic basis of uncharacterized bacterial phenotypes.

On the Impact of the Pangenome and Annotation Discrepancies While Building Protein Sequence Databases for Bacteria Proteogenomics

In proteomics, peptide information within mass spectrometry (MS) data from a specific organism sample is routinely matched against a protein sequence database that best represent such organism. However, if the species/strain in the sample is unknown or genetically poorly characterized, it becomes challenging to determine a database which can represent such sample. Building customized protein sequence databases merging multiple strains for a given species has become a strategy to overcome such restrictions. However, as more genetic information is publicly available and interesting genetic features such as the existence of pan- and core genes within a species are revealed, we questioned how efficient such merging strategies are to report relevant information. To test this assumption, we constructed databases containing conserved and unique sequences for 10 different species. Features that are relevant for probabilistic-based protein identification by proteomics were then monitored. As expected, increase in database complexity correlates with pangenomic complexity. However, Mycobacterium tuberculosis and Bordetella pertussis generated very complex databases even having low pangenomic complexity. We further tested database performance by using MS data from eight clinical strains from M. tuberculosis, and from two published datasets from Staphylococcus aureus. We show that by using an approach where database size is controlled by removing repeated identical tryptic sequences across strains/species, computational time can be reduced drastically as database complexity increases.

DNA microarray technology and bioinformatic web services

The pan-genomic microarray technique is used for environmental and/or clinical studies. Although microarray is an accurate and sharp diagnostic tool, the expertized bioinformaticians were able to minimize the outcome biases and maximize the flexibility and accuracy of the technique. The knowledge of bioinformatics plays a key role in association with probe designing and the utilization of correct probe sets and platforms. This technique is divided into two parts as dry lab (in silico studies) and wet lab (in vitro studies). Each part covers the other and are known as complementary divisions. In the case of microarray probe designing, a wide range of software, tools, and databases are necessary. Obviously, the application of right databases, software, and tools decreases the probable biases in the outcomes. Due to the importance of suitable probe designing, this article has focused its look onto a variety of online/offline databases, software, and tools.

A human gut phage catalog correlates the gut phageome with type 2 diabetes

BACKGROUND

Substantial efforts have been made to link the gut bacterial community to many complex human diseases. Nevertheless, the gut phages are often neglected.

RESULTS

In this study, we used multiple bioinformatic methods to catalog gut phages from whole-community metagenomic sequencing data of fecal samples collected from both type II diabetes (T2D) patients (n = 71) and normal Chinese adults (n = 74). The definition of phage operational taxonomic units (pOTUs) and identification of large phage scaffolds (n = 2567, ≥ 10 k) revealed a comprehensive human gut phageome with a substantial number of novel sequences encoding genes that were unrelated to those in known phages. Interestingly, we observed a significant increase in the number of gut phages in the T2D group and, in particular, identified 7 pOTUs specific to T2D. This finding was further validated in an independent dataset of 116 T2D and 109 control samples. Co-occurrence/exclusion analysis of the bacterial genera and pOTUs identified a complex core interaction between bacteria and phages in the human gut ecosystem, suggesting that the significant alterations of the gut phageome cannot be explained simply by co-variation with the altered bacterial hosts.

CONCLUSIONS

Alterations in the gut bacterial community have been linked to the chronic disease T2D, but the role of gut phages therein is not well understood. This is the first study to identify a T2D-specific gut phageome, indicating the existence of other mechanisms that might govern the gut phageome in T2D patients. These findings suggest the importance of the phageome in T2D risk, which warrants further investigation.

Transcriptomic profiling of Burkholderia phymatum STM815, Cupriavidus taiwanensis LMG19424 and Rhizobium mesoamericanum STM3625 in response to Mimosa pudica root exudates illuminates the molecular basis of their nodulation competitiveness and symbiotic evolutionary history

BACKGROUND

Rhizobial symbionts belong to the classes Alphaproteobacteria and Betaproteobacteria (called "alpha" and "beta"-rhizobia). Most knowledge on the genetic basis of symbiosis is based on model strains belonging to alpha-rhizobia. Mimosa pudica is a legume that offers an excellent opportunity to study the adaptation toward symbiotic nitrogen fixation in beta-rhizobia compared to alpha-rhizobia. In a previous study (Melkonian et al., Environ Microbiol 16:2099-111, 2014) we described the symbiotic competitiveness of M. pudica symbionts belonging to Burkholderia, Cupriavidus and Rhizobium species.

RESULTS

In this article we present a comparative analysis of the transcriptomes (by RNAseq) of B. phymatum STM815 (BP), C. taiwanensis LMG19424 (CT) and R. mesoamericanum STM3625 (RM) in conditions mimicking the early steps of symbiosis (i.e. perception of root exudates). BP exhibited the strongest transcriptome shift both quantitatively and qualitatively, which mirrors its high competitiveness in the early steps of symbiosis and its ancient evolutionary history as a symbiont, while CT had a minimal response which correlates with its status as a younger symbiont (probably via acquisition of symbiotic genes from a Burkholderia ancestor) and RM had a typical response of Alphaproteobacterial rhizospheric bacteria. Interestingly, the upregulation of nodulation genes was the only common response among the three strains; the exception was an up-regulated gene encoding a putative fatty acid hydroxylase, which appears to be a novel symbiotic gene specific to Mimosa symbionts.

CONCLUSION

The transcriptional response to root exudates was correlated to each strain nodulation competitiveness, with Burkholderia phymatum appearing as the best specialised symbiont of Mimosa pudica.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

The N-Acetylmuramic Acid 6-Phosphate Phosphatase MupP Completes the Pseudomonas Peptidoglycan Recycling Pathway Leading to Intrinsic Fosfomycin Resistance

Bacterial cells are encased in and stabilized by a netlike peptidoglycan (PGN) cell wall that undergoes turnover during bacterial growth. PGN turnover fragments are frequently salvaged by the cells via a pathway referred to as PGN recycling. Two different routes for the recycling of the cell wall sugar N-acetylmuramic acid (MurNAc) have been recognized in bacteria. In Escherichia coli and related enterobacteria, as well as in most Gram-positive bacteria, MurNAc is recovered via a catabolic route requiring a MurNAc 6-phosphate etherase (MurQ in E. coli) enzyme. However, many Gram-negative bacteria, including Pseudomonas species, lack a MurQ ortholog and use an alternative, anabolic recycling route that bypasses the de novo biosynthesis of uridyldiphosphate (UDP)-MurNAc, the first committed precursor of PGN. Bacteria featuring the latter pathway become intrinsically resistant to the antibiotic fosfomycin, which targets the de novo biosynthesis of UDP-MurNAc. We report here the identification and characterization of a phosphatase enzyme, named MupP, that had been predicted to complete the anabolic recycling pathway of Pseudomonas species but has remained unknown so far. It belongs to the large haloacid dehalogenase family of phosphatases and specifically converts MurNAc 6-phosphate to MurNAc. A ΔmupP mutant of Pseudomonas putida was highly susceptible to fosfomycin, accumulated large amounts of MurNAc 6-phosphate, and showed lower levels of UDP-MurNAc than wild-type cells, altogether consistent with a role for MupP in the anabolic PGN recycling route and as a determinant of intrinsic resistance to fosfomycin.

Many Gram-negative bacteria, but not E. coli, make use of a cell wall salvage pathway that contributes to the pool of UDP-MurNAc, the first committed precursor of cell wall synthesis in bacteria. This salvage pathway is of particular interest because it confers intrinsic resistance to the antibiotic fosfomycin, which blocks de novo UDP-MurNAc biosynthesis. Here we identified and characterized a previously missing enzyme within the salvage pathway, the MurNAc 6-phosphate phosphatase MupP of P. putida. MupP, together with the other enzymes of the anabolic recycling pathway, AnmK, AmgK, and MurU, yields UDP-MurNAc, renders bacteria intrinsically resistant to fosfomycin, and thus may serve as a novel drug target for antimicrobial therapy.

Functional characterization of the Helicobacter pylori chaperone protein HP0795

Trigger factor (TF) is one of the multiple bacterial chaperone proteins interacting with nascent peptides and facilitating their folding in bacteria. While TF is well-characterized in E. coli, HP0795, a TF-like homologue gene identified earlier in the pathogenic Helicobacter pylori (H. pylori), has not been studied biochemically to date. To characterize its function as a chaperone, we performed 3D-modeling, cross-linking and in vitro enzyme assays to HP0795 in vitro. Our results show that HP0795 possesses peptidyl prolyl cis-trans isomerase activity and exhibits a dimeric structure in solution. In addition, stable expression of HP0795 in a series of well-characterized E. coli chaperone-deficient strains rescued the growth defects in these mutants. Furthermore, we showed that the presence of HP0795 greatly reduced protein aggregation caused by deficiencies of chaperones in these strains. In contrast to other chaperone genes in H. pylori, gene expression of HP0795 displays little induction under acidic pH conditions. Together, our results suggest that HP0795 is a constitutively expressed TF-like protein of the prokaryotic chaperone family that may not play a major role in acid response. Given the pathogenic properties of H. pylori, our insights might provide new avenues for potential future medical intervention for H. pylori-related conditions.

AlgU Controls Expression of Virulence Genes in Pseudomonas syringae pv. tomato DC3000

Plant-pathogenic bacteria are able to integrate information about their environment and adjust gene expression to provide adaptive functions. AlgU, an extracytoplasmic function (ECF) sigma factor encoded by Pseudomonas syringae, controls expression of genes for alginate biosynthesis and genes involved with resisting osmotic and oxidative stress. AlgU is active while these bacteria are associated with plants, where its presence supports bacterial growth and disease symptoms. We found that AlgU is an important virulence factor for P. syringae pv. tomato DC3000 but that alginate production is dispensable for disease in host plants. This implies that AlgU regulates additional genes that facilitate bacterial pathogenesis. We used transcriptome sequencing (RNA-seq) to characterize the AlgU regulon and chromatin immunoprecipitation sequencing (ChIP-seq) to identify AlgU-regulated promoters associated with genes directly controlled by this sigma factor. We found that in addition to genes involved with alginate and osmotic and oxidative stress responses, AlgU regulates genes with known virulence functions, including components of the Hrp type III secretion system, virulence effectors, and the hrpL and hrpRS transcription regulators. These data suggest that P. syringae pv. tomato DC3000 has adapted to use signals that activate AlgU to induce expression of important virulence functions that facilitate survival and disease in plants.

IMPORTANCE

Plant immune systems produce antimicrobial and bacteriostatic conditions in response to bacterial infection. Plant-pathogenic bacteria are adapted to suppress and/or tolerate these conditions; however, the mechanisms controlling these bacterial systems are largely uncharacterized. The work presented here provides a mechanistic explanation for how P. syringae pv. tomato DC3000 coordinates expression of multiple genetic systems, including those dedicated to pathogenicity, in response to environmental conditions. This work demonstrates the scope of AlgU regulation in P. syringae pv. tomato DC3000 and characterizes the promoter sequence regulated by AlgU in these bacteria.

Recovering complete and draft population genomes from metagenome datasets

Assembly of metagenomic sequence data into microbial genomes is of fundamental value to improving our understanding of microbial ecology and metabolism by elucidating the functional potential of hard-to-culture microorganisms. Here, we provide a synthesis of available methods to bin metagenomic contigs into species-level groups and highlight how genetic diversity, sequencing depth, and coverage influence binning success. Despite the computational cost on application to deeply sequenced complex metagenomes (e.g., soil), covarying patterns of contig coverage across multiple datasets significantly improves the binning process. We also discuss and compare current genome validation methods and reveal how these methods tackle the problem of chimeric genome bins i.e., sequences from multiple species. Finally, we explore how population genome assembly can be used to uncover biogeographic trends and to characterize the effect of in situ functional constraints on the genome-wide evolution.

Bacterial Interactomes: Interacting Protein Partners Share Similar Function and Are Validated in Independent Assays More Frequently Than Previously Reported

Numerous affinity purification-mass spectrometry (AP-MS) and yeast two-hybrid screens have each defined thousands of pairwise protein-protein interactions (PPIs), most of which are between functionally unrelated proteins. The accuracy of these networks, however, is under debate. Here, we present an AP-MS survey of the bacterium Desulfovibrio vulgaris together with a critical reanalysis of nine published bacterial yeast two-hybrid and AP-MS screens. We have identified 459 high confidence PPIs from D. vulgaris and 391 from Escherichia coli Compared with the nine published interactomes, our two networks are smaller, are much less highly connected, and have significantly lower false discovery rates. In addition, our interactomes are much more enriched in protein pairs that are encoded in the same operon, have similar functions, and are reproducibly detected in other physical interaction assays than the pairs reported in prior studies. Our work establishes more stringent benchmarks for the properties of protein interactomes and suggests that bona fide PPIs much more frequently involve protein partners that are annotated with similar functions or that can be validated in independent assays than earlier studies suggested.

What Makes a Bacterial Species Pathogenic?:Comparative Genomic Analysis of the Genus Leptospira

Leptospirosis, caused by spirochetes of the genus Leptospira, is a globally widespread, neglected and emerging zoonotic disease. While whole genome analysis of individual pathogenic, intermediately pathogenic and saprophytic Leptospira species has been reported, comprehensive cross-species genomic comparison of all known species of infectious and non-infectious Leptospira, with the goal of identifying genes related to pathogenesis and mammalian host adaptation, remains a key gap in the field. Infectious Leptospira, comprised of pathogenic and intermediately pathogenic Leptospira, evolutionarily diverged from non-infectious, saprophytic Leptospira, as demonstrated by the following computational biology analyses: 1) the definitive taxonomy and evolutionary relatedness among all known Leptospira species; 2) genomically-predicted metabolic reconstructions that indicate novel adaptation of infectious Leptospira to mammals, including sialic acid biosynthesis, pathogen-specific porphyrin metabolism and the first-time demonstration of cobalamin (B12) autotrophy as a bacterial virulence factor; 3) CRISPR/Cas systems demonstrated only to be present in pathogenic Leptospira, suggesting a potential mechanism for this clade's refractoriness to gene targeting; 4) finding Leptospira pathogen-specific specialized protein secretion systems; 5) novel virulence-related genes/gene families such as the Virulence Modifying (VM) (PF07598 paralogs) proteins and pathogen-specific adhesins; 6) discovery of novel, pathogen-specific protein modification and secretion mechanisms including unique lipoprotein signal peptide motifs, Sec-independent twin arginine protein secretion motifs, and the absence of certain canonical signal recognition particle proteins from all Leptospira; and 7) and demonstration of infectious Leptospira-specific signal-responsive gene expression, motility and chemotaxis systems. By identifying large scale changes in infectious (pathogenic and intermediately pathogenic) vs. non-infectious Leptospira, this work provides new insights into the evolution of a genus of bacterial pathogens. This work will be a comprehensive roadmap for understanding leptospirosis pathogenesis. More generally, it provides new insights into mechanisms by which bacterial pathogens adapt to mammalian hosts.

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph--and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration--specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.

A novel method of consensus pan-chromosome assembly and large-scale comparative analysis reveal the highly flexible pan-genome of Acinetobacter baumannii

Background

Infections by pan-drug resistant Acinetobacter baumannii plague military and civilian healthcare systems. Previous A. baumannii pan-genomic studies used modest sample sizes of low diversity and comparisons to a single reference genome, limiting our understanding of gene order and content. A consensus representation of multiple genomes will provide a better framework for comparison. A large-scale comparative study will identify genomic determinants associated with their diversity and adaptation as a successful pathogen.

Results

We determine draft-level genomic sequence of 50 diverse military isolates and conduct the largest bacterial pan-genome analysis of 249 genomes. The pan-genome of A. baumannii is open when the input genomes are normalized for diversity with 1867 core proteins and a paralog-collapsed pan-genome size of 11,694 proteins. We developed a novel graph-based algorithm and use it to assemble the first consensus pan-chromosome, identifying both the order and orientation of core genes and flexible genomic regions. Comparative genome analyses demonstrate the existence of novel resistance islands and isolates with increased numbers of resistance island insertions over time, from single insertions in the 1950s to triple insertions in 2011. Gene clusters responsible for carbon utilization, siderophore production, and pilus assembly demonstrate frequent gain or loss among isolates.

Conclusions

The highly variable and dynamic nature of the A. baumannii genome may be the result of its success in rapidly adapting to both abiotic and biotic environments through the gain and loss of gene clusters controlling fitness. Importantly, some archaic adaptation mechanisms appear to have reemerged among recent isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0701-6) contains supplementary material, which is available to authorized users.

Databases for Microbiologists

Databases play an increasingly important role in biology. They archive, store, maintain, and share information on genes, genomes, expression data, protein sequences and structures, metabolites and reactions, interactions, and pathways. All these data are critically important to microbiologists. Furthermore, microbiology has its own databases that deal with model microorganisms, microbial diversity, physiology, and pathogenesis. Thousands of biological databases are currently available, and it becomes increasingly difficult to keep up with their development. The purpose of this minireview is to provide a brief survey of current databases that are of interest to microbiologists.

Genomic and transcriptomic analyses of colistin-resistant clinical isolates of Klebsiella pneumoniae reveal multiple pathways of resistance

The emergence of multidrug-resistant (MDR) Klebsiella pneumoniae has resulted in a more frequent reliance on treatment using colistin. However, resistance to colistin (Col(r)) is increasingly reported from clinical settings. The genetic mechanisms that lead to Col(r) in K. pneumoniae are not fully characterized. Using a combination of genome sequencing and transcriptional profiling by RNA sequencing (RNA-Seq) analysis, distinct genetic mechanisms were found among nine Col(r) clinical isolates. Col(r) was related to mutations in three different genes in K. pneumoniae strains, with distinct impacts on gene expression. Upregulation of the pmrH operon encoding 4-amino-4-deoxy-L-arabinose (Ara4N) modification of lipid A was found in all Col(r) strains. Alteration of the mgrB gene was observed in six strains. One strain had a mutation in phoQ. Common among these seven strains was elevated expression of phoPQ and unaltered expression of pmrCAB, which is involved in phosphoethanolamine addition to lipopolysaccharide (LPS). In two strains, separate mutations were found in a previously uncharacterized histidine kinase gene that is part of a two-component regulatory system (TCRS) now designated crrAB. In these strains, expression of pmrCAB, crrAB, and an adjacent glycosyltransferase gene, but not that of phoPQ, was elevated. Complementation with the wild-type allele restored colistin susceptibility in both strains. The crrAB genes are present in most K. pneumoniae genomes, but not in Escherichia coli. Additional upregulated genes in all strains include those involved in cation transport and maintenance of membrane integrity. Because the crrAB genes are present in only some strains, Col(r) mechanisms may be dependent on the genetic background.

Genome-wide analyses in bacteria show small-RNA enrichment for long and conserved intergenic regions

Interest in finding small RNAs (sRNAs) in bacteria has significantly increased in recent years due to their regulatory functions. Development of high-throughput methods and more sophisticated computational algorithms has allowed rapid identification of sRNA candidates in different species. However, given their various sizes (50 to 500 nucleotides [nt]) and their potential genomic locations in the 5' and 3' untranslated regions as well as in intergenic regions, identification and validation of true sRNAs have been challenging. In addition, the evolution of bacterial sRNAs across different species continues to be puzzling, given that they can exert similar functions with various sequences and structures. In this study, we analyzed the enrichment patterns of sRNAs in 13 well-annotated bacterial species using existing transcriptome and experimental data. All intergenic regions were analyzed by WU-BLAST to examine conservation levels relative to species within or outside their genus. In total, more than 900 validated bacterial sRNAs and 23,000 intergenic regions were analyzed. The results indicate that sRNAs are enriched in intergenic regions, which are longer and more conserved than the average intergenic regions in the corresponding bacterial genome. We also found that sRNA-coding regions have different conservation levels relative to their flanking regions. This work provides a way to analyze how noncoding RNAs are distributed in bacterial genomes and also shows conserved features of intergenic regions that encode sRNAs. These results also provide insight into the functions of regions surrounding sRNAs and into optimization of RNA search algorithms.

Insyght: navigating amongst abundant homologues, syntenies and gene functional annotations in bacteria, it's that symbol!

High-throughput techniques have considerably increased the potential of comparative genomics whilst simultaneously posing many new challenges. One of those challenges involves efficiently mining the large amount of data produced and exploring the landscape of both conserved and idiosyncratic genomic regions across multiple genomes. Domains of application of these analyses are diverse: identification of evolutionary events, inference of gene functions, detection of niche-specific genes or phylogenetic profiling. Insyght is a comparative genomic visualization tool that combines three complementary displays: (i) a table for thoroughly browsing amongst homologues, (ii) a comparator of orthologue functional annotations and (iii) a genomic organization view designed to improve the legibility of rearrangements and distinctive loci. The latter display combines symbolic and proportional graphical paradigms. Synchronized navigation across multiple species and interoperability between the views are core features of Insyght. A gene filter mechanism is provided that helps the user to build a biologically relevant gene set according to multiple criteria such as presence/absence of homologues and/or various annotations. We illustrate the use of Insyght with scenarios. Currently, only Bacteria and Archaea are supported. A public instance is available at http://genome.jouy.inra.fr/Insyght. The tool is freely downloadable for private data set analysis.

Inflammation-associated adherent-invasive Escherichia coli are enriched in pathways for use of propanediol and iron and M-cell translocation

BACKGROUND

Perturbations of the intestinal microbiome, termed dysbiosis, are linked to intestinal inflammation. Isolation of adherent-invasive Escherichia coli (AIEC) from intestines of patients with Crohn's disease (CD), dogs with granulomatous colitis, and mice with acute ileitis suggests these bacteria share pathoadaptive virulence factors that promote inflammation.

METHODS

To identify genes associated with AIEC, we sequenced the genomes of phylogenetically diverse AIEC strains isolated from people with CD (4), dogs with granulomatous colitis (2), and mice with ileitis (2) and 1 non-AIEC strain from CD ileum and compared them with 38 genome sequences of E. coli and Shigella. We then determined the prevalence of AIEC-associated genes in 49 E. coli strains from patients with CD and controls and correlated genotype with invasion of intestinal epithelial cells, persistence within macrophages, AIEC pathotype, and growth in standardized conditions.

RESULTS

Genes encoding propanediol utilization (pdu operon) and iron acquisition (yersiniabactin, chu operon) were overrepresented in AIEC relative to nonpathogenic E. coli. PduC (propanediol dehydratase) was enriched in CD-derived AIEC, correlated with increased cellular invasion, and persistence in vitro and was increasingly expressed in fucose-containing media. Growth of AIEC required iron, and the presence of chuA (heme acquisition) correlated with persistence in macrophages. CD-associated AIEC with lpfA 154 (long polar fimbriae) demonstrated increased invasion of epithelial cells and translocation across M cells.

CONCLUSIONS

Our findings provide novel insights into the genetic basis of the AIEC pathotype, supporting the concept that AIEC are equipped to exploit and promote intestinal inflammation and reveal potential targets for intervention against AIEC and inflammation-associated dysbiosis.

Transcriptome analysis of Escherichia coli O157:H7 grown in vitro in the sterile-filtrated cecal content of human gut microbiota associated rats reveals an adaptive expression of metabolic and virulence genes

In developed countries, enterohemorrhagic Escherichia coli (EHEC) O157:H7 is a leading cause of bloody diarrhea and renal failures in human. Understanding strategies employed by EHEC to colonize the intestine is of major importance since to date no cure exists to eradicate the pathogen. In this study, the adaptive response of EHEC to the intestinal milieu conditioned by a human microbiota was examined. A transcriptomic analysis was performed on the EHEC strain EDL933 incubated in vitro in the sterile-filtrated cecal content of human microbiota-associated rats (HMC) compared with EDL933 incubated in the sterile-filtrated cecal content of germ-free rat (GFC). EDL933 switches from a glycolytic metabolic profile in the GFC to an anaplerotic metabolic profile in HMC. The expression of several catabolism genes was strongly affected such as those involved in the utilization of sugars, glycerol, N-acetylneuraminic acid, amino acids and secondary metabolites. Interestingly, expression level of critical EHEC O157:H7 virulence genes including genes from the locus of enterocyte effacement was reduced in HMC. Altogether, these results contribute to the understanding of EHEC adaptive response to a digestive content and highlight the ability of the microbiota to repress EHEC virulence gene expression.

BB0238, a presumed tetratricopeptide repeat-containing protein, is required during Borrelia burgdorferi mammalian infection

The Lyme disease spirochete, Borrelia burgdorferi, occupies both a tick vector and mammalian host in nature. Considering the unique enzootic life cycle of B. burgdorferi, it is not surprising that a large proportion of its genome is composed of hypothetical proteins not found in other bacterial pathogens. bb0238 encodes a conserved hypothetical protein of unknown function that is predicted to contain a tetratricopeptide repeat (TPR) domain, a structural motif responsible for mediating protein-protein interactions. To evaluate the role of bb0238 during mammalian infection, a bb0238-deficient mutant was constructed. The bb0238 mutant was attenuated in mice infected via needle inoculation, and complementation of bb0238 expression restored infectivity to wild-type levels. bb0238 expression does not change in response to varying culture conditions, and thus, it appears to be constitutively expressed under in vitro conditions. bb0238 is expressed in murine tissues during infection, though there was no significant change in expression levels among different tissue types. Localization studies indicate that BB0238 is associated with the inner membrane of the spirochete and is therefore unlikely to promote interaction with host ligands during infection. B. burgdorferi clones containing point mutations in conserved residues of the putative TPR motif of BB0238 demonstrated attenuation in mice that was comparable to that in the bb0238 deletion mutant, suggesting that BB0238 may contain a functional TPR domain.

Variations in metabolic pathways create challenges for automated metabolic reconstructions: Examples from the tetrahydrofolate synthesis pathway

The availability of thousands of sequenced genomes has revealed the diversity of biochemical solutions to similar chemical problems. Even for molecules at the heart of metabolism, such as cofactors, the pathway enzymes first discovered in model organisms like Escherichia coli or Saccharomyces cerevisiae are often not universally conserved. Tetrahydrofolate (THF) (or its close relative tetrahydromethanopterin) is a universal and essential C1-carrier that most microbes and plants synthesize de novo. The THF biosynthesis pathway and enzymes are, however, not universal and alternate solutions are found for most steps, making this pathway a challenge to annotate automatically in many genomes. Comparing THF pathway reconstructions and functional annotations of a chosen set of folate synthesis genes in specific prokaryotes revealed the strengths and weaknesses of different microbial annotation platforms. This analysis revealed that most current platforms fail in metabolic reconstruction of variant pathways. However, all the pieces are in place to quickly correct these deficiencies if the different databases were built on each other's strengths.

Comparative functional genomics of Lactobacillus spp. reveals possible mechanisms for specialization of vaginal lactobacilli to their environment

Lactobacilli are found in a wide variety of habitats. Four species, Lactobacillus crispatus, L. gasseri, L. iners, and L. jensenii, are common and abundant in the human vagina and absent from other habitats. These may be adapted to the vagina and possess characteristics enabling them to thrive in that environment. Furthermore, stable codominance of multiple Lactobacillus species in a single community is infrequently observed. Thus, it is possible that individual vaginal Lactobacillus species possess unique characteristics that confer to them host-specific competitive advantages. We performed comparative functional genomic analyses of representatives of 25 species of Lactobacillus, searching for habitat-specific traits in the genomes of the vaginal lactobacilli. We found that the genomes of the vaginal species were significantly smaller and had significantly lower GC content than those of the nonvaginal species. No protein families were found to be specific to the vaginal species analyzed, but some were either over- or underrepresented relative to nonvaginal species. We also found that within the vaginal species, each genome coded for species-specific protein families. Our results suggest that even though the vaginal species show no general signatures of adaptation to the vaginal environment, each species has specific and perhaps unique ways of interacting with its environment, be it the host or other microbes in the community. These findings will serve as a foundation for further exploring the role of lactobacilli in the ecological dynamics of vaginal microbial communities and their ultimate impact on host health.

A second-generation protein-protein interaction network of Helicobacter pylori

Helicobacter pylori infections cause gastric ulcers and play a major role in the development of gastric cancer. In 2001, the first protein interactome was published for this species, revealing over 1500 binary protein interactions resulting from 261 yeast two-hybrid screens. Here we roughly double the number of previously published interactions using an ORFeome-based, proteome-wide yeast two-hybrid screening strategy. We identified a total of 1515 protein-protein interactions, of which 1461 are new. The integration of all the interactions reported in H. pylori results in 3004 unique interactions that connect about 70% of its proteome. Excluding interactions of promiscuous proteins we derived from our new data a core network consisting of 908 interactions. We compared our data set to several other bacterial interactomes and experimentally benchmarked the conservation of interactions using 365 protein pairs (interologs) of E. coli of which one third turned out to be conserved in both species.

New insights into dissemination and variation of the health care-associated pathogen Acinetobacter baumannii from genomic analysis

Acinetobacter baumannii is a globally important nosocomial pathogen characterized by an increasing incidence of multidrug resistance. Routes of dissemination and gene flow among health care facilities are poorly resolved and are important for understanding the epidemiology of A. baumannii, minimizing disease transmission, and improving patient outcomes. We used whole-genome sequencing to assess diversity and genome dynamics in 49 isolates from one United States hospital system during one year from 2007 to 2008. Core single-nucleotide-variant-based phylogenetic analysis revealed multiple founder strains and multiple independent strains recovered from the same patient yet was insufficient to fully resolve strain relationships, where gene content and insertion sequence patterns added additional discriminatory power. Gene content comparisons illustrated extensive and redundant antibiotic resistance gene carriage and direct evidence of gene transfer, recombination, gene loss, and mutation. Evidence of barriers to gene flow among hospital components was not found, suggesting complex mixing of strains and a large reservoir of A. baumannii strains capable of colonizing patients.

Genome sequencing was used to characterize multidrug-resistant Acinetobacter baumannii strains from one United States hospital system during a 1-year period to better understand how A. baumannii strains that cause infection are related to one another. Extensive variation in gene content was found, even among strains that were very closely related phylogenetically and epidemiologically. Several mechanisms contributed to this diversity, including transfer of mobile genetic elements, mobilization of insertion sequences, insertion sequence-mediated deletions, and genome-wide homologous recombination. Variation in gene content, however, lacked clear spatial or temporal patterns, suggesting a diverse pool of circulating strains with considerable interaction between strains and hospital locations. Widespread genetic variation among strains from the same hospital and even the same patient, particularly involving antibiotic resistance genes, reinforces the need for molecular diagnostic testing and genomic analysis to determine resistance profiles, rather than a reliance primarily on strain typing and antimicrobial resistance phenotypes for epidemiological studies.

A taxonomic signature of obesity in the microbiome? Getting to the guts of the matter

Obesity is an important and intractable public health problem. In addition to the well-known risk factors of behavior, diet, and genetics, gut microbial communities were recently identified as another possible source of risk and a potential therapeutic target. However, human and animal-model studies have yielded conflicting results about the precise nature of associations between microbiome composition and obesity. In this paper, we use publicly available data from the Human Microbiome Project (HMP) and MetaHIT, both surveys of healthy adults that include obese individuals, plus two smaller studies that specifically examined lean versus obese adults. We find that inter-study variability in the taxonomic composition of stool microbiomes far exceeds differences between lean and obese individuals within studies. Our analyses further reveal a high degree of variability in stool microbiome composition and diversity across individuals. While we confirm the previously published small, but statistically significant, differences in phylum-level taxonomic composition between lean and obese individuals in several cohorts, we find no association between BMI and taxonomic composition of stool microbiomes in the larger HMP and MetaHIT datasets. We explore a range of different statistical techniques and show that this result is robust to the choice of methodology. Differences between studies are likely due to a combination of technical and clinical factors. We conclude that there is no simple taxonomic signature of obesity in the microbiota of the human gut.

Network portal: a database for storage, analysis and visualization of biological networks

The ease of generating high-throughput data has enabled investigations into organismal complexity at the systems level through the inference of networks of interactions among the various cellular components (genes, RNAs, proteins and metabolites). The wider scientific community, however, currently has limited access to tools for network inference, visualization and analysis because these tasks often require advanced computational knowledge and expensive computing resources. We have designed the network portal (http://networks.systemsbiology.net) to serve as a modular database for the integration of user uploaded and public data, with inference algorithms and tools for the storage, visualization and analysis of biological networks. The portal is fully integrated into the Gaggle framework to seamlessly exchange data with desktop and web applications and to allow the user to create, save and modify workspaces, and it includes social networking capabilities for collaborative projects. While the current release of the database contains networks for 13 prokaryotic organisms from diverse phylogenetic clades (4678 co-regulated gene modules, 3466 regulators and 9291 cis-regulatory motifs), it will be rapidly populated with prokaryotic and eukaryotic organisms as relevant data become available in public repositories and through user input. The modular architecture, simple data formats and open API support community development of the portal.

Analysis of Phakopsora pachyrhizi transcript abundance in critical pathways at four time-points during infection of a susceptible soybean cultivar using deep sequencing

BACKGROUND

Phakopsora pachyrhizi, the causal agent responsible for soybean rust, is among the top hundred most virulent plant pathogens and can cause soybean yield losses of up to 80% when appropriate conditions are met. We used mRNA-Seq by Illumina to analyze pathogen transcript abundance at 15 seconds (s), 7 hours (h), 48 h, and 10 days (d) after inoculation (ai) of susceptible soybean leaves with P. pachyrhizi to gain new insights into transcript abundance in soybean and the pathogen at specific time-points during the infection including the uredinial stage.

RESULTS

Over three million five hundred thousand sequences were obtained for each time-point. Energy, nucleotide metabolism, and protein synthesis are major priorities for the fungus during infection and development as indicated by our transcript abundance studies. At all time-points, energy production is a necessity for P. pachyrhizi, as indicated by expression of many transcripts encoding enzymes involved in oxidative phosphorylation and carbohydrate metabolism (glycolysis, glyoxylate and dicarboxylate, pentose phosphate, pyruvate). However, at 15 sai, transcripts encoding enzymes involved in ATP production were highly abundant in order to provide enough energy for the spore to germinate, as observed by the expression of many transcripts encoding proteins involved in electron transport. At this early time-point, transcripts encoding proteins involved in RNA synthesis were also highly abundant, more so than transcripts encoding genes involved in DNA and protein synthesis. At 7 hai, shortly after germination during tube elongation and penetration, transcripts encoding enzymes involved in deoxyribonucleotide and DNA synthesis were highly abundant. At 48 hai, transcripts encoding enzymes involved in amino acid metabolism were highly abundant to provide for increased protein synthesis during haustoria maturation. During sporulation at 10 dai, the fungus still required carbohydrate metabolism, but there also was increased expression of transcripts encoding enzymes involved in fatty acid metabolism.

CONCLUSION

This information provides insight into molecular events and their timing throughout the life cycle of the P. pachyrhizi, and it may be useful in the development of new methods of broadening resistance of soybean to soybean rust.

Aureolib - a proteome signature library: towards an understanding of staphylococcus aureus pathophysiology

Gel-based proteomics is a powerful approach to study the physiology of Staphylococcus aureus under various growth restricting conditions. We analyzed 679 protein spots from a reference 2-dimensional gel of cytosolic proteins of S. aureus COL by mass spectrometry resulting in 521 different proteins. 4,692 time dependent protein synthesis profiles were generated by exposing S. aureus to nine infection-related stress and starvation stimuli (H2O2, diamide, paraquat, NO, fermentation, nitrate respiration, heat shock, puromycin, mupirocin). These expression profiles are stored in an online resource called Aureolib (http://www.aureolib.de). Moreover, information on target genes of 75 regulators and regulatory elements were included in the database. Cross-comparisons of this extensive data collection of protein synthesis profiles using the tools implemented in Aureolib lead to the identification of stress and starvation specific marker proteins. Altogether, 226 protein synthesis profiles showed induction ratios of 2.5-fold or higher under at least one of the tested conditions with 157 protein synthesis profiles specifically induced in response to a single stimulus. The respective proteins might serve as marker proteins for the corresponding stimulus. By contrast, proteins whose synthesis was increased or repressed in response to more than four stimuli are rather exceptional. The only protein that was induced by six stimuli is the universal stress protein SACOL1759. Most strikingly, cluster analyses of synthesis profiles of proteins differentially synthesized under at least one condition revealed only in rare cases a grouping that correlated with known regulon structures. The most prominent examples are the GapR, Rex, and CtsR regulon. In contrast, protein synthesis profiles of proteins belonging to the CodY and σ(B) regulon are widely distributed. In summary, Aureolib is by far the most comprehensive protein expression database for S. aureus and provides an essential tool to decipher more complex adaptation processes in S. aureus during host pathogen interaction.

JBioWH: an open-source Java framework for bioinformatics data integration

The Java BioWareHouse (JBioWH) project is an open-source platform-independent programming framework that allows a user to build his/her own integrated database from the most popular data sources. JBioWH can be used for intensive querying of multiple data sources and the creation of streamlined task-specific data sets on local PCs. JBioWH is based on a MySQL relational database scheme and includes JAVA API parser functions for retrieving data from 20 public databases (e.g. NCBI, KEGG, etc.). It also includes a client desktop application for (non-programmer) users to query data. In addition, JBioWH can be tailored for use in specific circumstances, including the handling of massive queries for high-throughput analyses or CPU intensive calculations. The framework is provided with complete documentation and application examples and it can be downloaded from the Project Web site at http://code.google.com/p/jbiowh. A MySQL server is available for demonstration purposes at hydrax.icgeb.trieste.it:3307. Database URL: http://code.google.com/p/jbiowh.

Candidate genes that may be responsible for the unusual resistances exhibited by Bacillus pumilus SAFR-032 spores

The spores of several Bacillus species, including Bacillus pumilus SAFR-032 and B. safensis FO-36b, which were isolated from the spacecraft assembly facility at NASA's Jet Propulsion Laboratory, are unusually resistant to UV radiation and hydrogen peroxide. In order to identify candidate genes that might be associated with these resistances, the whole genome of B. pumilus SAFR-032, and the draft genome of B. safensis FO-36b were compared in detail with the very closely related type strain B. pumilus ATCC7061(T). 170 genes are considered characteristic of SAFR-032, because they are absent from both FO-36b and ATCC7061(T). Forty of these SAFR-032 characteristic genes are entirely unique open reading frames. In addition, four genes are unique to the genomes of the resistant SAFR-032 and FO-36b. Fifty three genes involved in spore coat formation, regulation and germination, DNA repair, and peroxide resistance, are missing from all three genomes. The vast majority of these are cleanly deleted from their usual genomic context without any obvious replacement. Several DNA repair and peroxide resistance genes earlier reported to be unique to SAFR-032 are in fact shared with ATCC7061(T) and no longer considered to be promising candidates for association with the elevated resistances. Instead, several SAFR-032 characteristic genes were identified, which along with one or more of the unique SAFR-032 genes may be responsible for the elevated resistances. These new candidates include five genes associated with DNA repair, namely, BPUM_0608 a helicase, BPUM_0652 an ATP binding protein, BPUM_0653 an endonuclease, BPUM_0656 a DNA cytosine-5- methyltransferase, and BPUM_3674 a DNA helicase. Three of these candidate genes are in immediate proximity of two conserved hypothetical proteins, BPUM_0654 and BPUM_0655 that are also absent from both FO-36b and ATCC7061(T). This cluster of five genes is considered to be an especially promising target for future experimental work.

New insights on the reorganization of gene transcription in Pseudomonas putida KT2440 at elevated pressure

Background

Elevated pressure, elevated oxygen tension (DOT) and elevated carbon dioxide tension (DCT) are readily encountered at the bottom of large industrial bioreactors and during bioprocesses where pressure is applied for enhancing the oxygen transfer. Yet information about their effect on bacteria and on the gene expression thereof is scarce. To shed light on the cellular functions affected by these specific environmental conditions, the transcriptome of Pseudomonas putida KT2440, a bacterium of great relevance for the production of medium-chain-length polyhydroxyalkanoates, was thoroughly investigated using DNA microarrays.

Results

Very well defined chemostat cultivations were carried out with P. putida to produce high quality RNA samples and ensure that differential gene expression was caused exclusively by changes of pressure, DOT and/or DCT. Cellular stress was detected at 7 bar and elevated DCT in the form of heat shock and oxidative stress-like responses, and indicators of cell envelope perturbations were identified as well.

Globally, gene transcription was not considerably altered when DOT was increased from 40 ± 5 to 235 ± 20% at 7 bar and elevated DCT. Nevertheless, differential transcription was observed for a few genes linked to iron-sulfur cluster assembly, terminal oxidases, glutamate metabolism and arginine deiminase pathway, which shows their particular sensitivity to variations of DOT.

Conclusions

This study provides a comprehensive overview on the changes occurring in the transcriptome of P. putida upon mild variations of pressure, DOT and DCT. Interestingly, whereas the changes of gene transcription were widespread, the cell physiology was hardly affected, which illustrates how efficient reorganization of the gene transcription is for dealing with environmental changes that may otherwise be harmful. Several particularly sensitive cellular functions were identified, which will certainly contribute to the understanding of the mechanisms involved in stress sensing/response and to finding ways of enhancing the stress tolerance of microorganisms.

WholeCellKB: model organism databases for comprehensive whole-cell models

Whole-cell models promise to greatly facilitate the analysis of complex biological behaviors. Whole-cell model development requires comprehensive model organism databases. WholeCellKB (http://wholecellkb.stanford.edu) is an open-source web-based software program for constructing model organism databases. WholeCellKB provides an extensive and fully customizable data model that fully describes individual species including the structure and function of each gene, protein, reaction and pathway. We used WholeCellKB to create WholeCellKB-MG, a comprehensive database of the Gram-positive bacterium Mycoplasma genitalium using over 900 sources. WholeCellKB-MG is extensively cross-referenced to existing resources including BioCyc, KEGG and UniProt. WholeCellKB-MG is freely accessible through a web-based user interface as well as through a RESTful web service.

OrtholugeDB: a bacterial and archaeal orthology resource for improved comparative genomic analysis

Prediction of orthologs (homologous genes that diverged because of speciation) is an integral component of many comparative genomics methods. Although orthologs are more likely to have similar function versus paralogs (genes that diverged because of duplication), recent studies have shown that their degree of functional conservation is variable. Also, there are inherent problems with several large-scale ortholog prediction approaches. To address these issues, we previously developed Ortholuge, which uses phylogenetic distance ratios to provide more precise ortholog assessments for a set of predicted orthologs. However, the original version of Ortholuge required manual intervention and was not easily accessible; therefore, we now report the development of OrtholugeDB, available online at http://www.pathogenomics.sfu.ca/ortholugedb. OrtholugeDB provides ortholog predictions for completely sequenced bacterial and archaeal genomes from NCBI based on reciprocal best Basic Local Alignment Search Tool hits, supplemented with further evaluation by the more precise Ortholuge method. The OrtholugeDB web interface facilitates user-friendly and flexible ortholog analysis, from single genes to genomes, plus flexible data download options. We compare Ortholuge with similar methods, showing how it may more consistently identify orthologs with conserved features across a wide range of taxonomic distances. OrtholugeDB facilitates rapid, and more accurate, bacterial and archaeal comparative genomic analysis and large-scale ortholog predictions.

TIGRFAMs and Genome Properties in 2013

TIGRFAMs, available online at http://www.jcvi.org/tigrfams is a database of protein family definitions. Each entry features a seed alignment of trusted representative sequences, a hidden Markov model (HMM) built from that alignment, cutoff scores that let automated annotation pipelines decide which proteins are members, and annotations for transfer onto member proteins. Most TIGRFAMs models are designated equivalog, meaning they assign a specific name to proteins conserved in function from a common ancestral sequence. Models describing more functionally heterogeneous families are designated subfamily or domain, and assign less specific but more widely applicable annotations. The Genome Properties database, available at http://www.jcvi.org/genome-properties, specifies how computed evidence, including TIGRFAMs HMM results, should be used to judge whether an enzymatic pathway, a protein complex or another type of molecular subsystem is encoded in a genome. TIGRFAMs and Genome Properties content are developed in concert because subsystems reconstruction for large numbers of genomes guides selection of seed alignment sequences and cutoff values during protein family construction. Both databases specialize heavily in bacterial and archaeal subsystems. At present, 4284 models appear in TIGRFAMs, while 628 systems are described by Genome Properties. Content derives both from subsystem discovery work and from biocuration of the scientific literature.

Bioinformatics Resource Manager v2.3: an integrated software environment for systems biology with microRNA and cross-species analysis tools

BACKGROUND

MicroRNAs (miRNAs) are noncoding RNAs that direct post-transcriptional regulation of protein coding genes. Recent studies have shown miRNAs are important for controlling many biological processes, including nervous system development, and are highly conserved across species. Given their importance, computational tools are necessary for analysis, interpretation and integration of high-throughput (HTP) miRNA data in an increasing number of model species. The Bioinformatics Resource Manager (BRM) v2.3 is a software environment for data management, mining, integration and functional annotation of HTP biological data. In this study, we report recent updates to BRM for miRNA data analysis and cross-species comparisons across datasets.

RESULTS

BRM v2.3 has the capability to query predicted miRNA targets from multiple databases, retrieve potential regulatory miRNAs for known genes, integrate experimentally derived miRNA and mRNA datasets, perform ortholog mapping across species, and retrieve annotation and cross-reference identifiers for an expanded number of species. Here we use BRM to show that developmental exposure of zebrafish to 30 uM nicotine from 6-48 hours post fertilization (hpf) results in behavioral hyperactivity in larval zebrafish and alteration of putative miRNA gene targets in whole embryos at developmental stages that encompass early neurogenesis. We show typical workflows for using BRM to integrate experimental zebrafish miRNA and mRNA microarray datasets with example retrievals for zebrafish, including pathway annotation and mapping to human ortholog. Functional analysis of differentially regulated (p<0.05) gene targets in BRM indicates that nicotine exposure disrupts genes involved in neurogenesis, possibly through misregulation of nicotine-sensitive miRNAs.

CONCLUSIONS

BRM provides the ability to mine complex data for identification of candidate miRNAs or pathways that drive phenotypic outcome and, therefore, is a useful hypothesis generation tool for systems biology. The miRNA workflow in BRM allows for efficient processing of multiple miRNA and mRNA datasets in a single software environment with the added capability to interact with public data sources and visual analytic tools for HTP data analysis at a systems level. BRM is developed using Java™ and other open-source technologies for free distribution (http://www.sysbio.org/dataresources/brm.stm).

MBGD update 2013: the microbial genome database for exploring the diversity of microbial world

The microbial genome database for comparative analysis (MBGD, available at http://mbgd.genome.ad.jp/) is a platform for microbial genome comparison based on orthology analysis. As its unique feature, MBGD allows users to conduct orthology analysis among any specified set of organisms; this flexibility allows MBGD to adapt to a variety of microbial genomic study. Reflecting the huge diversity of microbial world, the number of microbial genome projects now becomes several thousands. To efficiently explore the diversity of the entire microbial genomic data, MBGD now provides summary pages for pre-calculated ortholog tables among various taxonomic groups. For some closely related taxa, MBGD also provides the conserved synteny information (core genome alignment) pre-calculated using the CoreAligner program. In addition, efficient incremental updating procedure can create extended ortholog table by adding additional genomes to the default ortholog table generated from the representative set of genomes. Combining with the functionalities of the dynamic orthology calculation of any specified set of organisms, MBGD is an efficient and flexible tool for exploring the microbial genome diversity.

Genetic structure of community acquired methicillin-resistant Staphylococcus aureus USA300

Background

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) is a significant bacterial pathogen that poses considerable clinical and public health challenges. The majority of the CA-MRSA disease burden consists of skin and soft tissue infections (SSTI) not associated with significant morbidity; however, CA-MRSA also causes severe, invasive infections resulting in significant morbidity and mortality. The broad range of disease severity may be influenced by bacterial genetic variation.

Results

We sequenced the complete genomes of 36 CA-MRSA clinical isolates from the predominant North American community acquired clonal type USA300 (18 SSTI and 18 severe infection-associated isolates). While all 36 isolates shared remarkable genetic similarity, we found greater overall time-dependent sequence diversity among SSTI isolates. In addition, pathway analysis of non-synonymous variations revealed increased sequence diversity in the putative virulence genes of SSTI isolates.

Conclusions

Here we report the first whole genome survey of diverse clinical isolates of the USA300 lineage and describe the evolution of the pathogen over time within a defined geographic area. The results demonstrate the close relatedness of clinically independent CA-MRSA isolates, which carry implications for understanding CA-MRSA epidemiology and combating its spread.

PanOCT: automated clustering of orthologs using conserved gene neighborhood for pan-genomic analysis of bacterial strains and closely related species

Pan-genome ortholog clustering tool (PanOCT) is a tool for pan-genomic analysis of closely related prokaryotic species or strains. PanOCT uses conserved gene neighborhood information to separate recently diverged paralogs into orthologous clusters where homology-only clustering methods cannot. The results from PanOCT and three commonly used graph-based ortholog-finding programs were compared using a set of four publicly available strains of the same bacterial species. All four methods agreed on ∼70% of the clusters and ∼86% of the proteins. The clusters that did not agree were inspected for evidence of correctness resulting in 85 high-confidence manually curated clusters that were used to compare all four methods.

MicrobeDB: a locally maintainable database of microbial genomic sequences

Summary: Analysis of microbial genomes often requires the general organization and comparison of tens to thousands of genomes both from public repositories and unpublished sources. MicrobeDB provides a foundation for such projects by the automation of downloading published, completed bacterial and archaeal genomes from key sources, parsing annotations of all genomes (both public and private) into a local database, and allowing interaction with the database through an easy to use programming interface. MicrobeDB creates a simple to use, easy to maintain, centralized local resource for various large-scale comparative genomic analyses and a back-end for future microbial application design.

Availability: MicrobeDB is freely available under the GNU-GPL at: http://github.com/mlangill/microbedb/

Contact:morgan.g.i.langille@gmail.com

A Gene-By-Gene Approach to Bacterial Population Genomics: Whole Genome MLST of Campylobacter

Campylobacteriosis remains a major human public health problem world-wide. Genetic analyses of Campylobacter isolates, and particularly molecular epidemiology, have been central to the study of this disease, particularly the characterization of Campylobacter genotypes isolated from human infection, farm animals, and retail food. These studies have demonstrated that Campylobacter populations are highly structured, with distinct genotypes associated with particular wild or domestic animal sources, and that chicken meat is the most likely source of most human infection in countries such as the UK. The availability of multiple whole genome sequences from Campylobacter isolates presents the prospect of identifying those genes or allelic variants responsible for host-association and increased human disease risk, but the diversity of Campylobacter genomes present challenges for such analyses. We present a gene-by-gene approach for investigating the genetic basis of phenotypes in diverse bacteria such as Campylobacter, implemented with the BIGSDB software on the pubMLST.org/campylobacter website.

Comparative genomic analysis of the genus Staphylococcus including Staphylococcus aureus and its newly described sister species Staphylococcus simiae

Background

Staphylococcus belongs to the Gram-positive low G + C content group of the Firmicutes division of bacteria. Staphylococcus aureus is an important human and veterinary pathogen that causes a broad spectrum of diseases, and has developed important multidrug resistant forms such as methicillin-resistant S. aureus (MRSA). Staphylococcus simiae was isolated from South American squirrel monkeys in 2000, and is a coagulase-negative bacterium, closely related, and possibly the sister group, to S. aureus. Comparative genomic analyses of closely related bacteria with different phenotypes can provide information relevant to understanding adaptation to host environment and mechanisms of pathogenicity.

Results

We determined a Roche/454 draft genome sequence for S. simiae and included it in comparative genomic analyses with 11 other Staphylococcus species including S. aureus. A genome based phylogeny of the genus confirms that S. simiae is the sister group to S. aureus and indicates that the most basal Staphylococcus lineage is Staphylococcus pseudintermedius, followed by Staphylococcus carnosus. Given the primary niche of these two latter taxa, compared to the other species in the genus, this phylogeny suggests that human adaptation evolved after the split of S. carnosus. The two coagulase-positive species (S. aureus and S. pseudintermedius) are not phylogenetically closest but share many virulence factors exclusively, suggesting that these genes were acquired by horizontal transfer. Enrichment in genes related to mobile elements such as prophage in S. aureus relative to S. simiae suggests that pathogenesis in the S. aureus group has developed by gene gain through horizontal transfer, after the split of S. aureus and S. simiae from their common ancestor.

Conclusions

Comparative genomic analyses across 12 Staphylococcus species provide hypotheses about lineages in which human adaptation has taken place and contributions of horizontal transfer in pathogenesis.

Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage

Bacteria in the 16S rRNA clade SAR86 are among the most abundant uncultivated constituents of microbial assemblages in the surface ocean for which little genomic information is currently available. Bioinformatic techniques were used to assemble two nearly complete genomes from marine metagenomes and single-cell sequencing provided two more partial genomes. Recruitment of metagenomic data shows that these SAR86 genomes substantially increase our knowledge of non-photosynthetic bacteria in the surface ocean. Phylogenomic analyses establish SAR86 as a basal and divergent lineage of γ-proteobacteria, and the individual genomes display a temperature-dependent distribution. Modestly sized at 1.25-1.7 Mbp, the SAR86 genomes lack several pathways for amino-acid and vitamin synthesis as well as sulfate reduction, trends commonly observed in other abundant marine microbes. SAR86 appears to be an aerobic chemoheterotroph with the potential for proteorhodopsin-based ATP generation, though the apparent lack of a retinal biosynthesis pathway may require it to scavenge exogenously-derived pigments to utilize proteorhodopsin. The genomes contain an expanded capacity for the degradation of lipids and carbohydrates acquired using a wealth of tonB-dependent outer membrane receptors. Like the abundant planktonic marine bacterial clade SAR11, SAR86 exhibits metabolic streamlining, but also a distinct carbon compound specialization, possibly avoiding competition.

Functional bias of positively selected genes in Streptococcus genomes

Rates of nonsynonymous substitution (dN) significantly higher than rates of synonymous substitution (dS) have been used as evidence of positive selection for the fixation of advantageous point mutations. It has been suggested that positive selection contributes to the evolution of virulence factors and certain functional categories in bacterial pathogens. The genus Streptococcus contains a number of important human and agricultural pathogens. Here we assessed positive selection across 13 Streptococcus species, and their relationship with virulence factors and functional categories. We found that known virulence genes were subject to positive selection pressure as much as other genes. After false discovery rate correction for multiple comparisons, no functional categories were significantly over- or under-represented in positively selected genes relative to other genes. Our results suggest that within the genus Streptococcus positive selection based on dN/dS ratios is not distributed with bias across biological functions.

CharProtDB: a database of experimentally characterized protein annotations

CharProtDB (http://www.jcvi.org/charprotdb/) is a curated database of biochemically characterized proteins. It provides a source of direct rather than transitive assignments of function, designed to support automated annotation pipelines. The initial data set in CharProtDB was collected through manual literature curation over the years by analysts at the J. Craig Venter Institute (JCVI) [formerly The Institute of Genomic Research (TIGR)] as part of their prokaryotic genome sequencing projects. The CharProtDB has been expanded by import of selected records from publicly available protein collections whose biocuration indicated direct rather than homology-based assignment of function. Annotations in CharProtDB include gene name, symbol and various controlled vocabulary terms, including Gene Ontology terms, Enzyme Commission number and TransportDB accession. Each annotation is referenced with the source; ideally a journal reference, or, if imported and lacking one, the original database source.

BacMap: an up-to-date electronic atlas of annotated bacterial genomes

Originally released in 2005, BacMap is an electronic, interactive atlas of fully sequenced bacterial genomes. It contains fully labeled, zoomable and searchable chromosome maps for essentially all sequenced prokaryotic (archaebacterial and eubacterial) species. Each map can be zoomed to the level of individual genes and each gene is hyperlinked to a richly annotated gene card. The latest release of BacMap (http://bacmap.wishartlab.com/) now contains data for more than 1700 bacterial species (~10× more than the 2005 release), corresponding to more than 2800 chromosome and plasmid maps. All bacterial genome maps are now supplemented with separate prophage genome maps as well as separate tRNA and rRNA maps. Each bacterial chromosome entry in BacMap also contains graphs and tables on a variety of gene and protein statistics. Likewise, every bacterial species entry contains a bacterial 'biography' card, with taxonomic details, phenotypic details, textual descriptions and images (when available). Improved data browsing and searching tools have also been added to allow more facile filtering, sorting and display of the chromosome maps and their contents.

The UniProt-GO Annotation database in 2011

The GO annotation dataset provided by the UniProt Consortium (GOA: http://www.ebi.ac.uk/GOA) is a comprehensive set of evidenced-based associations between terms from the Gene Ontology resource and UniProtKB proteins. Currently supplying over 100 million annotations to 11 million proteins in more than 360 000 taxa, this resource has increased 2-fold over the last 2 years and has benefited from a wealth of checks to improve annotation correctness and consistency as well as now supplying a greater information content enabled by GO Consortium annotation format developments. Detailed, manual GO annotations obtained from the curation of peer-reviewed papers are directly contributed by all UniProt curators and supplemented with manual and electronic annotations from 36 model organism and domain-focused scientific resources. The inclusion of high-quality, automatic annotation predictions ensures the UniProt GO annotation dataset supplies functional information to a wide range of proteins, including those from poorly characterized, non-model organism species. UniProt GO annotations are freely available in a range of formats accessible by both file downloads and web-based views. In addition, the introduction of a new, normalized file format in 2010 has made for easier handling of the complete UniProt-GOA data set.