Scalable methods for analyzing and visualizing phylogenetic placement of metagenomic samples

Microbial Community Shifts and Nitrogen Utilization in Peritidal Microbialites: The Role of Salinity and pH in Microbially Induced Carbonate Precipitation

Microbialites have the potential to record environmental changes and act as biosignatures of past geochemical conditions. As such, they could be used as indicators to decipher ancient rock records. Modern microbialites are primarily found in environments where competitors and destructors are absent or where biogeochemical conditions favor their continuous formation. Many previous studies have essentially focused on the role of photosynthetic microbes in controlling pH and carbonate speciation and potentially overlooked alternative non-photosynthetic pathways of carbonate precipitation. Given that microbial activity induces subtle geochemical changes, microbially induced carbonate precipitation (MICP) can involve several mechanisms, from extracellular polymeric substances (EPS), sulfate reduction, anaerobic oxidation of methane, to nitrogen cycling processes, such as ammonification, ureolysis, and denitrification. Moreover, the peritidal zone where temperate microbialites are mostly found today, is under the influence of both freshwater and seawater, arguing for successive biogeochemical processes leading to mineral saturation, and questioning interpretations of fossil records. This study investigates microbialites in three tide pools from the peritidal zone of Fongchueisha, Hengchun, Taiwan, to address the influence of salinity on microbial community composition and carbonate precipitation mechanisms. Microbial samples were collected across varying salinity gradients at multiple time points and analyzed using next-generation sequencing (NGS) of bacterial 16S and eukaryotic 18S rRNA genes. Our results indicate that dominant bacterial groups, including Cyanobacteria and Alphaproteobacteria, were largely influenced by salinity variations, albeit pH exhibited stronger correlation with community composition. Combining our results on geochemistry and taxonomic diversity over time, we inferred a shift in the trophic mode under high salinity conditions, during which the use of urea and amino acids as a nitrogen source outcompetes diazotrophy, ureolysis and ammonification of amino acids reinforcing carbonate precipitation dynamics by triggering an increase in both pH and dissolved inorganic carbon.

Absolute abundance calculation enhances the significance of microbiome data in antibiotic treatment studies

Background

The intestinal microbiota contributes to the colonization resistance of the gut towards bacterial pathogens. Antibiotic treatment often negatively affects the microbiome composition, rendering the host more susceptible for infections. However, a correct interpretation of such a perturbation requires quantitative microbiome profiling to reflect accurately the direction and magnitude of compositional changes within a microbiota. Standard 16S rRNA gene amplicon sequencing of microbiota samples offers compositional data in relative, but not absolute abundancies, and the presence of multiple copies of 16S rRNA genes in bacterial genomes introduces bias into compositional data. We explored whether improved sequencing data analysis influences the significance of the effect exerted by antibiotics on the faecal microbiota of young pigs using two veterinary antibiotics. Calculation of absolute abundances, either by flow cytometry-based bacterial cell counts or by spike-in of synthetic 16S rRNA genes, was employed and 16S rRNA gene copy numbers (GCN) were corrected.

Results

Cell number determination exhibited large interindividual variability in two pig studies, using either tylosin or tulathromycin. Following tylosin application, flow cytometry-based cell counting revealed decreased absolute abundances of five families and ten genera. These results were not detectable by standard 16S analysis based on relative abundances. Here, GCN correction additionally uncovered significant decreases of Lactobacillus and Faecalibacterium. In another experimental setting with tulathromycin treatment, bacterial abundance quantification by flow cytometry and by a spike-in method yielded similar results only on the phylum level. Even though the spike-in method identified the decrease of four genera, analysis by fluorescence-activated cell sorting (FACS) uncovered eight significantly reduced genera, such as Prevotella and Paraprevotella upon antibiotic treatment. In contrast, analysis of relative abundances only showed a decrease of Faecalibacterium and Rikenellaceae RC9 gut group and, thus, a much less detailed antibiotic effect.

Conclusion

Flow cytometry is a laborious method, but identified a higher number of significant microbiome changes in comparison to common compositional data analysis and even revealed to be superior to a spike-in method. Calculation of absolute abundances and GCN correction are valuable methods that should be standards in microbiome analyses in veterinary as well as human medicine.

Diversity of arbuscular mycorrhiza fungi in roots of giant miscanthus (Miscanthus × giganteus) and prairie cordgrass (Spartina pectinata) cultivated on heavy metal-contaminated areas

The use of indigenous AMF species from heavy metal contaminated areas can be a promising tool to support the phytostabilisation of such areas. The aim of the study was to evaluate the AMF species diversity in the roots of the perennial energy grasses Miscanthus × giganteus and Spartina pectinata grown in areas with different levels of heavy metal contamination with regard to the potential use of the dominant AMF species to support phytostabilisation of soils contaminated with Pb, Cd and Zn. Samples were taken from two sites with different levels of Pb, Cd and Zn contamination and from an uncontaminated site as a control. The AMF colonisation of the roots of Miscanthus × giganteus and Spartina pectinata was investigated. The composition of AMF species in the plant roots was determined by sequencing the D2 region of the LSU rDNA of Glomeromycota. Soil contamination had a significant effect on the composition of AMF communities in the roots. Diversispora and Claroideoglomus were the predominant genera in the communities in the heavily heavy metal contaminated area. The AMF communities at moderately contaminated and uncontaminated areas showed a similar structure, with Rhizoglomus as the dominant genus. Species such as Palaeospora spainiae, Rhizoglomus silesianum, Septoglomus sp., Septoglomus nigrum, Ambispora sp., Claroideoglomus etunicatum and Diversispora sp3. were identified exclusively in the roots of Miscanthus × giganteus and Spartina pectinata grown in contaminated areas. They could potentially be used to support phytostabilisation of areas contaminated with Pb, Cd and Zn, but further studies are needed.

Previse preterm birth in early pregnancy through vaginal microbiome signatures using metagenomics and dipstick assays

Annually, in India, 13% of all newborns are preterm, accounting for 23.4% of preterm birth (PTB) globally. The composition and diversity of the vaginal microbiome have a notable degree of ethnic inequality. For understanding differences in vaginal microbiome composition and functions between adverse and normal pregnancy, we have collected, processed and sequenced 600 high vaginal swab (HVS) samples across the three trimesters of pregnancy from 140 women who delivered at term and 60 women who delivered PTB, adopting a targeted metagenomics approach. The microbial signatures in HVS samples showed Lactobacillus genera to be highly abundant in term birth (TB), while in early pregnancy the abundances of Gardnerella, Atopobium, and Sneathia were found to be high in PTB. We further extended our analysis, identified specific microbial genomic signatures, and developed a dipstick assay for rapid identification of PTB-associated microbiota in HVS samples in low-resource settings.

Temporal dynamics of fecal microbiota community succession in broiler chickens, calves, and piglets under aerobic exposure

Researchers have extensively studied the effect of oxygen on the growth and survival of bacteria. However, the impact of oxygen on bacterial community structure, particularly its ability to select for taxa within the context of a complex microbial community, is still unclear. In a 21-day microcosm experiment, we investigated the effect of aerobic exposure on the fecal community structure and succession pattern in broiler, calf, and piglet feces (n = 10 for each feces type). Bacterial diversity decreased and community structure changed rapidly in the broiler microbiome (P < 0.001), while the fecal community of calves and piglets, which have higher initial diversity, was stable after initial exposure but decreased in diversity after 3 days (P < 0.001). The response to aerobic exposure was host animal specific, but in all three animals, the change in community structure was driven by a decrease in anaerobic species, primarily belonging to Firmicutes and Bacteroidetes (except in broilers where Bacteroidetes increased), along with an increase in aerobic species belonging to Proteobacteria and Actinobacteria. Using random forest regression, we identified microbial features that predict aerobic exposure. In all three animals, host-beneficial Prevotella-related ASVs decreased after exposure, while ASVs belonging to Acinetobacter, Corynbacterium, and Tissierella were increased. The decrease of Prevotella was rapid in broilers but delayed in calves and piglets. Knowing when these pathobionts increase in abundance after aerobic exposure could inform farm sanitation practices and could be important in designing animal experiments that modulate the microbiome.IMPORTANCEThe fecal microbial community is contained within a dynamic ecosystem of interacting microbes that varies in biotic and abiotic components across different animal species. Although oxygen affects bacterial growth, its specific impact on the structure of complex communities, such as those found in feces, and how these effects vary between different animal species are poorly understood. In this study, we demonstrate that the effect of aerobic exposure on the fecal microbiota was host-animal-specific, primarily driven by a decrease in Firmicutes and Bacteroidetes, but accompanied by an increase in Actinobacteria, Proteobacteria, and other pathobionts. Interestingly, we observed that more complex communities from pig and cattle exhibited initial resilience, while a less diverse community from broilers displayed a rapid response to aerobic exposure. Our findings offer insights that can inform farm sanitation practices, as well as experimental design, sample collection, and processing protocols for microbiome studies across various animal species.

