Update on RefSeq microbial genomes resources

The quest for environmental analytical microbiology: absolute quantitative microbiome using cellular internal standards

BACKGROUND

High-throughput sequencing has revolutionized environmental microbiome research, providing both quantitative and qualitative insights into nucleic acid targets in the environment. The resulting microbial composition (community structure) data are essential for environmental analytical microbiology, enabling characterization of community dynamics and assessing microbial pollutants for the development of intervention strategies. However, the relative abundances derived from sequencing impede comparisons across samples and studies.

RESULTS

This review systematically summarizes various absolute quantification (AQ) methods and their applications to obtain the absolute abundance of microbial cells and genetic elements. By critically comparing the strengths and limitations of AQ methods, we advocate the use of cellular internal standard-based high-throughput sequencing as an appropriate AQ approach for studying environmental microbiome originated from samples of complex matrices and high heterogeneity. To minimize ambiguity and facilitate cross-study comparisons, we outline essential reporting elements for technical considerations, and provide a checklist as a reference for environmental microbiome research.

CONCLUSIONS

In summary, we propose absolute microbiome quantification using cellular internal standards for environmental analytical microbiology, and we anticipate that this approach will greatly benefit future studies. Video Abstract.

Evaluation of 16S rRNA genes sequences and genome-based analysis for identification of non-pathogenic Yersinia

16S rRNA genes sequencing has been used for routine species identification and phylogenetic studies of bacteria. However, the high sequence similarity between some species and heterogeneity within copies at the intragenomic level could be a limiting factor of discriminatory ability. In this study, we aimed to compare 16S rRNA genes sequences and genome-based analysis (core SNPs and ANI) for identification of non-pathogenic Yersinia. We used complete and draft genomes of 373 Yersinia strains from the NCBI Genome database. The taxonomic affiliations of 34 genomes based on core SNPs and the ANI results did not match those specified in the GenBank database (NCBI). The intragenic homology of the 16S rRNA gene copies exceeded 99.5% in complete genomes, but above 50% of genomes have four or more variants of the 16S rRNA gene. Among 327 draft genomes of non-pathogenic Yersinia, 11% did not have a full-length 16S rRNA gene. Most of draft genomes has one copy of gene and it is not possible to define the intragenomic heterogenicity. The average homology of 16S rRNA gene was 98.76%, and the maximum variability was 2.85%. The low degree of genetic heterogenicity of the gene (0.36%) was determined in group Y. pekkanenii/Y. proxima/Y. aldovae/Y. intermedia/Y. kristensenii/Y. rochesterensis. The identical gene sequences were found in the genomes of the Y. intermedia and Y. rochesterensis strains identified using ANI and core SNPs analyses. The phylogenetic tree based on 16S rRNA genes differed from the tree based on core SNPs of the genomes and did not represent phylogenetic relationship between the Yersinia species. These findings will help to fill the data gaps in genome characteristics of deficiently studied non-pathogenic Yersinia.

Beav: a bacterial genome and mobile element annotation pipeline

Comprehensive and accurate genome annotation is crucial for inferring the predicted functions of an organism. Numerous tools exist to annotate genes, gene clusters, mobile genetic elements, and other diverse features. However, these tools and pipelines can be difficult to install and run, be specialized for a particular element or feature, or lack annotations for larger elements that provide important genomic context. Integrating results across analyses is also important for understanding gene function. To address these challenges, we present the Beav annotation pipeline. Beav is a command-line tool that automates the annotation of bacterial genome sequences, mobile genetic elements, molecular systems and gene clusters, key regulatory features, and other elements. Beav uses existing tools in addition to custom models, scripts, and databases to annotate diverse elements, systems, and sequence features. Custom databases for plant-associated microbes are incorporated to improve annotation of key virulence and symbiosis genes in agriculturally important pathogens and mutualists. Beav includes an optional Agrobacterium-specific pipeline that identifies and classifies oncogenic plasmids and annotates plasmid-specific features. Following the completion of all analyses, annotations are consolidated to produce a single comprehensive output. Finally, Beav generates publication-quality genome and plasmid maps. Beav is on Bioconda and is available for download at https://github.com/weisberglab/beav.

IMPORTANCE

Annotation of genome features, such as the presence of genes and their predicted function, or larger loci encoding secretion systems or biosynthetic gene clusters, is necessary for understanding the functions encoded by an organism. Genomes can also host diverse mobile genetic elements, such as integrative and conjugative elements and/or phages, that are often not annotated by existing pipelines. These elements can horizontally mobilize genes encoding for virulence, antimicrobial resistance, or other adaptive functions and alter the phenotype of an organism. We developed a software pipeline, called Beav, that combines new and existing tools for the comprehensive annotation of these and other major features. Existing pipelines often misannotate loci important for virulence or mutualism in plant-associated bacteria. Beav includes custom databases and optional workflows for the improved annotation of plant-associated bacteria. Beav is designed to be easy to install and run, making comprehensive genome annotation broadly available to the research community.

MIST: A microbial identification and source tracking system for next-generation sequencing data

