ABSTRACT Current supervised phylogeny-based methods fall short on recognizing species assembled from metagenomic datasets from under-investigated habitats, as they are often incomplete or lack closely known relatives. Here, we report an efficient software suite, “Genome Constellation”, that estimates similarities between genomes based on their k-mer matches, and subsequently uses these similarities for classification, clustering, and visualization. The clusters of reference genomes formed by Genome Constellation closely resemble known phylogenetic relationships while simultaneously revealing unexpected connections. In a dataset containing 1,693 draft genomes assembled from the Antarctic lake communities where only 40% could be placed in a phylogenetic tree, Genome Constellation improves taxa assignment to 61%. It revealed six clusters derived from new bacterial phyla and 63 new giant viruses, 3 of which missed by the traditional marker-based approach. In summary, we demonstrate that Genome Constellation can tackle the computational and algorithmic challenges in large-scale taxonomy analyses in metagenomics.

ABSTRACT The perennial ice-covered lakes of the Antarctic McMurdo Dry Valleys harbor oligotrophic microbial communities that are separated geographically from other aquatic systems. Their microbiomes include planktonic microbes as well as lift-off mat communities that emerge from the ice. We used ShortBRED to quantify the antibiotic resistance genes (ARGs) from metagenomes of lift-off mats emerging from ice, and from filtered water samples of Lake Fryxell and Lake Bonney. The overall proportion of ARG hits was similar to that found in temperate-zone rural water bodies with moderate human inputs. Specific ARGs showed distinct distributions for the two lakes, and for mat versus planktonic sources. Metagenomic taxa distributions showed that mat phototrophs consisted mainly of cyanobacteria or betaproteobacteria, whereas the water column phototrophs were mainly protists. An enrichment culture of the betaproteobacterium Rhodoferax antarcticus from a Lake Fryxell mat sample showed an unusual mat-forming phenotype not previously reported for this species. Its genome showed no ARGs associated with betaproteobacteria, but had ARGs consistent with a minor Pseudomonas component. The Antarctic lake mats and water showed specific ARGs distinctive to the mat and water sources, but overall ARG levels were similar to those of temperate water bodies with moderate human inputs.

Microbial communities of small freshwater bodies are poorly understood. Four ponds in Knox County, Ohio, were sampled over two years to investigate the relationship between the microbial taxa profiles, antibiotic resistance genes (ARGs), and environmental factors such as pH and tannin concentrations. For each site, microbial communities were collected by filtration and metagenomes were analyzed by short-read sequencing. Taxa profiles were predicted by the Kraken2/Bracken pipelines. Bacterial taxa with high abundance in these ponds included Betaproteobacteria (Polynucleobacter and Methylopumilus) and Actinobacteria (Planktophila, Nanopelagicus, and Mycolicibacterium). One pond, a former quarry with elevated pH, showed high prevalence of Cyanobacteria with a seasonal shift from Synechococcus to Planktothrix in the fall. Planktothrix increase was associated with acidification. ARGs were quantified using the ShortBRED pipeline to detect and quantify hits to a marker set derived from the Comprehensive Antibiotic Resistance Database (CARD). The top two ARGs with the largest marker hits encode components of a Stenotrophomonas drug efflux pump powered by proton-motive force (smeABC) and a mycobacterial global regulator that activates a drug pump and other cell defenses (mtrA). Pump function and global activation of transcription incur large energy expenditures, whose fitness cost may increase at high external pH where the cell's proton-motive force is diminished. The smeABC and mtrA prevalence showed a modest correlation with acidifying conditions (low pH and high tannins) which contribute a large transmembrane pH difference to the proton-motive force, thus increasing the cell's energy available for pump function and global gene expression. Compared to rivers and lakes, pond microbial ecosystems are understudied despite close contact with agriculture and recreation. Environmental microbes offer health benefits as well as hazards for human contact. Small water bodies may act as reservoirs for drug-resistant organisms and transfer of antibiotic resistance genes. Yet, the public is rarely aware of the potential for exposure to ARG-carrying organisms in recreational water bodies. Little is known about the capacity for freshwater microbial communities to remediate drug pollution and which biochemical factors may select against antibiotic resistance genes. This study analyzes the bacterial taxa composition and ARG prevalence including possible influence of factors such as pH and tannic acid levels.

ABSTRACT The entry of antibiotic resistance genes (ARGs) into aquatic systems has been documented for large municipal wastewater treatment plants, but there is less study of the impact of smaller plants that are situated on small rural rivers. We sampled water metagenomes for ARG and taxa composition from the Kokosing River, a small rural river in Knox County, Ohio, which has been designated an Ohio State Scenic River for retention of natural character. Samples were obtained 1.0 km upstream, 120 m downstream, and 6.4 km downstream from the effluent release of the Mount Vernon wastewater treatment plant (WWTP). ARGS were identified in metagenomes using ShortBRED markers from the CARD database screened against UniPROT. Through all seasons, the metagenome just downstream of the WWTP effluent showed a substantial elevation of at least 15 different ARGs, including 6 ARGs commonly associated with Acinetobacter baumannii such as msrE, mphE (macrolide resistance) and tet (39) (tetracycline resistance). The ARGs most prevalent near the effluent pipe persisted 6.4 km downriver. Using MetaPhlAn2 clade-specific marker genes, the taxa distribution near the effluent showed elevation of reads annotated as Acinetobacter species as well as gut-associated taxa, Bacteroides and Firmicutes. The ARG levels and taxa prevalence showed little dependence on seasonal chlorination of the effluent. Nitrogen and phosphorus were elevated near the effluent pipe but had no consistent correlation with ARG levels. We show that in a rural river microbiome, year-round wastewater effluent substantially elevates ARGs including those associated with multidrug-resistant A. baumanii . IMPORTANCE Antibiotic resistance is a growing problem worldwide, with frequent transmission between pathogens and environmental organisms. Rural rivers can support high levels of recreational use by people unaware of inputs from treated wastewater, while WWTPs can generate a small but significant portion of flow volume into a river surrounded by forest and agriculture. There is little information on the rural impacts of WWTP effluent on the delivery and transport of antibiotic resistance genes. In our study, the river water proximal to wastewater effluent shows evidence for the influx of multidrug-resistant Acinetobacter baumanii , an opportunistic pathogen of concern for hospitals but also widespread in natural environments. Our work highlights the importance of wastewater effluent in management of environmental antibiotic resistance, even in high quality, rural river systems.