Metaxa2 Diversity Tools: Easing microbial community analysis with Metaxa2

Similarity of Chinese and Pakistani oral microbiome

Oral microbiota is vital for human health and can be affected by various factors (i.e. diets, ethnicity). However, few studies have compared oral microbiota of individuals from different nationalities in the same environment. Here, we explored the assembly and interaction of oral microbial communities of Chinese and Pakistanis in one university. Firmicutes and Proteobacteria were the predominant microorganisms in the oral cavity of Chinese and Pakistanis. Streptococcus and Neisseria were the dominant genera of China, while Streptococcus and Haemophilus were the dominant genera of Pakistanis. In addition, the oral community membership and structure were not influenced by season, Chinese/Pakistani student and gender, reflecting the stability of the human oral microbiome. The beta diversity of oral microbiomes between Chinese and Pakistanis significantly differed in winter, but not in spring. The alpha diversity of Chinese students and Pakistani students was similar. Moreover, oral microbial community of both Chinese and Pakistani students was mainly driven by stochastic processes. The microbial network of Chinese was more complexity and stability than that of Pakistanis. Our study uncovers the characteristics of human oral microbiota, which is of great significance for oral and human health.

Microbial and Viral Genome and Proteome Nitrogen Demand Varies across Multiple Spatial Scales within a Marine Oxygen Minimum Zone

Nutrient availability can significantly influence microbial genomic and proteomic streamlining, for example, by selecting for lower nitrogen to carbon ratios. Oligotrophic open ocean microbes have streamlined genomic nitrogen requirements relative to those of their counterparts in nutrient-rich coastal waters. However, steep gradients in nutrient availability occur at meter-level, and even micron-level, spatial scales. It is unclear whether such gradients also structure genomic and proteomic stoichiometry. Focusing on the eastern tropical North Pacific oxygen minimum zone (OMZ), we use comparative metagenomics to examine how nitrogen availability shapes microbial and viral genome properties along the vertical gradient across the OMZ and between two size fractions, distinguishing free-living microbes versus particle-associated microbes. We find a substantial increase in the nitrogen content of encoded proteins in particle-associated over free-living bacteria and archaea across nitrogen availability regimes over depth. Within each size fraction, we find that bacterial and viral genomic nitrogen tends to increase with increasing nitrate concentrations with depth. In contrast to cellular genes, the nitrogen content of virus proteins does not differ between size fractions. We identified arginine as a key amino acid in the modulation of the C:N ratios of core genes for bacteria, archaea, and viruses. Functional analysis reveals that particle-associated bacterial metagenomes are enriched for genes that are involved in arginine metabolism and organic nitrogen compound catabolism. Our results are consistent with nitrogen streamlining in both cellular and viral genomes on spatial scales of meters to microns. These effects are similar in magnitude to those previously reported across scales of thousands of kilometers. IMPORTANCE The genomes of marine microbes can be shaped by nutrient cycles, with ocean-scale gradients in nitrogen availability being known to influence microbial amino acid usage. It is unclear, however, how genomic properties are shaped by nutrient changes over much smaller spatial scales, for example, along the vertical transition into oxygen minimum zones (OMZs) or from the exterior to the interior of detrital particles. Here, we measure protein nitrogen usage by marine bacteria, archaea, and viruses by using metagenomes from the nitracline of the eastern tropical North Pacific OMZ, including both particle-associated and nonassociated biomass. Our results show higher genomic and proteomic nitrogen content in particle-associated microbes and at depths with higher nitrogen availability for cellular and viral genomes. This discovery suggests that stoichiometry influences microbial and viral evolution across multiple scales, including the micrometer to millimeter scale associated with particle-associated versus free-living lifestyles.

Potato root-associated microbiomes adapt to combined water and nutrient limitation and have a plant genotype-specific role for plant stress mitigation

BACKGROUND

Due to climate change and reduced use of fertilizers combined stress scenarios are becoming increasingly frequent in crop production. In a field experiment we tested the effect of combined water and phosphorus limitation on the growth performance and plant traits of eight tetraploid and two diploid potato varieties as well as on root-associated microbiome diversity and functional potential. Microbiome and metagenome analysis targeted the diversity and potential functions of prokaryotes, fungi, plasmids, and bacteriophages and was linked to plant traits like tuber yield or timing of canopy closure.