EPIK: precise and scalable evolutionary placement with informative k-mers

MOTIVATION

Phylogenetic placement enables phylogenetic analysis of massive collections of newly sequenced DNA, when de novo tree inference is too unreliable or inefficient. Assuming that a high-quality reference tree is available, the idea is to seek the correct placement of the new sequences in that tree. Recently, alignment-free approaches to phylogenetic placement have emerged, both to circumvent the need to align the new sequences and to avoid the calculations that typically follow the alignment step. A promising approach is based on the inference of k-mers that can be potentially related to the reference sequences, also called phylo-k-mers. However, its usage is limited by the time and memory-consuming stage of reference data preprocessing and the large numbers of k-mers to consider.

RESULTS

We suggest a filtering method for selecting informative phylo-k-mers based on mutual information, which can significantly improve the efficiency of placement, at the cost of a small loss in placement accuracy. This method is implemented in IPK, a new tool for computing phylo-k-mers that significantly outperforms the software previously available. We also present EPIK, a new software for phylogenetic placement, supporting filtered phylo-k-mer databases. Our experiments on real-world data show that EPIK is the fastest phylogenetic placement tool available, when placing hundreds of thousands and millions of queries while still providing accurate placements.

AVAILABILITY AND IMPLEMENTATION

IPK and EPIK are freely available at https://github.com/phylo42/IPK and https://github.com/phylo42/EPIK. Both are implemented in C++ and Python and supported on Linux and MacOS.

Phylogenetic diversity only weakly mitigates climate-change-driven biodiversity loss in insect communities

To help address the underrepresentation of arthropods and Asian biodiversity from climate-change assessments, we carried out year-long, weekly sampling campaigns with Malaise traps at different elevations and latitudes in Gaoligongshan National Park in southwestern China. From these 623 samples, we barcoded 10,524 beetles and compared scenarios of climate-change-induced biodiversity loss, by designating seasonal, elevational, and latitudinal subsets of beetles as communities that plausibly could go extinct as a group, which we call "loss sets". The availability of a published mitochondrial-genome-based phylogeny of the Coleoptera allowed us to compare the loss of species diversity with and without accounting for phylogenetic relatedness. We hypothesised that phylogenetic relatedness would mitigate extinction, since the extinction of any loss set would result in the disappearance of all its species but only part of its evolutionary history, which is still extant in the remaining loss sets. We found different patterns of community clustering by season and latitude, depending on whether phylogenetic information was incorporated. However, accounting for phylogeny only slightly mitigated the amount of biodiversity loss under climate change scenarios, against our expectations: there is no phylogenetic "escape clause" for biodiversity conservation. We achieve the same results whether phylogenetic information was derived from the mitogenome phylogeny or from a de novo barcode-gene tree. We encourage interested researchers to use this data set to study lineage-specific community assembly patterns in conjunction with life-history traits and environmental covariates.

Saposhnikovia divaricata root and its major components ameliorate inflammation and altered gut microbial diversity and compositions in DSS-induced colitis

Background

The root of Saposhnikovia divaricata (Turcz.) Schischk is a well-known traditional medicinal plant, containing various bioactive compounds with anti-inflammatory, antioxidant, and analgesic properties. However, no scientific studies have validated its clinical use as an anti-inflammatory agent against inflammatory bowel disease (IBD). This study aimed to investigate whether the root extract of S. divaricata ameliorates IBD and induces gut microbial alteration, using a RAW 264.7 cell line and a DSS-induced colitis mouse model.

Methods

To investigate the anti-inflammatory effects and alleviation of IBD, using a methanol extract of Saposhnikovia divaricata (Turcz.) Schischk. root (MESD), RAW 264.7, murine macrophages and a dextran sodium sulfate (DSS)-induced colitis mouse model were employed. 16S rRNA gene sequencing was conducted to determine the alterations in the gut microbiota of mice with DSS-induced colitis.

Results

MESD significantly decreased nitric oxide (NO) and inflammatory cytokine levels in lipopolysaccharide (LPS)-induced RAW 264.7 cells in vitro. Oral administration of MESD reduced the expression of inflammatory cytokines in the colons of mice with DSS-induced colitis. Additionally, MESD inhibited the abundance of Clostridium sensu stricto 1 and enhanced the predicted functional pathways, including l-glutamate degradation VIII (to propanoic acid). Seven compounds with anti-inflammatory properties were isolated from the MESD. Among them, 3'-O-acetylhamaudol and 3'-O-angeloylhamaudol exhibited strong anti-inflammatory effects in vitro.

Conclusion

Overall, MESD may be a potential natural product for the treatment of IBD by lowering inflammatory cytokine levels and altering gut microbiota composition.

The photic-aphotic divide is a strong ecological and evolutionary force determining the distribution of ciliates (Alveolata, Ciliophora) in the ocean

