Soil microbiota represent one of the ancient evolutionary origins of antibiotic resistance and have been proposed as a reservoir of resistance genes available for exchange with clinical pathogens. Using a high-throughput functional metagenomic approach in conjunction with a pipeline for the de-novo assembly of short-read sequence data from functional selections (termed PARFuMS), we provide evidence for recent exchange of antibiotic resistance genes between environmental bacteria and clinical pathogens. We describe multidrug-resistant soil bacteria containing resistance cassettes against five classes of antibiotics (β-lactams, aminoglycosides, amphenicols, sulfonamides, and tetracyclines) which have perfect nucleotide identity to genes from diverse human pathogens. This identity encompasses non-coding regions as well as multiple mobilization sequences, offering not only evidence of lateral exchange, but also a mechanism by which antibiotic resistance disseminates.

Summary Ancient and diverse antibiotic resistance genes (ARGs) have previously been identified from soil1–3, including genes identical to those in human pathogens4. Despite the apparent overlap between soil and clinical resistomes4–6, factors influencing ARG composition in soil and their movement between genomes and habitats remain largely unknown3. General metagenome functions often correlate with the underlying structure of bacterial communities7–12. However, ARGs are hypothesized to be highly mobile4,5,13, prompting speculation that resistomes may not correlate with phylogenetic signatures or ecological divisions13,14. To investigate these relationships, we performed functional metagenomic selections for resistance to 18 antibiotics from 18 agricultural and grassland soils. The 2895 ARGs we discovered were predominantly novel, and represent all major resistance mechanisms15. We demonstrate that distinct soil types harbor distinct resistomes, and that nitrogen fertilizer amendments strongly influenced soil ARG content. Resistome composition also correlated with microbial phylogenetic and taxonomic structure, both across and within soil types. Consistent with this strong correlation, mobility elements syntenic with ARGs were rare in soil compared to sequenced pathogens, suggesting that ARGs in the soil may not transfer between bacteria as readily as is observed in the clinic. Together, our results indicate that bacterial community composition is the primary determinant of soil ARG content, challenging previous hypotheses that horizontal gene transfer effectively decouples resistomes from phylogeny13,14.

Sequencing-based DNA methylation profiling methods are comprehensive and, as accuracy and affordability improve, will increasingly supplant microarrays for genome-scale analyses. Here, four sequencing-based methodologies were applied to biological replicates of human embryonic stem cells to compare their CpG coverage genome-wide and in transposons, resolution, cost, concordance and its relationship with CpG density and genomic context. The two bisulfite methods reached concordance of 82% for CpG methylation levels and 99% for non-CpG cytosine methylation levels. Using binary methylation calls, two enrichment methods were 99% concordant, while regions assessed by all four methods were 97% concordant. To achieve comprehensive methylome coverage while reducing cost, an approach integrating two complementary methods was examined. The integrative methylome profile along with histone methylation, RNA, and SNP profiles derived from the sequence reads allowed genome-wide assessment of allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.

Antibiotic resistance is a dire clinical problem with important ecological dimensions. While antibiotic resistance in human pathogens continues to rise at alarming rates, the impact of environmental resistance on human health is still unclear. To investigate the relationship between human-associated and environmental resistomes, we analyzed functional metagenomic selections for resistance against 18 clinically relevant antibiotics from soil and human gut microbiota as well as a set of multidrug-resistant cultured soil isolates. These analyses were enabled by Resfams, a new curated database of protein families and associated highly precise and accurate profile hidden Markov models, confirmed for antibiotic resistance function and organized by ontology. We demonstrate that the antibiotic resistance functions that give rise to the resistance profiles observed in environmental and human-associated microbial communities significantly differ between ecologies. Antibiotic resistance functions that most discriminate between ecologies provide resistance to b-lactams and tetracyclines, two of the most widely used classes of antibiotics in the clinic and agriculture. We also analyzed the antibiotic resistance gene composition of over 6000 sequenced microbial genomes, revealing significant enrichment of resistance functions by both ecology and phylogeny. Together, our results indicate that environmental and human-associated microbial communities harbor distinct resistance genes, suggesting that antibiotic resistance functions are largely constrained by ecology.

Summary Enzymes capable of inactivating tetracycline are paradoxically rare compared with enzymes that inactivate other natural product antibiotics. We describe a family of flavoenzymes, previously unrecognizable as resistance genes, which are capable of degrading tetracycline antibiotics. From soil functional metagenomic selections, we discovered nine genes that confer high-level tetracycline resistance by enzymatic inactivation. We also demonstrate that a tenth enzyme, an uncharacterized homolog in the human pathogen Legionella longbeachae, similarly inactivates tetracycline. These enzymes catalyze the oxidation of tetracyclines in vitro both by known mechanisms and via previously undescribed activity. Tetracycline inactivation genes were identified in diverse soil types, encompass substantial sequence diversity, and are adjacent to genes implicated in horizontal gene transfer. Because tetracycline inactivation is scarcely observed in hospitals, these enzymes may fill an empty niche in pathogenic organisms, and should therefore be monitored for their dissemination potential into the clinic.