The Professional Committee of Microbiology of the National Pharmacopoeia Commission organized the drafting of the Technical Guidelines for Microbial Whole Genome Sequencing (WGS), aiming to standardize the method process and technical indicators of microbial WGS and ensure the accuracy of sequencing and identification. On the basis of the Guidelines, we developed an integrated microbial identification and source tracking (MIST) system, which could meet the needs of microbial identification and contamination investigation in food and drug quality control. MIST integrates three analysis pipelines: 16S/18S/internal transcribed spacer amplicon-based microbial identification, WGS-based microbial identification, and single-nucleotide polymorphism-based microbial source tracking. MIST can analyze sequence data in a variety of formats, such as Fasta, base call file, and FASTQ. It can be connected to a high-throughput sequencing instrument to acquire sequencing data directly. We also developed a publicly accessible web server for MIST (http://syj.i-sanger.cn).

Genomic and phenotypic features of Acinetobacter baumannii isolated from oil reservoirs reveal a novel subspecies specialized in degrading hazardous hydrocarbons

The genus Acinetobacter encompasses biotechnologically relevant species and nosocomial pathogens. In this study, nine isolates recovered from different oil reservoir samples showed the ability to grow with petroleum as the only carbon source and possessed the ability to emulsify kerosene. The whole genomes of the nine strains were sequenced and analyzed. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values of all strains were compared to the reference strains, and the results were below the reference values (<97.88 and 82, respectively), suggesting that the isolates belong to a new subspecies of Acinetobacter baumannii. The name Acinetobacter baumannii oleum ficedula is proposed. A comparison of the whole genome repertoire of 290 Acinetobacter species indicated that the strains in this study resemble non-pathogenic Acinetobacter strains. However, the new isolates resemble A. baumannii when comparing virulence factors. The isolates in this study carry many genes involved in hydrocarbon degradation, indicating the potential to degrade most toxic compounds listed by environmental regulatory agencies such as ATSDR, EPA, and CONAMA. In addition, despite the absence of known biosurfactant or bioemulsifier genes, the strains showed emulsifying activity, suggesting the presence of new pathways or genes related to this process. This study investigated the genomic, phenotypic, and biochemical features of the novel environmental subspecies A. baumannii oleum ficedula, revealing their potential to degrade hydrocarbons and to produce biosurfactants or bioemulsifiers. Applying these environmental subspecies in bioaugmentation strategies sheds light on future approaches to bioremediation. The study shows the importance of genomic analysis of environmental strains and their inclusion in metabolic pathways databases, highlighting unique enzymes/alternative pathways for consuming hazardous hydrocarbons.

Translational informatics for human microbiota: data resources, models and applications

With the rapid development of human intestinal microbiology and diverse microbiome-related studies and investigations, a large amount of data have been generated and accumulated. Meanwhile, different computational and bioinformatics models have been developed for pattern recognition and knowledge discovery using these data. Given the heterogeneity of these resources and models, we aimed to provide a landscape of the data resources, a comparison of the computational models and a summary of the translational informatics applied to microbiota data. We first review the existing databases, knowledge bases, knowledge graphs and standardizations of microbiome data. Then, the high-throughput sequencing techniques for the microbiome and the informatics tools for their analyses are compared. Finally, translational informatics for the microbiome, including biomarker discovery, personalized treatment and smart healthcare for complex diseases, are discussed.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Comparison of Metagenomics and Metatranscriptomics Tools: A Guide to Making the Right Choice

The study of microorganisms is a field of great interest due to their environmental (e.g., soil contamination) and biomedical (e.g., parasitic diseases, autism) importance. The advent of revolutionary next-generation sequencing techniques, and their application to the hypervariable regions of the 16S, 18S or 23S ribosomal subunits, have allowed the research of a large variety of organisms more in-depth, including bacteria, archaea, eukaryotes and fungi. Additionally, together with the development of analysis software, the creation of specific databases (e.g., SILVA or RDP) has boosted the enormous growth of these studies. As the cost of sequencing per sample has continuously decreased, new protocols have also emerged, such as shotgun sequencing, which allows the profiling of all taxonomic domains in a sample. The sequencing of hypervariable regions and shotgun sequencing are technologies that enable the taxonomic classification of microorganisms from the DNA present in microbial communities. However, they are not capable of measuring what is actively expressed. Conversely, we advocate that metatranscriptomics is a "new" technology that makes the identification of the mRNAs of a microbial community possible, quantifying gene expression levels and active biological pathways. Furthermore, it can be also used to characterise symbiotic interactions between the host and its microbiome. In this manuscript, we examine the three technologies above, and discuss the implementation of different software and databases, which greatly impact the obtaining of reliable results. Finally, we have developed two easy-to-use pipelines leveraging Nextflow technology. These aim to provide everything required for an average user to perform a metagenomic analysis of marker genes with QIMME2 and a metatranscriptomic study using Kraken2/Bracken.

Analysis of the Evolution of the MoxR ATPases

MoxR proteins comprise a family of ATPases Associated with diverse cellular Activities (AAA+). These proteins are widespread and found across the diversity of prokaryotic species. Despite their ubiquity, members of the group remain poorly characterized. Only a few examples of MoxR proteins have been associated with cellular roles, where they have been shown to perform chaperone-like functions. A characteristic feature of MoxR proteins is their association with proteins containing the von Willebrand factor type A (VWA) domain. In an effort to understand the spread and diversity of the MoxR family, an evolutionary approach was undertaken. Phylogenetic techniques were used to define nine major subfamilies within the MoxR family. A combination of phylogenetic and genomic approaches was utilized to explore the extent of the partnership between the MoxR and VWA domain containing proteins (VWA proteins). These analyses led to the clarification of genetic linkages between MoxR and VWA proteins. A significant partnership is described here, as seven of nine MoxR subfamilies were found to be linked to VWA proteins. Available genomic data were also used to assess the intraprotein diversification of MoxR and VWA protein sequences. Data clearly indicated that, in MoxR proteins, the ATPase domain is maintained with high conservation while the remaining protein sequence evolves at a faster rate; a similar pattern was observed for the VWA domain in VWA proteins. Overall, our data present insights into the modular evolution of MoxR ATPases.

Colonization of White-Tailed Deer (Odocoileus virginianus) from Urban and Suburban Environments with Cephalosporinase- and Carbapenemase-Producing Enterobacterales

Wildlife play a role in the acquisition, maintenance, and dissemination of antimicrobial resistance (AMR). This is especially true at the human-domestic animal-wildlife interface, like urbanized areas, where interactions occur that can promote the cross-over of AMR bacteria and genes. We conducted a 2-year fecal surveillance (n = 783) of a white-tailed deer (WTD) herd from an urban park system in Ohio to identify and characterize cephalosporin-resistant and carbapenemase-producing bacteria using selective enrichment. Using generalized linear mixed models we found that older (OR = 2.3, P < 0.001), male (OR = 1.8, P = 0.001) deer from urbanized habitats (OR = 1.4, P = 0.001) were more likely to harbor extended-spectrum cephalosporin-resistant Enterobacterales. In addition, we isolated two carbapenemase-producing Enterobacterales (CPE), a Klebsiella quasipneumoniae harboring bla_KPC-2 and an Escherichia coli harboring bla_NDM-5, from two deer from urban habitats. The genetic landscape of the plasmid carrying bla_KPC-2 was unique, not clustering with other reported plasmids encoding KPC-2, and only sharing 78% of its sequence with its nearest match. While the plasmid carrying bla_NDM-5 shared sequence similarity with other reported plasmids encoding NDM-5, the intact IS26 mobile genetic elements surrounding multiple drug resistant regions, including the bla_NDM-5, has been reported infrequently. Both carbapenemase genes were successfully conjugated to a J53 recipient conferring a carbapenem-resistant phenotype. Our findings highlight that urban environments play a significant role on the transmission of AMR bacteria and genes to wildlife and suggest WTD may play a role in the dissemination of clinically and epidemiologically relevant antimicrobial resistant bacteria. IMPORTANCE The role of wildlife in the spread of antimicrobial resistance is not fully characterized. Some wildlife, including white-tailed deer (WTD), can thrive in suburban and urban environments. This may result in the exchange of antimicrobial resistant bacteria and resistance genes between humans and wildlife, and lead to their spread in the environment. We found that WTD living in an urban park system carried antimicrobial resistant bacteria that were important to human health and resistant to antibiotics used to treat serious bacterial infections. This included two deer that carried bacteria resistant to carbapenem antibiotics which are critically important for treatment of life-threatening infections. These two bacteria had the ability to transfer their AMR resistance genes to other bacteria, making them a threat to public health. Our results suggest that WTD may contribute to the spread of antimicrobial resistant bacteria in the environment.

Genome-resolved metagenomics using environmental and clinical samples

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

A Win-Loss Interaction on Fe0 Between Methanogens and Acetogens From a Climate Lake

Diverse physiological groups congregate into environmental corrosive biofilms, yet the interspecies interactions between these corrosive physiological groups are seldom examined. We, therefore, explored Fe⁰-dependent cross-group interactions between acetogens and methanogens from lake sediments. On Fe⁰, acetogens were more corrosive and metabolically active when decoupled from methanogens, whereas methanogens were more metabolically active when coupled with acetogens. This suggests an opportunistic (win-loss) interaction on Fe⁰ between acetogens (loss) and methanogens (win). Clostridia and Methanobacterium were the major candidates doing acetogenesis and methanogenesis after four transfers (metagenome sequencing) and the only groups detected after 11 transfers (amplicon sequencing) on Fe⁰. Since abiotic H₂ failed to explain the high metabolic rates on Fe⁰, we examined whether cell exudates (spent media filtrate) promoted the H₂-evolving reaction on Fe⁰ above abiotic controls. Undeniably, spent media filtrate generated three- to four-fold more H₂ than abiotic controls, which could be partly explained by thermolabile enzymes and partly by non-thermolabile constituents released by cells. Next, we examined the metagenome for candidate enzymes/shuttles that could catalyze H₂ evolution from Fe⁰ and found candidate H₂-evolving hydrogenases and an almost complete pathway for flavin biosynthesis in Clostridium. Clostridial ferredoxin-dependent [FeFe]-hydrogenases may be catalyzing the H₂-evolving reaction on Fe⁰, explaining the significant H₂ evolved by spent media exposed to Fe⁰. It is typical of Clostridia to secrete enzymes and other small molecules for lytic purposes. Here, they may secrete such molecules to enhance their own electron uptake from extracellular electron donors but indirectly make their H₂-consuming neighbors-Methanobacterium-fare five times better in their presence. The particular enzymes and constituents promoting H₂ evolution from Fe⁰ remain to be determined. However, we postulate that in a static environment like corrosive crust biofilms in lake sediments, less corrosive methanogens like Methanobacterium could extend corrosion long after acetogenesis ceased, by exploiting the constituents secreted by acetogens.

Unraveling bacteria-mediated degradation of lignin-derived aromatic compounds in a freshwater environment

Terrestrial organic carbon-lignin plays a crucial role in the global carbon balance. However, limited studies presented the functional and ecological traits of lignin decomposers population in natural aquatic ecosystem. In this study, we performed a multi-omics analysis by deploying amplicon, metagenomic, and metatranscriptomic approaches to identify the key potential degraders and pathways involved lignin-derived aromatic compounds in the later stage of lignin degradation. By establishing microcosms with model lignin-derived aromatic compound (vanillic acid, VAN), based on the estimated absolute abundance (EAA) and the metagenome-assembled genomes (MAGs), novel potential lignin-derived aromatic compounds degraders were identified in the aquatic ecosystem. Furthermore, members of the phyla Proteobacteria and Actinobacteria were the potential major lignin-derived aromatic compounds degraders in the studied ecosystem. Our study demonstrated that genomes of the class Betaproteobacteria (Proteobacteria) possess a complete enzymatic system for the degradation of diarylpropanes, vanillate and protocatechuate, besides having the capacity to degrade other lignin-derived aromatic compounds. This study provides strong evidence for the ability of aquatic bacteria to degrade lignin-derived aromatic compounds and suggest that different microbes might occupy different niches in the later stage of lignin degradation.

Gut microbiota profiles of treatment-naïve adult acute myeloid leukemia patients with neutropenic fever during intensive chemotherapy

The intestinal bacterial flora of febrile neutropenic patients has been found to be significantly diverse. However, there are few reports of alterations of in adult acute myeloid leukemia (AML) patients. Stool samples of each treatment-naïve AML patient were collected the day before initiation of induction chemotherapy (pretreatment), on the first date of neutropenic fever and first date of bone marrow recovery. Bacterial DNA was extracted from stool samples and bacterial 16s ribosomal RNA genes were sequenced by next-generation sequencing. Relative abundance, overall richness, Shannon's diversity index and Simpson's diversity index were calculated. No antimicrobial prophylaxis was in placed in all participants. Ten cases of AML patients (4 male and 6 female) were included with a median age of 39 years (range: 19–49) and all of patients developed febrile neutropenia. Firmicutes dominated during the period of neutropenic fever, subsequently declining after bone marrow recovery a pattern in contrast to that shown by Bacteroidetes and Proteobacteria. Enterococcus was more abundant in the febrile neutropenia period compared to pretreatment (mean difference +20.2; p < 0.0001) while Escherichia notably declined during the same period (mean difference -11.2; p = 0.0064). At the operational taxonomic unit (OTU) level, there was a significantly higher level of overall richness in the pretreatment period than in the febrile neutropenic episode (mean OTU of 203.1 vs. 131.7; p = 0.012). Both of the diversity indexes of Shannon and Simpson showed a significant decrease during the febrile neutropenic period. Adult AML patients with a first episode of febrile neutropenia after initial intensive chemotherapy demonstrated a significant decrease in gut microbiota diversity and the level of diversity remained constant despite recovery of bone marrow.

Complete Genomic Analysis of VRE From a Cattle Feedlot: Focus on 2 Antibiotic Resistance

Practices in intensive animal farming such as the extensive use of antimicrobials have significant impacts on the genetic make-up of bacterial communities, especially on that of human/animal commensals. In this report, whole genome sequencing of two vancomycin-resistant enterococci (VRE) isolates from a cattle feedlot in the North West Province, South Africa, was used to highlight the threats that extensive antimicrobial usage in intensive animal rearing represents for environmental microbiomes and the food chain. The genomic DNA of the studied strains was extracted using a DNA extraction kit. Whole-genome sequencing was performed through next-generation sequencing. The genomes of Enterococcus durans strain NWUTAL1 and Enterococcus gallinarum strain S52016 consisted of 3,279,618 and 2,374,946 bp, respectively with G + C contents of 40.76 and 43.13%, respectively. Antibiotic resistance genes (ARG), plasmids and virulence factors (involved in biofilm formation, colonization and copper/silver efflux system), were detected in the genomes of both strains. The presence of these genetic determinants in the studied strains is a cause for concern as they may disseminate and find their way into the food chain via horizontal gene transfer amongst bacteria of the different ecological niches. Issues of this nature cannot be undermined and are relevant as far as food safety is concerned.

Metabolism of Soy Isoflavones by Intestinal Bacteria: Genome Analysis of an Adlercreutzia Equolifaciens Strain That Does Not Produce Equol

Isoflavones are transformed in the gut into more estrogen-like compounds or into inactive molecules. However, neither the intestinal microbes nor the pathways leading to the synthesis of isoflavone-derived metabolites are fully known. In the present work, 73 fecal isolates from three women with an equol-producing phenotype were considered to harbor equol-related genes by qPCR. After typing, 57 different strains of different taxa were tested for their ability to act on the isoflavones daidzein and genistein. Strains producing small to moderate amounts of dihydrodaidzein and/or O-desmethylangolensin (O-DMA) from daidzein and dihydrogenistein from genistein were recorded. However, either alone or in several strain combinations, equol producers were not found, even though one of the strains, W18.34a (also known as IPLA37004), was identified as Adlercreutzia equolifaciens, a well-described equol-producing species. Analysis and comparison of A. equolifaciens W18.34a and A. equolifaciens DSM19450^T (an equol producer bacterium) genome sequences suggested a deletion in the former involving a large part of the equol operon. Furthermore, genome comparison of A. equolifaciens and Asaccharobacter celatus (other equol-producing species) strains from databases indicated many of these also showed deletions within the equol operon. The present results contribute to our knowledge to the activity of gut bacteria on soy isoflavones.

Toxin-Antitoxin Gene Pairs Found in Tn3 Family Transposons Appear To Be an Integral Part of the Transposition Module

Transposable elements (TEs) are important in genetic diversification due to their recombination properties and their ability to promote horizontal gene transfer. Over the last decades, much effort has been made to understand TE transposition mechanisms and their impact on prokaryotic genomes. For example, the Tn3 family is ubiquitous in bacteria, molding their host genomes by the paste-and-copy mechanism. In addition to the transposition module, Tn3 members often carry additional passenger genes (e.g., conferring antibiotic or heavy metal resistance and virulence), and three were previously known to carry a toxin-antitoxin (TA) system often associated with plasmid maintenance; however, the role of TA systems within the Tn3 family is unknown. The genetic context of TA systems in Tn3 members suggests that they may play a regulatory role in ensuring stable invasion of these Tns during transposition.

Much of the diversity of prokaryotic genomes is contributed by the tightly controlled recombination activity of transposons (Tns). The Tn3 family is arguably one of the most widespread transposon families. Members carry a large range of passenger genes incorporated into their structures. Family members undergo replicative transposition using a DDE transposase to generate a cointegrate structure which is then resolved by site-specific recombination between specific DNA sequences (res) on each of the two Tn copies in the cointegrate. These sites also carry promoters controlling expression of the recombinase and transposase. We report here that a number of Tn3 members encode a type II toxin-antitoxin (TA) system, typically composed of a stable toxin and a labile antitoxin that binds the toxin and inhibits its lethal activity. This system serves to improve plasmid maintenance in a bacterial population and, until recently, was believed to be associated with bacterial persistence. At least six different TA gene pairs are associated with various Tn3 members. Our data suggest that several independent acquisition events have occurred. In contrast to most Tn3 family passenger genes, which are generally located away from the transposition module, the TA gene pairs abut the res site upstream of the resolvase genes. Although their role when part of Tn3 family transposons is unclear, this finding suggests a potential role for the embedded TA in stabilizing the associated transposon with the possibility that TA expression is coupled to expression of transposase and resolvase during the transposition process itself.

Functional cooperation of the glycine synthase-reductase and Wood-Ljungdahl pathways for autotrophic growth of Clostridium drakei

Despite sharing the first four reactions, coutilization of the Wood–Ljungdahl pathway (WLP) with the glycine synthase-reductase pathway (GSRP) and reductive glycine pathway (RGP) to fix C1 compounds has remained unknown. In this study, using Clostridium drakei, we elucidated the role of the GSRP and RGP in the presence of the WLP, via a genome-scale metabolic model, RNA-seq, ¹³C isotope-based metabolite-tracing experiments, biochemical assays, and heterologous expression. Overall, the data suggested the pathways are functional under autotrophic conditions. Along with the WLP, GSRP and RGP convert CO₂ to glycine and then to acetyl-phosphate and serine, which then obtain ATP by producing acetate and operate with limited reducing power. This is a unique coutilization of the pathways under autotrophic conditions in acetogens.

Among CO₂-fixing metabolic pathways in nature, the linear Wood–Ljungdahl pathway (WLP) in phylogenetically diverse acetate-forming acetogens comprises the most energetically efficient pathway, requires the least number of reactions, and converts CO₂ to formate and then into acetyl-CoA. Despite two genes encoding glycine synthase being well-conserved in WLP gene clusters, the functional role of glycine synthase under autotrophic growth conditions has remained uncertain. Here, using the reconstructed genome-scale metabolic model iSL771 based on the completed genome sequence, transcriptomics, ¹³C isotope-based metabolite-tracing experiments, biochemical assays, and heterologous expression of the pathway in another acetogen, we discovered that the WLP and the glycine synthase pathway are functionally interconnected to fix CO₂, subsequently converting CO₂ into acetyl-CoA, acetyl-phosphate, and serine. Moreover, the functional cooperation of the pathways enhances CO₂ consumption and cellular growth rates via bypassing reducing power required reactions for cellular metabolism during autotrophic growth of acetogens.

proGenomes2: an improved database for accurate and consistent habitat, taxonomic and functional annotations of prokaryotic genomes

Microbiology depends on the availability of annotated microbial genomes for many applications. Comparative genomics approaches have been a major advance, but consistent and accurate annotations of genomes can be hard to obtain. In addition, newer concepts such as the pan-genome concept are still being implemented to help answer biological questions. Hence, we present proGenomes2, which provides 87 920 high-quality genomes in a user-friendly and interactive manner. Genome sequences and annotations can be retrieved individually or by taxonomic clade. Every genome in the database has been assigned to a species cluster and most genomes could be accurately assigned to one or multiple habitats. In addition, general functional annotations and specific annotations of antibiotic resistance genes and single nucleotide variants are provided. In short, proGenomes2 provides threefold more genomes, enhanced habitat annotations, updated taxonomic and functional annotation and improved linkage to the NCBI BioSample database. The database is available at http://progenomes.embl.de/.

High-resolution mycobiota analysis reveals dynamic intestinal translocation preceding invasive candidiasis

The intestinal microbiota is a complex community of bacteria, archaea, viruses, protists and fungi^1,2. While the composition of bacterial constituents has been linked to immune homeostasis and to infectious susceptibility^3–7, the role of non-bacterial constituents and of cross-kingdom microbial interactions in these processes is poorly understood^2,8. Fungi represent a major cause of infectious morbidity and mortality in immune-compromised individuals, though the relationship of intestinal fungi (i.e., the mycobiota) with fungal bloodstream infections (BSI) remains undefined⁹. We integrated an optimized bioinformatics pipeline with high-resolution mycobiota sequencing and comparative genomic analyses of fecal and blood specimens from recipients of allogeneic hematopoietic cell transplant (allo-HCT). Patients with Candida BSI experienced a prior marked intestinal expansion of pathogenic Candida species; this expansion consisted of a complex dynamic between multiple species and subspecies with a stochastic translocation pattern into the bloodstream. The intestinal expansion of pathogenic Candida species was associated with a significant loss in bacterial burden and diversity, particularly in the anaerobes. Thus, simultaneous analysis of intestinal fungi and bacteria identifies dysbiosis states across kingdoms that may promote fungal translocation and facilitate invasive disease. These findings support microbiota-driven approaches to identify patients at risk for fungal BSI for pre-emptive therapeutic intervention.

Baltic Sea methanogens compete with acetogens for electrons from metallic iron

Microbially induced corrosion of metallic iron (Fe⁰)-containing structures is an environmental and economic hazard. Methanogens are abundant in low-sulfide environments and yet their specific role in Fe⁰ corrosion is poorly understood. In this study, Sporomusa and Methanosarcina dominated enrichments from Baltic Sea methanogenic sediments that were established with Fe⁰ as the sole electron donor and CO₂ as the electron acceptor. The Baltic-Sporomusa was phylogenetically affiliated to the electroactive acetogen S. silvacetica. Baltic-Sporomusa adjusted rapidly to growth on H₂. On Fe⁰, spent filtrate enhanced growth of this acetogen suggesting that it was using endogenous enzymes to retrieve electrons and produce acetate. Previous studies have proposed that acetate produced by acetogens can feed commensal acetoclastic methanogens such as Methanosarcina. However, Baltic-methanogens could not generate methane from acetate, plus the decrease or absence of acetogens stimulated their growth. The decrease in numbers of Sporomusa was concurrent with an upsurge in Methanosarcina and increased methane production, suggesting that methanogens compete with acetogens for electrons from Fe⁰. Furthermore, Baltic-methanogens were unable to use H₂ (1.5 atm) for methanogenesis and were inhibited by spent filtrate additions, indicating that enzymatically produced H₂ is not a favorable electron donor. We hypothesize that Baltic-methanogens retrieve electrons from Fe⁰ via a yet enigmatic direct electron uptake mechanism.

Machine Learning Approaches for Epidemiological Investigations of Food-Borne Disease Outbreaks

Foodborne diseases (FBDs) are infections of the gastrointestinal tract caused by foodborne pathogens (FBPs) such as bacteria [Salmonella, Listeria monocytogenes and Shiga toxin-producing E. coli (STEC)] and several viruses, but also parasites and some fungi. Artificial intelligence (AI) and its sub-discipline machine learning (ML) are re-emerging and gaining an ever increasing popularity in the scientific community and industry, and could lead to actionable knowledge in diverse ranges of sectors including epidemiological investigations of FBD outbreaks and antimicrobial resistance (AMR). As genotyping using whole-genome sequencing (WGS) is becoming more accessible and affordable, it is increasingly used as a routine tool for the detection of pathogens, and has the potential to differentiate between outbreak strains that are closely related, identify virulence/resistance genes and provide improved understanding of transmission events within hours to days. In most cases, the computational pipeline of WGS data analysis can be divided into four (though, not necessarily consecutive) major steps: de novo genome assembly, genome characterization, comparative genomics, and inference of phylogeny or phylogenomics. In each step, ML could be used to increase the speed and potentially the accuracy (provided increasing amounts of high-quality input data) of identification of the source of ongoing outbreaks, leading to more efficient treatment and prevention of additional cases. In this review, we explore whether ML or any other form of AI algorithms have already been proposed for the respective tasks and compare those with mechanistic model-based approaches.

A review of methods and databases for metagenomic classification and assembly

Microbiome research has grown rapidly over the past decade, with a proliferation of new methods that seek to make sense of large, complex data sets. Here, we survey two of the primary types of methods for analyzing microbiome data: read classification and metagenomic assembly, and we review some of the challenges facing these methods. All of the methods rely on public genome databases, and we also discuss the content of these databases and how their quality has a direct impact on our ability to interpret a microbiome sample.

A review of methods and databases for metagenomic classification and assembly

Microbiome research has grown rapidly over the past decade, with a proliferation of new methods that seek to make sense of large, complex data sets. Here, we survey two of the primary types of methods for analyzing microbiome data: read classification and metagenomic assembly, and we review some of the challenges facing these methods. All of the methods rely on public genome databases, and we also discuss the content of these databases and how their quality has a direct impact on our ability to interpret a microbiome sample.

Genome-wide correlation analysis suggests different roles of CRISPR-Cas systems in the acquisition of antibiotic resistance genes in diverse species

CRISPR-Cas systems are widespread in bacterial and archaeal genomes, and in their canonical role in phage defence they confer a fitness advantage. However, CRISPR-Cas may also hinder the uptake of potentially beneficial genes. This is particularly true under antibiotic selection, where preventing the uptake of antibiotic resistance genes could be detrimental. Newly discovered features within these evolutionary dynamics are anti-CRISPR genes, which inhibit specific CRISPR-Cas systems. We hypothesized that selection for antibiotic resistance might have resulted in an accumulation of anti-CRISPR genes in genomes that harbour CRISPR-Cas systems and horizontally acquired antibiotic resistance genes. To assess that question, we analysed correlations between the CRISPR-Cas, anti-CRISPR and antibiotic resistance gene content of 104 947 reference genomes, including 5677 different species. In most species, the presence of CRISPR-Cas systems did not correlate with the presence of antibiotic resistance genes. However, in some clinically important species, we observed either a positive or negative correlation of CRISPR-Cas with antibiotic resistance genes. Anti-CRISPR genes were common enough in four species to be analysed. In Pseudomonas aeruginosa, the presence of anti-CRISPRs was associated with antibiotic resistance genes. This analysis indicates that the role of CRISPR-Cas and anti-CRISPRs in the spread of antibiotic resistance is likely to be very different in particular pathogenic species and clinical environments. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Brachyspira catarrhinii sp. nov., an anaerobic intestinal spirochaete isolated from vervet monkeys may have been misidentified as Brachyspira aalborgi in previous studies

To date nine species of anaerobic intestinal spirochaetes have been validly assigned to the genus Brachyspira. These include both pathogenic and non-pathogenic species. In the current study a genomic analysis of a novel spirochaete isolate was undertaken to determine whether it is a distinct species that previously has been misidentified as Brachyspira aalborgi. The genome of spirochaete strain Z12 isolated from the faeces of a vervet monkey was sequenced and compared to the genomes of the type strains of the nine assigned Brachyspira species. Genome to Genome Distance (GGD) values and Average Nucleotide Identity (ANI) values were determined. Single nucleotide polymorphisms (SNP) were used to create a phylogenetic tree to assess relatedness. The 16S rRNA gene sequences of the strains were aligned and the similarity amongst the Brachyspira species was recorded. Multilocus sequence typing (MLST) using five loci was conducted on Z12 and results compared with those for other Brachyspira isolates. Assembly of the Z12 sequences revealed a 2,629,108 bp genome with an average G + C content of 31.3%. The GGD, ANI, 16S rRNA gene sequence comparisons and the MLST results all indicated that Z12 represents a distinct species within the genus Brachyspira, with its nearest neighbour being B. aalborgi. Spirochaete strain Z12^T was assigned as the type strain of a new species, Brachyspira catarrhinii sp. nov. The diagnostic PCR currently in use to detect B. aalborgi cross-reacts with Z12, but RFLP analysis of PCR product can be used to distinguish the two species. Previous reports of non-human primates being colonised by B. aalborgi based on PCR results may have been incorrect. The development of an improved diagnostic method will allow future studies on the distribution and possible clinical significance of these two anaerobic spirochaete species.

Curated BLAST for Genomes

Given a microbe’s genome sequence, we often want to predict what capabilities the organism has, such as which nutrients it requires or which energy sources it can use. Or, we know the organism has a capability and we want to find the genes involved. Scientists often use automated gene annotations to find relevant genes, but automated annotations are often vague or incorrect. Curated BLAST finds candidate genes for a capability without relying on automated annotations. First, Curated BLAST finds proteins (usually from other organisms) whose functions have been studied experimentally and whose curated descriptions match a query. Then, it searches the genome of interest for similar proteins and returns a list of candidates. Curated BLAST is fast and often finds relevant genes that are missed by automated annotation.

Curated BLAST for Genomes finds candidate genes for a process or an enzymatic activity within a genome of interest. In contrast to annotation tools, which usually predict a single activity for each protein, Curated BLAST asks if any of the proteins in the genome are similar to characterized proteins that are relevant. Given a query such as an enzyme’s name or an EC number, Curated BLAST searches the curated descriptions of over 100,000 characterized proteins, and it compares the relevant characterized proteins to the predicted proteins in the genome of interest. In case of errors in the gene models, Curated BLAST also searches the six-frame translation of the genome. Curated BLAST is available at http://papers.genomics.lbl.gov/curated.

IMPORTANCE Given a microbe’s genome sequence, we often want to predict what capabilities the organism has, such as which nutrients it requires or which energy sources it can use. Or, we know the organism has a capability and we want to find the genes involved. Scientists often use automated gene annotations to find relevant genes, but automated annotations are often vague or incorrect. Curated BLAST finds candidate genes for a capability without relying on automated annotations. First, Curated BLAST finds proteins (usually from other organisms) whose functions have been studied experimentally and whose curated descriptions match a query. Then, it searches the genome of interest for similar proteins and returns a list of candidates. Curated BLAST is fast and often finds relevant genes that are missed by automated annotation.

Corynebacterium alimapuense sp. nov., an obligate marine actinomycete isolated from sediment of Valparaíso bay, Chile

A novel Gram-positive, non-motile, non-spore-forming and aerobic bacterium, designated strain VA37-3^T, was isolated from a marine sediment sample collected at 19.2 m water depth from Valparaíso bay, Chile. Strain VA37-3^T exhibits 97.6 % 16S rRNA gene sequence similarity to Corynebacterium marinum D7015^T, 96.4 % to Corynebacterium humireducens MFC-5^T and 96 % to Corynebacterium testudinoris M935/96/4^T; and a rpoB gene sequence similarity of 85.1 % to Corynebacterium pollutisoli VMS11^T, both analyses suggesting that strain VA37-3^T represents a novel species of Corynebacterium. Physiological testing indicated that strain VA37-3^T requires artificial sea water or sodium-supplemented media for growth, representing the first obligate marine actinomycete of the genus Corynebacterium. The genome of the proposed new species, along with the type strains of its most closely related species were sequenced and characterized. In silico genome-based similarity analyses revealed an ANIb of 72.8 % (C. marinum D7015^T), ANIm of 85.0 % (Corynebacterium mustelae DSM 45274^T), tetra of 0.90 (Corynebacterium callunae DSM 20147^T) and ggdc of 24.7 % (Corynebacterium kutscheri DSM 20755^T) when compared with the closest related strains. The genomic DNA G+C content of strain VA37-3^T was 57.0 %. Chemotaxonomic assessment of strain VN6-2^T showed the major fatty acids were C18 : 1ω9c and C16 : 0. Menaquinones predominantly consisted of MK-8(II-H2). Polar lipids consisted of diphosphatidylglycerol, glycolipids, phosphatidylglycerol, phosphoglycolipid and phosphatidylinositol. Mycolic acids also were present. Overall, the results from phylogenetic, phenotypic and genomic analyses confirmed that strain VA37-3^T represents a novel species of the genus Corynebacterium, for which the name Corynebacterium alimapuense sp. nov. is proposed, with VA37-3^T as the type strain (=CCUG 69366^T=NCIMB 15118^T).

Proteotyping bacteria: Characterization, differentiation and identification of pneumococcus and other species within the Mitis Group of the genus Streptococcus by tandem mass spectrometry proteomics

A range of methodologies may be used for analyzing bacteria, depending on the purpose and the level of resolution needed. The capability for recognition of species distinctions within the complex spectrum of bacterial diversity is necessary for progress in microbiological research. In clinical settings, accurate, rapid and cost-effective methods are essential for early and efficient treatment of infections. Characterization and identification of microorganisms, using, bottom-up proteomics, or "proteotyping", relies on recognition of species-unique or associated peptides, by tandem mass spectrometry analyses, dependent upon an accurate and comprehensive foundation of genome sequence data, allowing for differentiation of species, at amino acid-level resolution. In this study, the high resolution and accuracy of MS/MS-based proteotyping was demonstrated, through analyses of the three phylogenetically and taxonomically most closely-related species of the Mitis Group of the genus Streptococcus: i.e., the pathogenic species, Streptococcus pneumoniae (pneumococcus), and the commensal species, Streptococcus pseudopneumoniae and Streptococcus mitis. To achieve high accuracy, a genome sequence database used for matching peptides was created and carefully curated. Here, MS-based, bottom-up proteotyping was observed and confirmed to attain the level of resolution necessary for differentiating and identifying the most-closely related bacterial species, as demonstrated by analyses of species of the Streptococcus Mitis Group, even when S. pneumoniae were mixed with S. pseudopneumoniae and S. mitis, by matching and identifying more than 200 unique peptides for each species.

Genome-scale analysis of syngas fermenting acetogenic bacteria reveals the translational regulation for its autotrophic growth

Background

Acetogenic bacteria constitute promising biocatalysts for the conversion of CO₂/H₂ or synthesis gas (H₂/CO/CO₂) into biofuels and value-added biochemicals. These microorganisms are naturally capable of autotrophic growth via unique acetogenesis metabolism. Despite their biosynthetic potential for commercial applications, a systemic understanding of the transcriptional and translational regulation of the acetogenesis metabolism remains unclear.

Results

By integrating genome-scale transcriptomic and translatomic data, we explored the regulatory logic of the acetogenesis to convert CO₂ into biomass and metabolites in Eubacterium limosum. The results indicate that majority of genes associated with autotrophic growth including the Wood-Ljungdahl pathway, the reduction of electron carriers, the energy conservation system, and gluconeogenesis were transcriptionally upregulated. The translation efficiency of genes in cellular respiration and electron bifurcation was also highly enhanced. In contrast, the transcriptionally abundant genes involved in the carbonyl branch of the Wood-Ljungdahl pathway, as well as the ion-translocating complex and ATP synthase complex in the energy conservation system, showed decreased translation efficiency. The translation efficiencies of genes were regulated by 5′UTR secondary structure under the autotrophic growth condition.

Conclusions

The results illustrated that the acetogenic bacteria reallocate protein synthesis, focusing more on the translation of genes for the generation of reduced electron carriers via electron bifurcation, rather than on those for carbon metabolism under autotrophic growth.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5238-0) contains supplementary material, which is available to authorized users.

