Interpretable metric learning in comparative metagenomics: The adaptive Haar-like distance

Random forests have emerged as a promising tool in comparative metagenomics because they can predict environmental characteristics based on microbial composition in datasets where β-diversity metrics fall short of revealing meaningful relationships between samples. Nevertheless, despite this efficacy, they lack biological insight in tandem with their predictions, potentially hindering scientific advancement. To overcome this limitation, we leverage a geometric characterization of random forests to introduce a data-driven phylogenetic β-diversity metric, the adaptive Haar-like distance. This new metric assigns a weight to each internal node (i.e., split or bifurcation) of a reference phylogeny, indicating the relative importance of that node in discerning environmental samples based on their microbial composition. Alongside this, a weighted nearest-neighbors classifier, constructed using the adaptive metric, can be used as a proxy for the random forest while maintaining accuracy on par with that of the original forest and another state-of-the-art classifier, CoDaCoRe. As shown in datasets from diverse microbial environments, however, the new metric and classifier significantly enhance the biological interpretability and visualization of high-dimensional metagenomic samples.

Microdiverse bacterial clades prevail across Antarctic wetlands

Antarctica's extreme environmental conditions impose selection pressures on microbial communities. Indeed, a previous study revealed that bacterial assemblages at the Cierva Point Wetland Complex (CPWC) are shaped by strong homogeneous selection. Yet which bacterial phylogenetic clades are shaped by selection processes and their ecological strategies to thrive in such extreme conditions remain unknown. Here, we applied the phyloscore and feature-level βNTI indexes coupled with phylofactorization to successfully detect bacterial monophyletic clades subjected to homogeneous (HoS) and heterogenous (HeS) selection. Remarkably, only the HoS clades showed high relative abundance across all samples and signs of putative microdiversity. The majority of the amplicon sequence variants (ASVs) within each HoS clade clustered into a unique 97% sequence similarity operational taxonomic unit (OTU) and inhabited a specific environment (lotic, lentic or terrestrial). Our findings suggest the existence of microdiversification leading to sub-taxa niche differentiation, with putative distinct ecotypes (consisting of groups of ASVs) adapted to a specific environment. We hypothesize that HoS clades thriving in the CPWC have phylogenetically conserved traits that accelerate their rate of evolution, enabling them to adapt to strong spatio-temporally variable selection pressures. Variable selection appears to operate within clades to cause very rapid microdiversification without losing key traits that lead to high abundance. Variable and homogeneous selection, therefore, operate simultaneously but on different aspects of organismal ecology. The result is an overall signal of homogeneous selection due to rapid within-clade microdiversification caused by variable selection. It is unknown whether other systems experience this dynamic, and we encourage future work evaluating the transferability of our results.

Viral genomic features predict orthopoxvirus reservoir hosts

Orthopoxviruses (OPVs), including the causative agents of smallpox and mpox have led to devastating outbreaks in human populations worldwide. However, the discontinuation of smallpox vaccination, which also provides cross-protection against related OPVs, has diminished global immunity to OPVs more broadly. We apply machine learning models incorporating both host ecological and viral genomic features to predict likely reservoirs of OPVs. We demonstrate that incorporating viral genomic features in addition to host ecological traits enhanced the accuracy of potential OPV host predictions, highlighting the importance of host-virus molecular interactions in predicting potential host species. We identify hotspots for geographic regions rich with potential OPV hosts in parts of southeast Asia, equatorial Africa, and the Amazon, revealing high overlap between regions predicted to have a high number of potential OPV host species and those with the lowest smallpox vaccination coverage, indicating a heightened risk for the emergence or establishment of zoonotic OPVs. Our findings can be used to target wildlife surveillance, particularly related to concerns about mpox establishment beyond its historical range.

Gill-associated bacteria are homogeneously selected in amphibious mangrove crabs to sustain host intertidal adaptation

BACKGROUND

The transition from water to air is a key event in the evolution of many marine organisms to access new food sources, escape water hypoxia, and exploit the higher and temperature-independent oxygen concentration of air. Despite the importance of microorganisms in host adaptation, their contribution to overcoming the challenges posed by the lifestyle changes from water to land is not well understood. To address this, we examined how microbial association with a key multifunctional organ, the gill, is involved in the intertidal adaptation of fiddler crabs, a dual-breathing organism.

RESULTS

Electron microscopy revealed a rod-shaped bacterial layer tightly connected to the gill lamellae of the five crab species sampled across a latitudinal gradient from the central Red Sea to the southern Indian Ocean. The gill bacterial community diversity assessed with 16S rRNA gene amplicon sequencing was consistently low across crab species, and the same actinobacterial group, namely Ilumatobacter, was dominant regardless of the geographic location of the host. Using metagenomics and metatranscriptomics, we detected that these members of actinobacteria are potentially able to convert ammonia to amino acids and may help eliminate toxic sulphur compounds and carbon monoxide to which crabs are constantly exposed.

CONCLUSIONS

These results indicate that bacteria selected on gills can play a role in the adaptation of animals in dynamic intertidal ecosystems. Hence, this relationship is likely to be important in the ecological and evolutionary processes of the transition from water to air and deserves further attention, including the ontogenetic onset of this association. Video Abstract.

Network embedding unveils the hidden interactions in the mammalian virome

Predicting host-virus interactions is fundamentally a network science problem. We develop a method for bipartite network prediction that combines a recommender system (linear filtering) with an imputation algorithm based on low-rank graph embedding. We test this method by applying it to a global database of mammal-virus interactions and thus show that it makes biologically plausible predictions that are robust to data biases. We find that the mammalian virome is under-characterized anywhere in the world. We suggest that future virus discovery efforts could prioritize the Amazon Basin (for its unique coevolutionary assemblages) and sub-Saharan Africa (for its poorly characterized zoonotic reservoirs). Graph embedding of the imputed network improves predictions of human infection from viral genome features, providing a shortlist of priorities for laboratory studies and surveillance. Overall, our study indicates that the global structure of the mammal-virus network contains a large amount of information that is recoverable, and this provides new insights into fundamental biology and disease emergence.

