How Much Do rRNA Gene Surveys Underestimate Extant Bacterial Diversity?

Potential routes of plastics biotransformation involving novel plastizymes revealed by global multi-omic analysis of plastic associated microbes

Despite increasing efforts across various disciplines, the fate, transport, and impact of synthetic plastics on the environment and public health remain poorly understood. To better elucidate the microbial ecology of plastic waste and its potential for biotransformation, we conducted a large-scale analysis of all publicly available meta-omic studies investigating plastics (n = 27) in the environment. Notably, we observed low prevalence of known plastic degraders throughout most environments, except for substantial enrichment in riverine systems. This indicates rivers may be a highly promising environment for discovery of novel plastic bioremediation products. Ocean samples associated with degrading plastics showed clear differentiation from non-degrading polymers, showing enrichment of novel putative biodegrading taxa in the degraded samples. Regarding plastisphere pathogenicity, we observed significant enrichment of antimicrobial resistance genes on plastics but not of virulence factors. Additionally, we report a co-occurrence network analysis of 10 + million proteins associated with the plastisphere. This analysis revealed a localized sub-region enriched with known and putative plastizymes-these may be useful for deeper investigation of nature's ability to biodegrade man-made plastics. Finally, the combined data from our meta-analysis was used to construct a publicly available database, the Plastics Meta-omic Database (PMDB)-accessible at plasticmdb.org. These data should aid in the integrated exploration of the microbial plastisphere and facilitate research efforts investigating the fate and bioremediation potential of environmental plastic waste.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Phenotypic and Genomic Characterization of Streptomyces pakalii sp. nov., a Novel Species with Anti-Biofilm and Anti-Quorum Sensing Activity in ESKAPE Bacteria

The increasing number of infections caused by antimicrobial multi-resistant microorganisms has led to the search for new microorganisms capable of producing novel antibiotics. This work proposes Streptomyces pakalii sp. nov. as a new member of the Streptomycetaceae family. The strain ENCB-J15 was isolated from the jungle soil in Palenque National Park, Chiapas, Mexico. The strain formed pale brown, dry, tough, and buried colonies in the agar with no diffusible pigment in GAE (glucose-asparagine-yeast extract) medium. Scanning electron micrographs showed typical mycelium with long chains of smooth and oval-shaped spores (3-10 m). The strain grew in all of the International Streptomyces Project (ISP)'s media at 28-37 °C with a pH of 6-9 and 0-10% NaCl. S. pakalii ENCB-J15 assimilated diverse carbon as well as organic and inorganic nitrogen sources. The strain also exhibited significant inhibitory activity against the prodigiosin synthesis of Serratia marcescens and the inhibition of the formation and destruction of biofilms of ESKAPE strains of Acinetobacter baumannii and Klebsiella pneumoniae. The draft genome sequencing of ENCB-J15 revealed a 7.6 Mb genome with a high G + C content (71.6%), 6833 total genes, and 6746 genes encoding putative proteins. A total of 26 accessory clusters of proteins associated with carbon sources and amino acid catabolism, DNA modification, and the antibiotic biosynthetic process were annotated. The 16S rRNA gene phylogeny, core-proteome phylogenomic tree, and virtual genome fingerprints support that S. pakalii ENCB-J15 is a new species related to Streptomyces badius and Streptomyces globisporus. Similarly, its average nucleotide identity (ANI) (96.4%), average amino acid identity (AAI) (96.06%), and virtual DNA-DNA hybridization (67.3%) provide evidence to recognize it as a new species. Comparative genomics revealed that S. pakalli and its closest related species maintain a well-conserved genomic synteny. This work proposes Streptomyces pakalii sp. nov. as a novel species that expresses anti-biofilm and anti-quorum sensing activities.

Microbial Diversity Biased Estimation Caused by Intragenomic Heterogeneity and Interspecific Conservation of 16S rRNA Genes

The 16S rRNA gene has been extensively used as a molecular marker to explore evolutionary relationships and profile microbial composition throughout various environments. Despite its convenience and prevalence, limitations are inevitable. Variable copy numbers, intragenomic heterogeneity, and low taxonomic resolution have caused biases in estimating microbial diversity. Here, analysis of 24,248 complete prokaryotic genomes indicated that the 16S rRNA gene copy number ranged from 1 to 37 in bacteria and 1 to 5 in archaea, and intragenomic heterogeneity was observed in 60% of prokaryotic genomes, most of which were below 1%. The overestimation of microbial diversity caused by intragenomic variation and the underestimation introduced by interspecific conservation were calculated when using full-length or partial 16S rRNA genes. Results showed that, at the 100% threshold, microbial diversity could be overestimated by as much as 156.5% when using the full-length gene. The V4 to V5 region-based analyses introduced the lowest overestimation rate (4.4%) but exhibited slightly lower species resolution than other variable regions under the 97% threshold. For different variable regions, appropriate thresholds rather than the canonical value 97% were proposed for minimizing the risk of splitting a single genome into multiple clusters and lumping together different species into the same cluster. This study has not only updated the 16S rRNA gene copy number and intragenomic variation information for the currently available prokaryotic genomes, but also elucidated the biases in estimating prokaryotic diversity with quantitative data, providing references for choosing amplified regions and clustering thresholds in microbial community surveys. IMPORTANCE Microbial diversity is typically analyzed using marker gene-based methods, of which 16S rRNA gene sequencing is the most widely used approach. However, obtaining an accurate estimation of microbial diversity remains a challenge, due to the intragenomic variation and low taxonomic resolution of 16S rRNA genes. Comprehensive examination of the bias in estimating such prokaryotic diversity using 16S rRNA genes within ever-increasing prokaryotic genomes highlights the importance of the choice of sequencing regions and clustering thresholds based on the specific research objectives.

Two novel species isolated from wheat rhizospheres in Serbia: Pseudomonas serbica sp. nov. and Pseudomonas serboccidentalis sp. nov