Synthetic microbe communities provide internal reference standards for metagenome sequencing and analysis

The complexity of microbial communities, combined with technical biases in next-generation sequencing, pose a challenge to metagenomic analysis. Here, we develop a set of internal DNA standards, termed “sequins” (sequencing spike-ins), that together constitute a synthetic community of artificial microbial genomes. Sequins are added to environmental DNA samples prior to library preparation, and undergo concurrent sequencing with the accompanying sample. We validate the performance of sequins by comparison to mock microbial communities, and demonstrate their use in the analysis of real metagenome samples. We show how sequins can be used to measure fold change differences in the size and structure of accompanying microbial communities, and perform quantitative normalization between samples. We further illustrate how sequins can be used to benchmark and optimize new methods, including nanopore long-read sequencing technology. We provide metagenome sequins, along with associated data sets, protocols, and an accompanying software toolkit, as reference standards to aid in metagenomic studies.

Complex microbial communities pose a challenge to metagenomic analysis. Here the authors develop ‘sequins’, internal DNA standards that represent a synthetic community of artificial genomes.

Using average nucleotide identity to improve taxonomic assignments in prokaryotic genomes at the NCBI

Average nucleotide identity analysis is a useful tool to verify taxonomic identities in prokaryotic genomes, for both complete and draft assemblies. Using optimum threshold ranges appropriate for different prokaryotic taxa, we have reviewed all prokaryotic genome assemblies in GenBank with regard to their taxonomic identity. We present the methods used to make such comparisons, the current status of GenBank verifications, and recent developments in confirming species assignments in new genome submissions.