Coronavirus sampling and surveillance in bats from 1996-2019: a systematic review and meta-analysis

The emergence of SARS-CoV-2 highlights a need for evidence-based strategies to monitor bat viruses. We performed a systematic review of coronavirus sampling (testing for RNA positivity) in bats globally. We identified 110 studies published between 2005 and 2020 that collectively reported positivity from 89,752 bat samples. We compiled 2,274 records of infection prevalence at the finest methodological, spatiotemporal and phylogenetic level of detail possible from public records into an open, static database named datacov, together with metadata on sampling and diagnostic methods. We found substantial heterogeneity in viral prevalence across studies, reflecting spatiotemporal variation in viral dynamics and methodological differences. Meta-analysis identified sample type and sampling design as the best predictors of prevalence, with virus detection maximized in rectal and faecal samples and by repeat sampling of the same site. Fewer than one in five studies collected and reported longitudinal data, and euthanasia did not improve virus detection. We show that bat sampling before the SARS-CoV-2 pandemic was concentrated in China, with research gaps in South Asia, the Americas and sub-Saharan Africa, and in subfamilies of phyllostomid bats. We propose that surveillance strategies should address these gaps to improve global health security and enable the origins of zoonotic coronaviruses to be identified.

Homogeneous Environmental Selection Structures the Bacterial Communities of Benthic Biofilms in Proglacial Floodplain Streams

In proglacial floodplains, glacier recession promotes biogeochemical and ecological gradients across relatively small spatial scales. The resulting environmental heterogeneity induces remarkable microbial biodiversity among proglacial stream biofilms. Yet the relative importance of environmental constraints in forming biofilm communities remains largely unknown. Extreme environmental conditions in proglacial streams may lead to the homogenizing selection of biofilm-forming microorganisms. However, environmental differences between proglacial streams may impose different selective forces, resulting in nested, spatially structured assembly processes. Here, we investigated bacterial community assembly processes by unraveling ecologically successful phylogenetic clades in two stream types (glacier-fed mainstems and non-glacier-fed tributaries) draining three proglacial floodplains in the Swiss Alps. Clades with low phylogenetic turnover rates were present in all stream types, including Gammaproteobacteria and Alphaproteobacteria, while the other clades were specific to one stream type. These clades constituted up to 34.8% and 31.1% of the community diversity and up to 61.3% and 50.9% of the relative abundances in mainstems and tributaries, respectively, highlighting their importance and success in these communities. Furthermore, the proportion of bacteria under homogeneous selection was inversely related to the abundance of photoautotrophs, and these clades may therefore decrease in abundance with the future "greening" of proglacial habitats. Finally, we found little effect of physical distance from the glacier on clades under selection in glacier-fed streams, probably due to the high hydrological connectivity of our study reaches. Overall, these findings shed new light on the mechanisms of microbial biofilm assembly in proglacial streams and help us to predict their future in a rapidly changing environment. IMPORTANCE Streams draining proglacial floodplains harbor benthic biofilms comprised of diverse microbial communities. These high-mountain ecosystems are rapidly changing with climate warming, and it is therefore critical to better understand the mechanisms underlying the assembly of their microbial communities. We found that homogeneous selection dominates the structuring of bacterial communities in benthic biofilms in both glacier-fed mainstems and nonglacier tributary streams within three proglacial floodplains in the Swiss Alps. However, differences between glacier-fed and tributary ecosystems may impose differential selective forces. Here, we uncovered nested, spatially structured assembly processes for proglacial floodplain communities. Our analyses additionally provided insights into linkages between aquatic photoautotrophs and the bacterial taxa under homogeneous selection, potentially by providing a labile source of carbon in these otherwise carbon-deprived systems. In the future, we expect a shift in the bacterial communities under homogeneous selection in glacier-fed streams as primary production becomes more important and streams become "greener".

Tropical bat ectoparasitism in continuous versus fragmented forests: A gap analysis and preliminary meta-analysis

Tropical regions are experiencing rapid rates of forest fragmentation, which can have several effects on wildlife, including altered parasite dynamics. Bats are a useful host group to consider the effects of fragmentation, because they are abundant in the tropics, serve important ecological roles, and harbor many parasites. Nevertheless, research on the effects of fragmentation on bat ectoparasites is still limited. To help guide ongoing and future research efforts, this study had two objectives: (1) conduct a gap analysis to characterize the state of currently available research on fragmentation effects on bat ectoparasites and (2) conduct a preliminary meta-analysis to identify current trends. We systematically highlighted several research gaps: Studies comparing the effects of fragmented versus continuous forests on ectoparasites are limited and have primarily been conducted in the Neotropics, with a focus on bats in the superfamily Noctilionidea (especially frugivorous phyllostomids). Our preliminary meta-analysis suggested that ectoparasite prevalence (but not the mean or variance in intensity) was higher in fragments than in continuous forests. Moreover, prevalence increased with increasing roost duration, and mean intensity was higher for bats with higher wing aspect ratios. Intensity variance was affected by an interaction between forest type and wing aspect ratio, such that variance increased for bats with high-wing aspect ratios in continuous forests but decreased in fragments. These results suggest that fragmentation can shape aspects of bat ectoparasitism and could have implications for the ecology, health, and conservation of bats in fragmented landscapes. However, existing research gaps could bias our current understanding of habitat change and bat health, and future research should thus investigate these effects in the Paleotropics and with other bat families.

The role of microbial ecology in improving the performance of anaerobic digestion of sewage sludge

