Whole-Genome Sequencing and Virulome Analysis of Escherichia coli Isolated from New Zealand Environments of Contrasting Observed Land Use.

Generic Escherichia coli is commonly used as an indicator of fecal contamination to assess water quality and human health risk. Where measured E. coli exceedances occur, the presence of other pathogenic microorganisms, such as Shiga toxin-producing E. coli (STEC), is assumed, but confirmatory data are lacking. Putative E. coli isolates (n = 709) were isolated from water, sediment, soil, periphyton, and feces samples (n = 189) from five sites representing native forest and agricultural environments. Ten E. coli isolates (1.41%) were stx2 positive, 19 (2.7%) were eae positive, and stx1-positive isolates were absent. At the sample level, stx2-positive E. coli (5 of 189, 2.6%) and eae-positive isolates (16 of 189, 8.5%) were rare. Using real-time PCR, these STEC-associated virulence factors were determined to be more prevalent in sample enrichments (stx1, 23.9%; stx2, 31.4%; eae, 53.7%) and positively correlated with generic E. coli isolate numbers (P < 0.05) determined using culture-based methods. Whole-genome sequencing (WGS) was undertaken on a subset of 238 isolates with assemblies representing seven E. coli phylogroups (A, B1, B2, C, D, E, and F), 22 Escherichia marmotae isolates, and 1 Escherichia ruysiae isolate. Virulence factors, including those from extraintestinal pathogenic E. coli, were extremely diverse in isolates from the different locations and were more common in phylogroup B2. Analysis of the virulome from WGS data permitted the identification of gene repertoires that may be involved in environmental fitness and broadly align with phylogroup. Although recovery of STEC isolates was low, our molecular data indicate that they are likely to be widely present in environmental samples containing diverse E. coli phylogroups. IMPORTANCE This study takes a systematic sampling approach to assess the public health risk of Escherichia coli recovered from freshwater sites within forest and farmland. The New Zealand landscape is dominated by livestock farming, and previous work has demonstrated that "recreational exposure to water" is a risk factor for human infection by Shiga toxin-producing Escherichia coli (STEC). Though STEC isolates were rarely isolated from water samples, STEC-associated virulence factors were identified more commonly from water sample culture enrichments and were associated with increased generic E. coli concentrations. Whole-genome sequencing data from both E. coli and newly described Escherichia spp. demonstrated the presence of virulence factors from E. coli pathotypes, including extraintestinal pathogenic E. coli. This has significance for understanding and interpreting the potential health risk from E. coli where water quality is poor and suggests a role of virulence factors in survival and persistence of E. coli and Escherichia spp.

Genomic insights into the physiology of Quinella, an iconic uncultured rumen bacterium.

Quinella is a genus of iconic rumen bacteria first reported in 1913. There are no cultures of these bacteria, and information on their physiology is scarce and contradictory. Increased abundance of Quinella was previously found in the rumens of some sheep that emit low amounts of methane (CH4) relative to their feed intake, but whether Quinella contributes to low CH4 emissions is not known. Here, we concentrate Quinella cells from sheep rumen contents, extract and sequence DNA, and reconstruct Quinella genomes that are >90% complete with as little as 0.20% contamination. Bioinformatic analyses of the encoded proteins indicate that lactate and propionate formation are major fermentation pathways. The presence of a gene encoding a potential uptake hydrogenase suggests that Quinella might be able to use free hydrogen (H2). None of the inferred metabolic pathways is predicted to produce H2, a major precursor of CH4, which is consistent with the lower CH4 emissions from those sheep with high abundances of this bacterium.

Occurrence and expression of genes encoding methyl-compound production in rumen bacteria.