The bulk of knowledge on marine ciliates is from shallow and/or sunlit waters. We studied ciliate diversity and distribution across epi- and mesopelagic oceanic waters, using DNA metabarcoding and phylogeny-based metrics. We analyzed sequences of the 18S rRNA gene (V4 region) from 369 samples collected at 12 depths (0-1000 m) at the Bermuda Atlantic Time-series Study site of the Sargasso Sea (North Atlantic) monthly for 3 years. The comprehensive depth and temporal resolutions analyzed led to three main findings. First, there was a gradual but significant decrease in alpha-diversity (based on Faith's phylogenetic diversity index) from surface to 1000-m waters. Second, multivariate analyses of beta-diversity (based on UniFrac distances) indicate that ciliate assemblages change significantly from photic to aphotic waters, with a switch from Oligotrichea to Oligohymenophorea prevalence. Third, phylogenetic placement of sequence variants and clade-level correlations (EPA-ng and GAPPA algorithms) show Oligotrichea, Litostomatea, Prostomatea, and Phyllopharyngea as anti-correlated with depth, while Oligohymenophorea (especially Apostomatia) have a direct relationship with depth. Two enigmatic environmental clades include either prevalent variants widely distributed in aphotic layers (the Oligohymenophorea OLIGO5) or subclades differentially distributed in photic versus aphotic waters (the Discotrichidae NASSO1). These results settle contradictory relationships between ciliate alpha-diversity and depth reported before, suggest functional changes in ciliate assemblages from photic to aphotic waters (with the prevalence of algivory and mixotrophy vs. omnivory and parasitism, respectively), and indicate that contemporary taxon distributions in the vertical profile have been strongly influenced by evolutionary processes. Integration of DNA sequences with organismal data (microscopy, functional experiments) and development of databases that link these sources of information remain as major tasks to better understand ciliate diversity, ecological roles, and evolution in the ocean.

Association of gut microbial dysbiosis with disease severity, response to therapy and disease outcomes in Indian patients with COVID-19

BACKGROUND

Severe coronavirus disease 2019 (COVID-19) is associated with systemic hyper-inflammation. An adaptive interaction between gut microbiota and host immune systems is important for intestinal homeostasis and systemic immune regulation. The association of gut microbial composition and functions with COVID-19 disease severity is sparse, especially in India. We analysed faecal microbial diversity and abundances in a cohort of Indian COVID-19 patients to identify key signatures in the gut microbial ecology in patients with severe COVID-19 disease as well as in response to different therapies. The composition of the gut microbiome was characterized using 16Sr RNA gene sequences of genomic DNA extracted from faecal samples of 52 COVID-19 patients. Metabolic pathways across the groups were predicted using PICRUSt2. All statistical analyses were done using Vegan in the R environment. Plasma cytokine abundance at recruitment was measured in a multiplex assay.

RESULTS

The gut microbiome composition of mild and severe patients was found to be significantly different. Immunomodulatory commensals, viz. Lachnospiraceae family members and Bifidobacteria producing butyrate and short-chain fatty acids (SCFAs), were under represented in patients with severe COVID-19, with an increased abundance of opportunistic pathogens like Eggerthella. The higher abundance of Lachnoclostridium in severe disease was reduced in response to convalescent plasma therapy. Specific microbial genera showed distinctive trends in enriched metabolic pathways, strong correlations with blood plasma cytokine levels, and associative link to disease outcomes.

CONCLUSION

Our study indicates that, along with SARS-CoV-2, a dysbiotic gut microbial community may also play an important role in COVID-19 severity through modulation of host immune responses.

Lactobacillus gasseri LG-G12 Restores Gut Microbiota and Intestinal Health in Obesity Mice on Ceftriaxone Therapy

Gut microbiota imbalance is associated with the occurrence of metabolic diseases such as obesity. Thus, its modulation is a promising strategy to restore gut microbiota and improve intestinal health in the obese. This paper examines the role of probiotics, antimicrobials, and diet in modulating gut microbiota and improving intestinal health. Accordingly, obesity was induced in C57BL/6J mice, after which they were redistributed and fed with an obesogenic diet (intervention A) or standard AIN-93 diet (intervention B). Concomitantly, all the groups underwent a treatment phase with Lactobacillus gasseri LG-G12, ceftriaxone, or ceftriaxone followed by L. gasseri LG-G12. At the end of the experimental period, the following analysis was conducted: metataxonomic analysis, functional profiling of gut microbiota, intestinal permeability, and caecal concentration of short-chain fatty acids. High-fat diet impaired bacterial diversity/richness, which was counteracted in association with L. gasseri LG-G12 and the AIN-93 diet. Additionally, SCFA-producing bacteria were negatively correlated with high intestinal permeability parameters, which was further confirmed via functional profile prediction of the gut microbiota. A novel perspective on anti-obesity probiotics is presented by these findings based on the improvement of intestinal health irrespective of undergoing antimicrobial therapy or not.

Network analysis of 16S rRNA sequences suggests microbial keystone taxa contribute to marine N2O cycling

The mechanisms by which large-scale microbial community function emerges from complex ecological interactions between individual taxa and functional groups remain obscure. We leveraged network analyses of 16S rRNA amplicon sequences obtained over a seven-month timeseries in seasonally anoxic Saanich Inlet (Vancouver Island, Canada) to investigate relationships between microbial community structure and water column N2O cycling. Taxa separately broadly into three discrete subnetworks with contrasting environmental distributions. Oxycline subnetworks were structured around keystone aerobic heterotrophs that correlated with nitrification rates and N2O supersaturations, linking N2O production and accumulation to taxa involved in organic matter remineralization. Keystone taxa implicated in anaerobic carbon, nitrogen, and sulfur cycling in anoxic environments clustered together in a low-oxygen subnetwork that correlated positively with nitrification N2O yields and N2O production from denitrification. Close coupling between N2O producers and consumers in the anoxic basin is indicated by strong correlations between the low-oxygen subnetwork, PICRUSt2-predicted nitrous oxide reductase (nosZ) gene abundances, and N2O undersaturation. This study implicates keystone taxa affiliated with common ODZ groups as a potential control on water column N2O cycling and provides a theoretical basis for further investigations into marine microbial interaction networks.

Phenol and Polyaromatic Hydrocarbons Are Stronger Drivers Than Host Plant Species in Shaping the Arbuscular Mycorrhizal Fungal Component of the Mycorrhizosphere

Changes in soil microbial communities in response to hydrocarbon pollution are critical indicators of disturbed ecosystem conditions. A core component of these communities that is functionally adjusted to the life-history traits of the host and environmental factors consists of arbuscular mycorrhizal fungi (AMF). AMF communities associated with Poa trivialis and Phragmites australis growing at a phenol and polynuclear aromatic hydrocarbon (PAH)-contaminated site and at an uncontaminated site were compared based on LSU rDNA sequencing. Dissimilarities in species composition and community structures indicated soil pollution as the main factor negatively affecting the AMF diversity. The AMF communities at the contaminated site were dominated by fungal generalists (Rhizophagus, Funneliformis, Claroideoglomus, Paraglomus) with wide ecological tolerance. At the control site, the AMF communities were characterized by higher taxonomic and functional diversity than those exposed to the contamination. The host plant identity was the main driver distinguishing the two AMF metacommunities. The AMF communities at the uncontaminated site were represented by Polonospora, Paraglomus, Oehlia, Nanoglomus, Rhizoglomus, Dominikia, and Microdominikia. Polonosporaceae and Paraglomeraceae were particularly dominant in the Ph. australis mycorrhizosphere. The high abundance of early diverging AMF could be due to the use of primers able to detect lineages such as Paraglomeracae that have not been recognized by previously used 18S rDNA primers.

Assembly of root-associated N2O-reducing communities of annual crops is governed by selection for nosZ clade I over clade II

The rhizosphere is a hotspot for denitrification. The nitrous oxide (N2O) reductase among denitrifiers and nondenitrifying N2O reducers is the only known N2O sink in the biosphere. We hypothesized that the composition of root-associated N2O-reducing communities when establishing on annual crops depend on soil type and plant species, but that assembly processes are independent of these factors and differ between nosZ clades I and II. Using a pot experiment with barley and sunflower and two soils, we analyzed the abundance, composition, and diversity of soil and root-associated N2O reducing communities by qPCR and amplicon sequencing of nosZ. Clade I was more abundant on roots compared to soil, while clade II showed the opposite. In barley, this pattern coincided with N2O availability, determined as potential N2O production rates, but for sunflower no N2O production was detected in the root compartment. Root and soil nosZ communities differed in composition and phylogeny-based community analyses indicated that assembly of root-associated N2O reducers was driven by the interaction between plant and soil type, with inferred competition being more influential than habitat selection. Selection between clades I and II in the root/soil interface is suggested, which may have functional consequences since most clade I microorganisms can produce N2O.

Cultivating marine bacteria under laboratory conditions: Overcoming the “unculturable” dogma

Underexplored seawater environments may contain biological resources with potential for new biotechnological applications. Metagenomic techniques revolutionized the study of bacterial communities but culture dependent methods will still be important to help the biodiscovery of new products and enzymes from marine bacteria. In this context, we promoted the growth of bacteria from a marine rock pond by culture dependent techniques and compared the results with culture independent methods. The total number of bacteria and diversity were studied in different agar plate media during 6 weeks. Agar plate counting was of the same order of magnitude of direct microscopy counts. The highest efficiency of cultivation was 45% attained in marine agar medium. Molecular analysis revealed 10 different phyla of which only four were isolated by the culture dependent method. On the other hand, four taxonomic orders were detected by cultivation but not by the molecular technique. These include bacteria from the phyla Bacillota and Actinomycetota. Our study shows that it is possible to grow more than the traditionally considered 1% of bacteria from a seawater sample using standard agar plate techniques and laboratorial conditions. The results also demonstrate the importance of culture methods to grow bacteria not detected by molecular approaches for future biotechnological applications.

Metagenomic Analysis Using Phylogenetic Placement-A Review of the First Decade

Phylogenetic placement refers to a family of tools and methods to analyze, visualize, and interpret the tsunami of metagenomic sequencing data generated by high-throughput sequencing. Compared to alternative (e. g., similarity-based) methods, it puts metabarcoding sequences into a phylogenetic context using a set of known reference sequences and taking evolutionary history into account. Thereby, one can increase the accuracy of metagenomic surveys and eliminate the requirement for having exact or close matches with existing sequence databases. Phylogenetic placement constitutes a valuable analysis tool per se, but also entails a plethora of downstream tools to interpret its results. A common use case is to analyze species communities obtained from metagenomic sequencing, for example via taxonomic assignment, diversity quantification, sample comparison, and identification of correlations with environmental variables. In this review, we provide an overview over the methods developed during the first 10 years. In particular, the goals of this review are 1) to motivate the usage of phylogenetic placement and illustrate some of its use cases, 2) to outline the full workflow, from raw sequences to publishable figures, including best practices, 3) to introduce the most common tools and methods and their capabilities, 4) to point out common placement pitfalls and misconceptions, 5) to showcase typical placement-based analyses, and how they can help to analyze, visualize, and interpret phylogenetic placement data.

Spatial and temporal dynamics of SAR11 marine bacteria across a nearshore to offshore transect in the tropical Pacific Ocean

Surveys of microbial communities across transitions coupled with contextual measures of the environment provide a useful approach to dissect the factors determining distributions of microorganisms across ecological niches. Here, monthly time-series samples of surface seawater along a transect spanning the nearshore coastal environment within Kāneʻohe Bay on the island of Oʻahu, Hawaiʻi, and the adjacent offshore environment were collected to investigate the diversity and abundance of SAR11 marine bacteria (order Pelagibacterales) over a 2-year time period. Using 16S ribosomal RNA gene amplicon sequencing, the spatiotemporal distributions of major SAR11 subclades and exact amplicon sequence variants (ASVs) were evaluated. Seven of eight SAR11 subclades detected in this study showed distinct subclade distributions across the coastal to offshore environments. The SAR11 community was dominated by seven (of 106 total) SAR11 ASVs that made up an average of 77% of total SAR11. These seven ASVs spanned five different SAR11 subclades (Ia, Ib, IIa, IV, and Va), and were recovered from all samples collected from either the coastal environment, the offshore, or both. SAR11 ASVs were more often restricted spatially to coastal or offshore environments (64 of 106 ASVs) than they were shared among coastal, transition, and offshore environments (39 of 106 ASVs). Overall, offshore SAR11 communities contained a higher diversity of SAR11 ASVs than their nearshore counterparts, with the highest diversity within the little-studied subclade IIa. This study reveals ecological differentiation of SAR11 marine bacteria across a short physiochemical gradient, further increasing our understanding of how SAR11 genetic diversity partitions into distinct ecological units.