RESULTS

The different potato genotypes responded differently to the combined stress and hosted distinct microbiota in the rhizosphere and the root endosphere. Proximity to the root, stress and potato genotype had significant effects on bacteria, whereas fungi were only mildly affected. To address the involvement of microbial functions, we investigated well and poorly performing potato genotypes (Stirling and Desirée, respectively) under stress conditions and executed a metagenome analysis of rhizosphere microbiomes subjected to stress and no stress conditions. Functions like ROS detoxification, aromatic amino acid and terpene metabolism were enriched and in synchrony with the metabolism of stressed plants. In Desirée, Pseudonocardiales had the genetic potential to take up assimilates produced in the fast-growing canopy and to reduce plant stress-sensing by degrading ethylene, but overall yield losses were high. In Stirling, Xanthomonadales had the genetic potential to reduce oxidative stress and to produce biofilms, potentially around roots. Biofilm formation could be involved in drought resilience and nutrient accessibility of Stirling and explain the recorded low yield losses. In the rhizosphere exposed to combined stress, the relative abundance of plasmids was reduced, and the diversity of phages was enriched. Moreover, mobile elements like plasmids and phages were affected by combined stresses in a genotype-specific manner.

CONCLUSION

Our study gives new insights into the interconnectedness of root-associated microbiota and plant stress responses in the field. Functional genes in the metagenome, phylogenetic composition and mobile elements play a role in potato stress adaption. In a poor and a well performing potato genotype grown under stress conditions, distinct functional genes pinpoint to a distinct stress sensing, water availability and compounds in the rhizospheres.

An improved workflow for accurate and robust healthcare environmental surveillance using metagenomics

BACKGROUND

Effective surveillance of microbial communities in the healthcare environment is increasingly important in infection prevention. Metagenomics-based techniques are promising due to their untargeted nature but are currently challenged by several limitations: (1) they are not powerful enough to extract valid signals out of the background noise for low-biomass samples, (2) they do not distinguish between viable and nonviable organisms, and (3) they do not reveal the microbial load quantitatively. An additional practical challenge towards a robust pipeline is the inability to efficiently allocate sequencing resources a priori. Assessment of sequencing depth is generally practiced post hoc, if at all, for most microbiome studies, regardless of the sample type. This practice is inefficient at best, and at worst, poor sequencing depth jeopardizes the interpretation of study results. To address these challenges, we present a workflow for metagenomics-based environmental surveillance that is appropriate for low-biomass samples, distinguishes viability, is quantitative, and estimates sequencing resources.

RESULTS

The workflow was developed using a representative microbiome sample, which was created by aggregating 120 surface swabs collected from a medical intensive care unit. Upon evaluating and optimizing techniques as well as developing new modules, we recommend best practices and introduce a well-structured workflow. We recommend adopting liquid-liquid extraction to improve DNA yield and only incorporating whole-cell filtration when the nonbacterial proportion is large. We suggest including propidium monoazide treatment coupled with internal standards and absolute abundance profiling for viability assessment and involving cultivation when demanding comprehensive profiling. We further recommend integrating internal standards for quantification and additionally qPCR when we expect poor taxonomic classification. We also introduce a machine learning-based model to predict required sequencing effort from accessible sample features. The model helps make full use of sequencing resources and achieve desired outcomes. Video Abstract CONCLUSIONS: This workflow will contribute to more accurate and robust environmental surveillance and infection prevention. Lessons gained from this study will also benefit the continuing development of methods in relevant fields.

Soil composition and rootstock genotype drive the root associated microbial communities in young grapevines