The use of next-generation diagnostic tools to optimise the anaerobic digestion of municipal sewage sludge has the potential to increase renewable natural gas recovery, improve the reuse of biosolid fertilisers and help operators expand circular economies globally. This review aims to provide perspectives on the role of microbial ecology in improving digester performance in wastewater treatment plants, highlighting that a systems biology approach is fundamental for monitoring mesophilic anaerobic sewage sludge in continuously stirred reactor tanks. We further highlight the potential applications arising from investigations into sludge ecology. The principal limitation for improvements in methane recoveries or in process stability of anaerobic digestion, especially after pre-treatment or during co-digestion, are ecological knowledge gaps related to the front-end metabolism (hydrolysis and fermentation). Operational problems such as stable biological foaming are a key problem, for which ecological markers are a suitable approach. However, no biomarkers exist yet to assist in monitoring and management of clade-specific foaming potentials along with other risks, such as pollutants and pathogens. Fundamental ecological principles apply to anaerobic digestion, which presents opportunities to predict and manipulate reactor functions. The path ahead for mapping ecological markers on process endpoints and risk factors of anaerobic digestion will involve numerical ecology, an expanding field that employs metrics derived from alpha, beta, phylogenetic, taxonomic, and functional diversity, as well as from phenotypes or life strategies derived from genetic potentials. In contrast to addressing operational issues (as noted above), which are effectively addressed by whole population or individual biomarkers, broad improvement and optimisation of function will require enhancement of hydrolysis and acidogenic processes. This will require a discovery-based approach, which will involve integrative research involving the proteome and metabolome. This will utilise, but overcome current limitations of DNA-centric approaches, and likely have broad application outside the specific field of anaerobic digestion.

Integrating phylogenetic and functional data in microbiome studies

MOTIVATION

Microbiome functional data are frequently analyzed to identify associations between microbial functions (e.g. genes) and sample groups of interest. However, it is challenging to distinguish between different possible explanations for variation in community-wide functional profiles by considering functions alone. To help address this problem, we have developed POMS, a package that implements multiple phylogeny-aware frameworks to more robustly identify enriched functions.

RESULTS

The key contribution is an extended balance-tree workflow that incorporates functional and taxonomic information to identify functions that are consistently enriched in sample groups across independent taxonomic lineages. Our package also includes a workflow for running phylogenetic regression. Based on simulated data we demonstrate that these approaches more accurately identify gene families that confer a selective advantage compared with commonly used tools. We also show that POMS in particular can identify enriched functions in real-world metagenomics datasets that are potential targets of strong selection on multiple members of the microbiome.

AVAILABILITY AND IMPLEMENTATION

These workflows are freely available in the POMS R package at https://github.com/gavinmdouglas/POMS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Optimising predictive models to prioritise viral discovery in zoonotic reservoirs

Despite the global investment in One Health disease surveillance, it remains difficult and costly to identify and monitor the wildlife reservoirs of novel zoonotic viruses. Statistical models can guide sampling target prioritisation, but the predictions from any given model might be highly uncertain; moreover, systematic model validation is rare, and the drivers of model performance are consequently under-documented. Here, we use the bat hosts of betacoronaviruses as a case study for the data-driven process of comparing and validating predictive models of probable reservoir hosts. In early 2020, we generated an ensemble of eight statistical models that predicted host-virus associations and developed priority sampling recommendations for potential bat reservoirs of betacoronaviruses and bridge hosts for SARS-CoV-2. During a time frame of more than a year, we tracked the discovery of 47 new bat hosts of betacoronaviruses, validated the initial predictions, and dynamically updated our analytical pipeline. We found that ecological trait-based models performed well at predicting these novel hosts, whereas network methods consistently performed approximately as well or worse than expected at random. These findings illustrate the importance of ensemble modelling as a buffer against mixed-model quality and highlight the value of including host ecology in predictive models. Our revised models showed an improved performance compared with the initial ensemble, and predicted more than 400 bat species globally that could be undetected betacoronavirus hosts. We show, through systematic validation, that machine learning models can help to optimise wildlife sampling for undiscovered viruses and illustrates how such approaches are best implemented through a dynamic process of prediction, data collection, validation, and updating.

Virus isolation data improve host predictions for New World rodent orthohantaviruses

Identifying reservoir host species is crucial for understanding the ecology of multi-host pathogens and predicting risks of pathogen spillover from wildlife to people. Predictive models are increasingly used for identifying ecological traits and prioritizing surveillance of likely zoonotic reservoirs, but these often employ different types of evidence for establishing host associations. Comparisons between models with different infection evidence are necessary to guide inferences about the trait profiles of likely hosts and identify which hosts and geographical regions are likely sources of spillover. Here, we use New World rodent-orthohantavirus associations to explore differences in the performance and predictions of models trained on two types of evidence for infection and onward transmission: RT-PCR and live virus isolation data, representing active infections versus host competence, respectively. Orthohantaviruses are primarily carried by muroid rodents and cause the diseases haemorrhagic fever with renal syndrome (HFRS) and hantavirus cardiopulmonary syndrome (HCPS) in humans. We show that although boosted regression tree (BRT) models trained on RT-PCR and live virus isolation data both performed well and capture generally similar trait profiles, rodent phylogeny influenced previously collected RT-PCR data, and BRTs using virus isolation data displayed a narrower list of predicted reservoirs than those using RT-PCR data. BRT models trained on RT-PCR data identified 138 undiscovered hosts and virus isolation models identified 92 undiscovered hosts, with 27 undiscovered hosts identified by both models. Distributions of predicted hosts were concentrated in several different regions for each model, with large discrepancies between evidence types. As a form of validation, virus isolation models independently predicted several orthohantavirus-rodent host associations that had been previously identified through empirical research using RT-PCR. Our model predictions provide a priority list of species and locations for future orthohantavirus sampling. More broadly, these results demonstrate the value of multiple data types for predicting zoonotic pathogen hosts. These methods can be applied across a range of systems to improve our understanding of pathogen maintenance and increase efficiency of pathogen surveillance.

Glacier shrinkage will accelerate downstream decomposition of organic matter and alters microbiome structure and function