Pseudomonas strains IT-194P, IT-215P, IT-P366^T and IT-P374^T were isolated from the rhizospheres of wheat grown in soils sampled from different fields (some of them known to be disease-suppressive) located near Mionica, Serbia. Phylogenetic analysis of the 16S rRNA genes and of whole genome sequences showed that these strains belong to two potentially new species, one containing strains IT-P366^T and IT-194P and clustering (whole genome analysis) next to P. umsongensis DSM16611^T, and another species containing strains IT-P374^T and IT-215P and clustering next to P. koreensis LMG21318^T. Genome analysis confirmed the proposition of novel species, as ANI was below the threshold of 95% and dDDH below 70% for strains IT-P366^T (compared with P. umsongensis DSM16611^T) and IT-P374^T (compared with P. koreensis LMG21318^T). Unlike P. umsongensis DSM16611^T, strains of P. serbica can grow on D-mannitol, but not on pectin, D-galacturonic acid, L-galactonic acid lactone and α-hydroxybutyric acid. In contrary to P. koreensis LMG21318^T, strains of P. serboccidentalis can use sucrose, inosine and α-ketoglutaric acid (but not L-histidine) as carbon sources. Altogether, these results indicate the existence of two novel species for which we propose the names Pseudomonas serbica sp. nov., with the type strain IT-P366^T (=CFBP 9060 ^T = LMG 32732 ^T = EML 1791 ^T) and Pseudomonas serboccidentalis sp. nov., with the type strain IT-P374^T (=CFBP 9061 ^T = LMG 32734 ^T = EML 1792 ^T). Strains from this study presented a set of phytobeneficial functions modulating plant hormonal balance, plant nutrition and plant protection, suggesting a potential as Plant Growth-Promoting Rhizobacteria (PGPR).

Microcystis pangenome reveals cryptic diversity within and across morphospecies

Microcystis, a common harmful algal bloom (HAB) taxon, threatens water supplies and human health, yet species delimitation is contentious in this taxon, leading to challenges in research and management of this threat. Historical and common morphology-based classifications recognize multiple morphospecies, most with variable and diverse ecologies, while DNA sequence-based classifications indicate a single species with multiple ecotypes. To better delimit Microcystis species, we conducted a pangenome analysis of 122 genomes. Core- and non-core gene phylogenetic analyses placed 113 genomes into 23 monophyletic clusters containing at least two genomes. Overall, genome-related indices revealed that Microcystis contains at least 16 putative genospecies. Fifteen genospecies included at least one Microcystis aeruginosa morphospecies, and 10 genospecies included two or more morphospecies. This classification system will enable consistent taxonomic identification of Microcystis and thereby aid in resolving some of the complexities and controversies that have long characterized eco-evolutionary research and management of this important HAB taxon.

Characterization of Mycobacterium salfingeri sp. nov.: A novel nontuberculous mycobacteria isolated from a human wound infection

Nontuberculous mycobacteria (NTM) are environmental bacteria commonly found in soil and water in almost every part of the world. While usually non-pathogenic, they can cause acute respiratory and cutaneous infections under certain circumstances or in patients with underlying medical conditions. Contrary to members of the Mycobacterium tuberculosis complex, documented human-to-human transmissions of NTM have been rarely reported and most cases result from direct environmental exposure. Here we describe the identification of a new NTM species isolated from a hand laceration of a New York State patient after a fall. This new NTM forms rough, orange pigmented colonies and is naturally resistant to doxycycline and tobramycin. Whole genome analysis reveal no close relatives present in public databases, and our findings are in accordance with the recognition of a new taxonomic species of NTM. We propose the name Mycobacterium salfingeri sp. nov. for this new NTM representative. The type strain is 20-157661T (DSM = 113368T, BCCM = ITM 501207T).

A culture-based and culture-independent approach to the study of landfill leachate bacterial and archaeal communities

The landfill is a convenient and affordable method of municipal solid waste (MSW) management. Landfill leachate contains a heavy load of pollutants and pathogens. Discharge of untreated leachate is the leading cause of surface and groundwater contamination and a threat to public and environmental health. To develop an efficient leachate treatment technology, an in-depth understanding of landfill chemistry and microbiology is essential. In the present manuscript, we conducted a comparative study of three different landfill leachate samples using cultivation-based and culture-independent molecular studies. We cultivated 85 species of aerobic, anaerobic bacteria and archaea from leachate represented by a total of 200 strains using extensive culturomics approaches. Twelve out of 200 cultivated strains of bacteria showed very low 16S rRNA gene sequence similarity (84-98.6%) with their closest relatives and could be the potential novel taxa, the first time cultivated from leachate. Members of the six genera only have 2-5 representatives from past studies from other habitats but first time cultivated from leachate. In addition to bacteria, we also cultivated and characterized different groups of methanogenic archaea. Our chemistry data indicate that leachate is a highly stressed ecosystem with an assemblage of many toxic wastes like sulfur, zinc, mercury, chromium, etc. 16S rRNA gene-based amplicon analysis showed the dominance of (30-55%) methanogens and haloarachaea. Our data suggest that archaea are the significant regulators of leachate ecology, and more in-depth studies with multiple leachate samples are required to understand their role in leachate nutrient cycling and the development of effective leachate treatment technology.

Comparative genome analysis of four Leuconostoc strains with a focus on carbohydrate-active enzymes and oligosaccharide utilization pathways

•

Comparative genomic analysis of four Leuconostoc strains was performed.

•

Leuconostoc spp. shared genomic similarity, but their genetic content differed.

•

Leuconostoc spp. showed different genes encoding CAZymes.

•

Oligosaccharide’s utilization and folate biosynthesis pathways were investigated.

Leuconostoc is mostly found in food, plants, and dairy products. Due to their innate genomic features, such as the presence of carbohydrate-active enzymes, bacteriocins, and plasmids, Leuconostoc spp. have great biotechnological potential. In this study, four strains were isolated and identified as Leuconostoc mesenteroides SG315 (LA), L. citreum SG255 (LB), L. lactis CCK940 (LC), and L. lactis SBC001 (LD). Comparative analysis was performed using their draft genome sequences. Differences among the four strains were analyzed using the average nucleotide identity, dot plot, and multiple alignments of conserved genomic sequences. Functional profiling revealed 2134, 1917, 1751, and 1816 open reading frames; 2023, 1823, 1655, and 1699 protein-coding genes; 60, 57, 83, and 82 RNA-coding genes; and GC content of 37.5 %, 38.8 %, 43.3 %, and 43.2 %, in LA, LB, LC, and LD, respectively. The total number of genes encoding carbohydrate-active enzymes was 76 (LA), 73 (LB), 57 (LC), and 67 (LD). These results indicate that the four strains shared a large number of genes, but their gene content is different. Furthermore, most genes with unknown functions were observed in the prophage regions of the genome. This study also elucidated the oligosaccharide utilization and folate biosynthesis pathways in Leuconostoc spp. Taken together, our findings provide useful information on the genomic diversity of CAZymes in the four Leuconostoc strains and suggest that these species could be used for potent exploitation.

Integrated Fecal Microbiome and Metabolomics Reveals a Novel Potential Biomarker for Predicting Tibial Dyschondroplasia in Chickens