Soil microbiota plays a significant role in plant development and health and appears to be a major component of certain forms of grapevine decline. A greenhouse experiment was conducted to study the impact of the microbiological quality of the soil and grapevine rootstock genotype on the root microbial community and development of young plants. Two rootstocks heterografted with the same scion were grown in two vineyard soils differing in microbial composition and activities. After 4 months, culture-dependent approaches and amplicon sequencing of bacterial 16S rRNA gene and fungal ITS were performed on roots, rhizosphere and bulk soil samples. The root mycorrhizal colonization and number of cultivable microorganisms in the rhizosphere compartment of both genotypes were clearly influenced by the soil status. The fungal diversity and richness were dependent on the soil status and the rootstock, whereas bacterial richness was affected by the genotype only. Fungal genera associated with grapevine diseases were more abundant in declining soil and related root samples. The rootstock affected the compartmentalization of microbial communities, underscoring its influence on microorganism selection. Fluorescence in situ hybridization (FISH) confirmed the presence of predominant root-associated bacteria. These results emphasized the importance of rootstock genotype and soil composition in shaping the microbiome of young vines.

Simulated global warming affects endophytic bacterial and fungal communities of Antarctic pearlwort leaves and some bacterial isolates support plant growth at low temperatures

Antarctica is one of the most stressful environments for plant life and the Antarctic pearlwort (Colobanthus quitensis) is adapted to the hostile conditions. Plant-associated microorganisms can contribute to plant survival in cold environments, but scarce information is available on the taxonomic structure and functional roles of C. quitensis-associated microbial communities. This study aimed at evaluating the possible impacts of climate warming on the taxonomic structure of C. quitensis endophytes and at investigating the contribution of culturable bacterial endophytes to plant growth at low temperatures. The culture-independent analysis revealed changes in the taxonomic structure of bacterial and fungal communities according to plant growth conditions, such as the collection site and the presence of open-top chambers (OTCs), which can simulate global warming. Plants grown inside OTCs showed lower microbial richness and higher relative abundances of biomarker bacterial genera (Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium, Aeromicrobium, Aureimonas, Hymenobacter, Novosphingobium, Pedobacter, Pseudomonas and Sphingomonas) and fungal genera (Alternaria, Cistella, and Vishniacozyma) compared to plants collected from open areas (OA), as a possible response to global warming simulated by OTCs. Culturable psychrotolerant bacteria of C. quitensis were able to endophytically colonize tomato seedlings and promote shoot growth at low temperatures, suggesting their potential contribution to plant tolerance to cold conditions.

Metabolite Production in Alkanna tinctoria Links Plant Development with the Recruitment of Individual Members of Microbiome Thriving at the Root-Soil Interface

Plants are naturally associated with diverse microbial communities, which play significant roles in plant performance, such as growth promotion or fending off pathogens. The roots of Alkanna tinctoria L. are rich in naphthoquinones, particularly the medicinally used enantiomers alkannin and shikonin and their derivatives. Former studies already have shown that microorganisms may modulate plant metabolism. To further investigate the potential interaction between A. tinctoria and associated microorganisms, we performed a greenhouse experiment in which A. tinctoria plants were grown in the presence of three distinct soil microbiomes. At four defined plant developmental stages, we made an in-depth assessment of bacterial and fungal root-associated microbiomes as well as all extracted primary and secondary metabolite content of root material. Our results showed that the plant developmental stage was the most important driver influencing the plant metabolite content, revealing peak contents of alkannin/shikonin derivatives at the fruiting stage. Plant root microbial diversity was influenced both by bulk soil origin and to a small extent by the developmental stage. The performed correlation analyses and cooccurrence networks on the measured metabolite content and the abundance of individual bacterial and fungal taxa suggested a dynamic and at times positive or negative relationship between root-associated microorganisms and root metabolism. In particular, the bacterial genera Labrys and Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium as well as four species of the fungal genus Penicillium were found to be positively correlated with higher content of alkannins. IMPORTANCE Previous studies have shown that individual, isolated microorganisms may influence secondary metabolism of plants and induce or stimulate the production of medicinally relevant secondary metabolism. Here, we analyzed the microbiome-metabolome linkage of the medicinal plant Alkanna tinctoria, which is known to produce valuable compounds, particularly the naphthoquinones alkannin and shikonin and their derivatives. A detailed bacterial and fungal microbiome and metabolome analysis of A. tinctoria roots revealed that the plant developmental stage influenced root metabolite production, whereas soil inoculants from three different geographical origins in which plants were grown shaped root-associated microbiota. Metabolomes of plant roots of the same developmental stage across different soils were highly similar, pinpointing to plant maturity as the primary driver of secondary metabolite production. Correlation and network analyses identified bacterial and fungal taxa showing a positive relationship between root-associated microorganisms and root metabolism. In particular, the bacterial genera Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium and Labrys as well as the fungal species of genus Penicillium were found to be positively correlated with higher content of alkannins.