BACKGROUND: Digestive processes in the rumen lead to the release of methyl-compounds, mainly methanol and methylamines, which are used by methyltrophic methanogens to form methane, an important agricultural greenhouse gas. Methylamines are produced from plant phosphatidylcholine degradation, by choline trimethylamine lyase, while methanol comes from demethoxylation of dietary pectins via pectin methylesterase activity. We have screened rumen metagenomic and metatranscriptomic datasets, metagenome assembled genomes, and the Hungate1000 genomes to identify organisms capable of producing methyl-compounds. We also describe the enrichment of pectin-degrading and methane-forming microbes from sheep rumen contents and the analysis of their genomes via metagenomic assembly. RESULTS: Screens of metagenomic data using the protein domains of choline trimethylamine lyase (CutC), and activator protein (CutD) found good matches only to Olsenella umbonata and to Caecibacter, while the Hungate1000 genomes and metagenome assembled genomes from the cattle rumen found bacteria within the phyla Actinobacteria, Firmicutes and Proteobacteria. The cutC and cutD genes clustered with genes that encode structural components of bacterial microcompartment proteins. Prevotella was the dominant genus encoding pectin methyl esterases, with smaller numbers of sequences identified from other fibre-degrading rumen bacteria. Some large pectin methyl esterases (> 2100 aa) were found to be encoded in Butyrivibrio genomes. The pectin-utilising, methane-producing consortium was composed of (i) a putative pectin-degrading bacterium (phylum Tenericutes, class Mollicutes), (ii) a galacturonate-using Sphaerochaeta sp. predicted to produce acetate, lactate, and ethanol, and (iii) a methylotrophic methanogen, Methanosphaera sp., with the ability to form methane via a primary ethanol-dependent, hydrogen-independent, methanogenesis pathway. CONCLUSIONS: The main bacteria that produce methyl-compounds have been identified in ruminants. Their enzymatic activities can now be targeted with the aim of finding ways to reduce the supply of methyl-compound substrates to methanogens, and thereby limit methylotrophic methanogenesis in the rumen.

Gene and transcript abundances of bacterial type III secretion systems from the rumen microbiome are correlated with methane yield in sheep.

BACKGROUND: Ruminants are important contributors to global methane emissions via microbial fermentation in their reticulo-rumens. This study is part of a larger program, characterising the rumen microbiomes of sheep which vary naturally in methane yield (g CH4/kg DM/day) and aims to define differences in microbial communities, and in gene and transcript abundances that can explain the animal methane phenotype. METHODS: Rumen microbiome metagenomic and metatranscriptomic data were analysed by Gene Set Enrichment, sparse partial least squares regression and the Wilcoxon Rank Sum test to estimate correlations between specific KEGG bacterial pathways/genes and high methane yield in sheep. KEGG genes enriched in high methane yield sheep were reassembled from raw reads and existing contigs and analysed by MEGAN to predict their phylogenetic origin. Protein coding sequences from Succinivibrio dextrinosolvens strains were analysed using Effective DB to predict bacterial type III secreted proteins. The effect of S. dextrinosolvens strain H5 growth on methane formation by rumen methanogens was explored using co-cultures. RESULTS: Detailed analysis of the rumen microbiomes of high methane yield sheep shows that gene and transcript abundances of bacterial type III secretion system genes are positively correlated with methane yield in sheep. Most of the bacterial type III secretion system genes could not be assigned to a particular bacterial group, but several genes were affiliated with the genus Succinivibrio, and searches of bacterial genome sequences found that strains of S. dextrinosolvens were part of a small group of rumen bacteria that encode this type of secretion system. In co-culture experiments, S. dextrinosolvens strain H5 showed a growth-enhancing effect on a methanogen belonging to the order Methanomassiliicoccales, and inhibition of a representative of the Methanobrevibacter gottschalkii clade. CONCLUSIONS: This is the first report of bacterial type III secretion system genes being associated with high methane emissions in ruminants, and identifies these secretions systems as potential new targets for methane mitigation research. The effects of S. dextrinosolvens on the growth of rumen methanogens in co-cultures indicate that bacteria-methanogen interactions are important modulators of methane production in ruminant animals.

Rumen metagenome and metatranscriptome analyses of low methane yield sheep reveals a Sharpea-enriched microbiome characterised by lactic acid formation and utilisation.