Tibial dyschondroplasia (TD) is a metabolic tibial-tarsal disorder occurring in fast-growing poultry, and its diagnosis is mainly based on an invasive method. Here, we profiled the fecal gut microbiome and metabolome of broilers with and without TD to identify potential non-invasive and non-stress biomarkers of TD. First, TD broilers with the most pronounced clinical signs during the experiment were screened and faecal samples were collected for integrated microbiome and metabolomics analysis. Moreover, the diagnostic potential of identified biomarkers was further validated throughout the experiment. It was noted that the microbial and metabolic signatures of TD broilers differed from those of normal broilers. TD broilers were characterized by enriched bacterial OTUs of the genus Klebsiella, and depleted genera [Ruminococcus], Dorea, Ruminococcus, Oscillospira, Ochrobactrum, and Sediminibacterium. In addition, a total of 189 fecal differential metabolites were identified, mainly enriched in the purine, vitamin and amino acid metabolism, which were closely associated with differential microbiota and tibia-related indicators. Furthermore, three fecal metabolites were screened, including 4-hydroxybenzaldehyde, which distinguished TD from normal broilers with extremely high specificity and was superior to serum bone markers. These results indicated that gut microbiota equilibrium might influence the pathogenesis of TD by modulating host metabolism, and the identified fecal metabolite 4-hydroxybenzaldehyde might be a potential and non-invasive biomarker for predicting TD in chickens.

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic strategies.

Genome Streamlining, Plasticity, and Metabolic Versatility Distinguish Co-occurring Toxic and Nontoxic Cyanobacterial Strains of Microcoleus

Harmful cyanobacterial bloom occurrences have increased worldwide due to climate change and eutrophication, causing nuisance and animal deaths. Species from the benthic cyanobacterial genus Microcoleus are ubiquitous and form thick mats in freshwater systems, such as rivers, that are sometimes toxic due to the production of potent neurotoxins (anatoxins). Anatoxin-producing (toxic) strains typically coexist with non-anatoxin-producing (nontoxic) strains in mats, although the reason for this is unclear. To determine the genetic mechanisms differentiating toxic and nontoxic Microcoleus, we sequenced and assembled genomes from 11 cultures and compared these to another 31 Microcoleus genomes. Average nucleotide identities (ANI) indicate that toxic and nontoxic strains are distinct species (ANI, <95%), and only 6% of genes are shared across all 42 genomes, suggesting a high level of genetic divergence among Microcoleus strains. Comparative genomics showed substantial genome streamlining in toxic strains and a potential dependency on external sources for thiamine and sucrose. Toxic and nontoxic strains are further differentiated by an additional set of putative nitrate transporter (nitrogen uptake) and cyanophycin (carbon and nitrogen storage) genes, respectively. These genes likely confer distinct competitive advantages based on nutrient availability and suggest nontoxic strains are more robust to nutrient fluctuations. Nontoxic strains also possess twice as many transposable elements, potentially facilitating greater genetic adaptation to environmental changes. Our results offer insights into the divergent evolution of Microcoleus strains and the potential for cooperative and competitive interactions that contribute to the co-occurrence of toxic and nontoxic species within mats. IMPORTANCE Microcoleus autumnalis, and closely related Microcoleus species, compose a geographically widespread group of freshwater benthic cyanobacteria. Canine deaths due to anatoxin-a poisoning, following exposure to toxic proliferations, have been reported globally. While Microcoleus proliferations are on the rise, the mechanisms underpinning competition between, or coexistence of, toxic and nontoxic strains are unknown. This study identifies substantial genetic differences between anatoxin-producing and non-anatoxin-producing strains, pointing to reduced metabolic flexibility in toxic strains, and potential dependence on cohabiting nontoxic strains. Results provide insights into the metabolic and evolutionary differences between toxic and nontoxic Microcoleus, which may assist in predicting and managing aquatic proliferations.

Amplicon Sequence Variants Artificially Split Bacterial Genomes into Separate Clusters

Amplicon sequencing variants (ASVs) have been proposed as an alternative to operational taxonomic units (OTUs) for analyzing microbial communities. ASVs have grown in popularity, in part because of a desire to reflect a more refined level of taxonomy since they do not cluster sequences based on a distance-based threshold. However, ASVs and the use of overly narrow thresholds to identify OTUs increase the risk of splitting a single genome into separate clusters. To assess this risk, I analyzed the intragenomic variation of 16S rRNA genes from the bacterial genomes represented in an rrn copy number database, which contained 20,427 genomes from 5,972 species. As the number of copies of the 16S rRNA gene increased in a genome, the number of ASVs also increased. There was an average of 0.58 ASVs per copy of the 16S rRNA gene for full-length 16S rRNA genes. It was necessary to use a distance threshold of 5.25% to cluster full-length ASVs from the same genome into a single OTU with 95% confidence for genomes with 7 copies of the 16S rRNA, such as Escherichia coli. This research highlights the risk of splitting a single bacterial genome into separate clusters when ASVs are used to analyze 16S rRNA gene sequence data. Although there is also a risk of clustering ASVs from different species into the same OTU when using broad distance thresholds, these risks are of less concern than artificially splitting a genome into separate ASVs and OTUs.

IMPORTANCE 16S rRNA gene sequencing has engendered significant interest in studying microbial communities. There has been tension between trying to classify 16S rRNA gene sequences to increasingly lower taxonomic levels and the reality that those levels were defined using more sequence and physiological information than is available from a fragment of the 16S rRNA gene. Furthermore, the naming of bacterial taxa reflects the biases of those who name them. One motivation for the recent push to adopt ASVs in place of OTUs in microbial community analyses is to allow researchers to perform their analyses at the finest possible level that reflects species-level taxonomy. The current research is significant because it quantifies the risk of artificially splitting bacterial genomes into separate clusters. Far from providing a better representation of bacterial taxonomy and biology, the ASV approach can lead to conflicting inferences about the ecology of different ASVs from the same genome.

Genome Annotation of Poly(lactic acid) Degrading Pseudomonas aeruginosa, Sphingobacterium sp. and Geobacillus sp

Pseudomonas aeruginosa and Sphingobacterium sp. are well known for their ability to decontaminate many environmental pollutants while Geobacillus sp. have been exploited for their thermostable enzymes. This study reports the annotation of genomes of P. aeruginosa S3, Sphingobacterium S2 and Geobacillus EC-3 that were isolated from compost, based on their ability to degrade poly(lactic acid), PLA. Draft genomes of the strains were assembled from Illumina reads, annotated and viewed with the aim of gaining insight into the genetic elements involved in degradation of PLA. The draft genome of Sphinogobacterium strain S2 (435 contigs) was estimated at 5,604,691 bp and the draft genome of P. aeruginosa strain S3 (303 contigs) was estimated at 6,631,638 bp. The draft genome of the thermophile Geobacillus strain EC-3 (111 contigs) was estimated at 3,397,712 bp. A total of 5385 (60% with annotation), 6437 (80% with annotation) and 3790 (74% with annotation) protein-coding genes were predicted for strains S2, S3 and EC-3, respectively. Catabolic genes for the biodegradation of xenobiotics, aromatic compounds and lactic acid as well as the genes attributable to the establishment and regulation of biofilm were identified in all three draft genomes. Our results reveal essential genetic elements that facilitate PLA metabolism at mesophilic and thermophilic temperatures in these three isolates.