Metabolic Diversity and Aero-Tolerance in Anammox Bacteria from Geochemically Distinct Aquifers

Anaerobic ammonium oxidation (anammox) is important for converting bioavailable nitrogen into dinitrogen gas, particularly in carbon-poor environments. However, the diversity and prevalence of anammox bacteria in the terrestrial subsurface-a typically oligotrophic environment-are little understood. To determine the distribution and activity of anammox bacteria across a range of aquifer lithologies and physicochemistries, we analyzed 16S rRNA genes and quantified hydrazine synthase genes and transcripts sampled from 59 groundwater wells and metagenomes and metatranscriptomes from an oxic-to-dysoxic subset. Data indicate that anammox and anammox-associated bacteria (class "Candidatus Brocadiae") are prevalent in the aquifers studied, and that anammox community composition is strongly differentiated by dissolved oxygen (DO), but not ammonia/nitrite. While "Candidatus Brocadiae" diversity decreased with increasing DO, "Candidatus Brocadiae" 16S rRNA genes and hydrazine synthase (hzsB) genes and transcripts were detected across a wide range of bulk groundwater DO concentrations (0 to 10 mg/L). Anammox genes and transcripts correlated significantly with those involved in aerobic ammonia oxidation (amoA), potentially representing a major source of nitrite for anammox. Eight "Candidatus Brocadiae" genomes (63 to 95% complete), representing 2 uncharacterized families and 6 novel species, were reconstructed. Six genomes have genes characteristic of anammox, all for chemolithoautotrophy. Anammox and aerotolerance genes of up to four "Candidatus Brocadiae" genomes were transcriptionally active under oxic and dysoxic conditions, although activity was highest in dysoxic groundwater. The coexpression of nrfAH nitrite reductase genes by "Candidatus Brocadiae" suggests active regeneration of ammonia for anammox. Our findings indicate that anammox bacteria contribute to loss of fixed N across diverse anoxic-to-oxic aquifer conditions, which is likely supported by nitrite from aerobic ammonia oxidation. IMPORTANCE Anammox is increasingly shown to play a major role in the aquatic nitrogen cycle and can outcompete heterotrophic denitrification in environments low in organic carbon. Given that aquifers are characteristically oligotrophic, anammox may represent a major route for the removal of fixed nitrogen in these environments, including agricultural nitrogen, a common groundwater contaminant. Our research confirms that anammox bacteria and the anammox process are prevalent in aquifers and occur across diverse lithologies (e.g., sandy gravel, sand-silt, and volcanic) and groundwater physicochemistries (e.g., various oxygen, carbon, nitrate, and ammonium concentrations). Results reveal niche differentiation among anammox bacteria largely driven by groundwater oxygen contents and provide evidence that anammox is supported by proximity to oxic niches and handoffs from aerobic ammonia oxidizers. We further show that this process, while anaerobic, is active in groundwater characterized as oxic, likely due to the availability of anoxic niches.

Effect of a Wood-Based Carrier of Trichoderma atroviride SC1 on the Microorganisms of the Soil

Wood pellets can sustain the growth of Trichoderma spp. in soil; however, little is known about their side effects on the microbiota. The aims of this study were to evaluate the effect of wood pellets on the growth of Trichoderma spp. in bulk soil and on the soil microbial population's composition and diversity. Trichoderma atroviride SC1 coated wood pellets and non-coated pellets were applied at the level of 10 g∙kg-1 of soil and at the final concentration of 5 × 103 conidia∙g-1 of soil and compared to a conidial suspension applied at the same concentration without the wood carrier. Untreated bulk soil served as a control. The non-coated wood pellets increased the total Trichoderma spp. population throughout the experiment (estimated as colony-forming unit g-1 of soil), while wood pellets coated with T. atroviride SC1 did not. The wood carrier increased the richness, and temporarily decreased the diversity, of the bacterial population, with Massilia being the most abundant bacterial genus, while it decreased both the richness and diversity of the fungal community. Wood pellets selectively increased fungal species having biocontrol potential, such as Mortierella, Cladorrhinum, and Stachybotrys, which confirms the suitability of such carriers of Trichoderma spp. for soil application.

