William L. Trimble

Intestinal Tissues Induce an SNP Mutation in Pseudomonas aeruginosa That Enhances Its Virulence: Possible Role in Anastomotic Leak

The most feared complication following intestinal resection is anastomotic leakage. In high risk areas (esophagus/rectum) where neoadjuvant chemoradiation is used, the incidence of anastomotic leaks remains unacceptably high (∼10%) even when performed by specialist surgeons in high volume centers. The aims of this study were to test the hypothesis that anastomotic leakage develops when pathogens colonizing anastomotic sites become in vivo transformed to express a tissue destroying phenotype. We developed a novel model of anastomotic leak in which rats were exposed to pre-operative radiation as in cancer surgery, underwent distal colon resection and then were intestinally inoculated with Pseudomonas aeruginosa, a common colonizer of the radiated intestine. Results demonstrated that intestinal tissues exposed to preoperative radiation developed a significant incidence of anastomotic leak (>60%; p<0.01) when colonized by P. aeruginosa compared to radiated tissues alone (0%). Phenotype analysis comparing the original inoculating strain (MPAO1- termed P1) and the strain retrieved from leaking anastomotic tissues (termed P2) demonstrated that P2 was altered in pyocyanin production and displayed enhanced collagenase activity, high swarming motility, and a destructive phenotype against cultured intestinal epithelial cells (i.e. apoptosis, barrier function, cytolysis). Comparative genotype analysis between P1 and P2 revealed a single nucleotide polymorphism (SNP) mutation in the mexT gene that led to a stop codon resulting in a non-functional truncated protein. Replacement of the mutated mexT gene in P2 with mexT from the original parental strain P1 led to reversion of P2 to the P1 phenotype. No spontaneous transformation was detected during 20 passages in TSB media. Use of a novel virulence suppressing compound PEG/Pi prevented P. aeruginosa transformation to the tissue destructive phenotype and prevented anastomotic leak in rats. This work demonstrates that in vivo transformation of microbial pathogens to a tissue destroying phenotype may have important implications in the pathogenesis of anastomotic leak.

A Platform-Independent Method for Detecting Errors in Metagenomic Sequencing Data: DRISEE

We provide a novel method, DRISEE (duplicate read inferred sequencing error estimation), to assess sequencing quality (alternatively referred to as “noise” or “error”) within and/or between sequencing samples. DRISEE provides positional error estimates that can be used to inform read trimming within a sample. It also provides global (whole sample) error estimates that can be used to identify samples with high or varying levels of sequencing error that may confound downstream analyses, particularly in the case of studies that utilize data from multiple sequencing samples. For shotgun metagenomic data, we believe that DRISEE provides estimates of sequencing error that are more accurate and less constrained by technical limitations than existing methods that rely on reference genomes or the use of scores (e.g. Phred). Here, DRISEE is applied to (non amplicon) data sets from both the 454 and Illumina platforms. The DRISEE error estimate is obtained by analyzing sets of artifactual duplicate reads (ADRs), a known by-product of both sequencing platforms. We present DRISEE as an open-source, platform-independent method to assess sequencing error in shotgun metagenomic data, and utilize it to discover previously uncharacterized error in de novo sequence data from the 454 and Illumina sequencing platforms.

The MG-RAST metagenomics database and portal in 2015