Opening a next-generation black box: Ecological trends for hundreds of species-like taxa uncovered within a single bacterial >99% 16S rRNA operational taxonomic unit

Current knowledge on environmental distribution and taxon richness of free-living bacteria is mainly based on cultivation-independent investigations employing 16S rRNA gene sequencing methods. Yet, 16S rRNA genes are evolutionarily rather conserved, resulting in limited taxonomic and ecological resolutions provided by this marker. The faster evolving protein-encoding gene priB was used to reveal ecological patterns hidden within a single operational taxonomic unit (OTU) defined by >99% 16S rRNA sequence similarity. The studied subcluster PnecC of the genus Polynucleobacter represents a ubiquitous group of abundant freshwater bacteria with cosmopolitan distribution, which is very frequently detected by diversity surveys of freshwater systems. Based on genome taxonomy and a large set of genome sequences, a sequence similarity threshold for delineation of species-like taxa could be established. In total, 600 species-like taxa were detected in 99 freshwater habitats scattered across three regions representing a latitudinal range of 3,400 km (42°N to 71°N) and a pH gradient of 4.2 to 8.6. In addition to the unexpectedly high richness, the increased taxonomic resolution revealed structuring of Polynucleobacter communities by a couple of macroecological trends, which was previously only demonstrated for phylogenetically much broader groups of bacteria. An unexpected pattern was the almost complete compositional separation of Polynucleobacter communities of Ca²⁺ -rich and Ca²⁺ -poor habitats. This compositional pattern strongly resembled the vicariance of plant species on silicate and limestone soils. The new cultivation-independent approach presented opened a window to an incredible, previously unseen diversity, and enables investigations aiming on deeper understanding of how environmental conditions shape bacterial communities and drive evolution of free-living bacteria.

All together now: Limitations and recommendations for the simultaneous analysis of all eukaryotic soil sequences

The soil environment contains a large, but historically underexplored, reservoir of biodiversity. Sequencing prokaryotic marker genes has become commonplace for the discovery and characterization of soil bacteria and archaea. Increasingly, this approach is also applied to eukaryotic marker genes to characterize the diversity and distribution of soil eukaryotes. However, understanding the properties and limitations of eukaryotic marker sequences is essential for correctly analysing, interpreting, and synthesizing the resulting data. Here, we illustrate several biases from sequencing data that affect measurements of biodiversity that arise from variation in morphology, taxonomy and phylogeny between organisms, as well as from sampling designs. We recommend analytical approaches to overcome these limitations, and outline how the benchmarking and standardization of sequencing protocols may improve the comparability of the data.

Pseudomonas allokribbensis sp. nov. and Pseudomonas gozinkensis sp. nov., Two New Species Isolated from a Volcanic Island, Izu Oshima, Japan

The genomes of two Pseudomonas strains, IzPS23^T and IzPS32d^T isolated from soil samples of Izu Oshima were compared to Pseudomonas type strains. Whole-genome sequence analysis revealed both belong to the Pseudomonas fluorescens lineage. The average nucleotide identity values of the whole-genome sequences of IzPS23^T and IzPS32d^T compared with other type strains showed high correlations with Pseudomonas kribbensis (93.1%) and Pseudomonas glycinae (93.5%), respectively. Genome-to-genome distances between the whole-genome sequences of IzPS23^T and IzPS32d^T showed correlations with Pseudomonas kribbensis (51.0%) and Pseudomonas glycinae (53.2%), respectively. Genotypic and phenotypic analysis indicated the two strains were novel species, and were named Pseudomonas allokribbensis (IzPS23^T = CECT 9961^T, = LMG 31525^T) and Pseudomonas gozinkensis (IzPS32d^T = CECT 9962^T, = LMG 31526^T), respectively.

Bioaugmentation with immobilized endophytic Penicillium restrictum to improve quorum quenching activity for biofouling control in an aerobic hollow-fiber membrane bioreactor treating antibiotic-containing wastewater

The effects of bioaugmentation with immobilized Penicillium restrictum on the removal efficiency of sulfamethoxazole (SMX), erythromycin (ERY) and tetracycline (TC) antibiotics as well as membrane biofouling was studied using hollow-fiber membrane bioreactor (HF-MBR). Bioaugmentation with P. restrictum led to a significant change in the antibiotic removal efficiency and relative abundance of aerobic microbial community, most probably as a result of its quorum quenching activity. Furthermore, in addition to its role in the increase of SMX and ERY removal efficiencies and the decrease of their sorption on solid phase, bioaugmentation significantly reduced the transmembrane pressure which in turn reduced membrane clogging. The most abundant phyla in sludge and biofilm samples in the presence of P. restrictum were observed to be Proteobacteria, Bacteroidetes and Firmicutes. Differences in bacterial compositions and their specificity in biodegradation of antibiotics in different reactors showed that bacteria were specifically selected under the pressure of antibiotics and growing fungus.

Microdiversity and phylogeographic diversification of bacterioplankton in pelagic freshwater systems revealed through long-read amplicon sequencing

BACKGROUND

Freshwater ecosystems are inhabited by members of cosmopolitan bacterioplankton lineages despite the disconnected nature of these habitats. The lineages are delineated based on > 97% 16S rRNA gene sequence similarity, but their intra-lineage microdiversity and phylogeography, which are key to understanding the eco-evolutional processes behind their ubiquity, remain unresolved. Here, we applied long-read amplicon sequencing targeting nearly full-length 16S rRNA genes and the adjacent ribosomal internal transcribed spacer sequences to reveal the intra-lineage diversities of pelagic bacterioplankton assemblages in 11 deep freshwater lakes in Japan and Europe.

RESULTS

Our single nucleotide-resolved analysis, which was validated using shotgun metagenomic sequencing, uncovered 7-101 amplicon sequence variants for each of the 11 predominant bacterial lineages and demonstrated sympatric, allopatric, and temporal microdiversities that could not be resolved through conventional approaches. Clusters of samples with similar intra-lineage population compositions were identified, which consistently supported genetic isolation between Japan and Europe. At a regional scale (up to hundreds of kilometers), dispersal between lakes was unlikely to be a limiting factor, and environmental factors or genetic drift were potential determinants of population composition. The extent of microdiversification varied among lineages, suggesting that highly diversified lineages (e.g., Iluma-A2 and acI-A1) achieve their ubiquity by containing a consortium of genotypes specific to each habitat, while less diversified lineages (e.g., CL500-11) may be ubiquitous due to a small number of widespread genotypes. The lowest extent of intra-lineage diversification was observed among the dominant hypolimnion-specific lineage (CL500-11), suggesting that their dispersal among lakes is not limited despite the hypolimnion being a more isolated habitat than the epilimnion.