Differences in Gut Microbiome Composition and Antibiotic Resistance Gene Distribution between Chinese and Pakistani University Students from a Common Peer Group

Gut microbiomes play important functional roles in human health and are also affected by many factors. However, few studies concentrate on gut microbiomes under exercise intervention. Additionally, antibiotic resistance genes (ARGs) carried by gut microbiomes may constantly pose a threat to human health. Here, ARGs and microbiomes of Chinese and Pakistanis participants were investigated using 16S rRNA gene sequencing and high-throughput quantitative PCR techniques. The exercise had no impact on gut microbiomes in the 12 individuals investigated during the observation period, while the different distribution of gut microbiomes was found in distinct nationalities. Overall, the dominant microbial phyla in the participants’ gut were Bacteroidota, Firmicutes and Proteobacteria. Some genera such as Prevotella and Dialister were more abundant in Pakistani participants and some other genera such as Bacteroides and Faecalibacterium were more abundant in Chinese participants. The microbial diversity in Chinese was higher than that in Pakistanis. Furthermore, microbial community structures were also different between Chinese and Pakistanis. For ARGs, the distribution of all detected ARGs is not distinct at each time point. Among these ARGs, floR was distributed differently in Chinese and Pakistani participants, and some ARGs such as tetQ and sul2 are positively correlated with several dominant microbiomes, particularly Bacteroidota and Firmicutes bacteria that did not fluctuate over time.

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan

Thermoflexus hugenholtzii JAD2^T, the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2^T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2^T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing “Candidatus Thermoflexus japonica,” “Candidatus Thermoflexus tengchongensis,” and “Candidatus Thermoflexus sinensis.” Genomics was integrated with targeted exometabolomics and ¹³C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2^T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2^T using ¹³C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific ¹³C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

The effect of antibiotics on the gut microbiome: a metagenomics analysis of microbial shift and gut antibiotic resistance in antibiotic treated mice

Background

Emergence of antibiotic resistance is a global public health concern. The relationships between antibiotic use, the gut community composition, normal physiology and metabolism, and individual and public health are still being defined. Shifts in composition of bacteria, antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) after antibiotic treatment are not well-understood.

Methods

This project used next-generation sequencing, custom-built metagenomics pipeline and differential abundance analysis to study the effect of antibiotic monotherapy on resistome and taxonomic composition in the gut of Balb/c mice infected with E. coli via transurethral catheterization to investigate the evolution and emergence of antibiotic resistance.

Results

There is a longitudinal decrease of gut microbiota diversity after antibiotic treatment. Various ARGs are enriched within the gut microbiota despite an overall reduction of the diversity and total amount of bacteria after antibiotic treatment. Sometimes treatment with a specific class of antibiotics selected for ARGs that resist antibiotics of a completely different class (e.g. treatment of ciprofloxacin or fosfomycin selected for cepA that resists ampicillin). Relative abundance of some MGEs increased substantially after antibiotic treatment (e.g. transposases in the ciprofloxacin group).

Conclusions

Antibiotic treatment caused a remarkable reduction in diversity of gut bacterial microbiota but enrichment of certain types of ARGs and MGEs. These results demonstrate an emergence of cross-resistance as well as a profound change in the gut resistome following oral treatment of antibiotics.

Mycobiome diversity: high-throughput sequencing and identification of fungi

Fungi are major ecological players in both terrestrial and aquatic environments by cycling organic matter and channelling nutrients across trophic levels. High-throughput sequencing (HTS) studies of fungal communities are redrawing the map of the fungal kingdom by hinting at its enormous - and largely uncharted - taxonomic and functional diversity. However, HTS approaches come with a range of pitfalls and potential biases, cautioning against unwary application and interpretation of HTS technologies and results. In this Review, we provide an overview and practical recommendations for aspects of HTS studies ranging from sampling and laboratory practices to data processing and analysis. We also discuss upcoming trends and techniques in the field and summarize recent and noteworthy results from HTS studies targeting fungal communities and guilds. Our Review highlights the need for reproducibility and public data availability in the study of fungal communities. If the associated challenges and conceptual barriers are overcome, HTS offers immense possibilities in mycology and elsewhere.