The shrinking of glaciers is among the most iconic consequences of climate change. Despite this, the downstream consequences for ecosystem processes and related microbiome structure and function remain poorly understood. Here, using a space-for-time substitution approach across 101 glacier-fed streams (GFSs) from six major regions worldwide, we investigated how glacier shrinkage is likely to impact the organic matter (OM) decomposition rates of benthic biofilms. To do this, we measured the activities of five common extracellular enzymes and estimated decomposition rates by using enzyme allocation equations based on stoichiometry. We found decomposition rates to average 0.0129 (% d^-1 ), and that decreases in glacier influence (estimated by percent glacier catchment coverage, turbidity, and a glacier index) accelerates decomposition rates. To explore mechanisms behind these relationships, we further compared decomposition rates with biofilm and stream water characteristics. We found that chlorophyll-a, temperature, and stream water N:P together explained 61% of the variability in decomposition. Algal biomass, which is also increasing with glacier shrinkage, showed a particularly strong relationship with decomposition, likely indicating their importance in contributing labile organic compounds to these carbon-poor habitats. We also found high relative abundances of chytrid fungi in GFS sediments, which putatively parasitize these algae, promoting decomposition through a fungal shunt. Exploring the biofilm microbiome, we then sought to identify bacterial phylogenetic clades significantly associated with decomposition, and found numerous positively (e.g., Saprospiraceae) and negatively (e.g., Nitrospira) related clades. Lastly, using metagenomics, we found evidence of different bacterial classes possessing different proportions of EEA-encoding genes, potentially informing some of the microbial associations with decomposition rates. Our results, therefore, present new mechanistic insights into OM decomposition in GFSs by demonstrating that an algal-based "green food web" is likely to increase in importance in the future and will promote important biogeochemical shifts in these streams as glaciers vanish.

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.

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

Axillary Microbiota Is Associated with Cognitive Impairment in Parkinson's Disease Patients

Cognitive impairment (CI) is among the most common non-motor symptoms of Parkinson’s disease (PD), with a substantially negative impact on patient management and outcome. The development and progression of CI exhibits high interindividual variability, which requires better diagnostic and monitoring strategies. PD patients often display sweating disorders resulting from autonomic dysfunction, which has been associated with CI. Because the axillary microbiota is known to change with humidity level and sweat composition, we hypothesized that the axillary microbiota of PD patients shifts in association with CI progression, and thus can be used as a proxy for classification of CI stages in PD. We compared the axillary microbiota compositions of 103 PD patients (55 PD patients with dementia [PDD] and 48 PD patients with mild cognitive impairment [PD-MCI]) and 26 cognitively normal healthy controls (HC). We found that axillary microbiota profiles differentiate HC, PD-MCI, and PDD groups based on differential ranking analysis, and detected an increasing trend in the log ratio of Corynebacterium to Anaerococcus in progression from HC to PDD. In addition, phylogenetic factorization revealed that the depletion of the Anaerococcus, Peptoniphilus, and W5053 genera is associated with PD-MCI and PDD. Moreover, functional predictions suggested significant increases in myo-inositol degradation, ergothioneine biosynthesis, propionate biosynthesis, menaquinone biosynthesis, and the proportion of aerobic bacteria and biofilm formation capacity, in parallel to increasing CI. Our results suggest that alterations in axillary microbiota are associated with CI in PD. Thus, axillary microbiota has the potential to be exploited as a noninvasive tool in the development of novel strategies.

IMPORTANCE Parkinson's disease (PD) is the second most common neurodegenerative disease. Cognitive impairment (CI) in PD has significant negative impacts on life quality of patients. The emergence and progression of cognitive impairment shows high variability among PD patients, and thus requires better diagnostic and monitoring strategies. Recent findings indicate a close link between autonomic dysfunction and cognitive impairment. Since thermoregulatory dysfunction and skin changes are among the main manifestations of autonomic dysfunction in PD, we hypothesized that alterations in the axillary microbiota may be useful for tracking cognitive impairment stages in PD. To our knowledge, this the first study characterizing the axillary microbiota of PD patients and exploring its association with cognitive impairment stages in PD. Future studies should include larger cohorts and multicenter studies to validate our results and investigate potential biological mechanisms.

Elevated atmospheric CO2 alters the microbial community composition and metabolic potential to mineralize organic phosphorus in the rhizosphere of wheat

BACKGROUND

Understanding how elevated atmospheric CO2 (eCO2) impacts on phosphorus (P) transformation in plant rhizosphere is critical for maintaining ecological sustainability in response to climate change, especially in agricultural systems where soil P availability is low.

METHODS

This study used rhizoboxes to physically separate rhizosphere regions (plant root-soil interface) into 1.5-mm segments. Wheat plants were grown in rhizoboxes under eCO2 (800 ppm) and ambient CO2 (400 ppm) in two farming soils, Chromosol and Vertosol, supplemented with phytate (organic P). Photosynthetic carbon flow in the plant-soil continuum was traced with 13CO2 labeling. Amplicon sequencing was performed on the rhizosphere-associated microbial community in the root-growth zone, and 1.5 mm and 3 mm away from the root.

RESULTS

Elevated CO2 accelerated the mineralization of phytate in the rhizosphere zones, which corresponded with increases in plant-derived 13C enrichment and the relative abundances of discreet phylogenetic clades containing Bacteroidetes and Gemmatimonadetes in the bacterial community, and Funneliformis affiliated to arbuscular mycorrhizas in the fungal community. Although the amplicon sequence variants (ASVs) associated the stimulation of phytate mineralization under eCO2 differed between the two soils, these ASVs belonged to the same phyla associated with phytase and phosphatase production. The symbiotic mycorrhizas in the rhizosphere of wheat under eCO2 benefited from increased plant C supply and increased P access from soil. Further supportive evidence was the eCO2-induced increase in the genetic pool expressing the pentose phosphate pathway, which is the central pathway for biosynthesis of RNA/DNA precursors.

CONCLUSIONS

The results suggested that an increased belowground carbon flow under eCO2 stimulated bacterial growth, changing community composition in favor of phylotypes capable of degrading aromatic P compounds. It is proposed that energy investments by bacteria into anabolic processes increase under eCO2 to level microbial P-use efficiencies and that synergies with symbiotic mycorrhizas further enhance the competition for and mineralization of organic P. Video Abstract.