16S rRNA Gene Sequencing for Deciphering the Colorectal Cancer Gut Microbiome: Current Protocols and Workflows

The human gut holds the densest microbiome ecosystem essential in maintaining a healthy host physiology, whereby disruption of this ecosystem has been linked to the development of colorectal cancer (CRC). The advent of next-generation sequencing technologies such as the 16S rRNA gene sequencing has enabled characterization of the CRC gut microbiome architecture in an affordable and culture-free approach. Nevertheless, the lack of standardization in handling and storage of biospecimens, nucleic acid extraction, 16S rRNA gene primer selection, length, and depth of sequencing and bioinformatics analyses have contributed to discrepancies found in various published studies of this field. Accurate characterization of the CRC microbiome found in different stages of CRC has the potential to be developed into a screening tool in the clinical setting. This mini review aims to concisely compile all available CRC microbiome studies performed till end of 2016 and to suggest standardized protocols that are crucial in developing a gut microbiome screening panel for CRC.

Functional Genomics and Phylogenetic Evidence Suggest Genus-Wide Cobalamin Production by the Globally Distributed Marine Nitrogen Fixer Trichodesmium

Only select prokaryotes can biosynthesize vitamin B₁₂ (i.e., cobalamins), but these organic co-enzymes are required by all microbial life and can be vanishingly scarce across extensive ocean biomes. Although global ocean genome data suggest cyanobacteria to be a major euphotic source of cobalamins, recent studies have highlighted that >95% of cyanobacteria can only produce a cobalamin analog, pseudo-B₁₂, due to the absence of the BluB protein that synthesizes the α ligand 5,6-dimethylbenzimidizole (DMB) required to biosynthesize cobalamins. Pseudo-B₁₂ is substantially less bioavailable to eukaryotic algae, as only certain taxa can intracellularly remodel it to one of the cobalamins. Here we present phylogenetic, metagenomic, transcriptomic, proteomic, and chemical analyses providing multiple lines of evidence that the nitrogen-fixing cyanobacterium Trichodesmium transcribes and translates the biosynthetic, cobalamin-requiring BluB enzyme. Phylogenetic evidence suggests that the Trichodesmium DMB biosynthesis gene, bluB, is of ancient origin, which could have aided in its ecological differentiation from other nitrogen-fixing cyanobacteria. Additionally, orthologue analyses reveal two genes encoding iron-dependent B₁₂ biosynthetic enzymes (cbiX and isiB), suggesting that iron availability may be linked not only to new nitrogen supplies from nitrogen fixation, but also to B₁₂ inputs by Trichodesmium. These analyses suggest that Trichodesmium contains the genus-wide genomic potential for a previously unrecognized role as a source of cobalamins, which may prove to considerably impact marine biogeochemical cycles.