Metaxa2 Database Builder: enabling taxonomic identification from metagenomic or metabarcoding data using any genetic marker

Motivation

Correct taxonomic identification of DNA sequences is central to studies of biodiversity using both shotgun metagenomic and metabarcoding approaches. However, no genetic marker gives sufficient performance across all the biological kingdoms, hampering studies of taxonomic diversity in many groups of organisms. This has led to the adoption of a range of genetic markers for DNA metabarcoding. While many taxonomic classification software tools can be re-trained on these genetic markers, they are often designed with assumptions that impair their utility on genes other than the SSU and LSU rRNA. Here, we present an update to Metaxa2 that enables the use of any genetic marker for taxonomic classification of metagenome and amplicon sequence data.

Results

We evaluated the Metaxa2 Database Builder on 11 commonly used barcoding regions and found that while there are wide differences in performance between different genetic markers, our software performs satisfactorily provided that the input taxonomy and sequence data are of high quality.

Availability and implementation

Freely available on the web as part of the Metaxa2 package at http://microbiology.se/software/metaxa2/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Industrial wastewater treatment plant enriches antibiotic resistance genes and alters the structure of microbial communities

Antibiotic resistance is an emerging global health crisis, driven largely by overuse and misuse of antibiotics. However, there are examples in which the production of these antimicrobial agents has polluted the environment with active antibiotic residues, selecting for antibiotic resistant bacteria and the genes they carry. In this work, we have used shotgun metagenomics to investigate the taxonomic structure and resistance gene composition of sludge communities in a treatment plant in Croatia receiving wastewater from production of the macrolide antibiotic azithromycin. We found that the total abundance of antibiotic resistance genes was three times higher in sludge from the treatment plant receiving wastewater from pharmaceutical production than in municipal sludge from a sewage treatment plant in Zagreb. Surprisingly, macrolide resistance genes did not have higher abundances in the industrial sludge, but genes associated with mobile genetic elements such as integrons had. We conclude that at high concentrations of antibiotics, selection may favor taxonomic shifts towards intrinsically resistant species or strains harboring chromosomal resistance mutations rather than acquisition of mobile resistance determinants. Our results underscore the need for regulatory action also within Europe to avoid release of antibiotics into the environment.

The dynamics of the antibiotic resistome in the feces of freshly weaned pigs following therapeutic administration of oxytetracycline

In this study, shotgun metagenomics was employed to monitor the effect of oxytetracycline, administered at a therapeutic dose, on the dynamics of the microbiota and resistome in the feces of weaned pigs. Sixteen weaning pigs were assigned to one of two treatments including standard starter diet for 21 days or antibiotic-supplemented diet (10 g oxytetracycline/100 kg body weight/day) for 7 days, followed by 14 days of standard starter diet. Feces were collected from the pigs on days 0, 8, and 21 for microbiota and resistome profiling. Pigs receiving oxytetracycline exhibited a significantly greater richness (ANOVA, P = 0.034) and diversity (ANOVA, P = 0.048) of antibiotic resistance genes (ARGs) than the control pigs. Antibiotic administration significantly enriched the abundances of 41 ARGs, mainly from the tetracycline, betalactam and multidrug resistance classes. Compositional shifts in the bacterial communities were observed following 7 days of antibiotic adminstration, with the medicated pigs showing an increase in Escherichia (Proteobacteria) and Prevotella (Bacteroidetes) populations compared with the nonmedicated pigs. This might be explained by the potential of these taxa to carry ARGs that may be transferred to other susceptible bacteria in the densely populated gut environment. These findings will help in the optimization of therapeutic schemes involving antibiotic usage in swine production.

The diversity of uncharacterized antibiotic resistance genes can be predicted from known gene variants-but not always

BACKGROUND