BACKGROUND: Enteric fermentation by farmed ruminant animals is a major source of methane and constitutes the second largest anthropogenic contributor to global warming. Reducing methane emissions from ruminants is needed to ensure sustainable animal production in the future. Methane yield varies naturally in sheep and is a heritable trait that can be used to select animals that yield less methane per unit of feed eaten. We previously demonstrated elevated expression of hydrogenotrophic methanogenesis pathway genes of methanogenic archaea in the rumens of high methane yield (HMY) sheep compared to their low methane yield (LMY) counterparts. Methane production in the rumen is strongly connected to microbial hydrogen production through fermentation processes. In this study, we investigate the contribution that rumen bacteria make to methane yield phenotypes in sheep. RESULTS: Using deep sequence metagenome and metatranscriptome datasets in combination with 16S rRNA gene amplicon sequencing from HMY and LMY sheep, we show enrichment of lactate-producing Sharpea spp. in LMY sheep bacterial communities. Increased gene and transcript abundances for sugar import and utilisation and production of lactate, propionate and butyrate were also observed in LMY animals. Sharpea azabuensis and Megasphaera spp. act as important drivers of lactate production and utilisation according to phylogenetic analysis and read mappings. CONCLUSIONS: Our findings show that the rumen microbiome in LMY animals supports a rapid heterofermentative growth, leading to lactate production. We postulate that lactate is subsequently metabolised mainly to butyrate in LMY animals, producing 2 mol of hydrogen and 0.5 mol of methane per mol hexose, which represents 24 % less than the 0.66 mol of methane formed from the 2.66 mol of hydrogen produced if hexose fermentation was directly to acetate and butyrate. These findings are consistent with the theory that a smaller rumen size with a higher turnover rate, where rapid heterofermentative growth would be an advantage, results in lower hydrogen production and lower methane formation. Together with previous methanogen gene expression data, this builds a strong concept of how animal traits and microbial communities shape the methane phenotype in sheep.

Quantification of bacterial and archaeal symbionts in high and low microbial abundance sponges using real-time PCR.

In spite of considerable insights into the microbial diversity of marine sponges, quantitative information on microbial abundances and community composition remains scarce. Here, we established qPCR assays for the specific quantification of four bacterial phyla of representative sponge symbionts as well as the kingdoms Eubacteria and Archaea. We could show that the 16S rRNA gene numbers of Archaea, Chloroflexi, and the candidate phylum Poribacteria were 4-6 orders of magnitude higher in high microbial abundance (HMA) than in low microbial abundance (LMA) sponges and that actinobacterial 16S rRNA gene numbers were 1-2 orders higher in HMA over LMA sponges, while those for Cyanobacteria were stable between HMA and LMA sponges. Fluorescence in situ hybridization of Aplysina aerophoba tissue sections confirmed the numerical dominance of Chloroflexi, which was followed by Poribacteria. Archaeal and actinobacterial cells were detected in much lower numbers. By use of fluorescence-activated cell sorting as a primer- and probe-independent approach, the dominance of Chloroflexi, Proteobacteria, and Poribacteria in A. aerophoba was confirmed. Our study provides new quantitative insights into the microbiology of sponges and contributes to a better understanding of the HMA/LMA dichotomy.

The candidate phylum Poribacteria by single-cell genomics: new insights into phylogeny, cell-compartmentation, eukaryote-like repeat proteins, and other genomic features.

The candidate phylum Poribacteria is one of the most dominant and widespread members of the microbial communities residing within marine sponges. Cell compartmentalization had been postulated along with their discovery about a decade ago and their phylogenetic association to the Planctomycetes, Verrucomicrobia, Chlamydiae superphylum was proposed soon thereafter. In the present study we revised these features based on genomic data obtained from six poribacterial single cells. We propose that Poribacteria form a distinct monophyletic phylum contiguous to the PVC superphylum together with other candidate phyla. Our genomic analyses supported the possibility of cell compartmentalization in form of bacterial microcompartments. Further analyses of eukaryote-like protein domains stressed the importance of such proteins with features including tetratricopeptide repeats, leucin rich repeats as well as low density lipoproteins receptor repeats, the latter of which are reported here for the first time from a sponge symbiont. Finally, examining the most abundant protein domain family on poribacterial genomes revealed diverse phyH family proteins, some of which may be related to dissolved organic posphorus uptake.