Rewired RNAi-mediated genome surveillance in house dust mites

House dust mites are common pests with an unusual evolutionary history, being descendants of a parasitic ancestor. Transition to parasitism is frequently accompanied by genome rearrangements, possibly to accommodate the genetic change needed to access new ecology. Transposable element (TE) activity is a source of genomic instability that can trigger large-scale genomic alterations. Eukaryotes have multiple transposon control mechanisms, one of which is RNA interference (RNAi). Investigation of the dust mite genome failed to identify a major RNAi pathway: the Piwi-associated RNA (piRNA) pathway, which has been replaced by a novel small-interfering RNA (siRNA)-like pathway. Co-opting of piRNA function by dust mite siRNAs is extensive, including establishment of TE control master loci that produce siRNAs. Interestingly, other members of the Acari have piRNAs indicating loss of this mechanism in dust mites is a recent event. Flux of RNAi-mediated control of TEs highlights the unusual arc of dust mite evolution.

Investigation of small RNA populations in dust mites revealed absence of the piwi-associated RNA (piRNA) pathway. Apart from several nematode and platyhelminths lineages, piRNAs are an essential component of animal genome surveillance, actively targeting and silencing transposable elements. In dust mites, expansion of Dicer produced small-interfering RNA (siRNA) biology compensates for loss of piRNAs. The dramatic difference we find in dust mites is likely a consequence of their evolutionary history, which is marked by descent from a parasite to the current free-living form. Our study highlights a correlation between perturbation of transposon surveillance and shifts in ecology.

