Structure and function of the global topsoil microbiome.

Soils harbour some of the most diverse microbiomes on Earth and are essential for both nutrient cycling and carbon storage. To understand soil functioning, it is necessary to model the global distribution patterns and functional gene repertoires of soil microorganisms, as well as the biotic and environmental associations between the diversity and structure of both bacterial and fungal soil communities1-4. Here we show, by leveraging metagenomics and metabarcoding of global topsoil samples (189 sites, 7,560 subsamples), that bacterial, but not fungal, genetic diversity is highest in temperate habitats and that microbial gene composition varies more strongly with environmental variables than with geographic distance. We demonstrate that fungi and bacteria show global niche differentiation that is associated with contrasting diversity responses to precipitation and soil pH. Furthermore, we provide evidence for strong bacterial-fungal antagonism, inferred from antibiotic-resistance genes, in topsoil and ocean habitats, indicating the substantial role of biotic interactions in shaping microbial communities. Our results suggest that both competition and environmental filtering affect the abundance, composition and encoded gene functions of bacterial and fungal communities, indicating that the relative contributions of these microorganisms to global nutrient cycling varies spatially.

CAFE: a software suite for analysis of paired-sample transposon insertion sequencing data.

SUMMARY: Sequencing of transposon insertion libraries is used to determine the relative fitness of individual mutants at a large scale. However, there is a lack of tools for specifically analyzing data from such experiments with paired sample designs. Here, we introduce CAFE-Coefficient-based Analysis of Fitness by read Enrichment-a software package that can analyze data from paired transposon mutant sequencing experiments, generate fitness coefficients for each gene and condition and perform appropriate statistical testing on these fitness coefficients. AVAILABILITY AND IMPLEMENTATION: CAFE is implemented in Perl and R. The source code is freely available for download under the MIT License from https://github.com/bengtssonpalme/cafe and http://microbiology.se/software/cafe/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

What Is the Role of the Environment in the Emergence of Novel Antibiotic Resistance Genes? A Modeling Approach.

It is generally accepted that intervention strategies to curb antibiotic resistance cannot solely focus on human and veterinary medicine but must also consider environmental settings. While the environment clearly has a role in transmission of resistant bacteria, its role in the emergence of novel antibiotic resistance genes (ARGs) is less clear. It has been suggested that the environment constitutes an enormous recruitment ground for ARGs to pathogens, but its extent is practically unknown. We have constructed a model framework for resistance emergence and used available quantitative data on relevant processes to identify limiting steps in the appearance of ARGs in human pathogens. We found that in a majority of possible scenarios, the environment would only play a minor role in the emergence of novel ARGs. However, the uncertainty is enormous, highlighting an urgent need for more quantitative data. Specifically, more data is most needed on the fitness costs of ARG carriage, the degree of dispersal of resistant bacteria from the environment to humans, and the rates of mobilization and horizontal transfer of ARGs. This type of data is instrumental to determine which processes should be targeted for interventions to curb development and transmission of ARGs in the environment.

The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications.