CONCLUSIONS

Our novel approach complemented the limited resolution of short-read amplicon sequencing and limited sensitivity of the metagenome assembly-based approach, and highlighted the complex ecological processes underlying the ubiquity of freshwater bacterioplankton lineages. To fully exploit the performance of the method, its relatively low read throughput is the major bottleneck to be overcome in the future. Video abstract.

Mycobacterium helveticum sp. nov., a novel slowly growing mycobacterial species associated with granulomatous lesions in adult swine

The occurrence of nontuberculous mycobacteria in different hosts and their implication as obligate or opportunistic pathogens remain mainly unclear. Mycobacteriosis in pigs is usually associated with members of the Mycobacterium avium complex and, in particular, with ‘Mycobacterium avium subsp. hominissuis’. Here we describe a novel slow-growing mycobacterial species isolated from lymph nodes obtained from two sows housed in different Swiss farms. The animals presented chronic inappetence and mild diarrhoea. Gross pathology revealed focal caseous lymphadenopathy of the mesenteric lymph nodes. Complete genome sequencing of the two isolates from the two sows was performed. The genomes comprised 5.76 Mb and an average nucleotide identity score of 99.97 %. Whole genome sequence, mycolic acid and matrix-assisted laser desorption ionization-time of flight mass spectrometry analyses revealed that the two isolates were not related to any previously described Mycobacterium species. The closest related species was Mycobacterium parmense, a slow-growing scotochromogenic mycobacterium first isolated from a cervical lymph node of a 3-year-old child. The name proposed for the new species is Mycobacterium helveticum sp. nov. and 16-83^T (=DSM 109965^T= LMG 2019-02457^T) is the type strain.

Nitrosophilus alvini gen. nov., sp. nov., a hydrogen-oxidizing chemolithoautotroph isolated from a deep-sea hydrothermal vent in the East Pacific Rise, inferred by a genome-based taxonomy of the phylum "Campylobacterota"

A novel bacterium, strain EPR55-1^T, was isolated from a deep-sea hydrothermal vent on the East Pacific Rise. The cells were motile rods. Growth was observed at temperatures between 50 and 60°C (optimum, 60°C), at pH values between 5.4 and 8.6 (optimum, pH 6.6) and in the presence of 2.4–3.2% (w/v) NaCl (optimum, 2.4%). The isolate used molecular hydrogen as its sole electron donor, carbon dioxide as its sole carbon source, ammonium as its sole nitrogen source, and thiosulfate, sulfite (0.01 to 0.001%, w/v) or elemental sulfur as its sole sulfur source. Nitrate, nitrous oxide (33%, v/v), thiosulfate, molecular oxygen (0.1%, v/v) or elemental sulfur could serve as the sole electron acceptor to support growth. Phylogenetic analyses based on both 16S rRNA gene sequences and whole genome sequences indicated that strain EPR55-1^T belonged to the family Nitratiruptoraceae of the class “Campylobacteria”, but it had the distinct phylogenetic relationship with the genus Nitratiruptor. On the basis of the physiological and molecular characteristics of the isolate, the name Nitrosophilus alvini gen. nov. sp. nov. is proposed, with EPR55-1^T as the type strain (= JCM 32893^T = KCTC 15925^T). In addition, it is shown that “Nitratiruptor labii” should be transferred to the genus Nitrtosophilus; the name Nitrosophilus labii comb. nov. (JCM 34002^T = DSM 111345^T) is proposed for this organism. Furthermore, 16S rRNA gene-based and genome-based analyses showed that Cetia pacifica is phylogenetically associated with Caminibacter species. We therefore propose the reclassification of Cetia pacifica as Caminibacter pacificus comb. nov. (DSM 27783^T = JCM 19563^T). Additionally, AAI thresholds for genus classification and the reclassification of subordinate taxa within “Campylobacteria” are also evaluated, based on the analyses using publicly available genomes of all the campylobacterial species.

The Limits and Avoidance of Biases in Metagenomic Analyses of Human Fecal Microbiota

An increasing body of evidence highlights the role of fecal microbiota in various human diseases. However, more than two-thirds of fecal bacteria cannot be cultivated by routine laboratory techniques. Thus, physicians and scientists use DNA sequencing and statistical tools to identify associations between bacterial subgroup abundances and disease. However, discrepancies between studies weaken these results. In the present study, we focus on biases that might account for these discrepancies. First, three different DNA extraction methods (G’NOME, QIAGEN, and PROMEGA) were compared with regard to their efficiency, i.e., the quality and quantity of DNA recovered from feces of 10 healthy volunteers. Then, the impact of the DNA extraction method on the bacteria identification and quantification was evaluated using our published cohort of sample subjected to both 16S rRNA sequencing and whole metagenome sequencing (WMS). WMS taxonomical assignation employed the universal marker genes profiler mOTU-v2, which is considered the gold standard. The three standard pipelines for 16S RNA analysis (MALT and MEGAN6, QIIME1, and DADA2) were applied for comparison. Taken together, our results indicate that the G’NOME-based method was optimal in terms of quantity and quality of DNA extracts. 16S rRNA sequence-based identification of abundant bacteria genera showed acceptable congruence with WMS sequencing, with the DADA2 pipeline yielding the highest congruent levels. However, for low abundance genera (<0.5% of the total abundance) two pipelines and/or validation by quantitative polymerase chain reaction (qPCR) or WMS are required. Hence, 16S rRNA sequencing for bacteria identification and quantification in clinical and translational studies should be limited to diagnostic purposes in well-characterized and abundant genera. Additional techniques are warranted for low abundant genera, such as WMS, qPCR, or the use of two bio-informatics pipelines.

Genome-based classification of Acidihalobacter prosperus F5 (=DSM 105917=JCM 32255) as Acidihalobacter yilgarnensis sp. nov

The genus Acidihalobacter has three validated species, Acidihalobacter ferrooxydans, Acidihalobacter prosperus and Acidihalobacter aeolinanus, all of which were isolated from Vulcano island, Italy. They are obligately chemolithotrophic, aerobic, acidophilic and halophilic in nature and use either ferrous iron or reduced sulphur as electron donors. Recently, a novel strain was isolated from an acidic, saline drain in the Yilgarn region of Western Australia. Strain F5T has an absolute requirement for sodium chloride (>5 mM) and is osmophilic, growing in elevated concentrations (>1 M) of magnesium sulphate. A defining feature of its physiology is its ability to catalyse the oxidative dissolution of the most abundant copper mineral, chalcopyrite, suggesting a potential role in biomining. Originally categorized as a strain of A. prosperus, 16S rRNA gene phylogeny and multiprotein phylogenies derived from clusters of orthologous proteins (COGS) of ribosomal protein families and universal protein families unambiguously demonstrate that strain F5T forms a well-supported separate branch as a sister clade to A. prosperus and is clearly distinguishable from A. ferrooxydans DSM 14175T and A. aeolinanus DSM14174T. Results of comparisons between strain F5T and the other Acidihalobacter species, using genome-based average nucleotide identity, average amino acid identity, correlation indices of tetra-nucleotide signatures (Tetra) and genome-to-genome distance (digital DNA-DNA hybridization), support the contention that strain F5T represents a novel species of the genus Acidihalobacter. It is proposed that strain F5T should be formally reclassified as Acidihalobacter yilgarnenesis F5T (=DSM 105917T=JCM 32255T).