Antibiotic resistance is considered one of the most urgent threats to modern healthcare, and the role of the environment in resistance development is increasingly recognized. It is often assumed that the abundance and diversity of known resistance genes are representative also for the non-characterized fraction of the resistome in a given environment, but this assumption has not been verified. In this study, this hypothesis is tested, using the resistance gene profiles of 1109 metagenomes from various environments.

RESULTS

This study shows that the diversity and abundance of known antibiotic resistance genes can generally predict the diversity and abundance of undescribed resistance genes. However, the extent of this predictability is dependent on the type of environment investigated. Furthermore, it is shown that carefully selected small sets of resistance genes can describe total resistance gene diversity remarkably well.

CONCLUSIONS

The results of this study suggest that knowledge gained from large-scale quantifications of known resistance genes can be utilized as a proxy for unknown resistance factors. This is important for current and proposed monitoring efforts for environmental antibiotic resistance and has implications for the design of risk ranking strategies and the choices of measures and methods for describing resistance gene abundance and diversity in the environment.

Occurrence of Traditional and Alternative Fecal Indicators in Tropical Urban Environments under Different Land Use Patterns

This study evaluated the geospatial distribution of fecal indicator bacteria (FIB) (i.e., Escherichia coli, Enterococcus spp.) and the alternative fecal indicator pepper mild mottle virus (PMMoV) in tropical freshwater environments under different land use patterns. Results show that the occurrence and concentration of microbial fecal indicators were higher for urban than for parkland-dominated areas, consistent with land use weightage. Significant positive correlations with traditional FIB indicate that PMMoV is a suitable indicator of fecal contamination in tropical catchments waters (0.549 ≤ rho ≤ 0.612; P < 0.01). PMMoV exhibited a strong significant correlation with land use weightage (rho = 0.728; P < 0.01) compared to traditional FIB (rho = 0.583; P < 0.01). In addition, chemical tracers were also added to evaluate the potential relationships with microbial fecal indicators. The relationships between diverse variables (e.g., environmental parameters, land use coverage, and chemical tracers) and the occurrence of FIB and PMMoV were evaluated. By using stepwise multiple linear regression (MLR), the empirical experimental models substantiate the impact of land use patterns and anthropogenic activities on microbial water quality, and the output results of the empirical models may be able to predict the sources and transportation of human fecal pollution or sewage contamination. In addition, the high correlation between PMMoV data obtained from quantitative real-time PCR (qPCR) and viral metagenomics data supports the possibility of using viral metagenomics to relatively quantify specific microbial indicators for monitoring microbial water quality (0.588 ≤ rho ≤ 0.879; P < 0.05).IMPORTANCE The results of this study may support the hypothesis of using PMMoV as an alternative indicator of human fecal contamination in tropical surface waters from the perspective of land use patterns. The predictive result of the occurrence of human fecal indicators with high accuracy may reflect the source and transportation of human fecal pollution, which are directly related to the risk to human health, and thereafter, steps can be taken to mitigate these risks.

Using metagenomics to investigate human and environmental resistomes

Antibiotic resistance is a global health concern declared by the WHO as one of the largest threats to modern healthcare. In recent years, metagenomic DNA sequencing has started to be applied as a tool to study antibiotic resistance in different environments, including the human microbiota. However, a multitude of methods exist for metagenomic data analysis, and not all methods are suitable for the investigation of resistance genes, particularly if the desired outcome is an assessment of risks to human health. In this review, we outline the current state of methods for sequence handling, mapping to databases of resistance genes, statistical analysis and metagenomic assembly. In addition, we provide an overview of important considerations related to the analysis of resistance genes, and recommend some of the currently used tools and methods that are best equipped to inform research and clinical practice related to antibiotic resistance.

The relation between Blastocystis and the intestinal microbiota in Swedish travellers

Background

Blastocystis sp. is a unicellular eukaryote that is commonly found in the human intestine. Its ability to cause disease is debated and a subject for ongoing research. In this study, faecal samples from 35 Swedish university students were examined through shotgun metagenomics before and after travel to the Indian peninsula or Central Africa. We aimed at assessing the impact of travel on Blastocystis carriage and seek associations between Blastocystis and the bacterial microbiota.

Results