App-SpaM: phylogenetic placement of short reads without sequence alignment

Motivation

Phylogenetic placement is the task of placing a query sequence of unknown taxonomic origin into a given phylogenetic tree of a set of reference sequences. A major field of application of such methods is, for example, the taxonomic identification of reads in metabarcoding or metagenomic studies. Several approaches to phylogenetic placement have been proposed in recent years. The most accurate of them requires a multiple sequence alignment of the references as input. However, calculating multiple alignments is not only time-consuming but also limits the applicability of these approaches.

Results

Herein, we propose Alignment-free phylogenetic placement algorithm based on Spaced-word Matches (App-SpaM), an efficient algorithm for the phylogenetic placement of short sequencing reads on a tree of a set of reference sequences. App-SpaM produces results of high quality that are on a par with the best available approaches to phylogenetic placement, while our software is two orders of magnitude faster than these existing methods. Our approach neither requires a multiple alignment of the reference sequences nor alignments of the queries to the references. This enables App-SpaM to perform phylogenetic placement on a broad variety of datasets.

Availability and implementation

The source code of App-SpaM is freely available on Github at https://github.com/matthiasblanke/App-SpaM together with detailed instructions for installation and settings. App-SpaM is furthermore available as a Conda-package on the Bioconda channel.

Contact

matthias.blanke@biologie.uni-goettingen.de.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Taxonomic and Functional Dysregulation in Salivary Microbiomes During Oral Carcinogenesis

Exploring microbial community compositions in humans with healthy versus diseased states is crucial to understand the microbe-host interplay associated with the disease progression. Although the relationship between oral cancer and microbiome was previously established, it remained controversial, and yet the ecological characteristics and their responses to oral carcinogenesis have not been well studied. Here, using the bacterial 16S rRNA gene amplicon sequencing along with the in silico function analysis by PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2), we systematically characterized the compositions and the ecological drivers of saliva microbiome in the cohorts of orally healthy, non-recurrent oral verrucous hyperplasia (a pre-cancer lesion), and oral verrucous hyperplasia–associated oral cancer at taxonomic and function levels, and compared them with the re-analysis of publicly available datasets. Diversity analyses showed that microbiome dysbiosis in saliva was significantly linked to oral health status. As oral health deteriorated, the number of core species declined, and metabolic pathways predicted by PICRUSt2 were dysregulated. Partitioned beta-diversity revealed an extremely high species turnover but low function turnover. Functional beta-diversity in saliva microbiome shifted from turnover to nestedness during oral carcinogenesis, which was not observed at taxonomic levels. Correspondingly, the quantitative analysis of stochasticity ratios showed that drivers of microbial composition and functional gene content of saliva microbiomes were primarily governed by the stochastic processes, yet the driver of functional gene content shifted toward deterministic processes as oral cancer developed. Re-analysis of publicly accessible datasets supported not only the distinctive family taxa of Veillonellaceae and Actinomycetaceae present in normal cohorts but also that Flavobacteriaceae and Peptostreptococcaceae as well as the dysregulated metabolic pathways of nucleotides, amino acids, fatty acids, and cell structure were related to oral cancer. Using predicted functional profiles to elucidate the correlations to the oral health status shows superior performance than using taxonomic data among different studies. These findings advance our understanding of the oral ecosystem in relation to oral carcinogenesis and provide a new direction to the development of microbiome-based tools to study the interplay of the oral microbiome, metabolites, and host health.

Ecological memory of recurrent drought modifies soil processes via changes in soil microbial community

Climate change is altering the frequency and severity of drought events. Recent evidence indicates that drought may produce legacy effects on soil microbial communities. However, it is unclear whether precedent drought events lead to ecological memory formation, i.e., the capacity of past events to influence current ecosystem response trajectories. Here, we utilize a long-term field experiment in a mountain grassland in central Austria with an experimental layout comparing 10 years of recurrent drought events to a single drought event and ambient conditions. We show that recurrent droughts increase the dissimilarity of microbial communities compared to control and single drought events, and enhance soil multifunctionality during drought (calculated via measurements of potential enzymatic activities, soil nutrients, microbial biomass stoichiometry and belowground net primary productivity). Our results indicate that soil microbial community composition changes in concert with its functioning, with consequences for soil processes. The formation of ecological memory in soil under recurrent drought may enhance the resilience of ecosystem functioning against future drought events.

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

Although ciliates are one of the most dominant microbial eukaryotic groups in many environments, there is a lack of updated global ciliate alignments and reference trees that can be used for phylogenetic placement methods to analyze environmental metabarcoding data. Here we fill this gap by providing reference alignments and trees for those ciliates taxa with available SSU-rDNA sequences derived from identified species. Each alignment contains 478 ciliate and six outgroup taxa, and they were made using different masking strategies for alignment positions (unmasked, masked and masked except the hypervariable V4 region). We constrained the monophyly of the major ciliate groups based on the recently updated classification of protists and based on phylogenomic data. Taxa of uncertain phylogenetic position were kept unconstrained, except for Mesodinium species that we constrained to form a clade with the Litostomatea. These ciliate reference alignments and trees can be used to perform taxonomic assignments of metabarcoding data, discover novel ciliate clades, estimate species richness, and overlay measured ecological parameters onto the phylogenetic placements.

PEWO: a collection of workflows to benchmark phylogenetic placement

MOTIVATION

Phylogenetic placement (PP) is a process of taxonomic identification for which several tools are now available. However, it remains difficult to assess which tool is more adapted to particular genomic data or a particular reference taxonomy. We developed Placement Evaluation WOrkflows (PEWO), the first benchmarking tool dedicated to PP assessment. Its automated workflows can evaluate PP at many levels, from parameter optimization for a particular tool, to the selection of the most appropriate genetic marker when PP-based species identifications are targeted. Our goal is that PEWO will become a community effort and a standard support for future developments and applications of PP.

AVAILABILITY AND IMPLEMENTATION

https://github.com/phylo42/PEWO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Microbiome Aggregated Traits and Assembly Are More Sensitive to Soil Management than Diversity

How soil is managed, particularly for agriculture, exerts stresses upon soil microbiomes, resulting in altered community structures and functional states. Understanding how soil microbiomes respond to combined stresses is important for predicting system performance under different land use scenarios, aids in identification of the most environmentally benign managements, and provides insight into how system function can be recovered in degraded soils. We use a long-established field experiment to study the effects of combined chronic (press) disturbance of the magnitude of organic carbon inputs with acute (pulse) effects of physical disturbance by tillage and chemical disturbance due to inorganic fertilization and pesticide application. We show that because of the variety of ways it can be assessed, biodiversity-here based on microbial small subunit rRNA gene phylotypes-does not provide a consistent view of community change. In contrast, aggregated traits associated with soil microbiomes indicate general loss of function, measured as a reduction of average genome lengths, associated with chronic reduction of organic inputs in arable or bare fallow soils and altered growth strategies associated with rRNA operon copy number in prokaryotes, as well as a switch to pathogenicity in fungal communities. In addition, pulse disturbance by soil tillage is associated with an increased influence of stochastic processes upon prokaryote community assembly, but fungicide used in arable soils results in niche assembly of fungal communities compared to untilled grassland. Overall, bacteria, archaea, and fungi do not share a common response to land management change, and estimates of biodiversity do not capture important facets of community adaptation to stresses adequately. IMPORTANCE Changes in soil microbiome diversity and function brought about by land management are predicted to influence a range of environmental services provided by soil, including provision of food and clean water. However, opportunities to compare the long-term effects of combinations of stresses imposed by different management approaches are limited. We exploit a globally unique 50-year field experiment, demonstrating that soil management practices alter microbiome diversity, community traits, and assembly. Grassland soil microbiomes are dominated by fewer-but phylogenetically more diverse-prokaryote phylotypes which sustain larger genomes than microbiomes in arable or bare fallow soil maintained free of plants. Dominant fungi in grassland soils are less phylogenetically diverse than those in arable or fallow soils. Soil tillage increases stochastic processes in microbiome assembly: this, combined with reduced plant biomass, presents opportunities for organisms with a capacity for pathogenesis to become established in stressed soils.