Iterative subtractive binning of freshwater chronoseries metagenomes identifies over 400 novel species and their ecologic preferences

Recent advances in sequencing technology and bioinformatic pipelines have allowed unprecedented access to the genomes of yet-uncultivated microorganisms from diverse environments. However, the catalogue of freshwater genomes remains limited, and most genome recovery attempts in freshwater ecosystems have only targeted specific taxa. Here, we present a genome recovery pipeline incorporating iterative subtractive binning, and apply it to a time series of 100 metagenomic datasets from seven connected lakes and estuaries along the Chattahoochee River (Southeastern USA). Our set of metagenome-assembled genomes (MAGs) represents >400 yet-unnamed genomospecies, substantially increasing the number of high-quality MAGs from freshwater lakes. We propose names for two novel species: 'Candidatus Elulimicrobium humile' ('Ca. Elulimicrobiota', 'Patescibacteria') and 'Candidatus Aquidulcis frankliniae' ('Chloroflexi'). Collectively, our MAGs represented about half of the total microbial community at any sampling point. To evaluate the prevalence of these genomospecies in the chronoseries, we introduce methodologies to estimate relative abundance and habitat preference that control for uneven genome quality and sample representation. We demonstrate high degrees of habitat-specialization and endemicity for most genomospecies in the Chattahoochee lakes. Wider ecological ranges characterized smaller genomes with higher coding densities, indicating an overall advantage of smaller, more compact genomes for cosmopolitan distributions.

Vast Differences in Strain-Level Diversity in the Gut Microbiota of Two Closely Related Honey Bee Species

Most bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.

Vast Differences in Strain-Level Diversity in the Gut Microbiota of Two Closely Related Honey Bee Species

Most bacterial species encompass strains with vastly different gene content. Strain diversity in microbial communities is therefore considered to be of functional importance. Yet little is known about the extent to which related microbial communities differ in diversity at this level and which underlying mechanisms may constrain and maintain strain-level diversity. Here, we used shotgun metagenomics to characterize and compare the gut microbiota of two honey bee species, Apis mellifera and Apis cerana, which diverged about 6 mya. Although the host species are colonized largely by the same bacterial 16S rRNA phylotypes, we find that their communities are host specific when analyzed with genomic resolution. Moreover, despite their similar ecology, A. mellifera displayed a much higher diversity of strains and functional gene content in the microbiota compared to A. cerana, both per colony and per individual bee. In particular, the gene repertoire for polysaccharide degradation was massively expanded in the microbiota of A. mellifera relative to A. cerana. Bee management practices, divergent ecological adaptation, or habitat size may have contributed to the observed differences in microbiota genomic diversity of these key pollinator species. Our results illustrate that the gut microbiota of closely related animal hosts can differ vastly in genomic diversity while displaying similar levels of diversity based on the 16S rRNA gene. Such differences are likely to have consequences for gut microbiota functioning and host-symbiont interactions, highlighting the need for metagenomic studies to understand the ecology and evolution of microbial communities.

•

Metagenomics reveals differences in gut microbiota diversity beyond the 16S rRNA gene

•

Apis cerana and Apis mellifera harbor distinct species and strains in their gut

•

Diversity is much higher in A. mellifera per individual bee and within colonies

•

Major differences in functions are related to polysaccharide degradation

Metabolic diversity and co-occurrence of multiple Ferrovum species at an acid mine drainage site

BACKGROUND

Ferrovum spp. are abundant in acid mine drainage sites globally where they play an important role in biogeochemical cycling. All known taxa in this genus are Fe(II) oxidizers. Thus, co-occurring members of the genus could be competitors within the same environment. However, we found multiple, co-occurring Ferrovum spp. in Cabin Branch, an acid mine drainage site in the Daniel Boone National Forest, KY.

RESULTS

Here we describe the distribution of Ferrovum spp. within the Cabin Branch communities and metagenome assembled genomes (MAGs) of two new Ferrovum spp. In contrast to previous studies, we recovered multiple 16S rRNA gene sequence variants suggesting the commonly used 97% cutoff may not be appropriate to differentiate Ferrovum spp. We also retrieved two nearly-complete Ferrovum spp. genomes from metagenomic data. The genomes of these taxa differ in several key ways relating to nutrient cycling, motility, and chemotaxis.

CONCLUSIONS

Previously reported Ferrovum genomes are also diverse with respect to these categories suggesting that the genus Ferrovum contains substantial metabolic diversity. This diversity likely explains how the members of this genus successfully co-occur in Cabin Branch and why Ferrovum spp. are abundant across geochemical gradients.

Temperature and Nutrient Levels Correspond with Lineage-Specific Microdiversification in the Ubiquitous and Abundant Freshwater Genus Limnohabitans

Most freshwater bacterial communities are characterized by a few dominant taxa that are often ubiquitous across freshwater biomes worldwide. Our understanding of the genomic diversity within these taxonomic groups is limited to a subset of taxa. Here, we investigated the genomic diversity that enables Limnohabitans, a freshwater genus key in funneling carbon from primary producers to higher trophic levels, to achieve abundance and ubiquity. We reconstructed eight putative Limnohabitans metagenome-assembled genomes (MAGs) from stations located along broad environmental gradients existing in Lake Michigan, part of Earth's largest surface freshwater system. De novo strain inference analysis resolved a total of 23 strains from these MAGs, which strongly partitioned into two habitat-specific clusters with cooccurring strains from different lineages. The largest number of strains belonged to the abundant LimB lineage, for which robust in situ strain delineation had not previously been achieved. Our data show that temperature and nutrient levels may be important environmental parameters associated with microdiversification within the Limnohabitans genus. In addition, strains predominant in low- and high-phosphorus conditions had larger genomic divergence than strains abundant under different temperatures. Comparative genomics and gene expression analysis yielded evidence for the ability of LimB populations to exhibit cellular motility and chemotaxis, a phenotype not yet associated with available Limnohabitans isolates. Our findings broaden historical marker gene-based surveys of Limnohabitans microdiversification and provide in situ evidence of genome diversity and its functional implications across freshwater gradients.IMPORTANCE Limnohabitans is an important bacterial taxonomic group for cycling carbon in freshwater ecosystems worldwide. Here, we examined the genomic diversity of different Limnohabitans lineages. We focused on the LimB lineage of this genus, which is globally distributed and often abundant, and its abundance has shown to be largely invariant to environmental change. Our data show that the LimB lineage is actually comprised of multiple cooccurring populations for which the composition and genomic characteristics are associated with variations in temperature and nutrient levels. The gene expression profiles of this lineage suggest the importance of chemotaxis and motility, traits that had not yet been associated with the Limnohabitans genus, in adapting to environmental conditions.