Nutrient-Colimited Trichodesmium as a Nitrogen Source or Sink in a Future Ocean

Nitrogen-fixing (N₂) cyanobacteria provide bioavailable nitrogen to vast ocean regions but are in turn limited by iron (Fe) and/or phosphorus (P), which may force them to employ alternative nitrogen acquisition strategies. The adaptive responses of nitrogen fixers to global-change drivers under nutrient-limited conditions could profoundly alter the current ocean nitrogen and carbon cycles. Here, we show that the globally important N₂ fixer Trichodesmium fundamentally shifts nitrogen metabolism toward organic-nitrogen scavenging following long-term high-CO₂ adaptation under iron and/or phosphorus (co)limitation. Global shifts in transcripts and proteins under high-CO₂/Fe-limited and/or P-limited conditions include decreases in the N₂-fixing nitrogenase enzyme, coupled with major increases in enzymes that oxidize trimethylamine (TMA). TMA is an abundant, biogeochemically important organic nitrogen compound that supports rapid Trichodesmium growth while inhibiting N₂ fixation. In a future high-CO₂ ocean, this whole-cell energetic reallocation toward organic nitrogen scavenging and away from N₂ fixation may reduce new-nitrogen inputs by Trichodesmium while simultaneously depleting the scarce fixed-nitrogen supplies of nitrogen-limited open-ocean ecosystems.IMPORTANCE Trichodesmium is among the most biogeochemically significant microorganisms in the ocean, since it supplies up to 50% of the new nitrogen supporting open-ocean food webs. We used Trichodesmium cultures adapted to high-CO₂ conditions for 7 years, followed by additional exposure to iron and/or phosphorus (co)limitation. We show that "future ocean" conditions of high CO₂ and concurrent nutrient limitation(s) fundamentally shift nitrogen metabolism away from nitrogen fixation and instead toward upregulation of organic nitrogen-scavenging pathways. We show that the responses of Trichodesmium to projected future ocean conditions include decreases in the nitrogen-fixing nitrogenase enzymes coupled with major increases in enzymes that oxidize the abundant organic nitrogen source trimethylamine (TMA). Such a shift toward organic nitrogen uptake and away from nitrogen fixation may substantially reduce new-nitrogen inputs by Trichodesmium to the rest of the microbial community in the future high-CO₂ ocean, with potential global implications for ocean carbon and nitrogen cycling.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