MG-RAST (http://metagenomics.anl.gov) is an open-submission data portal for processing, analyzing, sharing and disseminating metagenomic datasets. The system currently hosts over 200 000 datasets and is continuously updated. The volume of submissions has increased 4-fold over the past 24 months, now averaging 4 terabasepairs per month. In addition to several new features, we report changes to the analysis workflow and the technologies used to scale the pipeline up to the required throughput levels. To show possible uses for the data from MG-RAST, we present several examples integrating data and analyses from MG-RAST into popular third-party analysis tools or sequence alignment tools.

The complete genome sequence for putative H₂- and S-oxidizer Candidatus Sulfuricurvum sp., assembled de novo from an aquifer-derived metagenome.

We reconstructed the complete 2.4 Mb-long genome of a previously uncultivated epsilonproteobacterium, Candidatus Sulfuricurvum sp. RIFRC-1, via assembly of short-read shotgun metagenomic data using a complexity reduction approach. Genome-based comparisons indicate the bacterium is a novel species within the Sulfuricurvum genus, which contains one cultivated representative, S. kujiense. Divergence between the species appears due in part to extensive genomic rearrangements, gene loss and chromosomal versus plasmid encoding of certain (respiratory) genes by RIFRC-1. Deoxyribonucleic acid for the genome was obtained from terrestrial aquifer sediment, in which RIFRC-1 comprised ∼ 47% of the bacterial community. Genomic evidence suggests RIFRC-1 is a chemolithoautotrophic diazotroph capable of deriving energy for growth by microaerobic or nitrate-/nitric oxide-dependent oxidation of S°, sulfide or sulfite or H₂oxidation. Carbon may be fixed via the reductive tricarboxylic acid cycle. Consistent with these physiological attributes, the local aquifer was microoxic with small concentrations of available nitrate, small but elevated concentrations of reduced sulfur and NH(4)(+) /NH₃-limited. Additionally, various mechanisms for heavy metal and metalloid tolerance and virulence point to a lifestyle well-adapted for metal(loid)-rich environments and a shared evolutionary past with pathogenic Epsilonproteobacteria. Results expand upon recent findings highlighting the potential importance of sulfur and hydrogen metabolism in the terrestrial subsurface.

Short-read reading-frame predictors are not created equal: sequence error causes loss of signal

BackgroundGene prediction algorithms (or gene callers) are an essential tool for analyzing shotgun nucleic acid sequence data. Gene prediction is a ubiquitous step in sequence analysis pipelines; it reduces the volume of data by identifying the most likely reading frame for a fragment, permitting the out-of-frame translations to be ignored. In this study we evaluate five widely used ab initio gene-calling algorithms—FragGeneScan, MetaGeneAnnotator, MetaGeneMark, Orphelia, and Prodigal—for accuracy on short (75–1000 bp) fragments containing sequence error from previously published artificial data and “real” metagenomic datasets.ResultsWhile gene prediction tools have similar accuracies predicting genes on error-free fragments, in the presence of sequencing errors considerable differences between tools become evident. For error-containing short reads, FragGeneScan finds more prokaryotic coding regions than does MetaGeneAnnotator, MetaGeneMark, Orphelia, or Prodigal. This improved detection of genes in error-containing fragments, however, comes at the cost of much lower (50%) specificity and overprediction of genes in noncoding regions.ConclusionsAb initio gene callers offer a significant reduction in the computational burden of annotating individual nucleic acid reads and are used in many metagenomic annotation systems. For predicting reading frames on raw reads, we find the hidden Markov model approach in FragGeneScan is more sensitive than other gene prediction tools, while Prodigal, MGA, and MGM are better suited for higher-quality sequences such as assembled contigs.

Whole-Genome Sequence of the Human Pathogen Legionella pneumophila Serogroup 12 Strain 570-CO-H

We present the genomic sequence of the human pathogen Legionella pneumophila serogroup 12 strain 570-CO-H (ATCC 43290), a clinical isolate from the Colorado Department of Health, Denver, CO. This is the first example of a genome sequence of L. pneumophila from a serogroup other than serogroup 1. We highlight the similarities and differences relative to six genome sequences that have been reported for serogroup 1 strains.

A Metagenomics Portal for a Democratized Sequencing World

The democratized world of sequencing is leading to numerous data analysis challenges; MG-RAST addresses many of these challenges for diverse datasets, including amplicon datasets, shotgun metagenomes, and metatranscriptomes. The changes from version 2 to version 3 include the addition of a dedicated gene calling stage using FragGenescan, clustering of predicted proteins at 90% identity, and the use of BLAT for the computation of similarities. Together with changes in the underlying software infrastructure, this has enabled the dramatic scaling up of pipeline throughput while remaining on a limited hardware budget. The Web-based service allows upload, fully automated analysis, and visualization of results. As a result of the plummeting cost of sequencing and the readily available analytical power of MG-RAST, over 78,000 metagenomic datasets have been analyzed, with over 12,000 of them publicly available in MG-RAST.

Rapid quantification of sequence repeats to resolve the size, structure and contents of bacterial genomes.

A bstractBackgroundThe numerous classes of repeats often impede the assembly of genome sequences from the short reads provided by new sequencing technologies. We demonstrate a simple and rapid means to ascertain the repeat structure and total size of a bacterial or archaeal genome without the need for assembly by directly analyzing the abundances of distinct k-mers among reads.ResultsThe sensitivity of this procedure to resolve variation within a bacterial species is demonstrated: genome sizes and repeat structure of five environmental strains of E. coli from short Illumina reads were estimated by this method, and total genome sizes corresponded well with those obtained for the same strains by pulsed-field gel electrophoresis. In addition, this approach was applied to read-sets for completed genomes and shown to be accurate over a wide range of microbial genome sizes.ConclusionsApplication of these procedures, based solely on k-mer abundances in short read data sets, allows aspects of genome structure to be resolved that are not apparent from conventional short read assemblies. This knowledge of the repetitive content of genomes provides insights into genome evolution and diversity.

Eco-Evolutionary Dynamics of Episomes among Ecologically Cohesive Bacterial Populations

ABSTRACT Although plasmids and other episomes are recognized as key players in horizontal gene transfer among microbes, their diversity and dynamics among ecologically structured host populations in the wild remain poorly understood. Here, we show that natural populations of marine Vibrionaceae bacteria host large numbers of families of episomes, consisting of plasmids and a surprisingly high fraction of plasmid-like temperate phages. Episomes are unevenly distributed among host populations, and contrary to the notion that high-density communities in biofilms act as hot spots of gene transfer, we identified a strong bias for episomes to occur in free-living as opposed to particle-attached cells. Mapping of episomal families onto host phylogeny shows that, with the exception of all phage and a few plasmid families, most are of recent evolutionary origin and appear to have spread rapidly by horizontal transfer. Such high eco-evolutionary turnover is particularly surprising for plasmids that are, based on previously suggested categorization, putatively nontransmissible, indicating that this type of plasmid is indeed frequently transferred by currently unknown mechanisms. Finally, analysis of recent gene transfer among plasmids reveals a network of extensive exchange connecting nearly all episomes. Genes functioning in plasmid transfer and maintenance are frequently exchanged, suggesting that plasmids can be rapidly transformed from one category to another. The broad distribution of episomes among distantly related hosts and the observed promiscuous recombination patterns show how episomes can offer their hosts rapid assembly and dissemination of novel functions. IMPORTANCE Plasmids and other episomes are an integral part of bacterial biology in all environments, yet their study is heavily biased toward their role as vectors for antibiotic resistance genes. This study presents a comprehensive analysis of all episomes within several coexisting bacterial populations of Vibrionaceae from the coastal ocean and represents the largest-yet genomic survey of episomes from a single bacterial family. The host population framework allows analysis of the eco-evolutionary dynamics at unprecedented resolution, yielding several unexpected results. These include (i) discovery of novel, nonintegrative temperate phages, (ii) revision of a class of episomes, previously termed “nontransmissible,” as highly transmissible, and (iii) surprisingly high evolutionary turnover of episomes, manifest as frequent birth, spread, and loss. Plasmids and other episomes are an integral part of bacterial biology in all environments, yet their study is heavily biased toward their role as vectors for antibiotic resistance genes. This study presents a comprehensive analysis of all episomes within several coexisting bacterial populations of Vibrionaceae from the coastal ocean and represents the largest-yet genomic survey of episomes from a single bacterial family. The host population framework allows analysis of the eco-evolutionary dynamics at unprecedented resolution, yielding several unexpected results. These include (i) discovery of novel, nonintegrative temperate phages, (ii) revision of a class of episomes, previously termed “nontransmissible,” as highly transmissible, and (iii) surprisingly high evolutionary turnover of episomes, manifest as frequent birth, spread, and loss.

Draft Genome Sequence of Achromobacter piechaudii Strain HLE

Achromobacter piechaudii strain HLE is a betaproteobacterium (previously known as Alcaligenes faecalis) that was an early isolate with arsenite oxidase activity. This draft genome of 6.89 Mb is the second available genome for this species in the opportunistic pathogen Alcaligenaceae family.