The Ecobiomics project: Advancing metagenomics assessment of soil health and freshwater quality in Canada

Transformative advances in metagenomics are providing an unprecedented ability to characterize the enormous diversity of microorganisms and invertebrates sustaining soil health and water quality. These advances are enabling a better recognition of the ecological linkages between soil and water, and the biodiversity exchanges between these two reservoirs. They are also providing new perspectives for understanding microorganisms and invertebrates as part of interacting communities (i.e. microbiomes and zoobiomes), and considering plants, animals, and humans as holobionts comprised of their own cells as well as diverse microorganisms and invertebrates often acquired from soil and water. The Government of Canada's Genomics Research and Development Initiative (GRDI) launched the Ecobiomics Project to coordinate metagenomics capacity building across federal departments, and to apply metagenomics to better characterize microbial and invertebrate biodiversity for advancing environmental assessment, monitoring, and remediation activities. The Project has adopted standard methods for soil, water, and invertebrate sampling, collection and provenance of metadata, and nucleic acid extraction. High-throughput sequencing is located at a centralized sequencing facility. A centralized Bioinformatics Platform was established to enable a novel government-wide approach to harmonize metagenomics data collection, storage and bioinformatics analyses. Sixteen research projects were initiated under Soil Microbiome, Aquatic Microbiome, and Invertebrate Zoobiome Themes. Genomic observatories were established at long-term environmental monitoring sites for providing more comprehensive biodiversity reference points to assess environmental change.

phylogenize: correcting for phylogeny reveals genes associated with microbial distributions

Summary

Phylogenetic comparative methods are powerful but presently under-utilized ways to identify microbial genes underlying differences in community composition. These methods help to identify functionally important genes because they test for associations beyond those expected when related microbes occupy similar environments. We present phylogenize, a pipeline with web, QIIME 2 and R interfaces that allows researchers to perform phylogenetic regression on 16S amplicon and shotgun sequencing data and to visualize results. phylogenize applies broadly to both host-associated and environmental microbiomes. Using Human Microbiome Project and Earth Microbiome Project data, we show that phylogenize draws similar conclusions from 16S versus shotgun sequencing and reveals both known and candidate pathways associated with host colonization.

Availability and implementation

phylogenize is available at https://phylogenize.org and https://bitbucket.org/pbradz/phylogenize.

Supplementary information

Supplementary data are available at Bioinformatics online.

Next-generation sequencing as a screening tool for foodborne pathogens in fresh produce

Next generation sequencing (NGS) approaches are increasingly applied to tracing microbial contaminants entering the food chain due to NGS' untargeted nature and ability to investigate non-culturable (and/or difficult to culture) organisms while yielding genomic information about the microbiota. So far, a plethora of microbes has been shown to be associated with fresh produce, but few studies have utilised NGS to identify contamination with human pathogens. This study aims to establish the limit of detection (LoD) for Salmonella and phage MS2 (a Norovirus surrogate) contamination of fresh produce employing NGS approaches on the Illumina MiSeq: 16S amplicon-sequencing, and RNA-seq, using ScriptSeq (Illumina) and NEBNext (New England BioLabs) kits. ScriptSeq proved the most sensitive approach; delivering an LoD of 10⁴ CFU reaction^-1 (Colony Forming Units) for Salmonella and 10⁵ PFU reaction^-1 (Plaque Forming Units) for phage MS2. Use of the NEBNext kit resulted in detection of Salmonella at 10⁶ CFU reaction^-1 and phage MS2 at 10⁷ PFU reaction^-1. 16S amplicon-sequencing yielded a similar LoD of 10⁵ CFU reaction^-1 for Salmonella but could not detect MS2. The tested NGS methodologies, in combination with bioinformatics approaches applied, proved less sensitive than conventional microbial detection approaches.

The Microbial Genomes Atlas (MiGA) webserver: taxonomic and gene diversity analysis of Archaea and Bacteria at the whole genome level

The small subunit ribosomal RNA gene (16S rRNA) has been successfully used to catalogue and study the diversity of prokaryotic species and communities but it offers limited resolution at the species and finer levels, and cannot represent the whole-genome diversity and fluidity. To overcome these limitations, we introduced the Microbial Genomes Atlas (MiGA), a webserver that allows the classification of an unknown query genomic sequence, complete or partial, against all taxonomically classified taxa with available genome sequences, as well as comparisons to other related genomes including uncultivated ones, based on the genome-aggregate Average Nucleotide and Amino Acid Identity (ANI/AAI) concepts. MiGA integrates best practices in sequence quality trimming and assembly and allows input to be raw reads or assemblies from isolate genomes, single-cell sequences, and metagenome-assembled genomes (MAGs). Further, MiGA can take as input hundreds of closely related genomes of the same or closely related species (a so-called ‘Clade Project’) to assess their gene content diversity and evolutionary relationships, and calculate important clade properties such as the pangenome and core gene sets. Therefore, MiGA is expected to facilitate a range of genome-based taxonomic and diversity studies, and quality assessment across environmental and clinical settings. MiGA is available at http://microbial-genomes.org/.

Resource heterogeneity structures aquatic bacterial communities

Microorganisms are strongly influenced by the bottom-up effects of resource supply. While many species respond to fluctuations in the concentration of resources, microbial diversity may also be affected by the heterogeneity of the resource pool, which often reflects a mixture of distinct molecules. To test this hypothesis, we examined resource-diversity relationships for bacterioplankton in a set of north temperate lakes that varied in their concentration and composition of dissolved organic matter (DOM), which is an important resource for heterotrophic bacteria. Using 16S rRNA transcript sequencing and ecosystem metabolomics, we documented strong relationships between bacterial alpha-diversity (richness and evenness) and the bulk concentration and the number of molecules in the DOM pool. Similarly, bacterial community beta-diversity was related to both DOM concentration and composition. However, in some lakes the relative abundance of resource generalists, which was inversely related to the DOM concentration, may have reduced the effect of DOM heterogeneity on community composition. Together, our results demonstrate the potential metabolic interactions between bacteria and organic matter and suggest that changes in organic matter composition may alter the structure and function of bacterial communities.