The acetic acid bacterium Gluconobacter oxydans 621H is characterized by its exceptional ability to incompletely oxidize a great variety of carbohydrates in the periplasm. The metabolism of this α-proteobacterium has been characterized to some extent, yet little is known about its transcriptomes and related data. In this study, we applied two different RNAseq approaches. Primary transcriptomes enriched for 5′-ends of transcripts were sequenced to detect transcription start sites, which allow subsequent analysis of promoter motifs, ribosome binding sites, and 5´-UTRs. Whole transcriptomes were sequenced to identify expressed genes and operon structures.

Results

Sequencing of primary transcriptomes of G. oxydans revealed 2449 TSSs, which were classified according to their genomic context followed by identification of promoter and ribosome binding site motifs, analysis of 5´-UTRs including validation of predicted cis-regulatory elements and correction of start codons. 1144 (41%) of all genes were found to be expressed monocistronically, whereas 1634 genes were organized in 571 operons. Together, TSSs and whole transcriptome data were also used to identify novel intergenic (18), intragenic (328), and antisense transcripts (313).

Conclusions

This study provides deep insights into the transcriptional landscapes of G. oxydans. The comprehensive transcriptome data, which we made publicly available, facilitate further analysis of promoters and other regulatory elements. This will support future approaches for rational strain development and targeted gene expression in G. oxydans. The corrections of start codons further improve the high quality genome reference and support future proteome analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4415-x) contains supplementary material, which is available to authorized users.

First Report of a Whole-Genome Shotgun Sequence of a Clinical Enterococcus faecalis Sequence Type 6 Strain from South Africa

Enterococcus faecalis is a lactic acid-producing Gram-positive bacterium commonly found in the intestinal tract of humans and animals; it is implicated in multidrug-resistant nosocomial infections. The draft genome of this E. faecalis sequence type 6 (ST6) strain consists of 3,215,228 bp, with 37.20% GC content, 3,048 predicted coding sequences, and 61 RNA genes.

First insights into a type II toxin-antitoxin system from the clinical isolate Mycobacterium sp. MHSD3, similar to epsilon/zeta systems

A putative type II toxin-antitoxin (TA) system was found in the clinical isolate Mycobacterium sp. MHSD3, a strain closely related to Mycobacterium chelonae. Further analyses of the protein sequences of the two genes revealed the presence of domains related to a TA system. BLAST analyses indicated the presence of closely related proteins in the genomes of other recently published M. chelonae strains. The functionality of both elements of the TA system was demonstrated when expressed in Escherichia coli cells, and the predicted structure of the toxin is very similar to those of well-known zeta-toxins, leading to the definition of a type II TA system similar to epsilon/zeta TA systems in strains that are closely related to M. chelonae.

Draft Genome Sequence of a Clinical Enterococcus faecium Sequence Type 18 Strain from South Africa

We report the first draft genome sequence of an Enterococcus faecium sequence type 18 (ST18) strain isolated from a tuberculosis patient in Africa. The genome is comprised of 3,202,539 bp, 501 contigs, 37.70% GC content, 3,202 protein-encoding genes, and 61 RNA genes. The resistome and virulome of this important pathogen are presented herein.

AureoWiki ̵ The repository of the Staphylococcus aureus research and annotation community

In light of continuously accumulating data and knowledge on major human pathogens, comprehensive and up-to-date sources of easily accessible information are urgently required. The AureoWiki database (http://aureowiki.med.uni-greifswald.de) provides detailed information on the genes and proteins of clinically and experimentally relevant S. aureus strains, currently covering NCTC 8325, COL, Newman, USA300_FPR3757, and N315. By implementing a pan-genome approach, AureoWiki facilitates the transfer of knowledge gained in studies with different S. aureus strains, thus supporting functional annotation and better understanding of this organism. All data related to a given gene or gene product is compiled on a strain-specific gene page. The gene pages contain sequence-based information complemented by data on, for example, protein function and localization, transcriptional regulation, and gene expression. The information provided is connected via links to other databases and published literature. Importantly, orthologous genes of the individual strains, which are linked by a pan-genome gene identifier and a unified gene name, are presented side by side using strain-specific tabs. The respective pan-genome gene page contains an orthologue table for 32 S. aureus strains, a multiple-strain genome viewer, a protein sequence alignment as well as other comparative information. The data collected in AureoWiki is also accessible through various download options in order to support bioinformatics applications. In addition, based on two large-scale gene expression data sets, AureoWiki provides graphical representations of condition-dependent mRNA levels and protein profiles under various laboratory and infection-related conditions.

Biogeographic conservation of the cytosine epigenome in the globally important marine, nitrogen-fixing cyanobacterium Trichodesmium

Cytosine methylation has been shown to regulate essential cellular processes and impact biological adaptation. Despite its evolutionary importance, only a handful of bacterial, genome-wide cytosine studies have been conducted, with none for marine bacteria. Here, we examine the genome-wide, C⁵ -Methyl-cytosine (m5C) methylome and its correlation to global transcription in the marine nitrogen-fixing cyanobacterium Trichodesmium. We characterize genome-wide methylation and highlight conserved motifs across three Trichodesmium isolates and two Trichodesmium metagenomes, thereby identifying highly conserved, novel genomic signatures of potential gene regulation in Trichodesmium. Certain gene bodies with the highest methylation levels correlate with lower expression levels. Several methylated motifs were highly conserved across spatiotemporally separated Trichodesmium isolates, thereby elucidating biogeographically conserved methylation potential. These motifs were also highly conserved in Trichodesmium metagenomic samples from natural populations suggesting them to be potential in situ markers of m5C methylation. Using these data, we highlight predicted roles of cytosine methylation in global cellular metabolism providing evidence for a 'core' m5C methylome spanning different ocean regions. These results provide important insights into the m5C methylation landscape and its biogeochemical implications in an important marine N₂ -fixer, as well as advancing evolutionary theory examining methylation influences on adaptation.

Minimal and RNA-free RNase P in Aquifex aeolicus

RNase P is an essential tRNA-processing enzyme in all domains of life. We identified an unknown type of protein-only RNase P in the hyperthermophilic bacterium Aquifex aeolicus: Without an RNA subunit and the smallest of its kind, the 23-kDa polypeptide comprises a metallonuclease domain only. The protein has RNase P activity in vitro and rescued the growth of Escherichia coli and Saccharomyces cerevisiae strains with inactivations of their more complex and larger endogenous ribonucleoprotein RNase P. Homologs of Aquifex RNase P (HARP) were identified in many Archaea and some Bacteria, of which all Archaea and most Bacteria also encode an RNA-based RNase P; activity of both RNase P forms from the same bacterium or archaeon could be verified in two selected cases. Bioinformatic analyses suggest that A. aeolicus and related Aquificaceae likely acquired HARP by horizontal gene transfer from an archaeon.

Comparative Genomics of Ralstonia solanacearum Identifies Candidate Genes Associated with Cool Virulence

Strains of the Ralstonia solanacearum species complex in the phylotype IIB group are capable of causing Bacterial Wilt disease in potato and tomato at temperatures lower than 24°C. The capability of these strains to survive and to incite infection at temperatures colder than their normally tropical boundaries represents a threat to United States agriculture in temperate regions. In this work, we used a comparative genomics approach to identify orthologous genes linked to the lower temperature virulence phenotype. Six R. solanacearum cool virulent (CV) strains were compared to six strains non-pathogenic at low temperature (NPLT). CV strains can cause Bacterial Wilt symptoms at temperatures below 24°C, while NPLT cannot. Four R. solanacearum strains were sequenced for this work in order to complete the comparison. An orthologous genes comparison identified 44 genes present only in CV strains and 19 genes present only in NPLT strains. Gene annotation revealed a high percentage of genes compared with whole genomes in the transcriptional regulator and transport categories. A single nucleotide polymorphism (SNP) analysis identified 265 genes containing conserved non-synonymous SNPs in CV strains. Ten genes in the pathogenicity category were identified in this group. Comparisons of type 3 secretion system, type 6 secretion system (T6SS) clusters, and associated effectors did not indicate a correlation with the CV phenotype except for one T6SS VGR effector potentially associated with the CV phenotype. This is the first R. solanacearum genomic comparative analysis of multiple strains with different temperature related virulence. The candidate genes identified by this comparison are potential factors involved in virulence at low temperatures that need to be investigated. The high percentage of transcriptional regulators among the genes present only in CV strains supports the hypothesis that temperature dependent regulation of virulence genes explains the differential virulence phenotype at low temperatures. This comparison contributes to find new possible connections of temperature dependent virulence to the previously described complex regulatory system involving quorum-sensing, phenotype conversion (phcA), acyl-HSL production and responses to SA. It also added novel candidate T6SS effectors and useful detailed information about the T6SS in R. solanacearum.

Bioinformatics comparisons of RNA-binding proteins of pathogenic and non-pathogenic Escherichia coli strains reveal novel virulence factors

Background

Pathogenic bacteria have evolved various strategies to counteract host defences. They are also exposed to environments that are undergoing constant changes. Hence, in order to survive, bacteria must adapt themselves to the changing environmental conditions by performing regulations at the transcriptional and/or post-transcriptional levels. Roles of RNA-binding proteins (RBPs) as virulence factors have been very well studied. Here, we have used a sequence search-based method to compare and contrast the proteomes of 16 pathogenic and three non-pathogenic E. coli strains as well as to obtain a global picture of the RBP landscape (RBPome) in E. coli.

Results

Our results show that there are no significant differences in the percentage of RBPs encoded by the pathogenic and the non-pathogenic E. coli strains. The differences in the types of Pfam domains as well as Pfam RNA-binding domains, encoded by these two classes of E. coli strains, are also insignificant. The complete and distinct RBPome of E. coli has been established by studying all known E. coli strains till date. We have also identified RBPs that are exclusive to pathogenic strains, and most of them can be exploited as drug targets since they appear to be non-homologous to their human host proteins. Many of these pathogen-specific proteins were uncharacterised and their identities could be resolved on the basis of sequence homology searches with known proteins. Detailed structural modelling, molecular dynamics simulations and sequence comparisons have been pursued for selected examples to understand differences in stability and RNA-binding.

Conclusions

The approach used in this paper to cross-compare proteomes of pathogenic and non-pathogenic strains may also be extended to other bacterial or even eukaryotic proteomes to understand interesting differences in their RBPomes. The pathogen-specific RBPs reported in this study, may also be taken up further for clinical trials and/or experimental validations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4045-3) contains supplementary material, which is available to authorized users.