A multi-OMIC characterisation of biodegradation and microbial community succession within the PET plastisphere

BACKGROUND

Plastics now pollute marine environments across the globe. On entering these environments, plastics are rapidly colonised by a diverse community of microorganisms termed the plastisphere. Members of the plastisphere have a myriad of diverse functions typically found in any biofilm but, additionally, a number of marine plastisphere studies have claimed the presence of plastic-biodegrading organisms, although with little mechanistic verification. Here, we obtained a microbial community from marine plastic debris and analysed the community succession across 6 weeks of incubation with different polyethylene terephthalate (PET) products as the sole carbon source, and further characterised the mechanisms involved in PET degradation by two bacterial isolates from the plastisphere.

RESULTS

We found that all communities differed significantly from the inoculum and were dominated by Gammaproteobacteria, i.e. Alteromonadaceae and Thalassospiraceae at early time points, Alcanivoraceae at later time points and Vibrionaceae throughout. The large number of encoded enzymes involved in PET degradation found in predicted metagenomes and the observation of polymer oxidation by FTIR analyses both suggested PET degradation was occurring. However, we were unable to detect intermediates of PET hydrolysis with metabolomic analyses, which may be attributed to their rapid depletion by the complex community. To further confirm the PET biodegrading potential within the plastisphere of marine plastic debris, we used a combined proteogenomic and metabolomic approach to characterise amorphous PET degradation by two novel marine isolates, Thioclava sp. BHET1 and Bacillus sp. BHET2. The identification of PET hydrolytic intermediates by metabolomics confirmed that both isolates were able to degrade PET. High-throughput proteomics revealed that whilst Thioclava sp. BHET1 used the degradation pathway identified in terrestrial environment counterparts, these were absent in Bacillus sp. BHET2, indicating that either the enzymes used by this bacterium share little homology with those characterised previously, or that this bacterium uses a novel pathway for PET degradation.

CONCLUSIONS

Overall, the results of our multi-OMIC characterisation of PET degradation provide a significant step forwards in our understanding of marine plastic degradation by bacterial isolates and communities and evidences the biodegrading potential extant in the plastisphere of marine plastic debris. Video abstract.

Acute cyanotoxin poisoning reveals a marginal effect on mouse gut microbiome composition but indicates metabolic shifts related to liver and gut inflammation

Freshwater harmful algal blooms (HABs) are a major environmental health problem worldwide. HABs are caused by a predominance of cyanobacteria, some of which produce potent toxins. The most ubiquitous cyanotoxin is microcystin (MC) and the congener MC-LR is the most studied due to its toxicity. Short-term exposure to toxins can cause gut microbiome disturbances, but this has not been well described with MC-LR exposure. This study investigated the gut microbial communities of mice from a prior study, which identified significant liver toxicity from ingestion of MC-LR daily for 8 days. CD-1 mice were divided into three dosage groups: control, low exposure (sub-lethal MC-LR concentration), and high exposure (near-lethal MC-LR concentration). Fecal samples were analyzed using 16S rRNA sequencing. Results revealed that at population level, there were no significant shifts in bacterial diversity or the microbial community structure over the exposure period. However, there were significant differences between male and female mice. Predictive functional gene analysis indicated that several metabolic pathways were significantly different in the high dose group before exposure and following 7 doses of MC-LR, as well as between the control and high dose groups on Day 8. Significant differentially abundant taxa were also identified contributing to these pathways. Several pathways, including superpathway of N-acetylneuraminate degradation, were related to liver and gut inflammation. The outcome of this study suggests a need for in-depth investigation of metabolic activity and other functions in the gut in future studies, as well as potential consideration of the role of sex in MC-LR toxicity.

Crop type exerts greater influence upon rhizosphere phosphohydrolase gene abundance and phylogenetic diversity than phosphorus fertilization

Rock phosphate is an alternative form of phosphorus (P) fertilizer; however, there is no information regarding the influence of P fertilizer sources in Brazilian Cerrado soils upon microbial genes coding for phosphohydrolase enzymes in crop rhizospheres. Here, we analyze a field experiment comparing maize and sorghum grown under different P fertilization (rock phosphate and triple superphosphate) upon crop performance, phosphatase activity and rhizosphere microbiomes at three levels of diversity: small subunit rRNA marker genes of bacteria, archaea and fungi; a suite of alkaline and acid phosphatase and phytase genes; and ecotypes of individual genes. We found no significant difference in crop performance between the fertilizer sources, but the accumulation of fertilizer P into pools of organic soil P differed. Phosphatase activity was the only biological parameter influenced by P fertilization. Differences in rhizosphere microbiomes were observed at all levels of biodiversity due to crop type, but not fertilization. Inspection of phosphohydrolase gene ecotypes responsible for differences between the crops suggests a role for lateral genetic transfer in establishing ecotype distributions. Moreover, they were not reflected in microbial community composition, suggesting that they confer competitive advantage to individual cells rather than species in the sorghum rhizosphere.

A plasmid locus associated with Klebsiella clinical infections encodes a microbiome-dependent gut fitness factor

Klebsiella pneumoniae (Kp) is an important cause of healthcare-associated infections, which increases patient morbidity, mortality, and hospitalization costs. Gut colonization by Kp is consistently associated with subsequent Kp disease, and patients are predominantly infected with their colonizing strain. Our previous comparative genomics study, between disease-causing and asymptomatically colonizing Kp isolates, identified a plasmid-encoded tellurite (TeO3-2)-resistance (ter) operon as strongly associated with infection. However, TeO3-2 is extremely rare and toxic to humans. Thus, we used a multidisciplinary approach to determine the biological link between ter and Kp infection. First, we used a genomic and bioinformatic approach to extensively characterize Kp plasmids encoding the ter locus. These plasmids displayed substantial variation in plasmid incompatibility type and gene content. Moreover, the ter operon was genetically independent of other plasmid-encoded virulence and antibiotic resistance loci, both in our original patient cohort and in a large set (n = 88) of publicly available ter operon-encoding Kp plasmids, indicating that the ter operon is likely playing a direct, but yet undescribed role in Kp disease. Next, we employed multiple mouse models of infection and colonization to show that 1) the ter operon is dispensable during bacteremia, 2) the ter operon enhances fitness in the gut, 3) this phenotype is dependent on the colony of origin of mice, and 4) antibiotic disruption of the gut microbiota eliminates the requirement for ter. Furthermore, using 16S rRNA gene sequencing, we show that the ter operon enhances Kp fitness in the gut in the presence of specific indigenous microbiota, including those predicted to produce short chain fatty acids. Finally, administration of exogenous short-chain fatty acids in our mouse model of colonization was sufficient to reduce fitness of a ter mutant. These findings indicate that the ter operon, strongly associated with human infection, encodes factors that resist stress induced by the indigenous gut microbiota during colonization. This work represents a substantial advancement in our molecular understanding of Kp pathogenesis and gut colonization, directly relevant to Kp disease in healthcare settings.

An assessment of the persistence of putative pathogenic bacteria in chloraminated water distribution systems

This study investigated how a chloramine loss and nitrifying conditions influenced putative pathogenic bacterial diversity in bulk water and biofilm of a laboratory- and a full-scale chloraminated water distribution systems. Fifty-four reference databases containing full-length 16S rRNA gene sequences obtained from the National Centre for Biotechnology Information database were prepared to represent fifty-four pathogenic bacterial species listed in the World Health Organisation and Australian Drinking Water Quality Guidelines. When 16S rRNA gene sequences of all samples were screened against the fifty-four reference pathogenic databases, a total of thirty-one putative pathogenic bacteria were detected in both laboratory- and full-scale systems where total chlorine residuals ranged between 0.03 - 2.2 mg/L. Pathogenic bacterial species Mycolicibacterium fortuitum and Pseudomonas aeruginosa were noted in all laboratory (i.e. in bulk water and biofilm) and in bulk water of full-scale samples and Mycolicibacterium fortuitum dominated when chloramine residuals were high. Other different pathogenic bacterial species were observed dominant with decaying chloramine residuals. This study for the first time reports the diverse abundance of putative pathogenic bacteria resilient towards chloramine and highlights that metagenomics surveillance of drinking water can serve as a rapid assessment and an early warning of outbreaks of a large number of putative pathogenic bacteria.

Longitudinal assessment of the bovine ocular bacterial community dynamics in calves

BACKGROUND