Patterns of Exposure and Infection with Microparasites in Iberian Wild Carnivores: A Review and Meta-Analysis

We use a suite of meta-analytic and comparative methods to derive fundamental insights into how sampling effort, pathogen richness, infection prevalence, and seroprevalence vary across Carnivora taxa and Iberian geography. The red fox was the most studied species, the wolf and Iberian lynx were disproportionally studied, and the Arctoidea were understudied. Sampling effort was higher in Mediterranean areas, but central Spain showed the higher pathogen richness. Excluding studies analyzing fecal samples, 53 different pathogens have been detected in Iberian carnivores, including 16 viruses, 27 bacteria, and 10 protozoa but no fungi. Sampling effort and pathogen diversity were generally more similar among closely related carnivore species. Seropositivity to viruses was lower and higher in the Mustelinae and the Canidae, respectively, and seropositivity to protozoa was higher in both taxa. Canine distemper virus exposure was greatest in canids and mustelids. Carnivore protoparvovirus-1 exposure was greatest in the Atlantic regions, and the Felidae and the Musteloidea had lower infection prevalence. A subclade of the Mustelidae had a greater prevalence of Leishmania infection. We observed no relationships between host phylogenetic distance and pathogen sharing among species. Lastly, we identify important research pitfalls and future directions to improve the study of infectious disease in Iberian wild carnivore communities.

Microdiversity characterizes prevalent phylogenetic clades in the glacier-fed stream microbiome

Glacier-fed streams (GFSs) are extreme and rapidly vanishing ecosystems, and yet they harbor diverse microbial communities. Although our understanding of the GFS microbiome has recently increased, we do not know which microbial clades are ecologically successful in these ecosystems, nor do we understand potentially underlying mechanisms. Ecologically successful clades should be more prevalent across GFSs compared to other clades, which should be reflected as clade-wise distinctly low phylogenetic turnover. However, methods to assess such patterns are currently missing. Here we developed and applied a novel analytical framework, “phyloscore analysis”, to identify clades with lower spatial phylogenetic turnover than other clades in the sediment microbiome across twenty GFSs in New Zealand. These clades constituted up to 44% and 64% of community α-diversity and abundance, respectively. Furthermore, both their α-diversity and abundance increased as sediment chlorophyll a decreased, corroborating their ecological success in GFS habitats largely devoid of primary production. These clades also contained elevated levels of putative microdiversity than others, which could potentially explain their high prevalence in GFSs. This hitherto unknown microdiversity may be threatened as glaciers shrink, urging towards further genomic and functional exploration of the GFS microbiome.

Tree-aggregated predictive modeling of microbiome data

Modern high-throughput sequencing technologies provide low-cost microbiome survey data across all habitats of life at unprecedented scale. At the most granular level, the primary data consist of sparse counts of amplicon sequence variants or operational taxonomic units that are associated with taxonomic and phylogenetic group information. In this contribution, we leverage the hierarchical structure of amplicon data and propose a data-driven and scalable tree-guided aggregation framework to associate microbial subcompositions with response variables of interest. The excess number of zero or low count measurements at the read level forces traditional microbiome data analysis workflows to remove rare sequencing variants or group them by a fixed taxonomic rank, such as genus or phylum, or by phylogenetic similarity. By contrast, our framework, which we call trac (tree-aggregation of compositional data), learns data-adaptive taxon aggregation levels for predictive modeling, greatly reducing the need for user-defined aggregation in preprocessing while simultaneously integrating seamlessly into the compositional data analysis framework. We illustrate the versatility of our framework in the context of large-scale regression problems in human gut, soil, and marine microbial ecosystems. We posit that the inferred aggregation levels provide highly interpretable taxon groupings that can help microbiome researchers gain insights into the structure and functioning of the underlying ecosystem of interest.

treeclimbR pinpoints the data-dependent resolution of hierarchical hypotheses

treeclimbR is for analyzing hierarchical trees of entities, such as phylogenies or cell types, at different resolutions. It proposes multiple candidates that capture the latent signal and pinpoints branches or leaves that contain features of interest, in a data-driven way. It outperforms currently available methods on synthetic data, and we highlight the approach on various applications, including microbiome and microRNA surveys as well as single-cell cytometry and RNA-seq datasets. With the emergence of various multi-resolution genomic datasets, treeclimbR provides a thorough inspection on entities across resolutions and gives additional flexibility to uncover biological associations.

Taxonomic resolution affects host-parasite association model performance

Identifying the factors that structure host–parasite interactions is fundamental to understand the drivers of species distributions and to predict novel cross-species transmission events. More phylogenetically related host species tend to have more similar parasite associations, but parasite specificity may vary as a function of transmission mode, parasite taxonomy or life history. Accordingly, analyses that attempt to infer host−parasite associations using combined data on different parasite groups may perform quite differently relative to analyses on each parasite subset. In essence, are more data always better when predicting host−parasite associations, or does parasite taxonomic resolution matter? Here, we explore how taxonomic resolution affects predictive models of host−parasite associations using the London Natural History Museum's database of host–helminth interactions. Using boosted regression trees, we demonstrate that taxon-specific models (i.e. of Acanthocephalans, Nematodes and Platyhelminthes) consistently outperform full models in predicting mammal-helminth associations. At finer spatial resolutions, full and taxon-specific model performance does not vary, suggesting tradeoffs between phylogenetic and spatial scales of analysis. Although all models identify similar host and parasite covariates as important to such patterns, our results emphasize the importance of phylogenetic scale in the study of host–parasite interactions and suggest that using taxonomic subsets of data may improve predictions of parasite distributions and cross-species transmission. Predictive models of host–pathogen interactions should thus attempt to encompass the spatial resolution and phylogenetic scale desired for inference and prediction and potentially use model averaging or ensemble models to combine predictions from separately trained models.