The Phylogenomic Diversity of Herbivore-Associated Fibrobacter spp. Is Correlated to Lignocellulose-Degrading Potential

Members of the genus Fibrobacter are cellulose-degrading bacteria and common constituents of the gastrointestinal microbiota of herbivores. Although considerable phylogenetic diversity is observed among members of this group, few functional differences explaining the distinct ecological distributions of specific phylotypes have been described. In this study, we sequenced and performed a comparative analysis of whole genomes from 38 novel Fibrobacter strains against the type strains for the two formally described Fibrobacter species F. succinogenes strain S85 and F. intestinalis strain NR9. Significant differences in the number of genes encoding carbohydrate-active enzyme families involved in plant cell wall polysaccharide degradation were observed among Fibrobacter phylotypes. F. succinogenes genomes were consistently enriched in genes encoding carbohydrate-active enzymes compared to those of F. intestinalis strains. Moreover, genomes of F. succinogenes phylotypes that are dominant in the rumen had significantly more genes annotated to major families involved in hemicellulose degradation (e.g., CE6, GH10, and GH43) than did the genomes of F. succinogenes phylotypes typically observed in the lower gut of large hindgut-fermenting herbivores such as horses. Genes encoding a putative urease were also identified in 12 of the Fibrobacter genomes, which were primarily isolated from hindgut-fermenting hosts. Screening for growth on urea as the sole source of nitrogen provided strong evidence that the urease was active in these strains. These results represent the strongest evidence reported to date for specific functional differences contributing to the ecology of Fibrobacter spp. in the herbivore gut.IMPORTANCE The herbivore gut microbiome is incredibly diverse, and a functional understanding of this diversity is needed to more reliably manipulate this community for specific gain, such as increased production in ruminant livestock. Microbial degraders of plant cell wall polysaccharides in the herbivore gut, particularly Fibrobacter spp., are of fundamental importance to their hosts for digestion of a diet consisting primarily of recalcitrant plant fibers. Considerable phylogenetic diversity exists among members of the genus Fibrobacter, but much of this diversity remains cryptic. Here, we used comparative genomics, applied to a diverse collection of recently isolated Fibrobacter strains, to identify a robust association between carbohydrate-active enzyme gene content and the Fibrobacter phylogeny. Our results provide the strongest evidence reported to date for functional differences among Fibrobacter phylotypes associated with either the rumen or the hindgut and emphasize the general significance of carbohydrate-active enzymes in the evolution of fiber-degrading bacteria.

Quantifying the changes in genetic diversity within sequence-discrete bacterial populations across a spatial and temporal riverine gradient

Recent diversity studies have revealed that microbial communities of natural environments are dominated by species-like, sequence-discrete populations. However, how stable the sequence and gene-content diversity are within these populations and especially in highly dynamic lotic habitats remain unclear. Here we quantified the dynamics of intra-population diversity in samples spanning two years and five sites in the Kalamas River (Northwest Greece). A significant positive correlation was observed between higher intra-population sequence diversity and longer persistence over time, revealing that more diverse populations tended to represent more autochthonous (vs. allochthonous) community members. Assessment of intra-population gene-content changes caused by strain replacement or gene loss over time revealed different profiles with the majority of populations exhibiting gene-content changes close to 10% of the total genes, while one population exhibited ~21% change. The variable genes were enriched in hypothetical proteins and mobile elements, and thus, were probably functionally neutral or attributable to phage predation. A few notable exceptions to this pattern were also noted such as phototrophy-related proteins in summer vs. winter populations. Taken together, these results revealed that some freshwater genomes are remarkably dynamic, even across short time and spatial scales, and have implications for the bacterial species concept and microbial source tracking.

imGLAD: accurate detection and quantification of target organisms in metagenomes

Accurate detection of target microbial species in metagenomic datasets from environmental samples remains limited because the limit of detection of current methods is typically inaccessible and the frequency of false-positives, resulting from inadequate identification of regions of the genome that are either too highly conserved to be diagnostic (e.g., rRNA genes) or prone to frequent horizontal genetic exchange (e.g., mobile elements) remains unknown. To overcome these limitations, we introduce imGLAD, which aims to detect (target) genomic sequences in metagenomic datasets. imGLAD achieves high accuracy because it uses the sequence-discrete population concept for discriminating between metagenomic reads originating from the target organism compared to reads from co-occurring close relatives, masks regions of the genome that are not informative using the MyTaxa engine, and models both the sequencing breadth and depth to determine relative abundance and limit of detection. We validated imGLAD by analyzing metagenomic datasets derived from spinach leaves inoculated with the enteric pathogen Escherichia coli O157:H7 and showed that its limit of detection can be comparable to that of PCR-based approaches for these samples (∼1 cell/gram).

Nonpareil 3: Fast Estimation of Metagenomic Coverage and Sequence Diversity

Estimations of microbial community diversity based on metagenomic data sets are affected, often to an unknown degree, by biases derived from insufficient coverage and reference database-dependent estimations of diversity. For instance, the completeness of reference databases cannot be generally estimated since it depends on the extant diversity sampled to date, which, with the exception of a few habitats such as the human gut, remains severely undersampled. Further, estimation of the degree of coverage of a microbial community by a metagenomic data set is prohibitively time-consuming for large data sets, and coverage values may not be directly comparable between data sets obtained with different sequencing technologies. Here, we extend Nonpareil, a database-independent tool for the estimation of coverage in metagenomic data sets, to a high-performance computing implementation that scales up to hundreds of cores and includes, in addition, a k-mer-based estimation as sensitive as the original alignment-based version but about three hundred times as fast. Further, we propose a metric of sequence diversity (N_d ) derived directly from Nonpareil curves that correlates well with alpha diversity assessed by traditional metrics. We use this metric in different experiments demonstrating the correlation with the Shannon index estimated on 16S rRNA gene profiles and show that N_d additionally reveals seasonal patterns in marine samples that are not captured by the Shannon index and more precise rankings of the magnitude of diversity of microbial communities in different habitats. Therefore, the new version of Nonpareil, called Nonpareil 3, advances the toolbox for metagenomic analyses of microbiomes. IMPORTANCE Estimation of the coverage provided by a metagenomic data set, i.e., what fraction of the microbial community was sampled by DNA sequencing, represents an essential first step of every culture-independent genomic study that aims to robustly assess the sequence diversity present in a sample. However, estimation of coverage remains elusive because of several technical limitations associated with high computational requirements and limiting statistical approaches to quantify diversity. Here we described Nonpareil 3, a new bioinformatics algorithm that circumvents several of these limitations and thus can facilitate culture-independent studies in clinical or environmental settings, independent of the sequencing platform employed. In addition, we present a new metric of sequence diversity based on rarefied coverage and demonstrate its use in communities from diverse ecosystems.