- Open Access
A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes
Genome Biology volume 23, Article number: 212 (2022)
Earth’s environments harbor complex consortia of microbes that affect processes ranging from host health to biogeochemical cycles. Understanding their evolution and function is limited by an inability to isolate genomes in a high-throughput manner. Here, we present a workflow for bacterial whole-genome sequencing using open-source labware and the OpenTrons robotics platform, reducing costs to approximately $10 per genome. We assess genomic diversity within 45 gut bacterial species from wild-living chimpanzees and bonobos. We quantify intraspecific genomic diversity and reveal divergence of homologous plasmids between hosts. This enables population genetic analyses of bacterial strains not currently possible with metagenomic data alone.
Microbiota are complex mixtures of organisms, with dozens to hundreds of microbial species sharing genes both through ancestry with closely related strains, as well as through horizontal transfer to distantly related lineages [1,2,3]. Understanding how genetic variation arises and changes within these communities is critical if we hope to develop useful models of their evolution .
But despite the tremendous advances in sequencing technology in the past decades, the paired phenomena of within-species strain diversity and between-species horizontal gene transfer still present a challenge to assessing the genetic structure of populations within diverse metagenomes like the mammalian gut. Community metagenome sequencing can rapidly generate massive quantities of data from a microbiome, but with only limited ability to link genetic changes within the same genome or in populations of closely related cells [5, 6].
Mobile DNA elements, especially plasmids, are even more difficult to place in a metagenomic context [7,8,9]. In principle, cultivation offers a much more robust way to explore genomic variation within populations. Although cultivation necessarily introduces bias in the specific taxa that are recovered, and thus cannot replace metagenomic methods for understanding microbial communities, it offers a few clear advantages. By confidently drawing cellular bounds around genes, isolation represents a gold standard for describing genomic diversity and a necessary prerequisite for empirically demonstrating the functional consequences of such variation. For this reason, cultivation has seen renewed interest, with automation and screening techniques being employed to increase the breadth of diversity that can be reasonably assessed. However, such high-throughput approaches typically require enormous investments in capital equipment and labor [10,11,12,13], putting them out of reach for many researchers. This is especially true for those studying non-model systems where the bulk of unstudied microbial diversity is likely to be found. Advances in miniaturization and microfluidic technologies may one day permit rapid high-throughput cultivation from diverse environments [14, 15], though such approaches are not yet widely available. And while conventional isolation techniques using traditional solid media can easily generate thousands of isolates in a short period of time, generating genomic data from this many isolates is still a major barrier to most laboratories.
The recent availability of distributed, open-source laboratory automation and distributed manufacturing technologies suggests a potential solution: adapting high-throughput genome sequencing techniques to relatively inexpensive commercial and in-house-manufactured equipment. In combination with the extremely low per-base cost of modern sequencing, such an approach offers the potential to realize much of the benefits of capital-intensive conventional high-throughput culturing and sequencing pipelines at a fraction of the required investment.
Motivated by our desire to explore genomic evolution in the microbial populations associated with natural mammalian gut microbiomes, we set out to design an inexpensive end-to-end high-throughput genome sequencing protocol that could be easily replicated with a minimum of capital expenditure. While other high-throughput genomic protocols have been published that can reach low marginal costs per genome [16, 17], they typically rely on expensive high-precision robotics and other specialized equipment, thus making them more suited to well-funded laboratories or core facilities. We developed protocols, 3D-printed custom labware, and analysis pipelines to enable cost-effective high-throughput whole-genome sequencing of natural gut microbiota. These methods allowed us to circumvent traditional 16S rRNA-gene or mass spectrometry-based screening approaches, instead using full-genome sequencing to identify all cultivated isolates. Moreover, this approach enabled the generation and assembly of thousands of bacterial genomes from the hominid gut microbiota rapidly and at low cost relative to existing approaches. Results revealed substantial variation in the distribution of strain-level diversity among wild-living chimpanzees and bonobos (Pan), and, importantly, allowed us to link putative plasmids to their specific bacterial hosts across Pan individuals, populations, and subspecies.
For the purposes of validating the workflow, we carried every sample from DNA extraction through to sequencing. Even if, for example, an isolate failed to grow during liquid culture, we did not exclude it from downstream steps. This enabled us to determine appropriate exclusion criteria for future use.
In total, we picked, grew in liquid culture, extracted DNA from, and sequenced 1879 bacterial isolates (mean of 209, standard deviation of 143 per host individual). Of these, 1265 yielded extractions with DNA concentrations above 0.1 ng/µL; 1049 yielded library concentrations ≥ 0.5 ng/µL; 933 yielded ≥ 25 Mbp of sequence; and 715 yielded high-quality assemblies (> 90% complete and < 5% contaminated), 51 medium-quality assemblies (> 50% complete and < 5% contaminated), and 50 low-quality assemblies (≤ 50% complete and < 5% contaminated) (Fig. 1e). In total, 107 of the sequenced libraries gave assemblies that appeared to be contaminated based on CheckM results , indicating that around 10% of picked colonies may have not in fact been single clones.
The primary point of failure in the workflow appeared to be the liquid culture phase: only 67% of isolates yielded DNA concentrations above 0.1 ng/µL, and 31% above 1 ng/µL (Additional File 1: Fig. S1). Low turbidity of many tubes after incubation was consistent with either slow or no growth in liquid media for many of the colonies transferred from plated media. Initial DNA concentration was a good predictor of subsequent performance: 877 of the 1070 libraries with concentrations ≥ 0.5 ng/µL came from DNA extractions with concentrations above 0.1 ng/µL. Seven hundred four of the 715 high-quality assemblies (98%) came from samples with library concentrations above 0.5 ng/µL (Additional File 1: Fig. S2).
Isolate diversity and distribution
Of the 715 fully-assembled isolate genomes, 688 were classified successfully with GTDB-Tk . All 688 were classified as Firmicutes, with most (572) belonging to the Bacilli and 162 to the Clostridia. Together, these accounted for 9 unique taxonomic assignments at the order level, 13 at the level of family, and 30 at the genus level; all 688 genomes were assigned to a species. There was a mean of 10.1 unique GTDB taxon strings (standard deviation = 4.2) recovered per host individual. By far the most common genus among the assembled genomes was Streptococcus (308), followed by Enterococcus “D” group (76), Staphylococcus (53), Clostridium “P” group (44), and Blautia “A” group (42). Taxonomic classifications, assembly statistics, and other metadata for all isolates are presented in Additional File 3.
Sourmash, which assigns taxonomy based on kmer composition of reads rather than assemblies, was able to classify more of the isolates, with 828 being classified to at least the phylum level. 827 were classified to order, 812 to family, 797 to genus, and 778 to species level. These classifications were highly consistent with the full-genome taxonomies, with 99% matching at the phylum and class levels, 97% matching at the order and family levels, 97.0% matching at the genus level, and 94% matching at the species level. It should be noted that the phylogenetic placement using marker genes with the GTDB-Tk is likely to give more accurate results, so the additional taxonomic annotations estimated by Sourmash should be considered tentative. Phylogenetic reconstruction using concatenated marker gene sequences was also largely concordant with taxonomic assignment (Fig. 2).
Mapping metagenomic reads sequenced directly from wild chimpanzee and bonobo fecal samples against the assembled isolate genomes supported an origin from those samples. A mean of 2.97% (SD 1.12%) of metagenomic reads mapped to the isolate genome assemblies.
Isolate population genetics
Assembling individual isolate genomes also allowed us to explore the variation in within-species diversity that might have been hidden by 16S rRNA gene-based screening. Using dRep , we clustered the assembled isolates into 45 clusters sharing genome-wide estimated Average Nucleotide Identity (ANI) of > 95%. These clusters of strains correspond roughly to the convention for delineating bacterial species based on ANI divergence. All-by-all ANI comparisons within these clusters indicated differences in within-cluster diversity among clusters (Fig. 2, Additional File 1: Fig. S3). Clusters showed differences in similarity structure ranging from deep divisions with representatives recovered from multiple host individuals and species (e.g., cluster 7, Fig. 2); to clusters with more or less isogenic clones recovered from within individuals but which differed between individuals (e.g., cluster 16, Fig. 2c); to clusters that were entirely clonal, with identical genomes recovered from multiple individuals within the same host species (e.g., cluster 41, Fig. 2d).
A substantial amount of genomic diversity was observed within groups of genomes sharing identical full-length 16S rRNA gene haplotypes. Contigs containing the 16S rRNA gene were recovered from a total of 689 isolate genomes, with 594 genome assemblies containing only a single unique 16S rRNA gene haplotype and 95 containing more than one (see Additional File 3). Among the 39 unique 16S rRNA gene haplotype groups (i.e., groups of genomes sharing the same unique 16S rRNA gene haplotype), the average estimated genome-wide ANI ranged from 100% to 91.6% (mean = 99.15%, StdDev = 1.92%; Additional File 1: Fig. S4).
Putative plasmid diversity and distribution
One major advantage of bacterial genomes assembled from isolates relative to metagenome-assembled genomes is the ability to confidently associate plasmids with bacterial chromosomes. Here, we leveraged these data to assess the extent to which plasmid communities and sequences have diverged among bacterial lineages and among chimpanzee and bonobo subspecies sampled throughout equatorial Africa (Additional File 1: Fig. S5). We enumerated all plasmids within our isolate assemblies using a recently developed machine learning approach  which leverages gene content to predict plasmids. A scan of all isolate assemblies yielded a total of 516 putative plasmid contigs recovered from 258 individual genome assemblies. These belonged to 245 clusters as calculated by MobMess , with 64 of these clusters containing more than one contig. Only 31 of these putative plasmids matched existing plasmid sequences in the PLSDB plasmid database , indicating that more than 94% of the recovered plasmid sequences were novel.
Consistent with phylogenetic and geographic barriers to exchange of plasmids, we found that plasmids were most often shared between genomes of related taxa found within the same bacterial and Pan host species (Fig. 3A; Additional File 1: Fig. S6). However, in many cases, plasmids containing homologous stretches of DNA were found in genomes from multiple distantly related bacterial taxa (i.e., different families, Additional File 1: Fig. S6). Similarly, we observed several cases where similar plasmid sequences were detected in multiple Pan host individuals (Fig. 3B–D). Across all plasmids, we observed significantly greater nucleotide similarity (i.e., ANI) between homologous regions of plasmids recovered from different Pan individuals from the same Pan species than between homologous regions of plasmids recovered from different Pan individuals from different Pan species (Fig. 4; Mann–Whitney U test p-value = 0.036). Nucleotide similarity was greater still among homologous regions of different plasmids recovered from the same host individual (Fig. 4; Mann–Whitney U test p-value < 0.0001). These results indicate the divergence of homologous plasmid sequences between Pan host species.
Protocol cost estimates
Costs are difficult to estimate and communicate accurately, as purchasing prices and available equipment vary widely among laboratories. However, as one of the primary motivations of this manuscript is to make high-throughput isolate genome sequencing accessible to as many researchers as possible, we give our best estimates for both our required capital investment and per-sample consumable costs (Additional File 2: Table S1) as a point of reference.
Our laboratory already had basic molecular biology equipment, including PCR machines, centrifuges, manual pipettes, and access to a fluorescence plate reader and bead beater. Additional capital expenses required for this protocol included an OpenTrons OT-2 robot with 2 multi-channel pipettes and a magnetic plate expansion module, a strip tube bead beater adapter, and materials costs for the 3D-printed labware; in total, capital expenses amounted to approximately $13,000.
We estimate a per-sample consumables cost of around $10. Of this, liquid culture and DNA extraction account for about $1.50, library preparation around $3, and sequencing around $5. Even with the 36% success rate we observed here, with no culling of failed samples prior to sequencing, this equates to around $25 per high-quality genome assembly.
We developed a workflow for high-throughput bacterial genome sequencing from complex microbiota. Our workflow makes use of custom-designed 3D-printed labware, the relatively inexpensive OpenTrons liquid handling platform, and recently developed methods for Illumina library preparation using highly diluted reagents. Together, this combination of methods allowed the library preparation and whole-genome sequencing of hundreds of bacterial isolates in parallel for a marginal cost of ~ $10 per isolate. Importantly, by reducing per-isolate whole-genome sequencing costs substantially, our workflow alleviates the need for 16S rRNA gene- or mass spectrometry-based approaches for dereplicating bacterial strains prior to whole-genome sequencing. Of the bacterial isolates that grew in liquid culture and yielded appreciable DNA concentrations post-extraction (i.e., “Pass” in Columns 1 and 2 in Fig. 1g), > 80% yielded Hackflex libraries, nearly all of which yielded genome drafts upon sequencing (Fig. 1). Metagenomic data can be used to assemble contiguous sequences within bacterial chromosomes and plasmids, but these data alone struggle to capture fine-scale strain-level diversity and cannot fully determine the distributions of chromosomes and plasmids among bacterial cells. In contrast, whole-genome sequencing of bacterial isolates affords the opportunity to definitively associate plasmids with their bacterial hosts. We demonstrate the utility of this approach by isolating and profiling strain-level bacterial diversity in gut microbiota of wild chimpanzees and bonobos.
Machine-learning classification of assemblies discovered hundreds of previously undescribed plasmids in chimpanzee and bonobo gut bacterial isolates. Analyses of plasmid distributions among bacterial and Pan hosts revealed plasmids shared between distantly related bacterial lineages both within and between Pan individuals (Additional File 1: Fig. S6), consistent with horizontal gene transfer (HGT) within microbiota [13, 23, 24]. Although there were some observed instances of related plasmids being shared between bacterial taxa, the limited sample size of host individuals and restricted window into overall bacterial diversity limits our ability to make broad inferences about distribution. Interestingly, though, many homologous sequences were also shared between Pan-host species (Fig. 3). Of these sequences, sequence divergence between Pan species was significantly higher than that between conspecific Pan individuals (Fig. 4). These results indicate divergence of gut bacterial plasmid sequences between primate-host species lineages. Thus, isolate sequencing enables chromosome- and plasmid-resolved genomic analyses of bacterial species that remain difficult with metagenomic data alone.
Our workflow has several advantages and disadvantages relative to existing approaches for high-throughput bacterial isolation and whole-genome sequencing. One major advantage is its simplicity, as it relies on standard microbiological and molecular biology approaches and is fully automated on the OpenTrons platform. For example, relative to microfluidics-based isolation [25, 26] or single-cell genome sequencing approaches [27, 28], our method is readily applicable by labs without the need for capital-intensive specialized equipment. The equipment costs necessary to execute our full protocol are also dramatically lower than for a number of previously-developed high-throughput genome sequencing workflows that achieve low marginal costs using expensive robotics [16, 17, 29]. Similarly, while Hi-C-based approaches also have the ability to link plasmids with their host bacteria, these methods rely on labor-intensive protocols that crosslink chromatin with formaldehyde, then digested, and re-ligated to isolate covalently linked DNA fragments . Moreover, both droplet and Hi-C approaches typically capture only a fraction of the genome, and they in general do not allow for the retention of isolated cultures for further experimental study. In contrast, a weakness relative to single-cell and Hi-C approaches is that our workflow can only interrogate bacteria that can be cultured and isolated.
The data we report here represent the first two complete full-scale sequencing runs from this protocol, and there are still opportunities for improvement. First, the loss rate could likely be improved through further optimization of pipette accuracy and precision. We also note that approximately one third of bacterial colonies grown in isolation failed to yield appreciable concentrations of DNA (> 0.1 ng/uL), a failure rate that could likely be reduced by screening for isolate growth prior to extraction. Similarly, the rate of contamination (or sequence libraries containing DNA from multiple bacterial types) could be reduced by adding a secondary re-streaked plate culture step rather than picking directly into liquid culture; in this experiment, we chose to pick directly to maximize throughput. Second, although the protocols we provide can in principle be run with very little specific prior training or programming experience, some working knowledge of Python programming in general, and the Opentrons Python API in particular, is helpful. And third, the logistical challenges of moving from hundreds to potentially tens of thousands of samples—including storage, labeling, and in particular sample provenance validation and metadata tracking—are largely unaddressed here. We will be continuing to address each of these issues in future development of these protocols.
To ensure the greatest utility of our workflow for the research community, all protocols and hardware schematics are freely available for public use at https://github.com/tanaes/Moeller_Opentrons_protocol_library , https://github.com/tanaes/opentrons_functions , and https://github.com/CUMoellerLab/Labware . These repositories will be maintained and updated as we make further additions and improvements to the protocols in the future.
The isolates sequenced in this study represent, to our knowledge, the first large-scale compendium of cultured bacterial genomes from wild chimpanzee and bonobo gut microbiomes. The samples used to demonstrate the methods in this paper were far from ideal for purposes of generic cultivation: preserved in RNALater, and stored for many years at − 80 °C, a substantial portion of the diversity in the original samples was most likely no longer viable, resulting in a dearth of some common gut bacterial lineages in our dataset (Fig. 2). Although the entire process of cultivation, from initial inoculation of plates through colony isolation and regrowth in liquid media took place within an anaerobic chamber, we did not recover many expected anaerobes, leading to uncertainty of whether the preservation or growth conditions may have been responsible for their absence. Accordingly, only 3% of the metagenomic diversity detected by sequencing of the fecal samples used for cultivation was present in our dataset of isolate genomes. Fresh samples (or those collected into cryoprotectants specifically for purposes of later cultivation) will no doubt yield a greater diversity of original gut cultivars. However, the genome resources generated from this wild hominid gut bacterial isolate collection complement and enable comparative analyses with existing gut bacterial genome databases for a subset of chimpanzee and bonobo gut bacterial taxa. All isolates generated by this study have been preserved in glycerol stocks and are available upon request for research purposes.
The vast majority of global microbial genomic diversity remains unexplored. While centralized efforts to explore microbial diversity of particular significance to human health are generating enormous amounts of new data, exploration of most other environments most often occurs in a more decentralized fashion, often by researchers with less access to the capital equipment and economies of scale enjoyed by their medically-oriented peers. Many rare or endangered host species are represented in existing collections, often by samples collected for purposes other than microbial cultivation. Such collections represent a potentially vast resource for exploring naturally occurring host-associated microbes if cost-effective methods exist to access them. The protocols presented here expand the accessibility of high-throughput microbial genomics, thereby increasing the diversity of environments from which microbial isolates and reference genomes can be obtained. Given the interconnectedness of microbial genomic diversity in nature, expanding the breadth of such data will be of substantial benefit to researchers studying microbes from all sorts of environments.
To accomplish our goals of maximum isolate genome throughput with minimal capital and labor costs, we developed a workflow based around the OpenTrons OT-2 robotic liquid handling platform (Fig. 1). This instrument allows for repeatable automation of many protocols, while costing less than $10,000 as configured. Where possible, we took advantage of previously-published low-cost molecular biology protocols, adapting them for automation on the OpenTrons platform. All the protocols described here are available at https://github.com/tanaes/Moeller_Opentrons_protocol_library . In addition, we wrote extensions to the OpenTrons Protocol API to improve certain aspects of instrument behavior, especially relating to use with magnetic bead protocols. An installable library of these extensions is available at https://github.com/tanaes/opentrons_functions .
For some protocols, we found that there were key steps that would require laboratory apparatus that were either not available for commercial purchase, uncommon in a typical molecular biology lab, or would require substantial investment. For these steps, we designed our own versions suitable for rapid manufacture with 3D printers and/or laser cutters and using inexpensive commodity electronic components. These apparatus were designed using Fusion360 CAD software (Autodesk, Inc.). Full source files and component lists can be found at https://github.com/CUMoellerLab/Labware . All data in this paper were generated using versions of the 3D-printed apparatus described below, rather than commercially purchased alternatives. Commercial alternatives to all lab-built equipment are listed in Additional File 2: Table S2.
Sample collection, storage, and metagenome sequencing
To develop our high-throughput sequencing methodology, we used a set of 10 fecal specimens collected from wild chimpanzees and bonobos throughout equatorial Africa (Additional File 1: Fig. S5) between July 2003 and August 2014 (Additional File 2: Table S3). Samples were collected in the field at Pan nest sites, preserved and shipped in RNALater at room temperature, then frozen at − 80 °C for long-term storage [33,34,35,36,37,38,39]. These samples were expected to yield only a limited taxonomic fraction of the original microbial population: preservation in RNALater renders many types of microorganisms nonviable. However, our previous work has shown that some portions of the native microbiome from mammalian fecal samples remain viable even after long-term storage in this preservative , thereby allowing cultivation of a subset of the original community. Thus, this harsh selective filter against most bacteria was expected to facilitate the recovery of a common subset of gut bacterial species from multiple chimpanzee and bonobo hosts.
For metagenome sequencing, samples were centrifuged and approximately 50 mg of material removed from the pellet for DNA extraction. We extracted metagenomic DNA from pellets using the Qiagen PowerSoil extraction kit. Libraries were generated from metagenomic DNA using the “Illumina Equivalent” library prep method at the Cornell Biotechnology Research Center, and pooled libraries sequenced using an Illumina NovaSeq instrument at the UC Davis Sequencing Center.
Ten fecal samples—three from P. paniscus, three from P. t. troglodytes, and four from P. t. schweinfurthii—were selected for cultivation. To capture as much of the bacterial diversity that remained viable in RNAlater as possible, we used several different media for cultivation: Yeast Casitone Fatty Acids (YCFA), YCFA + Starch, Bifidobacterium selective media (BSM), Brain heart infusion-supplemented (BHIS), and Bacteroides Bile Esculin (BBE) (Additional File 2: Table S4). Recipes for all media were derived from . For each sample-by-medium combination, 100 μl of fecal material suspended in RNAlater was plated in an anaerobic chamber (Coy brand) on solid media. Plates were inoculated and incubated at 37 °C for five days in an anaerobic (5% hydrogen, 5% carbon dioxide, and 90% nitrogen) chamber (Coy Lab Products Inc). During each round of cultivation, blank control plates were kept in the anaerobic chamber along with the swabbed experimental plates to check for unintentional environmental contamination.
Liquid culture of picked colonies represented a potential throughput bottleneck, especially if isolates were cultured in conventional glass test tubes. To increase throughput, we instead grew colonies in 1.2 mL 96-place strip tube racks, which have the footprint and well spacing necessary for processing on the OpenTrons liquid handler. Individual colonies were picked by hand from plates into 900 µL of liquid media (Additional File 2: Table S4) using a sterile wooden toothpick without removal from the anaerobic chamber. Then, plates were incubated at 37 °C in the anaerobic chamber for 4 days.
To improve growth in liquid culture for cells that might benefit from increased waste gas diffusion or nutrient distribution, we designed small single-plate orbital shakers to fit inside our anaerobic incubator (Fig. 1b). Adapting an existing open-source design (https://learn.adafruit.com/crickit-lab-shaker/3d-printing), we simplified the electronic components, relocated all connections and controls to the front of the apparatus to facilitate use within the incubator, and changed it to use 5 V USB input for power, allowing us to use a single USB charger to power 7 individual shakers within the incubator.
Following anaerobic incubation, plates were removed from the chamber and 300 µL of media per tube was transferred to a clean deep-well plate and cells pelleted in a centrifuge at 16,000 g. After removal of supernatant, cells were resuspended in glycerol buffer and stored at -80 °C for future use.
Kit-based DNA extraction protocols typically cost between $3 and $5 per sample. For 16S rRNA gene amplicon-based screening, this step can sometimes be omitted with a chemical lysis prior to amplification. For whole-genome screening, we judged that the added complexity of a DNA extraction step was necessary. To reduce costs, we adapted the magnetic bead-based extraction methodology from Oberacker et al. , which uses laboratory-made reagents and either purchased or lab-made magnetic beads, for use on the OpenTrons platform.
For cell lysis, we chose to use beadbeating to ensure lysis of a broad range of bacterial cell types. We designed a 3D-printed and laser-cut loading system to precisely load 0.2 mm glass beads directly into the 96-well strip-tube plates (Fig. 1b) after pelleting cells and removing liquid media. After bead loading, 800 µL of guanidine HCL lysis buffer was added to the tubes, and they were capped and shaken on an Omni Bead Ruptor Elite at 6.5 m/s for 40 s. The tubes were then spun down on a centrifuge at 400 × g for 5 min, decapped, and then moved to the OpenTrons instrument for the remainder of the extraction. The detailed OpenTrons extraction protocol can be found in the project repository linked above. Briefly, the robot transfers 600 µL of lysate to a new plate, adds magnetic beads in a PEG-based binding buffer, and then goes through a series of magnetic binding and wash steps before eluting the extracted DNA in nuclease-free water.
We found extraction efficiency was greatly improved by gently agitating magnetic beads during the initial binding step. To accomplish this, we designed a 3D-printed rotator (Fig. 1c) with attachments for holding 96-well plates or microcentrifuge tubes. After transferring lysate and adding beads and binding buffer on the liquid handler, we programmed a pause to allow the user to remove the plate, seal it, and place it on the rotator for 10 min. Following this step, the plate was unsealed and returned to the liquid handler for the remainder of the protocol.
Extracted DNA was quantified in 384-well plates using a reduced-volume version of the QuantiFluor (Promega) fluorescence-based assay. Four 96-well plates (each the output from a single extraction protocol) were tested in each assay, using an OpenTrons protocol for sample transfer and a Tecan Infinite M200 plate reader for quantification.
Library prep and sequencing
To inexpensively generate sequencing libraries from thousands of DNA extractions, we adapted the Hackflex library prep protocol  to the OpenTrons liquid handler. Briefly, this protocol dilutes key reagents from the Illumina Library Prep protocol to stretch a single kit across more samples. Our adaptation of the protocol changes some reagent quantities to better fit the constraints of the OpenTrons format; for details, see the full protocol in the project repository linked above.
For the libraries presented here, we used barcoded library amplification primers purchased from the Cornell Biotechnology Resource Center. Initially, these shared a single i5 index per library plate, with unique i7 primers per sample. For later libraries, we switched to unique dual indexed (UDI) primers, with 96 unique i5 and i7 primers per plate. To facilitate multiplexing across library prep plates with UDIs, we created a version of the protocol to cycle column matches between i5 and i7 primer plates, allowing up to 12 library plates to be multiplexed without repeating an index combination. Libraries were amplified using 17 cycles of PCR prior to bead-based dual-sided size selection and final elution.
Final libraries were quantified by QuantiFluor (Promega) in 96 well plates, then pooled according to the following algorithm: the volume necessary to transfer 5 ng of library DNA was calculated; for samples requiring more than 10 µL to reach 5 ng transferred (likely failed libraries), 1 µL was transferred; for samples requiring less than 0.5 µL, 0.5 µL was transferred. Per-plate pools were combined and concentrated using magnetic beads and then provided to the Cornell Biotechnology Resource Center for sequencing on an Illumina NextSeq 500 instrument. Two separate sequence runs were performed, combining 13 and 10 library prep plates, respectively.
Isolate genomes derived from fecal samples preserved in RNAlater, which selects for a subset of the bacteria in the Pan microbiota , enabled analyses of intraspecific bacterial genomic diversity within and among Pan individuals and populations. To demonstrate the utility of isolate genomes for strain-level analyses of microbiota, we focused on two sets of analyses of genomic diversity that remain difficult or not possible with shotgun metagenomic data alone.
First, we characterized intraspecific bacterial genomic diversity both within and between individual Pan hosts. Isolate sequences were processed using the Bactopia pipeline . This pipeline does sequence trimming and QC with FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/), assembled sequences with Shovill (https://github.com/tseemann/shovill) and SKESA , and performs assembly quality checking with CheckM . Gene prediction and annotation was performed with Prokka . Taxonomy of each isolate sequence library was estimated from unassembled reads using Sourmash with the GTDB R06-RS202 LCA-formatted database and k = 31 . To create a phylogeny of isolates, assembled genomes predicted to be less than 5% contaminated with CheckM were processed using PhyloPhlAn2  using the Amphora2 marker set  and the “Fast / High Diversity” default settings. To estimate the relative abundances of isolates in original samples, we used CoverM (https://github.com/wwood/CoverM) to calculate coverage for each isolate genome in each of the available chimpanzee metagenomes. 16S rRNA gene sequences were recovered from isolate genome assemblies using PhyloFlash  as implemented in Bactopia. To identify unique haplotypes, the complete alignment of all 16S rRNA gene sequences from PhyloFlash was manually inspected for positions corresponding to the commonly-used 27F and 1492R primers and trimmed to just the portion internal to those priming sites.
Second, we identified plasmids within each isolate genome assembly to assess the distribution of and similarity of these mobile elements among bacterial and chimpanzee hosts. Putative plasmid contigs were identified using PlasX , which uses a machine learning algorithm to classify mobile elements based on their gene content. Open reading frames from contigs were annotated in Anvi’o  and supplied to the PlasX algorithm. Contigs scored 0.90 or higher by PlasX were considered putative plasmids. Circularity of putative plasmids was assessed based on read mapping according to the algorithm used by Yu, Fogarty, and Eren . To assess novelty relative to current known plasmid sequences, putative plasmid sequences were searched against the latest version of the PLSDB plasmid database (v. 2021_06_23_v2, ) using BLASTn , and considered as previously observed if they matched at least one database sequence across at least 90% of the query length at ≥ 60% sequence identity. Finally, putative plasmid similarity networks and gene alignment visualizations were generated using MobMess .
Availability of data and materials
All raw sequence data from this publication are available in the Qiita data repository, study number 14410 (https://qiita.ucsd.edu/study/description/14410), as well as at the EBI ENA repository with accession number ERP136830 . Opentrons protocols are available at https://github.com/CUMoellerLab/Moeller_Opentrons_protocol_library  and custom function library at https://github.com/CUMoellerLab/opentrons_functions . Printable labware files and assembly instructions are available at https://github.com/CUMoellerLab/Labware . Isolates are available upon request.
Lloyd-Price J, Mahurkar A, Rahnavard G, Crabtree J, Orvis J, Hall AB, et al. Strains, functions and dynamics in the expanded Human Microbiome Project. Nature. 2017;550:61–6.
Yan Y, Nguyen LH, Franzosa EA, Huttenhower C. Strain-level epidemiology of microbial communities and the human microbiome. Genome Med. 2020;12:71.
Van Rossum T, Ferretti P, Maistrenko OM, Bork P. Diversity within species: interpreting strains in microbiomes. Nat Rev Microbiol. 2020;18:491–506.
Segata N. On the Road to Strain-Resolved Comparative Metagenomics. mSystems. 2018;3. Available from: https://doi.org/10.1128/mSystems.00190-17
Sczyrba A, Hofmann P, Belmann P, Koslicki D, Janssen S, Dröge J, et al. Critical Assessment of Metagenome Interpretation-a benchmark of metagenomics software. Nat Methods. 2017;14:1063–71.
Meyer F, Fritz A, Deng Z-L, Koslicki D, Gurevich A, Robertson G, et al. Critical Assessment of Metagenome Interpretation - the second round of challenges. bioRxiv. bioRxiv; 2021. Available from: http://biorxiv.org/lookup/doi/10.1101/2021.07.12.451567.
Arredondo-Alonso S, Willems RJ, van Schaik W, Schürch AC. On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data. Microb Genom. 2017;3:e000128.
Burton JN, Liachko I, Dunham MJ, Shendure J. Species-Level Deconvolution of Metagenome Assemblies with Hi-C–Based Contact Probability Maps. G3 Genes|Genomes|Genetics. 2014. 1339–46. Available from: http://dx.doi.org/10.1534/g3.114.011825.
Marbouty M, Baudry L, Cournac A, Koszul R. Scaffolding bacterial genomes and probing host-virus interactions in gut microbiome by proximity ligation (chromosome capture) assay. Sci Adv. 2017;3:e1602105.
Lagier J-C, Dubourg G, Million M, Cadoret F, Bilen M, Fenollar F, et al. Culturing the human microbiota and culturomics. Nat Rev Microbiol. 2018;16:540–50.
Liu C, Du M-X, Abuduaini R, Yu H-Y, Li D-H, Wang Y-J, et al. Enlightening the taxonomy darkness of human gut microbiomes with a cultured biobank. Microbiome. 2021;9:119.
Poyet M, Groussin M, Gibbons SM, Avila-Pacheco J, Jiang X, Kearney SM, et al. A library of human gut bacterial isolates paired with longitudinal multiomics data enables mechanistic microbiome research. Nat Med. 2019;25:1442–52.
Groussin M, Poyet M, Sistiaga A, Kearney SM, Moniz K, Noel M, et al. Elevated rates of horizontal gene transfer in the industrialized human microbiome. Cell. 2021;184:2053-67.e18.
Fitzsimons MS, Novotny M, Lo C-C, Dichosa AEK, Yee-Greenbaum JL, Snook JP, et al. Nearly finished genomes produced using gel microdroplet culturing reveal substantial intraspecies genomic diversity within the human microbiome. Genome Res. 2013;23:878–88.
Jalili-Firoozinezhad S, Gazzaniga FS, Calamari EL, Camacho DM, Fadel CW, Nestor B, et al. Complex human gut microbiome cultured in anaerobic human intestine chips [Internet]. Available from: https://doi.org/10.1101/421404
Perez-Sepulveda BM, Heavens D, Pulford CV, Predeus AV, Low R, Webster H, et al. An accessible, efficient and global approach for the large-scale sequencing of bacterial genomes. Genome Biol [Internet]. Genome Biol. 2021;22. Available from: https://pubmed.ncbi.nlm.nih.gov/34930397/ [Cited 2022 Jun 22]
Sanders JG, Nurk S, Salido RA, Minich J, Xu ZZ, Zhu Q, et al. Optimizing sequencing protocols for leaderboard metagenomics by combining long and short reads. Genome Biol. 2019;20:226.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2019; Available from: https://doi.org/10.1093/bioinformatics/btz848
Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 2017;11:2864–8.
Yu MK, Fogarty EC, Eren AM. The genetic and ecological landscape of plasmids in the human gut. bioRxiv. bioRxiv; 2020. Available from: http://biorxiv.org/lookup/doi/10.1101/2020.11.01.361691.
Galata V, Fehlmann T, Backes C, Keller A. PLSDB: a resource of complete bacterial plasmids. Nucleic Acids Res. 2019;47:D195-202.
Brito IL. Examining horizontal gene transfer in microbial communities. Nat Rev Microbiol. 2021;19:442–53.
Kent AG, Vill AC, Shi Q, Satlin MJ, Brito IL. Widespread transfer of mobile antibiotic resistance genes within individual gut microbiomes revealed through bacterial Hi-C. Nat Commun. 2020;11:4379.
Ma L, Kim J, Hatzenpichler R, Karymov MA, Hubert N, Hanan IM, et al. Gene-targeted microfluidic cultivation validated by isolation of a gut bacterium listed in Human Microbiome Project’s Most Wanted taxa. Proc Natl Acad Sci. 2014:9768–73. Available from: https://doi.org/10.1073/pnas.1404753111
Villa MM, Bloom RJ, Silverman JD, Durand HK, Jiang S, Wu A, et al. Interindividual Variation in Dietary Carbohydrate Metabolism by Gut Bacteria Revealed with Droplet Microfluidic Culture. mSystems. 2020;5. Available from: https://doi.org/10.1128/mSystems.00864-19
Belton J-M, McCord RP, Gibcus JH, Naumova N, Zhan Y, Dekker J. Hi-C: a comprehensive technique to capture the conformation of genomes. Methods. 2012;58:268–76.
Stepanauskas R. Single cell genomics: an individual look at microbes. Curr Opin Microbiol. 2012:613–20. Available from: https://doi.org/10.1016/j.mib.2012.09.001
Mora-Castilla S, To C, Vaezeslami S, Morey R, Srinivasan S, Dumdie JN, et al. Miniaturization Technologies for Efficient Single-Cell Library Preparation for Next-Generation Sequencing. J Lab Autom. 2016;21:557–67.
Sanders J, MoellerLab. CUMoellerLab/Moeller_Opentrons_protocol_library: Sanders et al 2022. Zenodo; 2022. Available from: https://zenodo.org/record/7114630
Sanders J. CUMoellerLab/opentrons_functions: Sanders et al 2022. Zenodo; 2022. Available from: https://zenodo.org/record/7114632
Sanders J. CUMoellerLab/Labware: Sanders et al 2022. Zenodo; 2022. Available from: https://zenodo.org/record/7114625
Goldman SL, Sanders JG, Yan W, Denice A, Cornwall M, Ivey KN, et al. Culture‐enriched community profiling improves resolution of the vertebrate gut microbiota. Mol Ecol Res. 2022:122–36. Available from: https://doi.org/10.1111/1755-0998.13456
Keele BF, Van Heuverswyn F, Li Y, Bailes E, Takehisa J, Santiago ML, et al. Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science. 2006;313:523–6.
Keele BF, Jones JH, Terio KA, Estes JD, Rudicell RS, Wilson ML, et al. Increased mortality and AIDS-like immunopathology in wild chimpanzees infected with SIVcpz. Nature. 2009;460:515–9.
Rudicell RS, Holland Jones J, Wroblewski EE, Learn GH, Li Y, Robertson JD, et al. Impact of simian immunodeficiency virus infection on chimpanzee population dynamics. PLoS Pathog. 2010;6:e1001116.
Liu W, Li Y, Learn GH, Rudicell RS, Robertson JD, Keele BF, et al. Origin of the human malaria parasite Plasmodium falciparum in gorillas. Nature. 2010;467:420–5.
Liu W, Sherrill-Mix S, Learn GH, Scully EJ, Li Y, Avitto AN, et al. Wild bonobos host geographically restricted malaria parasites including a putative new Laverania species. Nat Commun. 2017;8:1635.
Li Y, Ndjango J-B, Learn GH, Ramirez MA, Keele BF, Bibollet-Ruche F, et al. Eastern chimpanzees, but not bonobos, represent a simian immunodeficiency virus reservoir. J Virol. 2012;86:10776–91.
Oberacker P, Stepper P, Bond DM, Höhn S, Focken J, Meyer V, et al. Bio-On-Magnetic-Beads (BOMB): Open platform for high-throughput nucleic acid extraction and manipulation. PLoS Biol. 2019;17:e3000107.
Gaio D, Anantanawat K, To J, Liu M, Monahan L, Darling AE. Hackflex: low-cost, high-throughput, Illumina Nextera Flex library construction. Microb Genom. 2022;8. Available from: https://doi.org/10.1099/mgen.0.000744
Petit RA 3rd, Read TD. Bactopia: a Flexible Pipeline for Complete Analysis of Bacterial Genomes. mSystems. 2020;5. Available from: https://doi.org/10.1128/mSystems.00190-20
Souvorov A, Agarwala R, Lipman DJ. SKESA: strategic k-mer extension for scrupulous assemblies. Genome Biol. 2018;19:153.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.
Brown CT, Titus Brown C, Irber L. sourmash: a library for MinHash sketching of DNA. J Open Source Software. 2016:27. Available from: https://doi.org/10.21105/joss.00027
Asnicar F, Thomas AM, Beghini F, Mengoni C, Manara S, Manghi P, et al. Precise phylogenetic analysis of microbial isolates and genomes from metagenomes using PhyloPhlAn 3.0. Nat Commun. 2020;11:2500.
Wu M, Scott AJ. Phylogenomic analysis of bacterial and archaeal sequences with AMPHORA2. Bioinformatics. 2012;28:1033–4.
Gruber-Vodicka HR, Seah BKB, Pruesse E. phyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes. mSystems. 2020;5. Available from: https://doi.org/10.1128/mSystems.00920-20
Eren AM, Esen ÖC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi’o: an advanced analysis and visualization platform for ’omics data. PeerJ. 2015;3:e1319.
Basic Local Alignment Search Tool (BLAST). Bioinformatics Funct Genom. 100–38. Available from: https://doi.org/10.1002/9780470451496.ch4
Sanders JG, Yan W, Mjungu D, Lonsdorf EV, Hart JA, Sanz CM, et al. Inexpensive high-throughput genome screening of isolate cultures. European Nucleotide Archive. 2022. Available from: https://www.ebi.ac.uk/ena/browser/view/PRJEB52146?show=reads [Cited 2022 Sep 26]
We would like to thank Ilana Brito and all the members of the Moeller Lab for helpful discussions; Yingying Li for quality control of fecal samples; and the Ithaca Generator makerspace, for contributing tools and expertise that aided in the design and production of the open-source labware in this study.
The review history is available as Additional file 4.
Peer review information
Kevin Pang was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
This work was supported by a laboratory start-up grant from Cornell University to AHM and grant R35 GM138284 from the National Institute of General Medical Sciences to AHM, and grant R01 AI 050529 to BHH. JGS was partially supported by grant T32 AI145821 to the Cornell Institute for Host–Microbe Interaction and Disease by the National Institutes of Health.
Ethics approval and consent to participate
Consent for publication
JGS is the founder and sole proprietor of Lightweight Labware LLC. All other authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1. Alluvial plot of protocol efficiency. Figure S2. Relationship between assembly quality and library concentration. Figure S3. Intraspecific nucleotide diversity. Figure S4. Genomic dissimilarity within 16S haplotypes.
Additional file 2: Table S1. Equipment sourcing and alternatives. Table S2. Sample information. Table S3. Culturing information. Table S4. Cost estimates.
About this article
Cite this article
Sanders, J.G., Yan, W., Mjungu, D. et al. A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes. Genome Biol 23, 212 (2022). https://doi.org/10.1186/s13059-022-02777-w