PaperBLAST: Text Mining Papers for Information about Homologs

With the recent explosion of genome sequencing data, there are now millions of uncharacterized proteins. If a scientist becomes interested in one of these proteins, it can be very difficult to find information as to its likely function. Often a protein whose sequence is similar, and which is likely to have a similar function, has been studied already, but this information is not available in any database. To help find articles about similar proteins, PaperBLAST searches the full text of scientific articles for protein identifiers or gene identifiers, and it links these articles to protein sequences. Then, given a protein of interest, it can quickly find similar proteins in its database by using standard software (BLAST), and it can show snippets of text from relevant papers. We hope that PaperBLAST will make it easier for biologists to predict proteins’ functions.

Large-scale genome sequencing has identified millions of protein-coding genes whose function is unknown. Many of these proteins are similar to characterized proteins from other organisms, but much of this information is missing from annotation databases and is hidden in the scientific literature. To make this information accessible, PaperBLAST uses EuropePMC to search the full text of scientific articles for references to genes. PaperBLAST also takes advantage of curated resources (Swiss-Prot, GeneRIF, and EcoCyc) that link protein sequences to scientific articles. PaperBLAST’s database includes over 700,000 scientific articles that mention over 400,000 different proteins. Given a protein of interest, PaperBLAST quickly finds similar proteins that are discussed in the literature and presents snippets of text from relevant articles or from the curators. PaperBLAST is available at http://papers.genomics.lbl.gov/.

IMPORTANCE With the recent explosion of genome sequencing data, there are now millions of uncharacterized proteins. If a scientist becomes interested in one of these proteins, it can be very difficult to find information as to its likely function. Often a protein whose sequence is similar, and which is likely to have a similar function, has been studied already, but this information is not available in any database. To help find articles about similar proteins, PaperBLAST searches the full text of scientific articles for protein identifiers or gene identifiers, and it links these articles to protein sequences. Then, given a protein of interest, it can quickly find similar proteins in its database by using standard software (BLAST), and it can show snippets of text from relevant papers. We hope that PaperBLAST will make it easier for biologists to predict proteins’ functions.

Draft Genome Sequence of Lawsonia intracellularis Strain E40504, Isolated from a Horse Diagnosed with Equine Proliferative Enteropathy

Reported herein is the draft genome sequence of equine-origin Lawsonia intracellularis strain E40504, an obligate intracellular bacterium and the etiological agent of equine proliferative enteropathy. The 1.69-Mb draft genome sequence includes 1,380 protein-coding genes and 49 RNA genes, and it lacks a genomic island reported in swine-origin L. intracellularis strain PHE/MN1-00.

High precision genome sequencing of engineered Gluconobacter oxydans 621H by combining long nanopore and short accurate Illumina reads

State of the art and novel high-throughput DNA sequencing technologies enable fascinating opportunities and applications in the life sciences including microbial genomics. Short high-quality read data already enable not only microbial genome sequencing, yet can be inadequately to solve problems in genome assemblies and for the analysis of structural variants, especially in engineered microbial cell factories. Single-molecule real-time sequencing technologies generating long reads promise to solve such assembly problems. In our study, we wanted to increase the average read length of long nanopore reads with R9 chemistry and conducted a hybrid approach for the analysis of structural variants to check the genome stability of a recombinant Gluconobacter oxydans 621H strain (IK003.1) engineered for improved growth. Therefore we combined accurate Illumina sequencing technology and low-cost single-molecule nanopore sequencing using the MinION^® device from Oxford Nanopore. In our hybrid approach with a modified library protocol we could increase the average size of nanopore 2D reads to about 18.9kb. Combining the long MinION nanopore reads with the high quality short Illumina reads enabled the assembly of the engineered chromosome into a single contig and comprehensive detection and clarification of 7 structural variants including all three known genetically engineered modifications. We found the genome of IK003.1 was stable over 70 generations of strain handling including 28h of process time in a bioreactor. The long read data revealed a novel 1420 bp transposon-flanked and ORF-containing sequence which was hitherto unknown in the G. oxydans 621H reference. Further analysis and genome sequencing showed that this region is already present in G. oxydans 621H wild-type strains. Our data of G. oxydans 621H wild-type DNA from different resources also revealed in 73 annotated coding sequences about 91 uniform nucleotide differences including InDels. Together, our results contribute to an improved high quality genome reference for G. oxydans 621H which is available via ENA accession PRJEB18739.

Typing and Characterization of Bacteria Using Bottom-up Tandem Mass Spectrometry Proteomics

Methods for rapid and reliable microbial identification are essential in modern healthcare. The ability to detect and correctly identify pathogenic species and their resistance phenotype is necessary for accurate diagnosis and efficient treatment of infectious diseases. Bottom-up tandem mass spectrometry (MS) proteomics enables rapid characterization of large parts of the expressed genes of microorganisms. However, the generated data are highly fragmented, making downstream analyses complex. Here we present TCUP, a new computational method for typing and characterizing bacteria using proteomics data from bottom-up tandem MS. TCUP compares the generated protein sequence data to reference databases and automatically finds peptides suitable for characterization of taxonomic composition and identification of expressed antimicrobial resistance genes. TCUP was evaluated using several clinically relevant bacterial species (Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, Streptococcus pneumoniae, Moraxella catarrhalis, and Haemophilus influenzae), using both simulated data generated by in silico peptide digestion and experimental proteomics data generated by liquid chromatography-tandem mass spectrometry (MS/MS). The results showed that TCUP performs correct peptide classifications at rates between 90.3 and 98.5% at the species level. The method was also able to estimate the relative abundances of individual species in mixed cultures. Furthermore, TCUP could identify expressed β-lactamases in an extended spectrum β-lactamase-producing (ESBL) E. coli strain, even when the strain was cultivated in the absence of antibiotics. Finally, TCUP is computationally efficient, easy to integrate in existing bioinformatics workflows, and freely available under an open source license for both Windows and Linux environments.