Disentangling interactions among mercury, immunity and infection in a Neotropical bat community

1. Contaminants such as mercury are pervasive and can have immunosuppressive effects on wildlife. Impaired immunity could be important for forecasting pathogen spillover, as many land-use changes that generate mercury contamination also bring wildlife into close contact with humans and domestic animals. However, the interactions among contaminants, immunity and infection are difficult to study in natural systems, and empirical tests of possible directional relationships remain rare. 2. We capitalized on extreme mercury variation in a diverse bat community in Belize to test association among contaminants, immunity and infection. By comparing a previous dataset of bats sampled in 2014 with new data from 2017, representing a period of rapid agricultural land conversion, we first confirmed bat species more reliant on aquatic prey had higher fur mercury. Bats in the agricultural habitat also had higher mercury in recent years. We then tested covariation between mercury and cellular immunity and determined if such relationships mediated associations between mercury and bacterial pathogens. As bat ecology can dictate exposure to mercury and pathogens, we also assessed species-specific patterns in mercury-infection relationships. 3. Across the bat community, individuals with higher mercury had fewer neutrophils but not lymphocytes, suggesting stronger associations with innate immunity. However, the odds of infection for haemoplasmas and Bartonella spp. were generally lowest in bats with high mercury, and relationships between mercury and immunity did not mediate infection patterns. Mercury also showed species- and clade-specific relationships with infection, being associated with especially low odds for haemoplasmas in Pteronotus mesoamericanus and Dermanura phaeotis. For Bartonella spp., mercury was associated with particularly low odds of infection in the genus Pteronotus but high odds in the subfamily Stenodermatinae. 4. Synthesis and application. Lower general infection risk in bats with high mercury despite weaker innate defense suggests contaminant-driven loss of pathogen habitat (i.e. anemia) or vector mortality as possible causes. Greater attention to these potential pathways could help disentangle relationships among contaminants, immunity and infection in anthropogenic habitats and help forecast disease risks. Our results also suggest that contaminants may increase infection risk in some taxa but not others, emphasizing the importance of considering surveillance and management at different phylogenetic scales.

Statistical and Machine Learning Techniques in Human Microbiome Studies: Contemporary Challenges and Solutions

The human microbiome has emerged as a central research topic in human biology and biomedicine. Current microbiome studies generate high-throughput omics data across different body sites, populations, and life stages. Many of the challenges in microbiome research are similar to other high-throughput studies, the quantitative analyses need to address the heterogeneity of data, specific statistical properties, and the remarkable variation in microbiome composition across individuals and body sites. This has led to a broad spectrum of statistical and machine learning challenges that range from study design, data processing, and standardization to analysis, modeling, cross-study comparison, prediction, data science ecosystems, and reproducible reporting. Nevertheless, although many statistics and machine learning approaches and tools have been developed, new techniques are needed to deal with emerging applications and the vast heterogeneity of microbiome data. We review and discuss emerging applications of statistical and machine learning techniques in human microbiome studies and introduce the COST Action CA18131 "ML4Microbiome" that brings together microbiome researchers and machine learning experts to address current challenges such as standardization of analysis pipelines for reproducibility of data analysis results, benchmarking, improvement, or development of existing and new tools and ontologies.

Highly decomposed organic carbon mediates the assembly of soil communities with traits for the biodegradation of chlorinated pollutants

To improve biodegradation strategies for chlorinated pollutants, the roles of soil organic matter and microbial function need to be clarified. It was hypothesised that microbial degradation of specific organic fractions in soils enhance community metabolic capability to degrade chlorinated pollutants. This field study used historic records of dieldrin concentrations since 1988 and established relationships between dieldrin dissipation and soil carbon fractions together with bacterial and fungal diversity in surface soils of Kurosol and Chromosol. Sparse partial least squares analysis linked dieldrin dissipation to metabolic activities associated with the highly decomposed carbon fraction. Dieldrin dissipation, after three decades of natural attenuation, was associated with increased bacterial species fitness for the decomposition of recalcitrant carbon substrates including synthetic chlorinated pollutants. These metabolic capabilities were linked to the decomposed carbon fraction, an important driver for the microbial community and function. Common bacterial traits among taxonomic groups enriched in samples with high dieldrin dissipation included their slow growth, large genome and complex metabolism which supported the notion that metabolic strategies for dieldrin degradation evolved in an energy-low soil environment. The findings provide new perspectives for bioremediation strategies and suggest that soil management should aim at stimulating metabolism at the decomposed, fine carbon fraction.

Analysis of microbial compositions: a review of normalization and differential abundance analysis

Increasingly, researchers are discovering associations between microbiome and a wide range of human diseases such as obesity, inflammatory bowel diseases, HIV, and so on. The first step towards microbiome wide association studies is the characterization of the composition of human microbiome under different conditions. Determination of differentially abundant microbes between two or more environments, known as differential abundance (DA) analysis, is a challenging and an important problem that has received considerable interest during the past decade. It is well documented in the literature that the observed microbiome data (OTU/SV table) are relative abundances with an excess of zeros. Since relative abundances sum to a constant, these data are necessarily compositional. In this article we review some recent methods for DA analysis and describe their strengths and weaknesses.

Worldwide occurrence of haemoplasmas in wildlife: Insights into the patterns of infection, transmission, pathology and zoonotic potential