Infectious Bovine Keratoconjunctivitis (IBK), commonly known as pinkeye, is one of the most significant diseases of beef cattle. As such, IBK costs the US beef industry at least 150 million annually. However, strategies to prevent IBK are limited, with most cases resulting in treatment with antibiotics once the disease has developed. Longitudinal studies evaluating establishment of the ocular microbiota may identify critical risk periods for IBK outbreaks or changes in the microbiota that may predispose animals to IBK.

RESULTS

In an attempt to characterize the establishment and colonization patterns of the bovine ocular microbiota, we conducted a longitudinal study consisting of 227 calves and evaluated the microbiota composition over time using amplicon sequence variants (ASVs) based on 16S rRNA sequencing data and culture-based approaches. Beef calves on trial consisted of both male (intact) and females. Breeds were composed of purebred Angus and composites with varying percentages of Simmental, Angus, and Red Angus breeds. Average age at the start of the trial was 65 days ±15.02 and all calves remained nursing on their dam until weaning (day 139 of the study). The trial consisted of 139 days with four sampling time points on day 0, 21, 41, and 139. The experimental population received three different vaccination treatments (autogenous, commercial (both inactivated bacteria), and adjuvant placebo), to assess the effectiveness of different vaccines for IBK prevention. A significant change in bacterial community composition was observed across time periods sampled compared to the baseline (p < 0.001). However, no treatment effect of vaccine was detected within the ocular bacterial community. The bacterial community composition with the greatest time span between sampling time periods (98d span) was most similar to the baseline sample collected, suggesting re-establishment of the ocular microbiota to baseline levels over time after perturbation. The effect of IgA levels on the microbial community was investigated in a subset of cattle within the study. However, no significant effect of IgA was observed. Significant changes in the ocular microbiota were identified when comparing communities pre- and post-clinical signs of IBK. Additionally, dynamic changes in opportunistic pathogens Moraxella spp. were observed and confirmed using culture based methods.

CONCLUSIONS

Our results indicate that the bovine ocular microbiota is well represented by opportunistic pathogens such as Moraxella and Mycoplasma. Furthermore, this study characterizes the diversity of the ocular microbiota in calves and demonstrates the plasticity of the ocular microbiota to change. Additionally, we demonstrate the ocular microbiome in calves is similar between the eyes and the perturbation of one eye results in similar changes in the other eye. We also demonstrate the bovine ocular microbiota is slow to recover post perturbation and as a result provide opportunistic pathogens a chance to establish within the eye leading to IBK and other diseases. Characterizing the dynamic nature of the ocular microbiota provides novel opportunities to develop potential probiotic intervention to reduce IBK outbreaks in cattle.

Metataxonomic profiling of bacterial communities and their predictive functional profiles in traditionally preserved meat products of Sikkim state in India

Traditionally preserved meat products are common food items in Sikkim state of India. We studied the high-throughput sequencing of four traditionally preserved meat products viz. beef kargyong, pork kargyong, yak satchu and khyopeh to profile the bacterial communities and also inferred their predictive functional profiles. Overall abundant OTUs in samples showed that Firmicutes was the abundant phylum followed by Proteobacteria and Bacteroidetes. Abundant species detected in each product were Psychrobacter pulmonis in beef kargyong, Lactobacillus sakei in pork kargyong, Bdellovibrio bacteriovorus and Ignatzschinera sp. in yak satchu and Lactobacillus sakei and Enterococcus sp. in khyopeh. Several genera unique to each product, based on analysis of shared OTUs contents, were observed among the samples except in khyopeh. Goods coverage recorded to 1.0 was observed, which reflected the maximum bacterial diversity in the samples. Alpha diversity metrics showed a maximum bacterial diversity in khyopeh and lowest in pork kargyong Community dissimilarities in the products were observed by PCoA plot. A total of 133 KEGG predictive functional pathways was observed in beef kargyong, 131 in pork kargyong, 125 in yak satchu and 101 in khyopeh. Metagenome contribution of the OTUs was computed using PICTRUSt2 and visualized by BURRITO software to predict the metabolic pathways. Several predictive functional profiles were contributed by abundant OTUs represented by Enterococcus, Acinetobacter, Agrobacterium, Bdellovibrio, Chryseobacterium, Lactococcus, Leuconostoc, Psychrobacter, and Staphylococcus.

Differential longitudinal establishment of human fecal bacterial communities in germ-free porcine and murine models

The majority of microbiome studies focused on understanding mechanistic relationships between the host and the microbiota have used mice and other rodents as the model of choice. However, the domestic pig is a relevant model that is currently underutilized for human microbiome investigations. In this study, we performed a direct comparison of the engraftment of fecal bacterial communities from human donors between human microbiota-associated (HMA) piglet and mouse models under identical dietary conditions. Analysis of 16S rRNA genes using amplicon sequence variants (ASVs) revealed that with the exception of early microbiota from infants, the more mature microbiotas tested established better in the HMA piglets compared to HMA mice. Of interest was the greater transplantation success of members belonging to phylum Firmicutes in the HMA piglets compared to the HMA mice. Together, these results provide evidence for the HMA piglet model potentially being more broadly applicable for donors with more mature microbiotas while the HMA mouse model might be more relevant for developing microbiotas such as those of infants. This study also emphasizes the necessity to exercise caution in extrapolating findings from HMA animals to humans, since up to 28% of taxa from some donors failed to colonize either model.

Genome Sequences of Bacteria Associated with the Diatom Cyclotella cryptica Strain CCMP332

We report draft genomes of two bacterial strains in the genera Hyphobacterium and Reichenbachiella, which are associated with the diatom Cyclotella cryptica strain CCMP332. Genomes from these strains were 2,691,501 bp and 3,325,829 bp in size, respectively, and will be useful for understanding interactions between diatoms and bacteria.

We report draft genomes of two bacterial strains in the genera Hyphobacterium and Reichenbachiella, which are associated with the diatom Cyclotella cryptica strain CCMP332. Genomes from these strains were 2,691,501 bp and 3,325,829 bp, respectively, and will be useful for understanding interactions between diatoms and bacteria.

Gut microbiota and metabolic health among overweight and obese individuals

Although obesity is associated with numerous diseases, the risks of disease may depend on metabolic health. Associations between the gut microbiota, obesity, and metabolic syndrome have been reported, but differences in microbiomes according to metabolic health in the obese population have not been explored in previous studies. Here, we investigated the composition of gut microbiota according to metabolic health status in obese and overweight subjects. A total of 747 overweight or obese adults were categorized by metabolic health status, and their fecal microbiota were profiled using 16S ribosomal RNA gene sequencing. We classified these adults into a metabolically healthy group (MH, N = 317) without any components of metabolic syndrome or a metabolically unhealthy group (MU, N = 430) defined as having at least one metabolic abnormality. The phylogenetic and non-phylogenetic alpha diversity for gut microbiota were lower in the MU group than the MH group, and there were significant differences in gut microbiota bacterial composition between the two groups. We found that the genus Oscillospira and the family Coriobacteriaceae were associated with good metabolic health in the overweight and obese populations. This is the first report to describe gut microbial diversity and composition in metabolically healthy and unhealthy overweight and obese individuals. Modulation of the gut microbiome may help prevent metabolic abnormalities in the obese population.

High-throughput sequence analysis of bacterial communities and their predictive functionalities in traditionally preserved fish products of Sikkim, India

Traditionally preserved fish products viz. suka ko maccha, a smoked fish product, sidra and sukuti, sun-dried fish products are commonly consumed in Sikkim state in India. Bacterial communities in these fish products were analysed by high-throughput sequence (HTS) method supported by bioinformatics tool. Metataxonomic of the overall bacterial communities in samples revealed the abundance of phylum Firmicutes followed by Proteobacteria. Psychrobacter was abundant genus in all traditionally preserved fish products of Sikkim, followed by Bacillus, Staphylococcus, Serratia, Clostridium, Enterobacter, Pseudomonas, Rummeliibacillus, Enterococcus, Photobacterium, Myroides, Peptostreptococcus, Plesiomonas and Achromobacter. Product-wise distribution showed that Bacillus was abundant in suka ko maacha and sidra samples, whereas Psychrobacter was abundant in sukuti samples. Unique genus to each product was observed on the basis of analysis of shared operational-taxonomic-unit (OTU) contents, Alpha diversity indices showed significantly differences among the samples, and also showed maximum coverage as per Good's coverage (0.99). Beta diversity showed clustering of bacterial compositions between suka ko maacha and sidra, whereas sukuti showed scattering pattern among the other samples, indicating a diverse population in suka ko maacha and sidra samples. Non-parametric analysis of abundant genera and predictive functionalities showed the complex bacterial inter-dependencies with predictive functionalities mostly in metabolism (79.88%).

Unraveling biogeochemical phosphorus dynamics in hyperarid Mars-analogue soils using stable oxygen isotopes in phosphate