Bacterial community profiles in low microbial abundance sponges.

It has long been recognized that sponges differ in the abundance of associated microorganisms, and they are therefore termed either 'low microbial abundance' (LMA) or 'high microbial abundance' (HMA) sponges. Many previous studies concentrated on the dense microbial communities in HMA sponges, whereas little is known about microorganisms in LMA sponges. Here, two LMA sponges from the Red Sea, two from the Caribbean and one from the South Pacific were investigated. With up to only five bacterial phyla per sponge, all LMA sponges showed lower phylum-level diversity than typical HMA sponges. Interestingly, each LMA sponge was dominated by a large clade within either Cyanobacteria or different classes of Proteobacteria. The overall similarity of bacterial communities among LMA sponges determined by operational taxonomic unit and UniFrac analysis was low. Also the number of sponge-specific clusters, which indicate bacteria specifically associated with sponges and which are numerous in HMA sponges, was low. A biogeographical or host-dependent distribution pattern was not observed. In conclusion, bacterial community profiles of LMA sponges are clearly different from profiles of HMA sponges and, remarkably, each LMA sponge seems to harbour its own unique bacterial community.

Single-cell genomics reveals complex carbohydrate degradation patterns in poribacterial symbionts of marine sponges.

Many marine sponges are hosts to dense and phylogenetically diverse microbial communities that are located in the extracellular matrix of the animal. The candidate phylum Poribacteria is a predominant member of the sponge microbiome and its representatives are nearly exclusively found in sponges. Here we used single-cell genomics to obtain comprehensive insights into the metabolic potential of individual poribacterial cells representing three distinct phylogenetic groups within Poribacteria. Genome sizes were up to 5.4 Mbp and genome coverage was as high as 98.5%. Common features of the poribacterial genomes indicated that heterotrophy is likely to be of importance for this bacterial candidate phylum. Carbohydrate-active enzyme database screening and further detailed analysis of carbohydrate metabolism suggested the ability to degrade diverse carbohydrate sources likely originating from seawater and from the host itself. The presence of uronic acid degradation pathways as well as several specific sulfatases provides strong support that Poribacteria degrade glycosaminoglycan chains of proteoglycans, which are important components of the sponge host matrix. Dominant glycoside hydrolase families further suggest degradation of other glycoproteins in the host matrix. We therefore propose that Poribacteria are well adapted to an existence in the sponge extracellular matrix. Poribacteria may be viewed as efficient scavengers and recyclers of a particular suite of carbon compounds that are unique to sponges as microbial ecosystems.

Exploring symbioses by single-cell genomics.

Single-cell genomics has advanced the field of microbiology from the analysis of microbial metagenomes where information is "drowning in a sea of sequences," to recognizing each microbial cell as a separate and unique entity. Single-cell genomics employs Phi29 polymerase-mediated whole-genome amplification to yield microgram-range genomic DNA from single microbial cells. This method has now been applied to a handful of symbiotic systems, including bacterial symbionts of marine sponges, insects (grasshoppers, termites), and vertebrates (mouse, human). In each case, novel insights were obtained into the functional genomic repertoire of the bacterial partner, which, in turn, led to an improved understanding of the corresponding host. Single-cell genomics is particularly valuable when dealing with uncultivated microorganisms, as is still the case for many bacterial symbionts. In this review, we explore the power of single-cell genomics for symbiosis research and highlight recent insights into the symbiotic systems that were obtained by this approach.