We found a prevalence of Blastocystis of 16/35 (46%) before travel and 15/35 (43%) after travel. The two most commonly Blastocystis subtypes (STs) found were ST3 and ST4, accounting for 20 of the 31 samples positive for Blastocystis. No mixed subtype carriage was detected. All ten individuals with a typable ST before and after travel maintained their initial ST. The composition of the gut bacterial community was not significantly different between Blastocystis-carriers and non-carriers. Interestingly, the presence of Blastocystis was accompanied with higher abundances of the bacterial genera Sporolactobacillus and Candidatus Carsonella. Blastocystis carriage was positively associated with high bacterial genus richness, and negatively correlated to the Bacteroides-driven enterotype. These associations were both largely dependent on ST4 – a subtype commonly described from Europe – while the globally prevalent ST3 did not show such significant relationships.

Conclusions

The high rate of Blastocystis subtype persistence found during travel indicates that long-term carriage of Blastocystis is common. The associations between Blastocystis and the bacterial microbiota found in this study could imply a link between Blastocystis and a healthy microbiota as well as with diets high in vegetables. Whether the associations between Blastocystis and the microbiota are resulting from the presence of Blastocystis, or are a prerequisite for colonization with Blastocystis, are interesting questions for further studies.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1139-7) contains supplementary material, which is available to authorized users.

In Vivo Analysis of the Viable Microbiota and Helicobacter pylori Transcriptome in Gastric Infection and Early Stages of Carcinogenesis

Emerging evidence shows that the human microbiota plays a larger role in disease progression and health than previously anticipated. Helicobacter pylori, the causative agent of gastric cancer and duodenal and gastric ulcers, was early associated with gastric disease, but it has also been proposed that the accompanying microbiota in Helicobacter pylori-infected individuals might affect disease progression and gastric cancer development. In this study, the composition of the transcriptionally active microbial community and H. pylori gene expression were determined using metatranscriptomic RNA sequencing of stomach biopsy specimens from individuals with different H. pylori infection statuses and premalignant tissue changes. The results show that H. pylori completely dominates the microbiota not only in infected individuals but also in most individuals classified as H. pylori uninfected using conventional methods. Furthermore, H. pylori abundance is positively correlated with the presence of Campylobacter, Deinococcus, and Sulfurospirillum Finally, we quantified the expression of a large number of Helicobacter pylori genes and found high expression of genes involved in pH regulation and nickel transport. Our study is the first to dissect the viable microbiota of the human stomach by metatranscriptomic analysis, and it shows that metatranscriptomic analysis of the gastric microbiota is feasible and can provide new insights into how bacteria respond in vivo to variations in the stomach microenvironment and at different stages of disease progression.

Elucidating selection processes for antibiotic resistance in sewage treatment plants using metagenomics

Sewage treatment plants (STPs) have repeatedly been suggested as "hotspots" for the emergence and dissemination of antibiotic-resistant bacteria. A critical question still unanswered is if selection pressures within STPs, caused by residual antibiotics or other co-selective agents, are sufficient to specifically promote resistance. To address this, we employed shotgun metagenomic sequencing of samples from different steps of the treatment process in three Swedish STPs. In parallel, concentrations of selected antibiotics, biocides and metals were analyzed. We found that concentrations of tetracycline and ciprofloxacin in the influent were above predicted concentrations for resistance selection, however, there was no consistent enrichment of resistance genes to any particular class of antibiotics in the STPs, neither for biocide and metal resistance genes. The most substantial change of the bacterial communities compared to human feces occurred already in the sewage pipes, manifested by a strong shift from obligate to facultative anaerobes. Through the treatment process, resistance genes against antibiotics, biocides and metals were not reduced to the same extent as fecal bacteria. The OXA-48 gene was consistently enriched in surplus and digested sludge. We find this worrying as OXA-48, still rare in Swedish clinical isolates, provides resistance to carbapenems, one of our most critically important classes of antibiotics. Taken together, metagenomics analyses did not provide clear support for specific antibiotic resistance selection. However, stronger selective forces affecting gross taxonomic composition, and with that resistance gene abundances, limit interpretability. Comprehensive analyses of resistant/non-resistant strains within relevant species are therefore warranted.