With annual precipitation less than 20 mm and extreme UV intensity, the Atacama Desert in northern Chile has long been utilized as an analogue for recent Mars. In these hyperarid environments, water and biomass are extremely limited, and thus, it becomes difficult to generate a full picture of biogeochemical phosphate-water dynamics. To address this problem, we sampled soils from five Atacama study sites and conducted three main analyses-stable oxygen isotopes in phosphate, enzyme pathway predictions, and cell culture experiments. We found that high sedimentation rates decrease the relative size of the organic phosphorus pool, which appears to hinder extremophiles. Phosphoenzyme and pathway prediction analyses imply that inorganic pyrophosphatase is the most likely catalytic agent to cycle P in these environments, and this process will rapidly overtake other P utilization strategies. In these soils, the biogenic δ¹⁸ O signatures of the soil phosphate (δ¹⁸ O_PO4 ) can slowly overprint lithogenic δ¹⁸ O_PO4 values over a timescale of tens to hundreds of millions of years when annual precipitation is more than 10 mm. The δ¹⁸ O_PO4 of calcium-bound phosphate minerals seems to preserve the δ¹⁸ O signature of the water used for biogeochemical P cycling, pointing toward sporadic rainfall and gypsum hydration water as key moisture sources. Where precipitation is less than 2 mm, biological cycling is restricted and bedrock δ¹⁸ O_PO4 values are preserved. This study demonstrates the utility of δ¹⁸ O_PO4 values as indicative of biogeochemical cycling and hydrodynamics in an extremely dry Mars-analogue environment.

SCRAPP: A tool to assess the diversity of microbial samples from phylogenetic placements

Microbial ecology research is currently driven by the continuously decreasing cost of DNA sequencing and the improving accuracy of data analysis methods. One such analysis method is phylogenetic placement, which establishes the phylogenetic identity of the anonymous environmental sequences in a sample by means of a given phylogenetic reference tree. However, assessing the diversity of a sample remains challenging, as traditional methods do not scale well with the increasing data volumes and/or do not leverage the phylogenetic placement information. Here, we present scrapp, a highly parallel and scalable tool that uses a molecular species delimitation algorithm to quantify the diversity distribution over the reference phylogeny for a given phylogenetic placement of the sample. scrapp employs a novel approach to cluster phylogenetic placements, called placement space clustering, to efficiently perform dimensionality reduction, so as to scale on large data volumes. Furthermore, it uses the phylogeny‐aware molecular species delimitation method mPTP to quantify diversity. We evaluated scrapp using both, simulated and empirical data sets. We use simulated data to verify our approach. Tests on an empirical data set show that scrapp‐derived metrics can classify samples by their diversity‐correlated features equally well or better than existing, commonly used approaches. scrapp is available at https://github.com/pbdas/scrapp.

PhyloMagnet: fast and accurate screening of short-read meta-omics data using gene-centric phylogenetics

Motivation

Metagenomic and metatranscriptomic sequencing have become increasingly popular tools for producing massive amounts of short-read data, often used for the reconstruction of draft genomes or the detection of (active) genes in microbial communities. Unfortunately, sequence assemblies of such datasets generally remain a computationally challenging task. Frequently, researchers are only interested in a specific group of organisms or genes; yet, the assembly of multiple datasets only to identify candidate sequences for a specific question is sometimes prohibitively slow, forcing researchers to select a subset of available datasets to address their question. Here, we present PhyloMagnet, a workflow to screen meta-omics datasets for taxa and genes of interest using gene-centric assembly and phylogenetic placement of sequences.

Results

Using PhyloMagnet, we could identify up to 87% of the genera in an in vitro mock community with variable abundances, while the false positive predictions per single gene tree ranged from 0 to 23%. When applied to a group of metagenomes for which a set of metagenome assembled genomes (MAGs) have been published, we could detect the majority of the taxonomic labels that the MAGs had been annotated with. In a metatranscriptomic setting, the phylogenetic placement of assembled contigs corresponds to that of transcripts obtained from transcriptome assembly.

Availability and implementation

PhyloMagnet is built using Nextflow, available at github.com/maxemil/PhyloMagnet and is developed and tested on Linux. It is released under the open source GNU GPL licence and documentation is available at phylomagnet.readthedocs.io. Version 0.5 of PhyloMagnet was used for all benchmarking experiments.

Supplementary information

Supplementary data are available at Bioinformatics online.

Microbial insights from Antarctic and Mediterranean shallow-water bone-eating worms

Bone-eating worms of the genus Osedax (Annelida, Siboglinidae) form unique holobionts (functional entity comprising host and associated microbiota), highly adapted to inhabit bone tissue of marine vertebrates. These gutless worms have developed nutritional symbioses housing intracellular, horizontally acquired, heterotrophic bacteria hypothesised to harness nutrients from organic compounds, sequestered within the bone. Despite previous efforts, critical mechanisms mediating activity and acquisition of diverse bacterial assemblages remain unclear. Using 16S rRNA amplicon sequencing, we performed detailed taxonomic and predicted functional analyses shedding light on the microbial communities of two shallow-water Osedax species (Osedax deceptionensis and Osedax ‘mediterranea’) from contrasting habitats (Antarctic and Mediterranean Sea), in two tissue types (roots and palps). Comparative assessments between host species revealed distinct microbial assemblages whilst, within host species and body tissue, relative symbiont frequencies retained high variability. We reported relatively high abundances of microbes previously classified as primary endosymbionts, Ribotype 1 (order Oceanospirillales), and diverse likely secondary epibionts warranting further exploration as recurrent Osedax associates. Surprisingly, O. ‘mediterranea’ exhibited relatively low abundance of Oceanospirillales, but increased abundance of other potentially hydrocarbon degrading bacteria from the family Alteromonadaceae. We hypothesise the presence of functionally similar, non-Oceanospirillales primary endosymbionts within O. ‘mediterranea’. Functional metagenomic profiling (using 16S rRNA sequences) predicted broad metabolic capabilities, encompassing relatively large abundances of genes associated with amino acid metabolism. Comparative analyses between host body tissue communities highlighted several genes potentially providing critical functions to the Osedax host or that confer adaptations for intracellular life, housed within bone embedded host root tissues.

Diverse variola virus (smallpox) strains were widespread in northern Europe in the Viking Age

Smallpox, one of the most devastating human diseases, killed between 300 million and 500 million people in the 20th century alone. We recovered viral sequences from 13 northern European individuals, including 11 dated to ~600-1050 CE, overlapping the Viking Age, and reconstructed near-complete variola virus genomes for four of them. The samples predate the earliest confirmed smallpox cases by ~1000 years, and the sequences reveal a now-extinct sister clade of the modern variola viruses that were in circulation before the eradication of smallpox. We date the most recent common ancestor of variola virus to ~1700 years ago. Distinct patterns of gene inactivation in the four near-complete sequences show that different evolutionary paths of genotypic host adaptation resulted in variola viruses that circulated widely among humans.

Reproductive hormones mediate changes in the gut microbiome during pregnancy and lactation in Phayre's leaf monkeys

Studies in multiple host species have shown that gut microbial diversity and composition change during pregnancy and lactation. However, the specific mechanisms underlying these shifts are not well understood. Here, we use longitudinal data from wild Phayre's leaf monkeys to test the hypothesis that fluctuations in reproductive hormone concentrations contribute to gut microbial shifts during pregnancy. We described the microbial taxonomic composition of 91 fecal samples from 15 females (n = 16 cycling, n = 36 pregnant, n = 39 lactating) using 16S rRNA gene amplicon sequencing and assessed whether the resulting data were better explained by overall reproductive stage or by fecal estrogen (fE) and progesterone (fP) concentrations. Our results indicate that while overall reproductive stage affected gut microbiome composition, the observed patterns were driven by reproductive hormones. Females had lower gut microbial diversity during pregnancy and fP concentrations were negatively correlated with diversity. Additionally, fP concentrations predicted both unweighted and weighted UniFrac distances, while reproductive state only predicted unweighted UniFrac distances. Seasonality (rainfall and periods of phytoprogestin consumption) additionally influenced gut microbial diversity and composition. Our results indicate that reproductive hormones, specifically progestagens, contribute to the shifts in the gut microbiome during pregnancy and lactation.

Meta-Analysis of Gut Dysbiosis in Parkinson's Disease

BACKGROUND

PD may begin with the intestinal accumulation of α-synuclein fibrils, which can be causally associated with gut dysbiosis. The variability of gut microbiota across countries prevented us from identifying shared gut dysbiosis in PD.

OBJECTIVES

To identify gut dysbiosis in PD across countries.