Haemotropic mycoplasmas (haemoplasmas) have increasingly attracted the attention of wildlife disease researchers due to a combination of wide host range, high prevalence and genetic diversity. A systematic review identified 75 articles that investigated haemoplasma infection in wildlife by molecular methods (chiefly targeting partial 16S rRNA gene sequences), which included 131 host genera across six orders. Studies were less common in the Eastern Hemisphere (especially Africa and Asia) and more frequent in the Artiodactyla and Carnivora. Meta-analysis showed that infection prevalence did not vary by geographic region nor host order, but wild hosts showed significantly higher prevalence than captive hosts. Using a taxonomically flexible machine learning algorithm, we also found vampire bats and cervids to have greater prevalence, whereas mink, a subclade of vesper bats, and true foxes all had lower prevalence compared to the remaining sampled mammal phylogeny. Haemoplasma genotype and nucleotide diversity varied little among wild mammals but were marginally lower in primates and bats. Coinfection with more than one haemoplasma species or genotype was always confirmed when assessed. Risk factors of infection identified were sociality, age, males and high trophic levels, and both prevalence and diversity were often higher in undisturbed environments. Haemoplasmas likely use different and concurrent transmission routes and typically display enzootic dynamics when wild populations are studied longitudinally. Haemoplasma pathology is poorly known in wildlife but appears subclinical. Candidatus Mycoplasma haematohominis, which causes disease in humans, probably has it natural host in bats. Haemoplasmas can serve as a model system in ecological and evolutionary studies, and future research on these pathogens in wildlife must focus on increasing the geographic range and taxa of studies and elucidating pathology, transmission and zoonotic potential. To facilitate such work, we recommend using universal PCR primers or NGS protocols to detect novel haemoplasmas and other genetic markers to differentiate among species and infer cross-species transmission.

Fast-lived Hosts and Zoonotic Risk

Because most emerging human pathogens originate in mammals, many studies aim to identify host traits that determine the risk of sourcing zoonotic outbreaks. Studies regularly assert that 'fast-lived' mammal species exhibiting greater fecundity and shorter lifespans tend to host more zoonoses; however, the causes of this association remain poorly understood and they cover a range of immune and nonimmune mechanisms. We discuss these drivers in the context of evolutionary ecology and wildlife-human interactions. Ultimately, differentiating these mechanisms will require linking interspecific variation in life history with immunity, pathogen diversity, transmissibility, and zoonotic risk, and critical data gaps currently limit our ability to do so. We highlight sampling and analytical frameworks to address this gap and to better inform zoonotic reservoir prediction.

Identifying Suspect Bat Reservoirs of Emerging Infections

Bats host a number of pathogens that cause severe disease and onward transmission in humans and domestic animals. Some of these pathogens, including henipaviruses and filoviruses, are considered a concern for future pandemics. There has been substantial effort to identify these viruses in bats. However, the reservoir hosts for Ebola virus are still unknown and henipaviruses are largely uncharacterized across their distribution. Identifying reservoir species is critical in understanding the viral ecology within these hosts and the conditions that lead to spillover. We collated surveillance data to identify taxonomic patterns in prevalence and seroprevalence and to assess sampling efforts across species. We systematically collected data on filovirus and henipavirus detections and used a machine-learning algorithm, phylofactorization, in order to search the bat phylogeny for cladistic patterns in filovirus and henipavirus infection, accounting for sampling efforts. Across sampled bat species, evidence for filovirus infection was widely dispersed across the sampled phylogeny. We found major gaps in filovirus sampling in bats, especially in Western Hemisphere species. Evidence for henipavirus infection was clustered within the Pteropodidae; however, no other clades have been as intensely sampled. The major predictor of filovirus and henipavirus exposure or infection was sampling effort. Based on these results, we recommend expanding surveillance for these pathogens across the bat phylogenetic tree.

Ecological and evolutionary drivers of haemoplasma infection and bacterial genotype sharing in a Neotropical bat community

Most emerging pathogens can infect multiple species, underlining the importance of understanding the ecological and evolutionary factors that allow some hosts to harbour greater infection prevalence and share pathogens with other species. However, our understanding of pathogen jumps is based primarily around viruses, despite bacteria accounting for the greatest proportion of zoonoses. Because bacterial pathogens in bats (order Chiroptera) can have conservation and human health consequences, studies that examine the ecological and evolutionary drivers of bacterial prevalence and barriers to pathogen sharing are crucially needed. Here were studied haemotropic Mycoplasma spp. (i.e., haemoplasmas) across a species-rich bat community in Belize over two years. Across 469 bats spanning 33 species, half of individuals and two-thirds of species were haemoplasma positive. Infection prevalence was higher for males and for species with larger body mass and colony sizes. Haemoplasmas displayed high genetic diversity (21 novel genotypes) and strong host specificity. Evolutionary patterns supported codivergence of bats and bacterial genotypes alongside phylogenetically constrained host shifts. Bat species centrality to the network of shared haemoplasma genotypes was phylogenetically clustered and unrelated to prevalence, further suggesting rare-but detectable-bacterial sharing between species. Our study highlights the importance of using fine phylogenetic scales when assessing host specificity and suggests phylogenetic similarity may play a key role in host shifts not only for viruses but also for bacteria. Such work more broadly contributes to increasing efforts to understand cross-species transmission and the epidemiological consequences of bacterial pathogens.

Response of the human gut and saliva microbiome to urbanization in Cameroon

Urban populations from highly industrialized countries are characterized by a lower gut bacterial diversity as well as by changes in composition compared to rural populations from less industrialized countries. To unveil the mechanisms and factors leading to this diversity loss, it is necessary to identify the factors associated with urbanization-induced shifts at a smaller geographical scale, especially in less industrialized countries. To do so, we investigated potential associations between a variety of dietary, medical, parasitological and socio-cultural factors and the gut and saliva microbiomes of 147 individuals from three populations along an urbanization gradient in Cameroon. We found that the presence of Entamoeba sp., a commensal gut protozoan, followed by stool consistency, were major determinants of the gut microbiome diversity and composition. Interestingly, urban individuals have retained most of their gut eukaryotic and bacterial diversity despite significant changes in diet compared to the rural areas, suggesting that the loss of bacterial microbiome diversity observed in industrialized areas is likely associated with medication. Finally, we observed a weak positive correlation between the gut and the saliva microbiome diversity and composition, even though the saliva microbiome is mainly shaped by habitat-related factors.

Rare microbial taxa emerge when communities collide: freshwater and marine microbiome responses to experimental mixing