UNITE (https://unite.ut.ee/) is a web-based database and sequence management environment for the molecular identification of fungi. It targets the formal fungal barcode-the nuclear ribosomal internal transcribed spacer (ITS) region-and offers all ∼1 000 000 public fungal ITS sequences for reference. These are clustered into ∼459 000 species hypotheses and assigned digital object identifiers (DOIs) to promote unambiguous reference across studies. In-house and web-based third-party sequence curation and annotation have resulted in more than 275 000 improvements to the data over the past 15 years. UNITE serves as a data provider for a range of metabarcoding software pipelines and regularly exchanges data with all major fungal sequence databases and other community resources. Recent improvements include redesigned handling of unclassifiable species hypotheses, integration with the taxonomic backbone of the Global Biodiversity Information Facility, and support for an unlimited number of parallel taxonomic classification systems.

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.

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.

Computational discovery and functional validation of novel fluoroquinolone resistance genes in public metagenomic data sets.

BACKGROUND: Fluoroquinolones are broad-spectrum antibiotics used to prevent and treat a wide range of bacterial infections. Plasmid-mediated qnr genes provide resistance to fluoroquinolones in many bacterial species and are increasingly encountered in clinical settings. Over the last decade, several families of qnr genes have been discovered and characterized, but their true prevalence and diversity still remain unclear. In particular, environmental and host-associated bacterial communities have been hypothesized to maintain a large and unknown collection of qnr genes that could be mobilized into pathogens. RESULTS: In this study we used computational methods to screen genomes and metagenomes for novel qnr genes. In contrast to previous studies, we analyzed an almost 20-fold larger dataset comprising almost 13 terabases of sequence data. In total, 362,843 potential qnr gene fragments were identified, from which 611 putative qnr genes were reconstructed. These gene sequences included all previously described plasmid-mediated qnr gene families. Fifty-two of the 611 identified qnr genes were reconstructed from metagenomes, and 20 of these were previously undescribed. All of the novel qnr genes were assembled from metagenomes associated with aquatic environments. Nine of the novel genes were selected for validation, and six of the tested genes conferred consistently decreased susceptibility to ciprofloxacin when expressed in Escherichia coli. CONCLUSIONS: The results presented in this study provide additional evidence for the ubiquitous presence of qnr genes in environmental microbial communities, expand the number of known qnr gene variants and further elucidate the diversity of this class of resistance genes. This study also strengthens the hypothesis that environmental bacterial communities act as sources of previously uncharacterized qnr genes.

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.

FARAO: the flexible all-round annotation organizer.

With decreasing costs of generating DNA sequence data, genome and metagenome projects have become accessible to a wider scientific community. However, to extract meaningful information and visualize the data remain challenging. We here introduce FARAO, a highly scalable software for organization, visualization and integration of annotation and read coverage data that can also combine output data from several bioinformatics tools. The capabilities of FARAO can greatly aid analyses of genomic and metagenomic datasets. AVAILABILITY AND IMPLEMENTATION: FARAO is implemented in Perl and is supported under Unix-like operative systems, including Linux and macOS. The Perl source code is freely available for download under the MIT License from http://microbiology.se/software/farao/ CONTACT: johan.bengtsson-palme@microbiology.seSupplementary information: Supplementary data are available at Bioinformatics online.

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.

Strategies to improve usability and preserve accuracy in biological sequence databases.

Biology is increasingly dependent on large-scale analysis, such as proteomics, creating a requirement for efficient bioinformatics. Bioinformatic predictions of biological functions rely upon correctly annotated database sequences, and the presence of inaccurately annotated or otherwise poorly described sequences introduces noise and bias to biological analyses. Accurate annotations are, for example, pivotal for correct identification of polypeptide fragments. However, standards for how sequence databases are organized and presented are currently insufficient. Here, we propose five strategies to address fundamental issues in the annotation of sequence databases: (i) to clearly separate experimentally verified and unverified sequence entries; (ii) to enable a system for tracing the origins of annotations; (iii) to separate entries with high-quality, informative annotation from less useful ones; (iv) to integrate automated quality-control software whenever such tools exist; and (v) to facilitate postsubmission editing of annotations and metadata associated with sequences. We believe that implementation of these strategies, for example as requirements for publication of database papers, would enable biology to better take advantage of large-scale data.

BacMet: antibacterial biocide and metal resistance genes database.

Antibiotic resistance has become a major human health concern due to widespread use, misuse and overuse of antibiotics. In addition to antibiotics, antibacterial biocides and metals can contribute to the development and maintenance of antibiotic resistance in bacterial communities through co-selection. Information on metal and biocide resistance genes, including their sequences and molecular functions, is, however, scattered. Here, we introduce BacMet (http://bacmet.biomedicine.gu.se)--a manually curated database of antibacterial biocide- and metal-resistance genes based on an in-depth review of the scientific literature. The BacMet database contains 470 experimentally verified resistance genes. In addition, the database also contains 25 477 potential resistance genes collected from public sequence repositories. All resistance genes in the BacMet database have been organized according to their molecular function and induced resistance phenotype.

Co-occurrence of resistance genes to antibiotics, biocides and metals reveals novel insights into their co-selection potential.

BACKGROUND: Antibacterial biocides and metals can co-select for antibiotic resistance when bacteria harbour resistance or tolerance genes towards both types of compounds. Despite numerous case studies, systematic and quantitative data on co-occurrence of such genes on plasmids and chromosomes is lacking, as is knowledge on environments and bacterial taxa that tend to carry resistance genes to such compounds. This effectively prevents identification of risk scenarios. Therefore, we aimed to identify general patterns for which biocide/metal resistance genes (BMRGs) and antibiotic resistance genes (ARGs) that tend to occur together. We also aimed to quantify co-occurrence of resistance genes in different environments and taxa, and investigate to what extent plasmids carrying both types of genes are conjugative and/or are carrying toxin-antitoxin systems. RESULTS: Co-occurrence patterns of resistance genes were derived from publicly available, fully sequenced bacterial genomes (n = 2522) and plasmids (n = 4582). The only BMRGs commonly co-occurring with ARGs on plasmids were mercury resistance genes and the qacE∆1 gene that provides low-level resistance to quaternary ammonium compounds. Novel connections between cadmium/zinc and macrolide/aminoglycoside resistance genes were also uncovered. Several clinically important bacterial taxa were particularly prone to carry both BMRGs and ARGs. Bacteria carrying BMRGs more often carried ARGs compared to bacteria without (p < 0.0001). BMRGs were found in 86 % of bacterial genomes, and co-occurred with ARGs in 17 % of the cases. In contrast, co-occurrences of BMRGs and ARGs were rare on plasmids from all external environments (<0.7 %) but more common on those of human and domestic animal origin (5 % and 7 %, respectively). Finally, plasmids with both BMRGs and ARGs were more likely to be conjugative (p < 0.0001) and carry toxin-antitoxin systems (p < 0.0001) than plasmids without resistance genes. CONCLUSIONS: This is the first large-scale identification of compounds, taxa and environments of particular concern for co-selection of resistance against antibiotics, biocides and metals. Genetic co-occurrences suggest that plasmids provide limited opportunities for biocides and metals to promote horizontal transfer of antibiotic resistance through co-selection, whereas ample possibilities exist for indirect selection via chromosomal BMRGs. Taken together, the derived patterns improve our understanding of co-selection potential between biocides, metals and antibiotics, and thereby provide guidance for risk-reducing actions.

The Human Gut Microbiome as a Transporter of Antibiotic Resistance Genes between Continents.

Previous studies of antibiotic resistance dissemination by travel have, by targeting only a select number of cultivable bacterial species, omitted most of the human microbiome. Here, we used explorative shotgun metagenomic sequencing to address the abundance of >300 antibiotic resistance genes in fecal specimens from 35 Swedish students taken before and after exchange programs on the Indian peninsula or in Central Africa. All specimens were additionally cultured for extended-spectrum beta-lactamase (ESBL)-producing enterobacteria, and the isolates obtained were genome sequenced. The overall taxonomic diversity and composition of the gut microbiome remained stable before and after travel, but there was an increasing abundance of Proteobacteria in 25/35 students. The relative abundance of antibiotic resistance genes increased, most prominently for genes encoding resistance to sulfonamide (2.6-fold increase), trimethoprim (7.7-fold), and beta-lactams (2.6-fold). Importantly, the increase observed occurred without any antibiotic intake. Of 18 students visiting the Indian peninsula, 12 acquired ESBL-producing Escherichia coli, while none returning from Africa were positive. Despite deep sequencing efforts, the sensitivity of metagenomics was not sufficient to detect acquisition of the low-abundant genes responsible for the observed ESBL phenotype. In conclusion, metagenomic sequencing of the intestinal microbiome of Swedish students returning from exchange programs in Central Africa or the Indian peninsula showed increased abundance of genes encoding resistance to widely used antibiotics.

Metagenomics reveals that detoxification systems are underrepresented in marine bacterial communities.

BACKGROUND: Environmental shotgun sequencing (metagenomics) provides a new way to study communities in microbial ecology. We here use sequence data from the Global Ocean Sampling (GOS) expedition to investigate toxicant selection pressures revealed by the presence of detoxification genes in marine bacteria. To capture a broad range of potential toxicants we selected detoxification protein families representing systems protecting microorganisms from a variety of stressors, such as metals, organic compounds, antibiotics and oxygen radicals. RESULTS: Using a bioinformatics procedure based on comparative analysis to finished bacterial genomes we found that the amount of detoxification genes present in marine microorganisms seems surprisingly small. The underrepresentation is particularly evident for toxicant transporters and proteins involved in detoxifying metals. Exceptions are enzymes involved in oxidative stress defense where peroxidase enzymes are more abundant in marine bacteria compared to bacteria in general. In contrast, catalases are almost completely absent from the open ocean environment, suggesting that peroxidases and peroxiredoxins constitute a core line of defense against reactive oxygen species (ROS) in the marine milieu. CONCLUSIONS: We found no indication that detoxification systems would be generally more abundant close to the coast compared to the open ocean. On the contrary, for several of the protein families that displayed a significant geographical distribution, like peroxidase, penicillin binding transpeptidase and divalent ion transport protein, the open ocean samples showed the highest abundance. Along the same lines, the abundance of most detoxification proteins did not increase with estimated pollution. The low level of detoxification systems in marine bacteria indicate that the majority of marine bacteria have a low capacity to adapt to increased pollution. Our study exemplifies the use of metagenomics data in ecotoxicology, and in particular how anthropogenic consequences on life in the sea can be examined.

Low coverage of species constrains the use of DNA barcoding to assess mosquito biodiversity.

Mosquitoes (Culicidae) represent the main vector insects globally, and they also inhabit many of the terrestrial and aquatic habitats of the world. DNA barcoding and metabarcoding are now widely used in both research and routine practices involving mosquitoes. However, these methodologies rely on information available in databases consisting of barcode sequences representing taxonomically identified voucher specimens. In this study, we assess the availability of public data for mosquitoes in the main online databases, focusing specifically on the two most widely used DNA barcoding markers in Culicidae: COI and ITS2. In addition, we test hypotheses on possible factors affecting species coverage (i.e., the percentage of species covered in the online databases) for COI in different countries and the occurrence of the DNA barcode gap for COI. Our findings showed differences in the data publicly available in the repositories, with a taxonomic or species coverage of 28.4-30.11% for COI in BOLD + GenBank, and 12.32% for ITS2 in GenBank. Afrotropical, Australian and Oriental biogeographic regions had the lowest coverages, while Nearctic, Palearctic and Oceanian had the highest. The Neotropical region had an intermediate coverage. In general, countries with a higher diversity of mosquitoes and higher numbers of medically important species had lower coverage. Moreover, countries with a higher number of endemic species tended to have a higher coverage. Although our DNA barcode gap analyses suggested that the species boundaries need to be revised in half of the mosquito species available in the databases, additional data must be gathered to confirm these results and to allow explaining the occurrence of the DNA barcode gap. We hope this study can help guide regional species inventories of mosquitoes and the completion of a publicly available reference library of DNA barcodes for all mosquito species.

Aberrant microbiomes are associated with increased antibiotic resistance gene load in hybrid mice.

Antibiotic resistance is a priority public health problem resulting from eco-evolutionary dynamics within microbial communities and their interaction at a mammalian host interface or geographical scale. The links between mammalian host genetics, bacterial gut community, and antimicrobial resistance gene (ARG) content must be better understood in natural populations inhabiting heterogeneous environments. Hybridization, the interbreeding of genetically divergent populations, influences different components of the gut microbial communities. However, its impact on bacterial traits such as antibiotic resistance is unknown. Here, we present that hybridization might shape bacterial communities and ARG occurrence. We used amplicon sequencing to study the gut microbiome and to predict ARG composition in natural populations of house mice (Mus musculus). We compared gastrointestinal bacterial and ARG diversity, composition, and abundance across a gradient of pure and hybrid genotypes in the European House Mouse Hybrid Zone. We observed an increased overall predicted richness of ARG in hybrid mice. We found bacteria-ARG interactions by their co-abundance and detected phenotypes of extreme abundances in hybrid mice at the level of specific bacterial taxa and ARGs, mainly multidrug resistance genes. Our work suggests that mammalian host genetic variation impacts the gut microbiome and chromosomal ARGs. However, it raises further questions on how the mammalian host genetics impact ARGs via microbiome dynamics or environmental covariates.

Towards a unified paradigm for sequence-based identification of fungi.

The nuclear ribosomal internal transcribed spacer (ITS) region is the formal fungal barcode and in most cases the marker of choice for the exploration of fungal diversity in environmental samples. Two problems are particularly acute in the pursuit of satisfactory taxonomic assignment of newly generated ITS sequences: (i) the lack of an inclusive, reliable public reference data set and (ii) the lack of means to refer to fungal species, for which no Latin name is available in a standardized stable way. Here, we report on progress in these regards through further development of the UNITE database (http://unite.ut.ee) for molecular identification of fungi. All fungal species represented by at least two ITS sequences in the international nucleotide sequence databases are now given a unique, stable name of the accession number type (e.g. Hymenoscyphus pseudoalbidus|GU586904|SH133781.05FU), and their taxonomic and ecological annotations were corrected as far as possible through a distributed, third-party annotation effort. We introduce the term 'species hypothesis' (SH) for the taxa discovered in clustering on different similarity thresholds (97-99%). An automatically or manually designated sequence is chosen to represent each such SH. These reference sequences are released (http://unite.ut.ee/repository.php) for use by the scientific community in, for example, local sequence similarity searches and in the QIIME pipeline. The system and the data will be updated automatically as the number of public fungal ITS sequences grows. We invite everybody in the position to improve the annotation or metadata associated with their particular fungal lineages of expertise to do so through the new Web-based sequence management system in UNITE.

A global baseline for qPCR-determined antimicrobial resistance gene prevalence across environments.

The environment is an important component in the emergence and transmission of antimicrobial resistance (AMR). Despite that, little effort has been made to monitor AMR outside of clinical and veterinary settings. Partially, this is caused by a lack of comprehensive reference data for the vast majority of environments. To enable monitoring to detect deviations from the normal background resistance levels in the environment, it is necessary to establish a baseline of AMR in a variety of settings. In an attempt to establish this baseline level, we here performed a comprehensive literature survey, identifying 150 scientific papers containing relevant qPCR data on antimicrobial resistance genes (ARGs) in environments associated with potential routes for AMR dissemination. The collected data included 1594 samples distributed across 30 different countries and 12 sample types, in a time span from 2001 to 2020. We found that for most ARGs, the typically reported abundances in human impacted environments fell in an interval from 10-5 to 10-3 copies per 16S rRNA, roughly corresponding to one ARG copy in a thousand bacteria. Altogether these data represent a comprehensive overview of the occurrence and levels of ARGs in different environments, providing background data for risk assessment models within current and future AMR monitoring frameworks.

Latent antibiotic resistance genes are abundant, diverse, and mobile in human, animal, and environmental microbiomes.

BACKGROUND: Bacterial communities in humans, animals, and the external environment maintain a large collection of antibiotic resistance genes (ARGs). However, few of these ARGs are well-characterized and thus established in existing resistance gene databases. In contrast, the remaining latent ARGs are typically unknown and overlooked in most sequencing-based studies. Our view of the resistome and its diversity is therefore incomplete, which hampers our ability to assess risk for promotion and spread of yet undiscovered resistance determinants. RESULTS: A reference database consisting of both established and latent ARGs (ARGs not present in current resistance gene repositories) was created. By analyzing more than 10,000 metagenomic samples, we showed that latent ARGs were more abundant and diverse than established ARGs in all studied environments, including the human- and animal-associated microbiomes. The pan-resistomes, i.e., all ARGs present in an environment, were heavily dominated by latent ARGs. In comparison, the core-resistome, i.e., ARGs that were commonly encountered, comprised both latent and established ARGs. We identified several latent ARGs shared between environments and/or present in human pathogens. Context analysis of these genes showed that they were located on mobile genetic elements, including conjugative elements. We, furthermore, identified that wastewater microbiomes had a surprisingly large pan- and core-resistome, which makes it a potentially high-risk environment for the mobilization and promotion of latent ARGs. CONCLUSIONS: Our results show that latent ARGs are ubiquitously present in all environments and constitute a diverse reservoir from which new resistance determinants can be recruited to pathogens. Several latent ARGs already had high mobile potential and were present in human pathogens, suggesting that they may constitute emerging threats to human health. We conclude that the full resistome-including both latent and established ARGs-needs to be considered to properly assess the risks associated with antibiotic selection pressures. Video Abstract.