Transformers enable accurate prediction of acute and chronic chemical toxicity in aquatic organisms.

Environmental hazard assessments are reliant on toxicity data that cover multiple organism groups. Generating experimental toxicity data is, however, resource-intensive and time-consuming. Computational methods are fast and cost-efficient alternatives, but the low accuracy and narrow applicability domains have made their adaptation slow. Here, we present a AI-based model for predicting chemical toxicity. The model uses transformers to capture toxicity-specific features directly from the chemical structures and deep neural networks to predict effect concentrations. The model showed high predictive performance for all tested organism groups-algae, aquatic invertebrates and fish-and has, in comparison to commonly used QSAR methods, a larger applicability domain and a considerably lower error. When the model was trained on data with multiple effect concentrations (EC50/EC10), the performance was further improved. We conclude that deep learning and transformers have the potential to markedly advance computational prediction of chemical toxicity.

The UNITE database for molecular identification and taxonomic communication of fungi and other eukaryotes: sequences, taxa and classifications reconsidered.

UNITE (https://unite.ut.ee) is a web-based database and sequence management environment for molecular identification of eukaryotes. It targets the nuclear ribosomal internal transcribed spacer (ITS) region and offers nearly 10 million such sequences for reference. These are clustered into ∼2.4M species hypotheses (SHs), each assigned a unique digital object identifier (DOI) to promote unambiguous referencing across studies. UNITE users have contributed over 600 000 third-party sequence annotations, which are shared with a range of databases and other community resources. Recent improvements facilitate the detection of cross-kingdom biological associations and the integration of undescribed groups of organisms into everyday biological pursuits. Serving as a digital twin for eukaryotic biodiversity and communities worldwide, the latest release of UNITE offers improved avenues for biodiversity discovery, precise taxonomic communication and integration of biological knowledge across platforms.

Dataset on aquatic ecotoxicity predictions of 2697 chemicals, using three quantitative structure-activity relationship platforms.

Empirical and in silico data on the aquatic ecotoxicology of 2697 organic chemicals were collected in order to compile a dataset for assessing the predictive power of current Quantitative Structure Activity Relationship (QSAR) models and software platforms. This document presents the dataset and the data pipeline for its creation. Empirical data were collected from the US EPA ECOTOX Knowledgebase (ECOTOX) and the EFSA (European Food Safety Authority) report "Completion of data entry of pesticide ecotoxicology Tier 1 study endpoints in a XML schema - database". Only data for OECD recommended algae, daphnia and fish species were retained. QSAR toxicity predictions were calculated for each chemical and each of six endpoints using ECOSAR, VEGA and the Toxicity Estimation Software Tool (T.E.S.T.) platforms. Finally, the dataset was amended with SMILES, InChIKey, pKa and logP collected from webchem and PubChem.

Extensive screening reveals previously undiscovered aminoglycoside resistance genes in human pathogens.

Antibiotic resistance is a growing threat to human health, caused in part by pathogens accumulating antibiotic resistance genes (ARGs) through horizontal gene transfer. New ARGs are typically not recognized until they have become widely disseminated, which limits our ability to reduce their spread. In this study, we use large-scale computational screening of bacterial genomes to identify previously undiscovered mobile ARGs in pathogens. From ~1 million genomes, we predict 1,071,815 genes encoding 34,053 unique aminoglycoside-modifying enzymes (AMEs). These cluster into 7,612 families (<70% amino acid identity) of which 88 are previously described. Fifty new AME families are associated with mobile genetic elements and pathogenic hosts. From these, 24 of 28 experimentally tested AMEs confer resistance to aminoglycoside(s) in Escherichia coli, with 17 providing resistance above clinical breakpoints. This study greatly expands the range of clinically relevant aminoglycoside resistance determinants and demonstrates that computational methods enable early discovery of potentially emerging ARGs.

How, not if, is the question mycologists should be asking about DNA-based typification.

Fungal metabarcoding of substrates such as soil, wood, and water is uncovering an unprecedented number of fungal species that do not seem to produce tangible morphological structures and that defy our best attempts at cultivation, thus falling outside the scope of the International Code of Nomenclature for algae, fungi, and plants. The present study uses the new, ninth release of the species hypotheses of the UNITE database to show that species discovery through environmental sequencing vastly outpaces traditional, Sanger sequencing-based efforts in a strongly increasing trend over the last five years. Our findings challenge the present stance of some in the mycological community - that the current situation is satisfactory and that no change is needed to "the code" - and suggest that we should be discussing not whether to allow DNA-based descriptions (typifications) of species and by extension higher ranks of fungi, but what the precise requirements for such DNA-based typifications should be. We submit a tentative list of such criteria for further discussion. The present authors hope for a revitalized and deepened discussion on DNA-based typification, because to us it seems harmful and counter-productive to intentionally deny the overwhelming majority of extant fungi a formal standing under the International Code of Nomenclature for algae, fungi, and plants.

Phenotypic and transcriptomic acclimation of the green microalga Raphidocelis subcapitata to high environmental levels of the herbicide diflufenican.

Herbicide pollution poses a worldwide threat to plants and freshwater ecosystems. However, the understanding of how organisms develop tolerance to these chemicals and the associated trade-off expenses are largely unknown. This study aims to investigate the physiological and transcriptional mechanisms underlying the acclimation of the green microalgal model species Raphidocelis subcapitata (Selenastraceae) towards the herbicide diflufenican, and the fitness costs associated with tolerance development. Algae were exposed for 12 weeks (corresponding to 100 generations) to diflufenican at the two environmental concentrations 10 and 310 ng/L. The monitoring of growth, pigment composition, and photosynthetic performance throughout the experiment revealed an initial dose-dependent stress phase (week 1) with an EC50 of 397 ng/L, followed by a time-dependent recovery phase during weeks 2 to 4. After week 4, R. subcapitata was acclimated to diflufenican exposure with a similar growth rate, content of carotenoids, and photosynthetic performance as the unexposed control algae. This acclimation state of the algae was explored in terms of tolerance acquisition, changes in the fatty acids composition, diflufenican removal rate, cell size, and changes in mRNA gene expression profile, revealing potential fitness costs associated with acclimation, such as up-regulation of genes related to cell division, structure, morphology, and reduction of cell size. Overall, this study demonstrates that R. subcapitata can quickly acclimate to environmental but toxic levels of diflufenican; however, the acclimation is associated with trade-off expenses that result in smaller cell size.

Evidence for wastewaters as environments where mobile antibiotic resistance genes emerge.

The emergence and spread of mobile antibiotic resistance genes (ARGs) in pathogens have become a serious threat to global health. Still little is known about where ARGs gain mobility in the first place. Here, we aimed to collect evidence indicating where such initial mobilization events of clinically relevant ARGs may have occurred. We found that the majority of previously identified origin species did not carry the mobilizing elements that likely enabled intracellular mobility of the ARGs, suggesting a necessary interplay between different bacteria. Analyses of a broad range of metagenomes revealed that wastewaters and wastewater-impacted environments had by far the highest abundance of both origin species and corresponding mobilizing elements. Most origin species were only occasionally detected in other environments. Co-occurrence of origin species and corresponding mobilizing elements were rare in human microbiota. Our results identify wastewaters and wastewater-impacted environments as plausible arenas for the initial mobilization of resistance genes.

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.

Strain-level bacterial typing directly from patient samples using optical DNA mapping.

BACKGROUND: Identification of pathogens is crucial to efficiently treat and prevent bacterial infections. However, existing diagnostic techniques are slow or have a too low resolution for well-informed clinical decisions. METHODS: In this study, we have developed an optical DNA mapping-based method for strain-level bacterial typing and simultaneous plasmid characterisation. For the typing, different taxonomical resolutions were examined and cultivated pure Escherichia coli and Klebsiella pneumoniae samples were used for parameter optimization. Finally, the method was applied to mixed bacterial samples and uncultured urine samples from patients with urinary tract infections. RESULTS: We demonstrate that optical DNA mapping of single DNA molecules can identify Escherichia coli and Klebsiella pneumoniae at the strain level directly from patient samples. At a taxonomic resolution corresponding to E. coli sequence type 131 and K. pneumoniae clonal complex 258 forming distinct groups, the average true positive prediction rates are 94% and 89%, respectively. The single-molecule aspect of the method enables us to identify multiple E. coli strains in polymicrobial samples. Furthermore, by targeting plasmid-borne antibiotic resistance genes with Cas9 restriction, we simultaneously identify the strain or subtype and characterize the corresponding plasmids. CONCLUSION: The optical DNA mapping method is accurate and directly applicable to polymicrobial and clinical samples without cultivation. Hence, it has the potential to rapidly provide comprehensive diagnostics information, thereby optimizing early antibiotic treatment and opening up for future precision medicine management.

For bacterial infections, it is important to rapidly and accurately identify and characterize the type of bacteria involved so that optimal antibiotic treatment can be given quickly to the patient. However, current diagnostic methods are sometimes slow and cannot be used for mixtures of bacteria. We have, therefore, developed a method to identify bacteria directly from patient samples. The method was tested on two common species of disease-causing bacteria ­ Escherichia coli and Klebsiella pneumoniae ­ and it could correctly identify the bacterial strain or subtype in both urine samples and mixtures. Hence, the method has the potential to provide fast diagnostic information for choosing the most suited antibiotic, thereby reducing the risk of death and suffering.

Economic interests cloud hazard reductions in the European regulation of substances of very high concern.

Here we investigate how the conflicts between hazard reduction and economic interests have shaped the regulation of substances of very high concern (SVHCs) under the Authorization program of the European chemical regulation Registration, Evaluation, Authorization, and Restriction of Chemicals (REACH). Since regulation starts with listing SVHCs on the Candidate List, we analyze the relative importance of toxicological properties, economic motivations, and available scientific knowledge on the probability of inclusion on the Candidate List. We find that the most important factor in whether a substance is listed is whether it is being produced in, or imported into, the European Economic Area (EEA), with the regulators less likely to place a substance on the list if it is currently being produced or imported in the EEA. This evidence suggests that regulators have listed chemicals of secondary importance leading to lesser than anticipated hazard reductions, either because production and imports had already ceased before the listing, or because the compound has never been produced or imported in the EEA at all.

The impact of the ST131 clone on recurrent ESBL-producing E. coli urinary tract infection: a prospective comparative study.

The global emergence of extended-spectrum beta-lactamase-producing Escherichia coli (ESBL-E. coli), mainly causing urinary tract infections (UTI), is of great concern. Almost one third of patients with UTI, develop recurrent UTI (RUTI). We followed 297 patients for one year after their first episode of UTI due to ESBL-E. coli. Our aim was to evaluate the impact of the globally dominant sequence type (ST)131 clone and its clades, on the risk of subsequent recurrences with ESBL-E. coli. Isolates from patients developing RUTI (68/297) were compared with those from patients with sporadic UTI (SUTI, 229/297). No association was found between RUTI and the two most prevalent phylogroups B2 and D, blaCTX-M genes, or resistance profile. Half of the patients with RUTI were infected with ST131 isolates. Clade C2 were in dominance (50/119) among ST131 isolates. They were more common in patients with RUTI than SUTI (28% vs 13%) and multivariate analysis showed an increased odds-ratio (OR = 2.21, p = 0.033) for recurrences in patients infected with these isolates as compared to non-ST131 isolates. Detecting specific biomarkers, as ST131 clade C2, in ESBL-E. coli UTI isolates may aid in prediction of RUTI and improve diagnostics and care of patients with a risk of ESBL-E. coli recurrences.

The curse of the uncultured fungus.

The international DNA sequence databases abound in fungal sequences not annotated beyond the kingdom level, typically bearing names such as "uncultured fungus". These sequences beget low-resolution mycological results and invite further deposition of similarly poorly annotated entries. What do these sequences represent? This study uses a 767,918-sequence corpus of public full-length fungal ITS sequences to estimate what proportion of the 95,055 "uncultured fungus" sequences that represent truly unidentifiable fungal taxa - and what proportion of them that would have been straightforward to annotate to some more meaningful taxonomic level at the time of sequence deposition. Our results suggest that more than 70% of these sequences would have been trivial to identify to at least the order/family level at the time of sequence deposition, hinting that factors other than poor availability of relevant reference sequences explain the low-resolution names. We speculate that researchers' perceived lack of time and lack of insight into the ramifications of this problem are the main explanations for the low-resolution names. We were surprised to find that more than a fifth of these sequences seem to have been deposited by mycologists rather than researchers unfamiliar with the consequences of poorly annotated fungal sequences in molecular repositories. The proportion of these needlessly poorly annotated sequences does not decline over time, suggesting that this problem must not be left unchecked.

Large-scale characterization of the macrolide resistome reveals high diversity and several new pathogen-associated genes.

Macrolides are broad-spectrum antibiotics used to treat a range of infections. Resistance to macrolides is often conferred by mobile resistance genes encoding Erm methyltransferases or Mph phosphotransferases. New erm and mph genes keep being discovered in clinical settings but their origins remain unknown, as is the type of macrolide resistance genes that will appear in the future. In this study, we used optimized hidden Markov models to characterize the macrolide resistome. Over 16 terabases of genomic and metagenomic data, representing a large taxonomic diversity (11 030 species) and diverse environments (1944 metagenomic samples), were searched for the presence of erm and mph genes. From this data, we predicted 28 340 macrolide resistance genes encoding 2892 unique protein sequences, which were clustered into 663 gene families (<70 % amino acid identity), of which 619 (94 %) were previously uncharacterized. This included six new resistance gene families, which were located on mobile genetic elements in pathogens. The function of ten predicted new resistance genes were experimentally validated in Escherichia coli using a growth assay. Among the ten tested genes, seven conferred increased resistance to erythromycin, with five genes additionally conferring increased resistance to azithromycin, showing that our models can be used to predict new functional resistance genes. Our analysis also showed that macrolide resistance genes have diverse origins and have transferred horizontally over large phylogenetic distances into human pathogens. This study expands the known macrolide resistome more than ten-fold, provides insights into its evolution, and demonstrates how computational screening can identify new resistance genes before they become a significant clinical problem.

GEnView: a gene-centric, phylogeny-based comparative genomics pipeline for bacterial genomes and plasmids.

SUMMARY: Comparing genomic loci of a given bacterial gene across strains and species can provide insights into their evolution, including information on e.g. acquired mobility, the degree of conservation between different taxa or indications of horizontal gene transfer events. While thousands of bacterial genomes are available to date, there is no software that facilitates comparisons of individual gene loci for a large number of genomes. GEnView (Genetic Environment View) is a Python-based pipeline for the comparative analysis of gene-loci in a large number of bacterial genomes, providing users with automated, taxon-selective access to the >800.000 genomes and plasmids currently available in the NCBI Assembly and RefSeq databases, and is able to process local genomes that are not deposited at NCBI, enabling searches for genomic sequences and to analyze their genetic environments through the interactive visualization and extensive metadata files created by GEnView. AVAILABILITY AND IMPLEMENTATION: GEnView is implemented in Python 3. Instructions for download and usage can be found at https://github.com/EbmeyerSt/GEnView under GLP3. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Detection of structural variations in densely-labelled optical DNA barcodes: A hidden Markov model approach.

Large-scale genomic alterations play an important role in disease, gene expression, and chromosome evolution. Optical DNA mapping (ODM), commonly categorized into sparsely-labelled ODM and densely-labelled ODM, provides sequence-specific continuous intensity profiles (DNA barcodes) along single DNA molecules and is a technique well-suited for detecting such alterations. For sparsely-labelled barcodes, the possibility to detect large genomic alterations has been investigated extensively, while densely-labelled barcodes have not received as much attention. In this work, we introduce HMMSV, a hidden Markov model (HMM) based algorithm for detecting structural variations (SVs) directly in densely-labelled barcodes without access to sequence information. We evaluate our approach using simulated data-sets with 5 different types of SVs, and combinations thereof, and demonstrate that the method reaches a true positive rate greater than 80% for randomly generated barcodes with single variations of size 25 kilobases (kb). Increasing the length of the SV further leads to larger true positive rates. For a real data-set with experimental barcodes on bacterial plasmids, we successfully detect matching barcode pairs and SVs without any particular assumption of the types of SVs present. Instead, our method effectively goes through all possible combinations of SVs. Since ODM works on length scales typically not reachable with other techniques, our methodology is a promising tool for identifying arbitrary combinations of genomic alterations.

A framework for identifying the recent origins of mobile antibiotic resistance genes.

Since the introduction of antibiotics as therapeutic agents, many bacterial pathogens have developed resistance to antibiotics. Mobile resistance genes, acquired through horizontal gene transfer, play an important role in this process. Understanding from which bacterial taxa these genes were mobilized, and whether their origin taxa share common traits, is critical for predicting which environments and conditions contribute to the emergence of novel resistance genes. This knowledge may prove valuable for limiting or delaying future transfer of novel resistance genes into pathogens. The literature on the origins of mobile resistance genes is scattered and based on evidence of variable quality. Here, we summarize, amend and scrutinize the evidence for 37 proposed origins of mobile resistance genes. Using state-of-the-art genomic analyses, we supplement and evaluate the evidence based on well-defined criteria. Nineteen percent of reported origins did not fulfill the criteria to confidently assign the respective origin. Of the curated origin taxa, >90% have been associated with infection in humans or domestic animals, some taxa being the origin of several different resistance genes. The clinical emergence of these resistance genes appears to be a consequence of antibiotic selection pressure on taxa that are permanently or transiently associated with the human/domestic animal microbiome.

SMAD4 haploinsufficiency in small intestinal neuroendocrine tumors.

BACKGROUND: Patients with small intestinal neuroendocrine tumors (SINETs) frequently present with lymph node and liver metastases at the time of diagnosis, but the molecular changes that lead to the progression of these tumors are largely unknown. Sequencing studies have only identified recurrent point mutations at low frequencies with CDKN1B being the most common harboring heterozygous mutations in less than 10% of all tumors. Although SINETs are genetically stable tumors with a low frequency of point mutations and indels, they often harbor recurrent hemizygous copy number alterations (CNAs) yet the functional implications of these CNA are unclear. METHODS: Utilizing comparative genomic hybridization (CGH) arrays we analyzed the CNA profile of 131 SINETs from 117 patients. Two tumor suppressor genes and corresponding proteins i.e. SMAD4, and CDKN1B, were further characterized using a tissue microarray (TMA) with 846 SINETs. Immunohistochemistry (IHC) was used to quantify protein expression in TMA samples and this was correlated with chromosome number evaluated with fluorescent in-situ hybridization (FISH). Intestinal tissue from a Smad4+/- mouse model was used to detect entero-endocrine cell hyperplasia with IHC. RESULTS: Analyzing the CGH arrays we found loss of chromosome 18q and SMAD4 in 71% of SINETs and that focal loss of chromosome 12 affecting the CDKN1B was present in 9.4% of SINETs. No homozygous loss of chromosome 18 was detected. Hemizygous loss of SMAD4, but not CDKN1B, significantly correlated with reduced protein levels but hemizygous loss of SMAD4 did not induce entero-endocrine cell hyperplasia in the Smad4+/- mouse model. In addition, patients with low SMAD4 protein expression in primary tumors more often presented with metastatic disease. CONCLUSIONS: Hemizygous loss of chromosome 18q and the SMAD4 gene is the most common genetic event in SINETs and our results suggests that this could influence SMAD4 protein expression and spread of metastases. Although SMAD4 haploinsufficiency alone did not induce tumor initiation, loss of chromosome 18 could represent an evolutionary advantage in SINETs explaining the high prevalence of this aberration. Functional consequences of reduced SMAD4 protein levels could hypothetically be a potential mechanism as to why loss of chromosome 18 appears to be clonally selected in SINETs.

An updated phylogeny of the metallo-ß-lactamases.

BACKGROUND: Metallo-ß-lactamases (MBLs) are enzymes that use zinc-dependent hydrolysis to confer resistance to almost all available ß-lactam antibiotics. They are hypothesized to originate from commensal and environmental bacteria, from where some have mobilized and transferred horizontally to pathogens. The current phylogeny of MBLs, however, is biased as it is founded largely on genes encountered in pathogenic bacteria. This incompleteness is emphasized by recent findings of environmental MBLs with new forms of zinc binding sites and atypical functional profiles. OBJECTIVES: To expand the phylogeny of MBLs to provide a more accurate view of their evolutionary history. METHODS: We searched more than 16 terabases of genomic and metagenomic data for MBLs of the three subclasses B1, B2 and B3 using the validated fARGene method. Predicted genes, together with the previously known ones, were used to infer phylogenetic trees. RESULTS: We identified 2290 unique MBL genes forming 817 gene families, of which 741 were previously uncharacterized. MBLs from subclasses B1 and B3 separated into distinct monophyletic groups, in agreement with their taxonomic and functional properties. We present evidence that clinically associated MBLs were mobilized from Proteobacteria. Additionally, we identified three new variants of the zinc binding sites, indicating that the functional repertoire is broader than previously reported. CONCLUSIONS: Based on our results, we recommend that the nomenclature of MBLs is refined into the phylogenetic groups B1.1-B1.5 and B3.1-B3.4 that more accurately describe their molecular and functional characteristics. Our results will also facilitate the annotation of novel MBLs, reflecting their taxonomic organization and evolutionary origin.

Predicting clinical resistance prevalence using sewage metagenomic data.

Antibiotic resistance surveillance through regional and up-to-date testing of clinical isolates is a foundation for implementing effective empirical treatment. Surveillance data also provides an overview of geographical and temporal changes that are invaluable for guiding interventions. Still, due to limited infrastructure and resources, clinical surveillance data is lacking in many parts of the world. Given that sewage is largely made up of human fecal bacteria from many people, sewage epidemiology could provide a cost-efficient strategy to partly fill the current gap in clinical surveillance of antibiotic resistance. Here we explored the potential of sewage metagenomic data to assess clinical antibiotic resistance prevalence using environmental and clinical surveillance data from across the world. The sewage resistome correlated to clinical surveillance data of invasive Escherichia coli isolates, but none of several tested approaches provided a sufficient resolution for clear discrimination between resistance towards different classes of antibiotics. However, in combination with socioeconomic data, the overall clinical resistance situation could be predicted with good precision. We conclude that analyses of bacterial genes in sewage could contribute to informing management of antibiotic resistance.

Comprehensive screening of genomic and metagenomic data reveals a large diversity of tetracycline resistance genes.

Tetracyclines are broad-spectrum antibiotics used to prevent or treat a variety of bacterial infections. Resistance is often mediated through mobile resistance genes, which encode one of the three main mechanisms: active efflux, ribosomal target protection or enzymatic degradation. In the last few decades, a large number of new tetracycline-resistance genes have been discovered in clinical settings. These genes are hypothesized to originate from environmental and commensal bacteria, but the diversity of tetracycline-resistance determinants that have not yet been mobilized into pathogens is unknown. In this study, we aimed to characterize the potential tetracycline resistome by screening genomic and metagenomic data for novel resistance genes. By using probabilistic models, we predicted 1254 unique putative tetracycline resistance genes, representing 195 gene families (<70 % amino acid sequence identity), whereof 164 families had not been described previously. Out of 17 predicted genes selected for experimental verification, 7 induced a resistance phenotype in an Escherichia coli host. Several of the predicted genes were located on mobile genetic elements or in regions that indicated mobility, suggesting that they easily can be shared between bacteria. Furthermore, phylogenetic analysis indicated several events of horizontal gene transfer between bacterial phyla. Our results also suggested that acquired efflux pumps originate from proteobacterial species, while ribosomal protection genes have been mobilized from Firmicutes and Actinobacteria. This study significantly expands the knowledge of known and putatively novel tetracycline resistance genes, their mobility and evolutionary history. The study also provides insights into the unknown resistome and genes that may be encountered in clinical settings in the future.