Whole microbial communities regularly merge with one another, often in tandem with their environments, in a process called community coalescence. Such events impose substantial changes: abiotic perturbation from environmental blending and biotic perturbation of community merging. We used an aquatic mixing experiment to unravel the effects of these perturbations on the whole microbiome response and on the success of individual taxa when distinct freshwater and marine communities coalesce. We found that an equal mix of freshwater and marine habitats and blended microbiomes resulted in strong convergence of the community structure toward that of the marine microbiome. The enzymatic potential of these blended microbiomes in mixed media also converged toward that of the marine, with strong correlations between the multivariate response patterns of the enzymes and of community structure. Exposing each endmember inocula to an axenic equal mix of their freshwater and marine source waters led to a 96% loss of taxa from our freshwater microbiomes and a 66% loss from our marine microbiomes. When both inocula were added together to this mixed environment, interactions amongst the communities led to a further loss of 29% and 49% of freshwater and marine taxa, respectively. Under both the axenic and competitive scenarios, the diversity lost was somewhat counterbalanced by increased abundance of microbial taxa that were too rare to detect in the initial inocula. Our study emphasizes the importance of the rare biosphere as a critical component of microbial community responses to community coalescence.

Widespread prevalence but contrasting patterns of intragenomic rRNA polymorphisms in nematodes: Implications for phylogeny, species delimitation and life history inference

Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free-living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population-level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.

The problem of scale in the prediction and management of pathogen spillover

Disease emergence events, epidemics and pandemics all underscore the need to predict zoonotic pathogen spillover. Because cross-species transmission is inherently hierarchical, involving processes that occur at varying levels of biological organization, such predictive efforts can be complicated by the many scales and vastness of data potentially required for forecasting. A wide range of approaches are currently used to forecast spillover risk (e.g. macroecology, pathogen discovery, surveillance of human populations, among others), each of which is bound within particular phylogenetic, spatial and temporal scales of prediction. Here, we contextualize these diverse approaches within their forecasting goals and resulting scales of prediction to illustrate critical areas of conceptual and pragmatic overlap. Specifically, we focus on an ecological perspective to envision a research pipeline that connects these different scales of data and predictions from the aims of discovery to intervention. Pathogen discovery and predictions focused at the phylogenetic scale can first provide coarse and pattern-based guidance for which reservoirs, vectors and pathogens are likely to be involved in spillover, thereby narrowing surveillance targets and where such efforts should be conducted. Next, these predictions can be followed with ecologically driven spatio-temporal studies of reservoirs and vectors to quantify spatio-temporal fluctuations in infection and to mechanistically understand how pathogens circulate and are transmitted to humans. This approach can also help identify general regions and periods for which spillover is most likely. We illustrate this point by highlighting several case studies where long-term, ecologically focused studies (e.g. Lyme disease in the northeast USA, Hendra virus in eastern Australia, Plasmodium knowlesi in Southeast Asia) have facilitated predicting spillover in space and time and facilitated the design of possible intervention strategies. Such studies can in turn help narrow human surveillance efforts and help refine and improve future large-scale, phylogenetic predictions. We conclude by discussing how greater integration and exchange between data and predictions generated across these varying scales could ultimately help generate more actionable forecasts and interventions. This article is part of the theme issue 'Dynamic and integrative approaches to understanding pathogen spillover'.

Percolation models of pathogen spillover

Predicting pathogen spillover requires counting spillover events and aligning such counts with process-related covariates for each spillover event. How can we connect our analysis of spillover counts to simple, mechanistic models of pathogens jumping from reservoir hosts to recipient hosts? We illustrate how the pathways to pathogen spillover can be represented as a directed graph connecting reservoir hosts and recipient hosts and the number of spillover events modelled as a percolation of infectious units along that graph. Percolation models of pathogen spillover formalize popular intuition and management concepts for pathogen spillover, such as the inextricably multilevel nature of cross-species transmission, the impact of covariance between processes such as pathogen shedding and human susceptibility on spillover risk, and the assumptions under which the effect of a management intervention targeting one process, such as persistence of vectors, will translate to an equal effect on the overall spillover risk. Percolation models also link statistical analysis of spillover event datasets with a mechanistic model of spillover. Linear models, one might construct for process-specific parameters, such as the log-rate of shedding from one of several alternative reservoirs, yield a nonlinear model of the log-rate of spillover. The resulting nonlinearity is approximately piecewise linear with major impacts on statistical inferences of the importance of process-specific covariates such as vector density. We recommend that statistical analysis of spillover datasets use piecewise linear models, such as generalized additive models, regression clustering or ensembles of linear models, to capture the piecewise linearity expected from percolation models. We discuss the implications of our findings for predictions of spillover risk beyond the range of observed covariates, a major challenge of forecasting spillover risk in the Anthropocene. This article is part of the theme issue 'Dynamic and integrative approaches to understanding pathogen spillover'.

Scalable methods for analyzing and visualizing phylogenetic placement of metagenomic samples

BACKGROUND

The exponential decrease in molecular sequencing cost generates unprecedented amounts of data. Hence, scalable methods to analyze these data are required. Phylogenetic (or Evolutionary) Placement methods identify the evolutionary provenance of anonymous sequences with respect to a given reference phylogeny. This increasingly popular method is deployed for scrutinizing metagenomic samples from environments such as water, soil, or the human gut.

NOVEL METHODS

Here, we present novel and, more importantly, highly scalable methods for analyzing phylogenetic placements of metagenomic samples. More specifically, we introduce methods for (a) visualizing differences between samples and their correlation with associated meta-data on the reference phylogeny, (b) clustering similar samples using a variant of the k-means method, and (c) finding phylogenetic factors using an adaptation of the Phylofactorization method. These methods enable to interpret metagenomic data in a phylogenetic context, to find patterns in the data, and to identify branches of the phylogeny that are driving these patterns.

RESULTS

To demonstrate the scalability and utility of our methods, as well as to provide exemplary interpretations of our methods, we applied them to 3 publicly available datasets comprising 9782 samples with a total of approximately 168 million sequences. The results indicate that new biological insights can be attained via our methods.