METHODS

We performed 16S ribosomal RNA gene sequencing analysis of gut microbiota in 223 patients with PD and 137 controls, and meta-analyzed gut dysbiosis by combining our dataset with four previously reported data sets from the United States, Finland, Russia, and Germany. We excluded uncommon taxa from our analyses. For pathway analysis, we developed the Kyoto Encyclopedia of Genes and Genomes orthology set enrichment analysis method.

RESULTS

After adjusting for confounding factors (body mass index, constipation, sex, age, and catechol-O-methyl transferase inhibitor), genera Akkermansia and Catabacter, as well as families Akkermansiaceae, were increased, whereas genera Roseburia, Faecalibacterium, and Lachnospiraceae ND3007 group were decreased in PD. Catechol-O-methyl transferase inhibitor intake markedly increased family Lactobacillaceae. Inspection of these bacteria in 12 datasets that were not included in the meta-analysis revealed that increased genus Akkermansia and decreased genera Roseburia and Faecalibacterium were frequently observed across countries. Kyoto Encyclopedia of Genes and Genomes orthology set enrichment analysis revealed changes in short-chain fatty acid metabolisms in our dataset.

CONCLUSIONS

We report that intestinal mucin layer-degrading Akkermansia is increased and that short-chain fatty acid-producing Roseburia and Faecalibacterium are decreased in PD across countries. © 2020 International Parkinson and Movement Disorder Society.

Major Lipids, Apolipoproteins, and Alterations of Gut Microbiota

The gut microbiota has been linked to blood lipids. However, the relationship between the gut microbiome and other lipid markers like apolipoproteins A1 (apoA1) and B (apoB) as well as classical lipid markers in Asians remain unclear. Here, we examined the associations between gut microbial diversity and taxonomic compositions with both apolipoproteins and lipid markers in a large number of Korean patients. The fecal 16S rRNA gene sequencing data from 1141 subjects were analyzed and subjects were categorized into control group (G0) or abnormal group (G1) according to blood lipid measurements. The microbial diversity and several taxa of the gut microbiota were significantly associated with triglyceride, apoA1, and apoB levels, but not with total cholesterol, low-density lipoprotein cholesterol, and high-density lipoprotein cholesterol levels. The alpha diversity of the gut microbiota was inversely associated with high triglyceride level. Interestingly, G1 of apoA1 showed increased microbial richness and distinct microbial community compared with G0 of apoA1. A high abundance of Fusobacteria and low abundance of Oscillospira were found in the hypertriglyceridemia group. In this large-scale study, we identified associations of gut microbiota with apolipoproteins and classical lipid markers, indicating that the gut microbiota may be an important target for regulating blood lipids.

Tissue-scale microbiota of the Mediterranean mussel (Mytilus galloprovincialis) and its relationship with the environment

In this study, we characterize the structural variation of the microbiota of Mytilus galloprovincialis at the tissue scale, also exploring the connection with the microbial ecosystem of the surrounding water. Mussels were sampled within a farm located in the North-Western Adriatic Sea and microbiota composition was analyzed in gills, hemolymph, digestive glands, stomach and foot by Next Generation Sequencing marker gene approach. Mussels showed a distinctive microbiota structure, with specific declinations at the tissue level. Indeed, each tissue is characterized by a distinct pattern of dominant families, reflecting a peculiar adaptation to the respective tissue niche. For instance, the microbiota of the digestive gland is characterized by Ruminococcaceae and Lachnospiraceae, being shaped to ferment complex polysaccharides of dietary origin into short-chain fatty acids, well matching the general asset of the animal gut microbiota. Conversely, the gill and hemolymph ecosystems are dominated by marine microorganisms with aerobic oxidative metabolism, consistent with the role played by these tissues as an interface with the external environment. Our findings highlight the putative importance of mussel microbiota for different aspects of host physiology, with ultimate repercussions on mussel health and productivity.

Genesis and Gappa: processing, analyzing and visualizing phylogenetic (placement) data

SUMMARY

We present genesis, a library for working with phylogenetic data, and gappa, an accompanying command-line tool for conducting typical analyses on such data. The tools target phylogenetic trees and phylogenetic placements, sequences, taxonomies and other relevant data types, offer high-level simplicity as well as low-level customizability, and are computationally efficient, well-tested and field-proven.

AVAILABILITY AND IMPLEMENTATION

Both genesis and gappa are written in modern C++11, and are freely available under GPLv3 at http://github.com/lczech/genesis and http://github.com/lczech/gappa.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Vitamin B and vitamin D as modulators of gut microbiota in overweight individuals

Growing evidence indicates that vitamin B (VB) and vitamin D (VD) are associated with the development of obesity. The effects of supplemental dietary VB and VD on the gut microbiota of obese individuals are still largely unknown. In the present study, 997 normoweight (NW) and 773 overweight (OW) samples were screened from the American Gut Project. The microbial α-diversity was not affected by VB, VD or VB combined with VD (VBD) supplementation in the normoweight population or the overweigh population. VD significantly modulated 20 gut microbial metabolic pathways in NW. Moreover, VBD supplementation significantly affected 3 phyla and 3 families and 22 pathways in the OW gut microbiota. In summary, our results highlight that VD and VBD might be necessary for reshaping the gut microbiota in NW and OW, respectively.

Trimethylamine-N-Oxide Postprandial Response in Plasma and Urine Is Lower After Fermented Compared to Non-Fermented Dairy Consumption in Healthy Adults

Trimethylamine-N-oxide (TMAO) can be produced by the gut microbiota from dietary substrates and is associated with cardiovascular disease. While dairy products contain TMAO precursors, the effect of fermented dairy on TMAO metabolism remains unclear. We used plasma and urine samples collected for two randomised cross-over studies to evaluate the effects of fermented dairy consumption on TMAO metabolism. In Study 1, thirteen healthy young men tested a yogurt and an acidified milk during postprandial tests and a two-week daily intervention. In Study 2, ten healthy adults tested milk and cheese during postprandial tests. TMAO and five related metabolites were measured in plasma and urine by LC-MS/MS and NMR. Faecal microbiota composition was assessed in Study 1 (16S rRNA metagenomics sequencing). Fermented milk products were associated with lower postprandial TMAO responses than non-fermented milks in urine (Study 1, p = 0.01; Study 2, p = 0.02) and in plasma, comparing yogurt and acidified milk (Study 1, p = 0.04). Daily consumption of dairy products did not differentially affect fasting TMAO metabolites. Significant correlations were observed between microbiota taxa and circulating or urinary TMAO concentrations. Fermentation of dairy products appear, at least transiently, to affect associations between dairy products and circulating TMAO levels.

Prokaryotic and eukaryotic microbiomes associated with blooms of the ichthyotoxic dinoflagellate Cochlodinium (Margalefidinium) polykrikoides in New York, USA, estuaries

While harmful algal blooms caused by the ichthyotoxic dinoflagellate, Cochlodinium (Margalefidinium) polykrikoides, are allelopathic and may have unique associations with bacteria, a comprehensive assessment of the planktonic communities associated with these blooms has been lacking. Here, we used high-throughput amplicon sequencing to assess size fractionated (0.2 and 5 μm) bacterial (16S) and phytoplankton assemblages (18S) associated with blooms of C. polykrikoides during recurrent blooms in NY, USA. Over a three-year period, samples were collected inside (‘patch’) and outside (‘non-patch’) dense accumulations of C. polykrikoides to assess the microbiome associated with these blooms. Eukaryotic plankton communities of blooms had significantly lower diversity than non-bloom samples, and non-bloom samples hosted 30 eukaryotic operational taxonomic units (OTUs) not found within blooms, suggesting they may have been allelopathically excluded from blooms. Differential abundance analyses revealed that C. polykrikoides blooms were significantly enriched in dinoflagellates (p<0.001) and the experimental enrichment of C. polykrikoides led to a significant increase in the relative abundance of eight genera of dinoflagellates but a significant decline in other eukaryotic plankton. Amoebophrya co-dominated both within- and near- C. polykrikoides blooms and was more abundant in bloom patches. The core bacterial microbiome of the >0.2μm fraction of blooms was dominated by an uncultured bacterium from the SAR11 clade, while the >5μm size fraction was co-dominated by an uncultured bacterium from Rhodobacteraceae and Coraliomargarita. Two bacterial lineages within the >0.2μm fraction, as well as the Gammaproteobacterium, Halioglobus, from the >5μm fraction were unique to the microbiome of blooms, while there were 154 bacterial OTUs only found in non-bloom waters. Collectively, these findings reveal the unique composition and potential function of eukaryotic and prokaryotic communities associated with C. polykrikoides blooms.