Single-cell genomics reveals the lifestyle of Poribacteria, a candidate phylum symbiotically associated with marine sponges.

In this study, we present a single-cell genomics approach for the functional characterization of the candidate phylum Poribacteria, members of which are nearly exclusively found in marine sponges. The microbial consortia of the Mediterranean sponge Aplysina aerophoba were singularized by fluorescence-activated cell sorting, and individual microbial cells were subjected to phi29 polymerase-mediated 'whole-genome amplification'. Pyrosequencing of a single amplified genome (SAG) derived from a member of the Poribacteria resulted in nearly 1.6 Mb of genomic information distributed among 554 contigs analyzed in this study. Approximately two-third of the poribacterial genome was sequenced. Our findings shed light on the functional properties and lifestyle of a possibly ancient bacterial symbiont of marine sponges. The Poribacteria are mixotrophic bacteria with autotrophic CO(2)-fixation capacities through the Wood-Ljungdahl pathway. The cell wall is of Gram-negative origin. The Poribacteria produce at least two polyketide synthases (PKSs), one of which is the sponge-specific Sup-type PKS. Several putative symbiosis factors such as adhesins (bacterial Ig-like domains, lamininin G domain proteins), adhesin-related proteins (ankyrin, fibronectin type III) and tetratrico peptide repeat domain-encoding proteins were identified, which might be involved in mediating sponge-microbe interactions. The discovery of genes coding for 24-isopropyl steroids implies that certain fossil biomarkers used to date the origins of metazoan life on earth may possibly be of poribacterial origin. Single-cell genomic approaches, such as those shown herein, contribute to a better understanding of beneficial microbial consortia, of which most members are, because of the lack of cultivation, inaccessible by conventional techniques.

The pathology of sponge orange band disease affecting the Caribbean barrel sponge Xestospongia muta.

The aim of this study was to examine sponge orange band (SOB) disease affecting the prominent Caribbean sponge Xestospongia muta. Scanning and transmission electron microscopy revealed that SOB is accompanied by the massive destruction of the pinacoderm. Chlorophyll a content and the main secondary metabolites, tetrahydrofurans, characteristic of X. muta, were significantly lower in bleached than in healthy tissues. Denaturing gradient gel electrophoresis using cyanobacteria-specific 16S rRNA gene primers revealed a distinct shift from the Synechococcus/Prochlorococcus clade of sponge symbionts towards several clades of unspecific cyanobacteria, including lineages associated with coral disease (i.e. Leptolyngbya sp.). Underwater infection experiments were conducted by transplanting bleached cores into healthy individuals, but revealed no signs of SOB development. This study provided no evidence for the involvement of a specific microbial pathogen as an etiologic agent of disease; hence, the cause of SOB disease in X. muta remains unidentified.

Activity profiles for marine sponge-associated bacteria obtained by 16S rRNA vs 16S rRNA gene comparisons.

The phylogenetic diversity of microorganisms in marine sponges is becoming increasingly well described, yet relatively little is known about the activities of these symbionts. Given the seemingly favourable environment provided to microbes by their sponge hosts, as indicated by the extraordinarily high abundance of sponge symbionts, we hypothesized that the majority of sponge-associated bacteria are active in situ. To test this hypothesis we compared, for the first time in sponges, 16S rRNA gene- vs 16S rRNA-derived bacterial community profiles to gain insights into symbiont composition and activity, respectively. Clone libraries revealed a highly diverse bacterial community in Ancorina alata, and a much lower diversity in Polymastia sp., which were identified by electron microscopy as a high- and a low-microbial abundance sponge, respectively. Substantial overlap between DNA and RNA libraries was evident at both phylum and phylotype levels, indicating in situ activity for a large fraction of sponge-associated bacteria. This active fraction included uncultivated, sponge-specific lineages within, for example, Actinobacteria, Chloroflexi and Gemmatimonadetes. This study shows the potential of RNA vs DNA comparisons based on the 16S rRNA gene to provide insights into the activity of sponge-associated microorganisms.